ENH: Regression test for film_gls (will fail until film!2 is merged)

unit_tests/film_gls/design.con

+/ContrastName1	Generation 
+/ContrastName2	Shadowing 
+/ContrastName3	Mean 
+/ContrastName4	Shad > Gen 
+/ContrastName5	Gen > Shad 
+/NumWaves	4
+/NumContrasts	5
+/PPheights		2.498861e-01	2.406885e-01	1.254726e-01	2.252375e-01	2.252375e-01
+/RequiredEffect		1.395	1.303	1.165	1.333	1.333
+
+/Matrix
+1.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 
+0.000000e+00 0.000000e+00 1.000000e+00 0.000000e+00 
+1.000000e+00 0.000000e+00 1.000000e+00 0.000000e+00 
+-1.000000e+00 0.000000e+00 1.000000e+00 0.000000e+00 
+1.000000e+00 0.000000e+00 -1.000000e+00 0.000000e+00 

unit_tests/film_gls/design.fts

+/NumWaves	5
+/NumContrasts	1
+
+/Matrix
+1 1 0 0 0 

unit_tests/film_gls/design_even.mat

+/NumWaves 4
+/NumPoints 50
+/Matrix
+1.648834999999999884e-02 -3.726071999999999978e-03 -4.833967000000000125e-02 -1.963018999999999875e-03
+1.119840000000000071e-02 -3.639564000000000080e-03 -4.980388000000000176e-02 -1.761013999999999982e-03
+5.874693000000000276e-03 -3.561888999999999941e-03 -5.085856999999999878e-02 3.680828000000000138e-03
+5.319972999999999759e-04 -3.439620999999999911e-03 -4.143511999999999890e-02 1.017567000000000055e-01
+-4.756011000000000093e-03 -3.213473000000000111e-03 1.536263999999999963e-01 1.921240000000000100e-02
+-9.853109999999999993e-03 -2.989758999999999819e-03 -2.774908999999999883e-03 -9.302396999999999749e-02
+-1.489293000000000040e-02 -2.856495000000000128e-03 -3.159589000000000147e-02 6.295581999999999567e-02
+-1.992050000000000071e-02 -2.696741999999999802e-03 1.240711000000000036e-01 3.985261000000000359e-03
+-2.483733000000000121e-02 -2.466271999999999822e-03 -2.327842000000000106e-02 -5.903798999999999852e-02
+-2.959615000000000162e-02 -2.225098999999999914e-03 6.898188999999999883e-03 6.379331000000000584e-02
+-3.421665000000000123e-02 -1.957060999999999957e-03 1.050973000000000046e-01 -2.023018000000000038e-02
+-3.862000000000000155e-02 5.954985000000000139e-02 -3.314361999999999864e-02 -8.629162999999999417e-02
+7.960121000000000557e-02 6.756821999999999828e-02 -6.654564999999999786e-02 -1.346018000000000041e-02
+9.585565999999999531e-02 -6.029484999999999706e-02 -5.899767000000000206e-02 8.135866000000000278e-03
+-4.101391000000000064e-02 -8.127914999999999446e-02 -4.923609000000000319e-02 8.894826999999999595e-02
+-7.207805999999999935e-02 7.568362000000000353e-03 1.198999000000000037e-01 6.050200000000000022e-02
+-3.246686000000000016e-02 1.030601000000000017e-01 7.213057999999999970e-02 -8.124520000000000353e-02
+1.290009000000000017e-01 1.234730999999999980e-03 -4.204758000000000107e-02 -6.670086999999999544e-02
+-2.872722000000000123e-02 -6.461886000000000019e-02 -6.029729999999999812e-02 -4.615990999999999844e-03
+-5.131997000000000121e-03 6.184020999999999962e-02 -5.023683000000000315e-02 8.354531999999999584e-03
+9.098035000000000150e-02 8.121815000000000284e-02 -4.258347999999999994e-02 4.637561000000000114e-03
+1.570882999999999863e-01 5.537315000000000119e-03 -3.998582999999999998e-02 1.284144999999999991e-03
+1.044230000000000019e-01 -4.815475999999999784e-02 -3.904910999999999799e-02 2.360751999999999924e-04
+6.108828000000000175e-02 1.124179999999999971e-02 -3.855113999999999763e-02 -6.901035000000000242e-05
+1.255222000000000004e-01 -5.403618000000000310e-02 -3.822646999999999851e-02 8.235082999999999986e-02
+-4.584987000000000090e-02 -1.079753999999999992e-01 1.271100000000000008e-01 7.520477000000000412e-02
+-9.599312000000000122e-02 -1.721912999999999916e-02 1.125185999999999964e-01 1.300007000000000072e-02
+-8.781260000000000454e-02 1.405311000000000060e-02 1.535006999999999899e-01 -6.825562000000000296e-02
+-7.509161999999999781e-02 4.979760999999999910e-02 -2.375682000000000141e-02 -1.096218999999999943e-01
+5.075144999999999657e-03 8.936959999999999349e-02 -6.483840999999999910e-02 -1.650048000000000140e-02
+1.000737000000000015e-01 -2.654621999999999885e-02 -5.569791999999999799e-02 9.277593000000000645e-03
+-4.787783000000000316e-02 -8.976066000000000600e-02 -4.525785000000000224e-02 6.143817999999999675e-03
+-8.509134000000000142e-02 -1.175694000000000045e-02 -4.242430999999999980e-02 6.635983000000000463e-04
+-7.849005000000000565e-02 1.149177000000000011e-02 -4.295538000000000134e-02 -1.457778999999999976e-03
+-6.894811000000000689e-02 9.960926999999999629e-03 -4.436257999999999868e-02 -2.035997999999999801e-03
+-6.503551999999999944e-02 2.326276000000000038e-02 -4.604478000000000043e-02 -2.305683000000000211e-03
+-2.873697000000000057e-02 9.938784999999999981e-02 -4.798501000000000155e-02 -2.483525999999999789e-03
+1.288446000000000036e-01 9.734888999999999348e-02 -5.001557000000000225e-02 -2.632933999999999969e-03
+1.672249000000000096e-01 -2.095228999999999842e-02 -5.224695000000000034e-02 -2.758229000000000063e-03
+8.970435999999999688e-02 -3.740351000000000092e-02 -5.451968000000000092e-02 -2.798154000000000179e-03
+9.215199999999999780e-02 2.984886999999999957e-02 -5.682232999999999729e-02 -2.860592999999999816e-03
+1.492319000000000007e-01 -5.370757000000000309e-02 -5.921124999999999999e-02 5.828088000000000019e-02
+-1.320214999999999930e-02 -1.116783999999999971e-01 6.077806000000000181e-02 6.865992999999999402e-02
+-7.841309999999999947e-02 -2.620781000000000152e-02 7.869441000000000630e-02 -5.900455999999999757e-02
+-7.245508000000000504e-02 1.375199000000000039e-02 -5.671285999999999689e-02 -5.013501999999999759e-04
+-5.751353000000000021e-02 1.412625999999999986e-02 7.872091000000000505e-02 6.318824000000000662e-02
+-5.022292999999999896e-02 7.876944000000000307e-03 7.018169999999999964e-02 1.913267000000000093e-02
+-4.749503000000000058e-02 4.945217000000000016e-03 1.175366000000000050e-01 3.301392000000000215e-02
+-4.596125000000000210e-02 4.169728000000000323e-02 1.365811000000000108e-01 -6.322851000000000166e-02
+3.033071000000000031e-02 8.680093999999999310e-02 -8.620730999999999533e-03 -1.161920000000000036e-01

