-
Notifications
You must be signed in to change notification settings - Fork 0
/
lppParameters.m
81 lines (65 loc) · 1.69 KB
/
lppParameters.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
function EEG = lppParameters(EEG,timeOutputs,eegChan)
%% Calculate LPP parameters, like amplitude and magnitude.
% By John LaRocco
neuCase=EEG.neuEpochs;
alcCase=EEG.alcEpochs;
fodCase=EEG.fodEpochs;
fixCase=EEG.fixEpochs;
baselineWindow=timeOutputs.baselineWindow;
sRate=EEG.srate;
lBnd=.5;
uBnd=.75;
lBnd=baselineWindow+lBnd;
uBnd=baselineWindow+uBnd;
lBnd=round(lBnd*sRate);
uBnd=round(uBnd*sRate);
neuCase=squeeze(neuCase(eegChan,:,:));
alcCase=squeeze(alcCase(eegChan,:,:));
fodCase=squeeze(fodCase(eegChan,:,:));
fixCase=squeeze(fixCase(eegChan,:,:));
neuCase=mean(neuCase');
alcCase=mean(alcCase');
fodCase=mean(fodCase');
fixCase=mean(fixCase');
badCase=.5*(alcCase+fodCase);
badLpp=badCase(lBnd:uBnd);
neuLpp=neuCase(lBnd:uBnd);
alcLpp=alcCase(lBnd:uBnd);
fodLpp=fodCase(lBnd:uBnd);
fixLpp=fixCase(lBnd:uBnd);
ampBad=max(abs(badLpp));
magBad=sqrt(ampBad.*ampBad);
ampBadStd=std(badLpp);
ampNeu=max(abs(neuLpp));
magNeu=sqrt(ampNeu.*ampNeu);
ampNeuStd=max(abs(neuLpp));
ampAlc=max(abs(alcLpp));
magAlc=sqrt(ampAlc.*ampAlc);
ampAlcStd=std(alcLpp);
ampFod=max(abs(fodLpp));
magFod=sqrt(ampFod.*ampFod);
ampFodStd=std(fodLpp);
ampFix=max(abs(fixLpp));
magFix=sqrt(ampFix.*ampFix);
ampFixStd=std(fixLpp);
EEG.ampBad=ampBad;
EEG.ampNeu=ampNeu;
EEG.ampAlc=ampAlc;
EEG.ampFod=ampFod;
EEG.ampFix=ampFix;
EEG.magBad=magBad;
EEG.magNeu=magNeu;
EEG.magAlc=magAlc;
EEG.magFod=magFod;
EEG.magFix=magFix;
EEG.ampBadStd=ampBadStd;
EEG.ampNeuStd=ampNeuStd;
EEG.ampAlcStd=ampAlcStd;
EEG.ampFodStd=ampFodStd;
EEG.ampFixStd=ampFixStd;
EEG.neuCase=neuCase;
EEG.alcCase=alcCase;
EEG.fodCase=fodCase;
EEG.fixCase=fixCase;
EEG.badCase=badCase;
end