rapidtide 3.0a12__py3-none-any.whl → 3.0a14__py3-none-any.whl

cloud/gmscalc-HCPYA

@@ -15,7 +15,7 @@ fi
 
 echo "Subject="${SUBJ}", RUN="${RUN}", PEDIR="${PEDIR}
 
-mkdir -p /data_out/HCP_1200/minpreprocglm/${SUBJ}
+mkdir -p /data_out/HCP_1200/minpreprocslfofilt/${SUBJ}
 gmscalc \
     /data_out/HCP_1200/${SUBJ}/rapidtide_rfMRI_${RUN}_${PEDIR}_desc-lfofilterCleaned_bold.nii.gz
     /data_out/HCP_1200/fixglm/${SUBJ}/rapidtide_rfMRI_${RUN}_${PEDIR}_lfofilteredCleanedGMS \

cloud/rapidtide-HCPYA

@@ -15,11 +15,11 @@ fi
 
 echo "Subject="${SUBJ}", RUN="${RUN}", PEDIR="${PEDIR}
 
-#    --glmsourcefile=/data_in/HCP_1200/${SUBJ}/MNINonLinear/Results/rfMRI_${RUN}_${PEDIR}/rfMRI_${RUN}_${PEDIR}_hp2000_clean.nii.gz \
-mkdir -p /data_out/HCP_1200/minpreprocglm/${SUBJ}
+#    --denoisesourcefile=/data_in/HCP_1200/${SUBJ}/MNINonLinear/Results/rfMRI_${RUN}_${PEDIR}/rfMRI_${RUN}_${PEDIR}_hp2000_clean.nii.gz \
+mkdir -p /data_out/HCP_1200/minpreprocslfofilt/${SUBJ}
 rapidtide \
     /data_in/HCP_1200/${SUBJ}/MNINonLinear/Results/rfMRI_${RUN}_${PEDIR}/rfMRI_${RUN}_${PEDIR}.nii.gz \
-    /data_out/HCP_1200/minpreprocglm/${SUBJ}/rapidtide_rfMRI_${RUN}_${PEDIR} \
+    /data_out/HCP_1200/minpreprocslfofilt/${SUBJ}/rapidtide_rfMRI_${RUN}_${PEDIR} \
     --despecklepasses 4 \
     --filterfreqs 0.009 0.15 \
     --searchrange -7.5 15.0 \

rapidtide/Colortables.py

@@ -66,7 +66,7 @@ def gen_thermal_state():
     thegradient = Gradients["thermal"]
     thegradient["name"] = "thermal"
     return thegradient
-    #return Gradients["thermal"]
+    # return Gradients["thermal"]
 
 
 def gen_flame_state():
@@ -75,7 +75,7 @@ def gen_flame_state():
     thegradient = Gradients["flame"]
     thegradient["name"] = "flame"
     return thegradient
-    #return Gradients["flame"]
+    # return Gradients["flame"]
 
 
 def gen_yellowy_state():
@@ -84,7 +84,7 @@ def gen_yellowy_state():
     thegradient = Gradients["yellowy"]
     thegradient["name"] = "yellowy"
     return thegradient
-    #return Gradients["yellowy"]
+    # return Gradients["yellowy"]
 
 
 def gen_bipolar_state():
@@ -93,7 +93,7 @@ def gen_bipolar_state():
     thegradient = Gradients["bipolar"]
     thegradient["name"] = "bipolar"
     return thegradient
-    #return Gradients["bipolar"]
+    # return Gradients["bipolar"]
 
 
 def gen_spectrum_state():
@@ -102,7 +102,7 @@ def gen_spectrum_state():
     thegradient = Gradients["spectrum"]
     thegradient["name"] = "spectrum"
     return thegradient
-    #return Gradients["spectrum"]
+    # return Gradients["spectrum"]
 
 
 def gen_turbo_state():
@@ -111,7 +111,7 @@ def gen_turbo_state():
     thegradient = Gradients["turbo"]
     thegradient["name"] = "turbo"
     return thegradient
