-
Notifications
You must be signed in to change notification settings - Fork 2
/
erf_osc_analysis_performance.m
executable file
·32 lines (26 loc) · 1.46 KB
/
erf_osc_analysis_performance.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
function erf_osc_analysis_perfomance(subj)
load(sprintf('/project/3011085.02/analysis/behavior/sub-%03d/sub-%03d_trialinfo.mat', subj, subj));
trialinfo_all = trialinfo;
nTrials_all = size(trialinfo,1);
idxNoReversal = find(trialinfo(:,8)==0);
noReversal = length(idxNoReversal);
trialinfo(idxNoReversal,:)=[];
nTrials_reversal = size(trialinfo,1);
nTrials_validResp = sum(trialinfo(:,9)>0 & ((trialinfo(:,9)-trialinfo(:,8))/1200)<0.7);
performanceAll = nTrials_validResp/nTrials_reversal;
clear trialinfo
load(sprintf('/project/3011085.02/processed/sub-%03d/ses-meg01/sub-%03d_cleandata.mat', subj, subj), 'dataClean');
trialinfo_cleaned = dataClean.trialinfo;
trialinfo = trialinfo_cleaned;
nTrials_cleaned = size(trialinfo,1);
idxNoReversal_cleaned = find(trialinfo(:,5)==0);
noReversal_cleaned = length(idxNoReversal_cleaned);
trialinfo(idxNoReversal_cleaned,:)=[];
nTrials_reversal_cleaned = size(trialinfo,1);
nTrials_validResp_cleaned = sum(trialinfo(:,6)>0 & ((trialinfo(:,6)-trialinfo(:,5))/1200)<0.7);
performance_cleaned = nTrials_validResp_cleaned/nTrials_reversal_cleaned;
save(sprintf('/project/3011085.02/analysis/behavior/sub-%03d/sub-%03d_performance.mat', subj, subj), ...
'performanceAll', 'nTrials_validResp', 'nTrials_all', 'nTrials_reversal', 'noReversal', 'trialinfo_all', ...
'trialinfo_cleaned', 'nTrials_cleaned', 'idxNoReversal_cleaned', 'noReversal_cleaned', 'nTrials_reversal_cleaned', ...
'nTrials_validResp_cleaned', 'performance_cleaned');
end