small change to eeg_icartifact_eyes and proc_subcomp to accommodate n…

…on-square icawinv - i.e., when fewer comps than channels exist
sjburwell · Jun 5, 2019 · 0eabadc · 0eabadc
1 parent dcda725
commit 0eabadc
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/functions/basefunc/eeg_icartifact_eyes.m b/functions/basefunc/eeg_icartifact_eyes.m
@@ -101,7 +101,7 @@
 eyetime   = abs(atanh([blink_corrs(:,2); hem_corrs(:,2)])); 
 [timemeas ,timecrit ] = thresholder(stat, eyetime , crit);
 
-all_ics = [1:length(EEG.icawinv), 1:length(EEG.icawinv)];
+all_ics = [1:size(EEG.icawinv,2), 1:size(EEG.icawinv,2)];
 if ~isempty(intersect(find(spacemeas>spacecrit),find(timemeas>timecrit)));
     disp([EEG.filename '; threshold(s): [' num2str([THRESHOLD.crit]) ']'])
     ICs = all_ics(intersect(find(spacemeas>spacecrit),find(timemeas>timecrit)));

diff --git a/functions/procfunc/proc_subcomp.m b/functions/procfunc/proc_subcomp.m
@@ -127,7 +127,7 @@
 [spacemeas,spacecrit] = thresholder(args.threshold(1).stat, space, args.threshold(1).crit);
 [timemeas ,timecrit ] = thresholder(args.threshold(1).stat, time , args.threshold(1).crit);
 
-all_ics = repmat([1:length(EEG.icawinv)],[1,length(args.spatial)]);
+all_ics = repmat([1:size(EEG.icawinv,2)],[1,length(args.spatial)]);
 if ~isempty(intersect(find(spacemeas>spacecrit),find(timemeas>timecrit)));
     disp([EEG.filename '; threshold(s): [' num2str([THRESHOLD.crit]) ']'])
     ICs = all_ics(intersect(find(spacemeas>spacecrit),find(timemeas>timecrit)));