-    #return Gradients["turbo"]
+    # return Gradients["turbo"]
 
 
 def gen_gray_state():
@@ -142,7 +142,7 @@ def gen_viridis_state():
     thegradient = Gradients["viridis"]
     thegradient["name"] = "viridis"
     return thegradient
-    #return Gradients["viridis"]
+    # return Gradients["viridis"]
 
 
 def gen_inferno_state():
@@ -151,7 +151,7 @@ def gen_inferno_state():
     thegradient = Gradients["inferno"]
     thegradient["name"] = "inferno"
     return thegradient
-    #return Gradients["inferno"]
+    # return Gradients["inferno"]
 
 
 def gen_plasma_state():
@@ -160,7 +160,7 @@ def gen_plasma_state():
     thegradient = Gradients["plasma"]
     thegradient["name"] = "plasma"
     return thegradient
-    #return Gradients["plasma"]
+    # return Gradients["plasma"]
 
 
 def gen_magma_state():
@@ -169,7 +169,7 @@ def gen_magma_state():
     thegradient = Gradients["magma"]
     thegradient["name"] = "magma"
     return thegradient
-    #return Gradients["magma"]
+    # return Gradients["magma"]
 
 
 def gen_g2y2r_state():

rapidtide/DerivativeDelay.py

@@ -0,0 +1,213 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+#
+#   Copyright 2016-2024 Blaise Frederick
+#
+#   Licensed under the Apache License, Version 2.0 (the "License");
+#   you may not use this file except in compliance with the License.
+#   You may obtain a copy of the License at
+#
+#       http://www.apache.org/licenses/LICENSE-2.0
+#
+#   Unless required by applicable law or agreed to in writing, software
+#   distributed under the License is distributed on an "AS IS" BASIS,
+#   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+#   See the License for the specific language governing permissions and
+#   limitations under the License.
+#
+#
+"""
+A class to impmement regressor refinement
+"""
+import copy
+
+import numpy as np
+
+import rapidtide.fit as tide_fit
+import rapidtide.io as tide_io
+import rapidtide.miscmath as tide_math
+import rapidtide.refinedelay as tide_refinedelay
+import rapidtide.resample as tide_resample
+import rapidtide.stats as tide_stats
+import rapidtide.util as tide_util
+from rapidtide.tests.utils import mse
+
+
+class DerivativeDelay:
+    def __init__(
+        self,
+        internalvalidfmrishape,
+        internalvalidpaddedfmrishape,
+        pid,
+        outputname,
+        initial_fmri_x,
+        paddedinitial_fmri_x,
+        os_fmri_x,
+        genlagtc,
+        sharedmem=False,
+        gausssigma=-1,
+        numderivs=1,
+        lagminthresh=0.25,
+        lagmaxthresh=3.0,
+        sigmathresh=1000.0,
+        cleanrefined=False,
+        bipolar=False,
+        fixdelay=False,
+        includemask=None,
+        excludemask=None,
+        LGR=None,
+        nprocs=1,
+        detrendorder=1,
+        alwaysmultiproc=False,
+        showprogressbar=True,
+        chunksize=50000,
+        padtrs=10,
+        refineprenorm="var",
+        refineweighting=None,
+        refinetype="pca",
+        pcacomponents=0.8,
+        dodispersioncalc=False,
+        dispersioncalc_lower=-5.0,
+        dispersioncalc_upper=5.0,
+        dispersioncalc_step=0.5,
+        windowfunc="hamming",
+        passes=3,
+        maxpasses=15,
+        convergencethresh=None,
+        interptype="univariate",
+        usetmask=False,
+        tmask_y=None,
+        tmaskos_y=None,
+        fastresamplerpadtime=45.0,
+        debug=False,
+        rt_floattype="float64",
+        rt_floatset=np.float64,
+    ):
+        self.outputname = outputname
+        self.numderivs = numderivs
+        if gausssigma < 0.0:
+            # set gausssigma automatically
+            self.gausssigma = np.mean([xdim, ydim, slicedim]) / 2.0
+        else:
+            self.gausssigma = gausssigma
+        self.debug = debug
+        self.setgenlagtc(genlagtc)
+
+    def setgenlagtc(self, genlagtc):
+        self.genlagtc = genlagtc
+
+    def getderivratios(
+        self, fmri_data_valid, validvoxels, initial_fmri_x, lagtimes_valid, corrmask_valid
+    ):
+        print("Refinement calibration start")
+        regressderivratios = tide_refinedelay.getderivratios(
+            fmri_data_valid,
+            validvoxels,
+            initial_fmri_x,
+            lagtimes_valid,
+            corrmask_valid,
+            self.genlagtc,
+            "glm",
+            self.outputname,
+            oversamptr,
+            sLFOfitmean,
+            rvalue,
+            r2value,
+            fitNorm[:, : (self.numderivs + 1)],
+            fitcoeff[:, : (self.numderivs + 1)],
+            movingsignal,
+            lagtc,
+            filtereddata,
+            LGR,
+            TimingLGR,
+            therunoptions,
+            regressderivs=self.numderivs,
+            debug=self.debug,
+        )
+
+    def calibrate(self):
+        if self.numderivs == 1:
+            medfiltregressderivratios, filteredregressderivratios, delayoffsetMAD = (
+                tide_refinedelay.filterderivratios(
+                    regressderivratios,
+                    (xsize, ysize, numslices),
+                    validvoxels,
+                    (xdim, ydim, slicedim),
+                    gausssigma=args.delayoffsetgausssigma,
+                    patchthresh=args.delaypatchthresh,
+                    fileiscifti=False,
+                    textio=False,
+                    rt_floattype=rt_floattype,
+                    debug=args.debug,
+                )
+            )
+
+            # find the mapping of derivative ratios to delays
+            tide_refinedelay.trainratiotooffset(
+                self.genlagtc,
+                initial_fmri_x,
+                self.outputname,
+                args.outputlevel,
+                mindelay=args.mindelay,
+                maxdelay=args.maxdelay,
+                numpoints=args.numpoints,
+                debug=args.debug,
+            )
+            TimingLGR.info("Refinement calibration end")
+
+            # now calculate the delay offsets
+            TimingLGR.info("Calculating delay offsets")
+            delayoffset = np.zeros_like(filteredregressderivratios)
+            if args.focaldebug:
+                print(f"calculating delayoffsets for {filteredregressderivratios.shape[0]} voxels")
+            for i in range(filteredregressderivratios.shape[0]):
+                delayoffset[i], closestoffset = tide_refinedelay.ratiotodelay(
+                    filteredregressderivratios[i]
+                )
+            refinedvoxelstoreport = filteredregressderivratios.shape[0]
+        else:
+            medfiltregressderivratios = np.zeros_like(regressderivratios)
+            filteredregressderivratios = np.zeros_like(regressderivratios)
+            delayoffsetMAD = np.zeros(args.refineregressderivs, dtype=float)
+            for i in range(args.refineregressderivs):
+                (
+                    medfiltregressderivratios[i, :],
+                    filteredregressderivratios[i, :],
+                    delayoffsetMAD[i],
+                ) = tide_refinedelay.filterderivratios(
+                    regressderivratios[i, :],
+                    (xsize, ysize, numslices),
+                    validvoxels,
+                    (xdim, ydim, slicedim),
+                    gausssigma=args.delayoffsetgausssigma,
+                    patchthresh=args.delaypatchthresh,
+                    fileiscifti=False,
+                    textio=False,
+                    rt_floattype=rt_floattype,
+                    debug=args.debug,
+                )
+
+    def getdelays(self):
+        # now calculate the delay offsets
+        delayoffset = np.zeros_like(filteredregressderivratios[0, :])
+        if self.debug:
+            print(f"calculating delayoffsets for {filteredregressderivratios.shape[1]} voxels")
+        for i in range(filteredregressderivratios.shape[1]):
+            delayoffset[i] = tide_refinedelay.coffstodelay(
+                filteredregressderivratios[:, i],
+                mindelay=self.mindelay,
+                maxdelay=self.maxdelay,
+            )
+        refinedvoxelstoreport = filteredregressderivratios.shape[1]
+
+    def savestats(self):
+        namesuffix = "_desc-delayoffset_hist"
+        tide_stats.makeandsavehistogram(
+            delayoffset,
+            therunoptions["histlen"],
+            1,
+            self.outputname + namesuffix,
+            displaytitle="Histogram of delay offsets calculated from coefficient ratios",
+            dictvarname="delayoffsethist",
+            thedict=None,
+        )

