-
Notifications
You must be signed in to change notification settings - Fork 0
/
tryRestLoader.m
85 lines (60 loc) · 1.74 KB
/
tryRestLoader.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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
%T = readtable('myfile.csv');
%T = readtable('myfile.csv','NumHeaderLines',3);
mainWindow=30;
esuf1='a_rest.bdf';
esuf2='b_rest.bdf';
suf1='a-Cue Reactivity EEG.log';
suf2='b-Cue Reactivity EEG.log';
hrValues=[];
bpmValues=[];
%https://www.mathworks.com/matlabcentral/fileexchange/73049-calculate-heart-rate-from-electrocardiogram-data
%https://github.com/VisLab/EEG-Beats
homeDir='C:\Users\John\Documents\MATLAB\tryEeg\CUE_REST\';
channelLocationFile = 'C:\Users\John\Documents\MATLAB\eeglab2021.1\plugins\dipfit\standard_BESA\standard-10-5-cap385.elp';
oldSR=1024;
list= dir([homeDir '\TRY*']);
T = struct2table(list);
subNames=table2cell(T(:,1));
ki=2;
subNum='TRY001';
for ki=1:length(subNames)
subNum=subNames{ki};
%% derive values
f11=strcat(homeDir,subNum);
f11=strcat(f11,'\');
f11=strcat(f11,subNum);
fName1=strcat(f11,suf1);
fName2=strcat(f11,suf2);
efName1=strcat(f11,esuf1);
efName2=strcat(f11,esuf2);
%% main loop
try
%% loading tables
timeOutputs1 = getTryTimestampsMeasures(fName1);
timeOutputs2 = getTryTimestampsMeasures(fName2);
%% preproc
EEG = pop_biosig(efName1);
[EEG,params] = tryPreproc(EEG);
[ekgPeaks, params] = eeg_beats(EEG, params);
[bpm1,hr1] = calcBpm(EEG,ekgPeaks);
EEG = pop_biosig(efName2);
[EEG,params] = tryPreproc(EEG);
[ekgPeaks, params] = eeg_beats(EEG, params);
[bpm2,hr2] = calcBpm(EEG,ekgPeaks);
hr=[hr1 hr2];
hrValues(ki,:)=hr;
bpm=[bpm1 bpm2];
bpmValues(ki,:)=bpm;
catch
hr=[0 0];
hrValues(ki,:)=hr;
bpm=[0 0];
bpmValues(ki,:)=bpm;
end
end
tryBatch={};
tryBatch{1}=hrValues;
tryBatch{2}=bpmValues;
save('tryHrRest.mat',"tryBatch");
%[h1,p1]=ttest2(bpmValues(:,1),bpmValues(:,2));
%[h2,p2]=ttest2(hrValues(:,1),hrValues(:,2));