add hack for MIGP to work on tICA-cleaned data that has removed dilat…

…ed areas
Washington-University · May 22, 2023 · d73c7d9 · d73c7d9
1 parent 0a8416e
commit d73c7d9
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/tICA/scripts/Compiled_MIGP/MIGP b/tICA/scripts/Compiled_MIGP/MIGP
diff --git a/tICA/scripts/MIGP.m b/tICA/scripts/MIGP.m
@@ -55,7 +55,9 @@ function MIGP(StudyFolder, Subjlist, fMRINamesRaw, ProcSTRING, dPCAinternal, dPC
  vnsum = vn.cdata * 0;
  end
  dtseries = ciftiopen(dtseriesname, wbcommand);
- grot = [grot demean(dtseries.cdata, 2) ./ repmat(max(vn.cdata, 0.001), 1, size(dtseries.cdata, 2))];
+ %FIXME: tICA cleanup appears to regress out dilated locations (from FoV issues), for now repurpose this divide-by-zero prevention to act like a stdev threshold
+ %presumably the input data will always be scaled to grand mean 10,000
+ grot = [grot demean(dtseries.cdata, 2) ./ repmat(max(vn.cdata, 10), 1, size(dtseries.cdata, 2))];
  vnsum = vnsum + vn.cdata;
  c = c + 1;
  else