rapidtide/{Refiner.py → RegressorRefiner.py}

@@ -33,7 +33,7 @@ import rapidtide.util as tide_util
 from rapidtide.tests.utils import mse
 
 
-class Refiner:
+class RegressorRefiner:
 
     refinemaskvoxels = None
 

rapidtide/__init__.py

@@ -1,3 +1,4 @@
 # -*- coding: utf-8 -*-
 from . import _version
-__version__ = _version.get_versions()['version']
+
+__version__ = _version.get_versions()["version"]

rapidtide/_version.py

@@ -12,12 +12,12 @@
 """Git implementation of _version.py."""
 
 import errno
+import functools
 import os
 import re
 import subprocess
 import sys
 from typing import Any, Callable, Dict, List, Optional, Tuple
-import functools
 
 
 def get_keywords() -> Dict[str, str]:

rapidtide/data/examples/src/test_mlregressallt.py

@@ -3,10 +3,10 @@
 import numpy as np
 from sklearn.linear_model import LinearRegression
 
-from rapidtide.fit import mlregress_alt
+from rapidtide.fit import mlregress, olsregress
 
 
-def test_mlregress_alt():
+def test_olsregress():
     # Test with a simple dataset
     X = np.array([[1, 2], [3, 4]])
     y = np.array([5, 6])
