Hello!
I cannot figure out why this is breaking. It should not be out of range as we are just comparing values. Could you please help me get to the bottom of this?
Function Call in main script:
%% Regional Channel Break down:
% each channel
% each sub-region
% included in the function call
R1 = [1,2,3,17,5,18,6,11]; % right ventrolateral prefrontal cortex
R2 = [4, 9, 10]; % right ventrolateral prefrontal cortex
R3 = [7, 8, 21,22, 12, 13,25, 26]; % right frontopolar prefrontal cortex
R4 = [14, 15, 16, 29, 30]; % right oribitofrontal cortex
L1 = [19,20,33,34,35, 38,39,43]; % left ventrolateral prefrontal cortex
L2 = [40, 44, 45]; % left ventrolateral prefrontal cortex
L3 = [23, 24, 27, 28, 36, 37, 41, 42]; % left frontopolar prefrontal cortex
L4 = [31, 32, 46, 47, 48]; % left oribitofrontal cortex
regions = input(['Which Region? Enter one of the following: ' ...
'R1, R2, R3, R4, L1, L2, L3, L4: '])
%Looping through a struct:
epochMatrixRegions = zeros(24,4);
fn=fieldnames(dataStruct); %datastruct contains the timestamps
fn=fn(2:end);
i = 1;
%loop through the fields
for i = 1:numel(fn)
%access the data
stamps=dataStruct.(fn{i});
oxyhb= Participants.(fn{i});
%call function
[epochMatrixRegions(i,1), epochMatrixRegions(i,2),...
epochMatrixRegions(i,3), epochMatrixRegions(i,4)] = ...
thinkAloudAUC_Regions(stamps,oxyhb, regions);
end
Part inside the fcn that is breaking:
%% Regions
R1 = [1,2,3,17,5,18,6,11]; % right ventrolateral prefrontal cortex
R2 = [4, 9, 10]; % right ventrolateral prefrontal cortex
R3 = [7, 8, 21,22, 12, 13,25, 26]; % right frontopolar prefrontal cortex
R4 = [14, 15, 16, 29, 30]; % right oribitofrontal cortex
L1 = [19,20,33,34,35, 38,39,43]; % left ventrolateral prefrontal cortex
L2 = [40, 44, 45]; % left ventrolateral prefrontal cortex
L3 = [23, 24, 27, 28, 36, 37, 41, 42]; % left frontopolar prefrontal cortex
L4 = [31, 32, 46, 47, 48]; % left oribitofrontal cortex
%% Load Participant Data:
% The odd columns are HbO, the even are HbR
% for now, we only want to look at oxygenated data
HbO = oxyhb.data(:,2:2:end);
%filter the HbO data for the desired region:
channels = 1:1:48;
for j = 1:length(channels)
if ismember(j,regions) == false
HbO(:,j) = [];
end
end
The error I get is:
Matrix index is out of range for deletion.
Error in thinkAloudAUC_Regions (line 45)
HbO(:,j) = [];
Error in ThinkAloudProcessing (line 265)
thinkAloudAUC_Regions(stamps,oxyhb, regions);