unit_tests/film_gls/design_odd.mat

+/NumWaves 4
+/NumPoints 51
+/Matrix
+1.648834999999999884e-02 -3.726071999999999978e-03 -4.833967000000000125e-02 -1.963018999999999875e-03
+1.119840000000000071e-02 -3.639564000000000080e-03 -4.980388000000000176e-02 -1.761013999999999982e-03
+5.874693000000000276e-03 -3.561888999999999941e-03 -5.085856999999999878e-02 3.680828000000000138e-03
+5.319972999999999759e-04 -3.439620999999999911e-03 -4.143511999999999890e-02 1.017567000000000055e-01
+-4.756011000000000093e-03 -3.213473000000000111e-03 1.536263999999999963e-01 1.921240000000000100e-02
+-9.853109999999999993e-03 -2.989758999999999819e-03 -2.774908999999999883e-03 -9.302396999999999749e-02
+-1.489293000000000040e-02 -2.856495000000000128e-03 -3.159589000000000147e-02 6.295581999999999567e-02
+-1.992050000000000071e-02 -2.696741999999999802e-03 1.240711000000000036e-01 3.985261000000000359e-03
+-2.483733000000000121e-02 -2.466271999999999822e-03 -2.327842000000000106e-02 -5.903798999999999852e-02
+-2.959615000000000162e-02 -2.225098999999999914e-03 6.898188999999999883e-03 6.379331000000000584e-02
+-3.421665000000000123e-02 -1.957060999999999957e-03 1.050973000000000046e-01 -2.023018000000000038e-02
+-3.862000000000000155e-02 5.954985000000000139e-02 -3.314361999999999864e-02 -8.629162999999999417e-02
+7.960121000000000557e-02 6.756821999999999828e-02 -6.654564999999999786e-02 -1.346018000000000041e-02
+9.585565999999999531e-02 -6.029484999999999706e-02 -5.899767000000000206e-02 8.135866000000000278e-03
+-4.101391000000000064e-02 -8.127914999999999446e-02 -4.923609000000000319e-02 8.894826999999999595e-02
+-7.207805999999999935e-02 7.568362000000000353e-03 1.198999000000000037e-01 6.050200000000000022e-02
+-3.246686000000000016e-02 1.030601000000000017e-01 7.213057999999999970e-02 -8.124520000000000353e-02
+1.290009000000000017e-01 1.234730999999999980e-03 -4.204758000000000107e-02 -6.670086999999999544e-02
+-2.872722000000000123e-02 -6.461886000000000019e-02 -6.029729999999999812e-02 -4.615990999999999844e-03
+-5.131997000000000121e-03 6.184020999999999962e-02 -5.023683000000000315e-02 8.354531999999999584e-03
+9.098035000000000150e-02 8.121815000000000284e-02 -4.258347999999999994e-02 4.637561000000000114e-03
+1.570882999999999863e-01 5.537315000000000119e-03 -3.998582999999999998e-02 1.284144999999999991e-03
+1.044230000000000019e-01 -4.815475999999999784e-02 -3.904910999999999799e-02 2.360751999999999924e-04
+6.108828000000000175e-02 1.124179999999999971e-02 -3.855113999999999763e-02 -6.901035000000000242e-05
+1.255222000000000004e-01 -5.403618000000000310e-02 -3.822646999999999851e-02 8.235082999999999986e-02
+-4.584987000000000090e-02 -1.079753999999999992e-01 1.271100000000000008e-01 7.520477000000000412e-02
+-9.599312000000000122e-02 -1.721912999999999916e-02 1.125185999999999964e-01 1.300007000000000072e-02
+-8.781260000000000454e-02 1.405311000000000060e-02 1.535006999999999899e-01 -6.825562000000000296e-02
+-7.509161999999999781e-02 4.979760999999999910e-02 -2.375682000000000141e-02 -1.096218999999999943e-01
+5.075144999999999657e-03 8.936959999999999349e-02 -6.483840999999999910e-02 -1.650048000000000140e-02
+1.000737000000000015e-01 -2.654621999999999885e-02 -5.569791999999999799e-02 9.277593000000000645e-03
+-4.787783000000000316e-02 -8.976066000000000600e-02 -4.525785000000000224e-02 6.143817999999999675e-03
+-8.509134000000000142e-02 -1.175694000000000045e-02 -4.242430999999999980e-02 6.635983000000000463e-04
+-7.849005000000000565e-02 1.149177000000000011e-02 -4.295538000000000134e-02 -1.457778999999999976e-03
+-6.894811000000000689e-02 9.960926999999999629e-03 -4.436257999999999868e-02 -2.035997999999999801e-03
+-6.503551999999999944e-02 2.326276000000000038e-02 -4.604478000000000043e-02 -2.305683000000000211e-03
+-2.873697000000000057e-02 9.938784999999999981e-02 -4.798501000000000155e-02 -2.483525999999999789e-03
+1.288446000000000036e-01 9.734888999999999348e-02 -5.001557000000000225e-02 -2.632933999999999969e-03
+1.672249000000000096e-01 -2.095228999999999842e-02 -5.224695000000000034e-02 -2.758229000000000063e-03
+8.970435999999999688e-02 -3.740351000000000092e-02 -5.451968000000000092e-02 -2.798154000000000179e-03
+9.215199999999999780e-02 2.984886999999999957e-02 -5.682232999999999729e-02 -2.860592999999999816e-03
+1.492319000000000007e-01 -5.370757000000000309e-02 -5.921124999999999999e-02 5.828088000000000019e-02
+-1.320214999999999930e-02 -1.116783999999999971e-01 6.077806000000000181e-02 6.865992999999999402e-02
+-7.841309999999999947e-02 -2.620781000000000152e-02 7.869441000000000630e-02 -5.900455999999999757e-02
+-7.245508000000000504e-02 1.375199000000000039e-02 -5.671285999999999689e-02 -5.013501999999999759e-04
+-5.751353000000000021e-02 1.412625999999999986e-02 7.872091000000000505e-02 6.318824000000000662e-02
+-5.022292999999999896e-02 7.876944000000000307e-03 7.018169999999999964e-02 1.913267000000000093e-02
+-4.749503000000000058e-02 4.945217000000000016e-03 1.175366000000000050e-01 3.301392000000000215e-02
+-4.596125000000000210e-02 4.169728000000000323e-02 1.365811000000000108e-01 -6.322851000000000166e-02
+3.033071000000000031e-02 8.680093999999999310e-02 -8.620730999999999533e-03 -1.161920000000000036e-01
+1.250539999999999985e-01 -2.703016999999999920e-02 -9.495515999999999679e-02 5.906059999999999782e-02

unit_tests/film_gls/feedsRun.autoregressive

+#!/usr/bin/env fslpython
+#
+# Test film_gls with the --ar option on both even and odd
+# (#timepoints) data.
+
+
+import os
+import sys
+import tempfile
+import os.path as op
+import subprocess as sp
+
+import numpy as np
+
+from pyfeeds.common import fslbin
+
+outdir   = sys.argv[1]
+
+fslroi   = fslbin('fslroi')
+film_gls = fslbin('film_gls')
+v2t      = fslbin('Vest2Text')
+t2v      = fslbin('Text2Vest')
+
+dcon = 'design.con'
+dfts = 'design.fts'
+
+
+def run(cmd):
+    sp.run(cmd.split(), check=True)
+
+
+# one run with even # timepoints,
+# and one with odd # timepoints
+for suffix, npts in [('even', 50), ('odd', 51)]:
+
+    datafile  = f'roi_{suffix}.nii.gz'
+    dmat      = f'design_{suffix}.mat'
+    runoutdir = op.join(outdir, suffix)
+
+    run(f'{film_gls} --in={datafile} --ar --rn={runoutdir} '
+        f'--pd={dmat} --con={dcon} --fcon={dfts} '
+         '--thr=1000.0 --sa --ms=5')