@@ -15,13 +15,20 @@ def test_mlregress_alt():
     expected_intercept = 4
     expected_R = 0.707
 
-    coeffs, R = mlregress_alt(X, y, intercept)
+    mlcoeffs, mlR = mlregress(X, y, intercept=intercept)
+    print(f"{mlcoeffs=}, {mlR=}")
 
-    print(coeffs, expected_coeffs)
-    assert np.allclose(coeffs, expected_coeffs)
-    assert np.isclose(R, expected_R)
+    olscoeffs, olsR = olsregress(X, y, intercept=intercept)
+    print(f"{olscoeffs=}, {olsR=}")
 
-def test_mlregress_alt_intercept():
+    print(mlcoeffs, expected_coeffs)
+    print(olscoeffs, expected_coeffs)
+    assert np.allclose(mlcoeffs, expected_coeffs)
+    assert np.isclose(mlR, expected_R)
+    assert np.allclose(olscoeffs, expected_coeffs)
+    assert np.isclose(olsR, expected_R)
+
+def test_olsregress_intercept():
     # Test with a simple dataset and no intercept
     X = np.array([[1, 2], [3, 4]])
     y = np.array([5, 6])
@@ -30,26 +37,34 @@ def test_mlregress_alt_intercept():
     expected_intercept = None
     expected_R = 0.707
 
-    coeffs, R = mlregress_alt(X, y, intercept)
+    coeffs, R = mlregress(X, y, intercept)
 
     assert np.allclose(coeffs, expected_coeffs)
     assert np.isclose(R, expected_R)
 
-def test_mlregress_alt_no_intercept():
+def test_olsregress_nointercept():
     # Test with a simple dataset and no intercept
     X = np.array([[1, 2], [3, 4]])
     y = np.array([5, 6])
