rapidtide 3.0.7.1__py3-none-any.whl → 3.0.8__py3-none-any.whl

rapidtide/linfitfiltpass.py

@@ -92,6 +92,7 @@ def linfitfiltpass(
     nprocs=1,
     alwaysmultiproc=False,
     confoundregress=False,
+    ratiosonly=False,
     procbyvoxel=True,
     showprogressbar=True,
     mp_chunksize=1000,
@@ -207,6 +208,18 @@ def linfitfiltpass(
                     r2value[voxel[0]] = voxel[3]
                     filtereddata[voxel[0], :] = voxel[7]
                     itemstotal += 1
+            elif ratiosonly:
+                for voxel in data_out:
+                    meanvalue[voxel[0]] = voxel[1]
+                    rvalue[voxel[0]] = voxel[2]
+                    r2value[voxel[0]] = voxel[3]
+                    if theevs.ndim > 1:
+                        fitcoeff[voxel[0], :] = voxel[4]
+                        fitNorm[voxel[0], :] = voxel[5]
+                    else:
+                        fitcoeff[voxel[0]] = voxel[4]
+                        fitNorm[voxel[0]] = voxel[5]
+                    itemstotal += 1
             else:
                 for voxel in data_out:
                     meanvalue[voxel[0]] = voxel[1]
@@ -227,6 +240,18 @@ def linfitfiltpass(
                     r2value[timepoint[0]] = timepoint[3]
                     filtereddata[:, timepoint[0]] = timepoint[7]
                     itemstotal += 1
+            elif ratiosonly:
+                for timepoint in data_out:
+                    meanvalue[timepoint[0]] = timepoint[1]
+                    rvalue[timepoint[0]] = timepoint[2]
+                    r2value[timepoint[0]] = timepoint[3]
+                    if theevs.ndim > 1:
+                        fitcoeff[:, timepoint[0]] = timepoint[4]
+                        fitNorm[:, timepoint[0]] = timepoint[5]
+                    else:
+                        fitcoeff[timepoint[0]] = timepoint[4]
+                        fitNorm[timepoint[0]] = timepoint[5]
+                    itemstotal += 1
             else:
                 for timepoint in data_out:
                     meanvalue[timepoint[0]] = timepoint[1]
@@ -271,6 +296,23 @@ def linfitfiltpass(
                             rt_floatset=rt_floatset,
                             rt_floattype=rt_floattype,
                         )
+                    elif ratiosonly:
+                        (
+                            dummy,
+                            meanvalue[vox],
+                            rvalue[vox],
+                            r2value[vox],
+                            fitcoeff[vox],
+                            fitNorm[vox],
+                            dummy,
+                            dummy,
+                        ) = _procOneRegressionFitItem(
+                            vox,
+                            theevs[vox, :],
+                            thedata,
+                            rt_floatset=rt_floatset,
+                            rt_floattype=rt_floattype,
+                        )
                     else:
                         (
                             dummy,
@@ -315,6 +357,23 @@ def linfitfiltpass(
                             rt_floatset=rt_floatset,
                             rt_floattype=rt_floattype,
                         )
+                    elif ratiosonly:
+                        (
+                            dummy,
+                            meanvalue[timepoint],
+                            rvalue[timepoint],
+                            r2value[timepoint],
+                            fitcoeff[timepoint],
+                            fitNorm[timepoint],
+                            dummy,
+                            dummy,
+                        ) = _procOneRegressionFitItem(
+                            timepoint,
+                            theevs[:, timepoint],
+                            thedata,
+                            rt_floatset=rt_floatset,
+                            rt_floattype=rt_floattype,
+                        )
                     else:
                         (
                             dummy,

rapidtide/makelaggedtcs.py

@@ -71,7 +71,14 @@ def makelaggedtcs(
     chunksize=1000,
     rt_floatset=np.float64,
     rt_floattype="float64",
+    debug=False,
 ):
+    if debug:
+        print("makelaggedtcs: Starting")
+        print(f"\t{lagtc.shape=}")
+        print(f"\t{lagtimes.shape=}")
+        print(f"\t{timeaxis.shape=}")
+
     inputshape = lagtc.shape
     voxelargs = [
         lagtcgenerator,
@@ -110,4 +117,7 @@ def makelaggedtcs(
         if LGR is not None:
             LGR.info("garbage collected")
 
+    if debug:
+        print("makelaggedtcs: End\n\n")
+
     return volumetotal

rapidtide/refinedelay.py

@@ -282,18 +282,6 @@ def ratiotodelay(theratio, offset=0.0, debug=False):
             closestindex = offsetindex
     closestoffset = funcoffsets[closestindex]
     distance = np.fabs(funcoffsets[closestindex] - offset)
-    """if (maplimits[0] < theratio < maplimits[1]) and (
-        distance < (funcoffsets[1] - funcoffsets[0]) / 2
-    ):
-        return (
-            ratiotooffsetfunc[closestindex](theratio) + (offset - closestoffset),
-            closestoffset,
-        )
-    else:
-        return (
-            0.0,
-            closestoffset,
-        )"""
 
     if theratio < maplimits[0]:
         return (
@@ -302,7 +290,7 @@ def ratiotodelay(theratio, offset=0.0, debug=False):
         )
     elif theratio > maplimits[1]:
         return (
-            ratiotooffsetfunc[closestindex](maplimits[1]),
+            ratiotooffsetfunc[closestindex](maplimits[1]) - (offset - closestoffset),
             closestoffset,
         )
     else:
@@ -371,12 +359,12 @@ def getderivratios(
     if endtr is None:
         endtr = fmri_data_valid.shape[1]
     if debug:
-        print("getderivratios")
-        print(f"{fitNorm.shape=}")
-        print(f"{fitcoeff.shape=}")
-        print(f"{regressderivs=}")
-        print(f"{starttr=}")
-        print(f"{endtr=}")
+        print("getderivratios: Starting")
+        print(f"\t{fitNorm.shape=}")
+        print(f"\t{fitcoeff.shape=}")
+        print(f"\t{regressderivs=}")
+        print(f"\t{starttr=}")
+        print(f"\t{endtr=}")
 
     voxelsprocessed_regressionfilt, regressorset, evset = tide_regressfrommaps.regressfrommaps(
         fmri_data_valid[:, starttr:endtr],
@@ -406,6 +394,7 @@ def getderivratios(
         mp_chunksize=optiondict["mp_chunksize"],
         showprogressbar=optiondict["showprogressbar"],
         alwaysmultiproc=optiondict["alwaysmultiproc"],
+        ratiosonly=True,
         debug=debug,
     )
 
@@ -418,6 +407,8 @@ def getderivratios(
         for i in range(regressderivs):
             regressderivratios[i, :] = np.nan_to_num(fitcoeff[:, i + 1] / fitcoeff[:, 0])
 
+    if debug:
+        print("getderivratios: End\n\n")
     return regressderivratios, rvalue