echoPred/plottingResults.m at master · lhuntneuro/echoPred · GitHub

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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
function plottingResults(prediction,residual,eeg,stims,avec,avecR,param,videoFlag)

% videoFlag: if 1 it shows the video.
% all other inputs come from echoPred as results of the simulation; except for param and stims defined before echoPred.


supportdim = 2*param.refreshrate; %number of datapoints in filter

randTrial=ceil(param.trialnumber*rand(1));

times = [1:supportdim]*1000/param.refreshrate; %in ms
eegRandomTrial=squeeze(eeg(:,randTrial,:));

for ii=1:size(eegRandomTrial,2)
    eegRandomTrialnew(:,ii)=zscore(eegRandomTrial(:,ii));
end


if videoFlag==1
    figure

    for ii=1:10:size(eeg,1)

        subplot(2,2,1)
        cla
        hold on
        plot(eegRandomTrial)
        myLim=get(gca,'ylim');
        plot([ii ii],myLim)
        set(gca,'ylim',myLim)
        title('not zscored')
        xlabel('time [ms]')

        subplot(2,2,2)
        cla
        hold on
        plot(eegRandomTrialnew(:,1))
        myLim2=get(gca,'ylim');
        plot([ii ii],myLim2)
        set(gca,'ylim',myLim2)
        title('zscored')
        xlabel('time [ms]')

        subplot(2,2,3)
        imagesc(eegRandomTrialnew(ii,:))
        set(gca,'clim',myLim2)
        colorbar
    	xlabel('layers')
        title('2D representation')

        subplot(2,2,4)
        plot(eegRandomTrialnew(ii,:))
        set(gca,'ylim',myLim2)
    	xlabel('layers')
        title('1D representation')
        pause(.00001)
    end
end


%% additional figures

figure
subplot(2,1,1)
cla
hold on
plot(eegRandomTrial)
myLim=get(gca,'ylim');
plot([ii ii],myLim)
set(gca,'ylim',myLim)
title('eegAv')
for ii=1:size(avec,2)
    legendStr{ii}=['L' int2str(ii)];
end

subplot(2,2,3)
hold on
plot(times,avec); title('IRF EEG')
for ii=1:size(avec,2)
    legendStr{ii}=['L' int2str(ii)];
end
legend(legendStr)

subplot(2,2,4)
hold on
plot(times,avecR); title('IRF residual ')
legend(legendStr)


figure
for ii=1:size(avec,2)
    subplot(2,ceil(size(avec,2)/2),ii)
    hold on
    plot(times,avec(:,ii));
    title(['IRF - L' int2str(ii)])
    xlabel('time [ms]')
end

figure
for ii=1:size(prediction,3)
    subplot(2,ceil(size(prediction,3)/2),ii)
    hold on
    plot((prediction(:,randTrial,ii)),'r');
    plot((residual(:,randTrial,ii)),'b');
    if ii==1
        plot((stims(:,randTrial)),'--k')
    else
        plot((prediction(:,randTrial,ii-1)),'k');
    end
    title(['Layer ' int2str(ii)])
    legend('pred','residual','input')
    xlabel('time [ms]')
    xlim([3900 4300 ])

end


end