-    intercept = None
+    intercept = False
     expected_coeffs = np.array([[1, 2], [3, 4]])
-    expected_intercept = None
+    expected_intercept = False
     expected_R = 0.707
 
-    coeffs, R = mlregress_alt(X, y, intercept)
+    mlcoeffs, mlR = mlregress(X, y, intercept=intercept)
+    print(f"{mlcoeffs=}, {mlR=}")
 
-    assert np.allclose(coeffs, expected_coeffs)
-    assert np.isclose(R, expected_R)
+    olscoeffs, olsR = olsregress(X, y, intercept=intercept)
+    print(f"{olscoeffs=}, {olsR=}")
+
+    print(mlcoeffs, expected_coeffs)
+    print(olscoeffs, expected_coeffs)
+    assert np.allclose(mlcoeffs, expected_coeffs)
+    assert np.isclose(mlR, expected_R)
+    assert np.allclose(olscoeffs, expected_coeffs)
+    assert np.isclose(olsR, expected_R)
 
-test_mlregress_alt()
-test_mlregress_alt_intercept()
-test_mlregress_alt_nointercept()
+test_olsregress_nointercept()
+test_olsregress_intercept()
+test_olsregress()
 

rapidtide/data/examples/src/testalign

@@ -41,7 +41,7 @@ alignedvec, maxdelay, maxval, failreason = tide_corr.aligntcwithref(
 )
 print(f"{maxdelay=}, {maxval=}, {failreason=}")
 
-filteredvec, datatoremove, R, coffs = tide_fit.glmfilt(fixedvec, alignedvec, debug=True)
+filteredvec, datatoremove, R, coffs = tide_fit.linfitfilt(fixedvec, alignedvec, debug=True)
 
 offset = 0.0
 plt.plot(xaxis, fixedvec + offset)

rapidtide/data/examples/src/testboth

@@ -10,7 +10,7 @@ rapidtide \
     sub-RAPIDTIDETEST.nii.gz \
     ../dst/sub-RAPIDTIDETEST
 
-retroglm \
+retroregress \
     --nprocs -1 \
     --alternateoutput ../dst/sub-RETROTEST_refined \
     --makepseudofile \

rapidtide/data/examples/src/testcifti

@@ -0,0 +1,11 @@
+#!/bin/bash
+
+rapidtide \
+	sub-RAPIDTIDETEST_cifti.ptseries.nii \
+        ../dst/sub-RAPIDTIDETEST5 \
+        --filterband lfo \
+        --searchrange -15 15 \
+        --passes 3 \
+        --globalmaskmethod variance \
+        --nprocs -1 \
+        --numnull 10000 

rapidtide/data/examples/src/testdelayvar

@@ -0,0 +1,14 @@
+#!/bin/bash
+
+delayvar \
+    sub-RAPIDTIDETEST.nii.gz \
+    ../dst/sub-RAPIDTIDETEST \
+    --alternateoutput ../dst/delayvartest \
+    --nprocs -1 \
+    --windowsize 30 \
+    --systemicfittype pca \
+    --trainrange -999 -999 \
+    --pcacomponents 0.8 \
+    --outputlevel max \
+    --focaldebug \
+    --windelayoffsetspatialfilt 3

rapidtide/data/examples/src/testfmri

@@ -2,6 +2,7 @@
 
 rapidtide \
     --spatialfilt -1 \
+    --delayoffsetspatialfilt -1 \
     --nprocs -1 \
     --searchrange -5 20 \
     --simcalcrange 50 -1 \

rapidtide/data/examples/src/testglmfilt

@@ -34,7 +34,9 @@ datavec = amp1 * evs[:, 0] + amp2 * evs[:, 1] + amp3 * evs[:, 2] + intercept
 noisevec = noiseamp * np.random.rand(arrlen)
 invec = datavec + noisevec
 
-filtered, datatoremove, R, outcoffs, outintercept = tide_fit.glmfilt(invec, evs, returnintercept=True, debug=True)
+filtered, datatoremove, R, outcoffs, outintercept = tide_fit.linfitfilt(
+    invec, evs, returnintercept=True, debug=True
+)
 
 incoffs = [amp1, amp2, amp3]
 coffdiffs = outcoffs - incoffs
@@ -58,7 +60,7 @@ plt.plot(datavec - datatoremove + offset)
 print(R)
 plt.show()
 
-filtered, datatoremove, R, coffs = tide_fit.glmfilt(invec, evs[:, 0], debug=True)
+filtered, datatoremove, R, coffs = tide_fit.linfitfilt(invec, evs[:, 0], debug=True)
 
 """
 offset = 0.0
@@ -80,7 +82,7 @@ invec2 = (
     + noiseamp * np.random.rand(arrlen)
 )
 
-filtered, datatoremove, R, coffs = tide_fit.glmfilt(invec2, evs, debug=True)
+filtered, datatoremove, R, coffs = tide_fit.linfitfilt(invec2, evs, debug=True)
 
 """
 offset = 0.0
@@ -95,7 +97,7 @@ print(R)
 plt.show()
 """
 
-filtered, thenewevs, datatoremove, R, coffs = tide_fit.expandedglmfilt(
+filtered, thenewevs, datatoremove, R, coffs = tide_fit.expandedlinfitfilt(
     invec2, evs, ncomps=3, debug=True
 )
 
@@ -116,7 +118,7 @@ print(R)
 plt.show()
 """
 
-filtered, datatoremove, R, coffs = tide_fit.glmfilt(invec2, thenewevs, debug=True)
+filtered, datatoremove, R, coffs = tide_fit.linfitfilt(invec2, thenewevs, debug=True)
 """
 offset = 0.0
 plt.plot(invec2 + offset)
@@ -138,7 +140,7 @@ alignedev, maxdelay, maxval, failreason = tide_corr.aligntcwithref(
     verbose=True,
 )
 print(f"{maxdelay=}, {maxval=}, {failreason=}")
-filtered, datatoremove, R, coffs = tide_fit.glmfilt(invec3, alignedev, debug=True)
+filtered, datatoremove, R, coffs = tide_fit.linfitfilt(invec3, alignedev, debug=True)
 offset = 0.0
 plt.plot(evs[:, 0])
 offset += 1.0

rapidtide/data/examples/src/testhappy

@@ -65,5 +65,5 @@
 happy \
     sub-HAPPYTEST.nii.gz \
     sub-HAPPYTEST.json \
-    --spatialglm \
+    --spatialregression \
     ../dst/happy

rapidtide/data/examples/src/testnewrefine

@@ -6,7 +6,7 @@ rapidtide \
     --searchrange -5 20 \
     --simcalcrange 50 -1 \
     --outputlevel max \
-    --refineglmderivs 1 \
+    --refineregressderivs 1 \
     sub-RAPIDTIDETEST.nii.gz \
     ../dst/sub-RAPIDTIDETEST_1deriv
 
@@ -16,34 +16,34 @@ rapidtide \
     --searchrange -5 20 \
     --simcalcrange 50 -1 \
     --outputlevel max \
-    --refineglmderivs 2 \
+    --refineregressderivs 2 \
     sub-RAPIDTIDETEST.nii.gz \
     ../dst/sub-RAPIDTIDETEST_2deriv
 
 
-retroglm \
+retroregress \
     sub-RAPIDTIDETEST.nii.gz \
     ../dst/sub-RAPIDTIDETEST_1deriv \
     --alternateoutput ../dst/retrotest_1to1 \
     --nprocs -1 \
-    --glmderivs 0 \
-    --refineglmderivs 1 \
+    --regressderivs 0 \
+    --refineregressderivs 1 \
     --outputlevel max
 
-retroglm \
+retroregress \
     sub-RAPIDTIDETEST.nii.gz \
     ../dst/sub-RAPIDTIDETEST_2deriv \
     --alternateoutput ../dst/retrotest_2to1 \
     --nprocs -1 \
-    --glmderivs 0 \
-    --refineglmderivs 1 \
+    --regressderivs 0 \
+    --refineregressderivs 1 \
     --outputlevel max
 
-retroglm \
+retroregress \
     sub-RAPIDTIDETEST.nii.gz \
     ../dst/sub-RAPIDTIDETEST_1deriv \
     --alternateoutput ../dst/retrotest_1to2 \
     --nprocs -1 \
-    --glmderivs 0 \
-    --refineglmderivs 2 \
+    --regressderivs 0 \
+    --refineregressderivs 2 \
     --outputlevel max

rapidtide/data/examples/src/testnoiseamp

@@ -1,9 +1,9 @@
 #!/usr/bin/env python
 
-import rapidtide.miscmath as tide_math
-import rapidtide.io as tide_io
 import matplotlib.pyplot as plt
 
+import rapidtide.io as tide_io
+import rapidtide.miscmath as tide_math
 
 windowsize = 40.0
 samplerate, starttime, inputdata = tide_io.readcolfrombidstsv("../dst/sub-RAPIDTIDETEST_desc-oversampledmovingregressor_timeseries.json", columnname="pass3")

rapidtide/data/examples/src/testretro

@@ -6,16 +6,25 @@
 #    --refinetype unweighted_average \
 #    --numnull 0 \
 #    --outputlevel max \
-#    --glmderivs 1 \
+#    --regressderivs 1 \
 #    sub-RAPIDTIDETEST.nii.gz \
 #    ../dst/sub-RAPIDTIDETEST
 
-##retroglm sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/0deriv --nprocs -1 
-retroglm sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/0deriv --nprocs -1 --glmderivs 0 --makepseudofile --outputlevel max
-retroglm sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/0deriv_refined --nprocs -1 --glmderivs 0 --makepseudofile --outputlevel max --refinedelay --nofilterwithrefineddelay
-retroglm sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/0deriv_refinedfit --nprocs -1 --glmderivs 0 --makepseudofile --outputlevel max --refinedelay
-##retroglm sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/2deriv --nprocs -1 --glmderivs 2
-##retroglm sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/3deriv --nprocs -1 --glmderivs 3
+##retroregress sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/0deriv --nprocs -1 
+
+retroregress \
+    sub-RAPIDTIDETEST.nii.gz \
+    ../dst/sub-RAPIDTIDETEST \
+    --alternateoutput ../dst/wintest \
+    --nprocs -1 \
+    --outputlevel max \
+    --focaldebug \
+    --delayoffsetspatialfilt -1 
+
+#retroregress sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/0deriv_refined --nprocs -1 --regressderivs 0 --makepseudofile --outputlevel max --refinedelay --nofilterwithrefineddelay
+#retroregress sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/0deriv_refinedfit --nprocs -1 --regressderivs 0 --makepseudofile --outputlevel max --refinedelay
+##retroregress sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/2deriv --nprocs -1 --regressderivs 2
+##retroregress sub-RAPIDTIDETEST.nii.gz ../dst/sub-RAPIDTIDETEST --alternateoutput ../dst/3deriv --nprocs -1 --regressderivs 3
 
 ALTNAME=1deriv
 

rapidtide/data/examples/src/testretrolagtcs

@@ -13,4 +13,4 @@ retrolagtcs \
     ../dst/sub-RAPIDTIDETEST_1deriv_desc-maxtime_map.nii.gz \
     ../dst/sub-RAPIDTIDETEST_1deriv_desc-lagtcgenerator_timeseries \
     ../dst/retrolag_1deriv \
-    --glmderivs 1 
+    --regressderivs 1

rapidtide/dlfilter.py

@@ -77,7 +77,6 @@ from tensorflow.keras.layers import (
     TimeDistributed,
     UpSampling1D,
 )
-
 from tensorflow.keras.models import Sequential, load_model
 from tensorflow.keras.optimizers import RMSprop