Percept-LFP-Toolbox/plotStimRunFS.m at main · C-ACT/Percept-LFP-Toolbox · 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
140
141
142
143
144
145
146
147
148
149
150
151
152
153
%S. N. Pitts: 06-08-2022
function [] = plotStimRunFS(C_uniform,S,k,logical)
close all


    %summary(C_uniform)

    % figure
    % hold on
    % Peak_Power = [C.R_peaks{:}];
    % histogram(Peak_Power)
    % xlabel('Power')
    % ylabel('Quantity')
    % title('All Beta Power Peaks, Right')
    % hold off

    % figure
    % Peak_Power = [C.L_peaks{:}];
    % histogram(Peak_Power)
    % xlabel('Power')
    % ylabel('Quantity')
    % title('All Beta Power Peaks, Left')
    % hold off

    % figure
    % hold on
    % Trough_Power = [C.R_trough{:}];
    % histogram(Trough_Power)
    % xlabel('Power')
    % ylabel('Quantity')
    % title('All Beta Power Troughs, Right')
    % hold off


    % figure
    % Trough_Power = [C.L_trough{:}];
    % histogram(Trough_Power)
    % xlabel('Power')
    % ylabel('Quantity')
    % title('All Beta Power Troughs, Left')
    % hold off

    %% Right Data
    figure
    subplot(2,1,1)
    hold on
    yyaxis left
    plot(C_uniform.Time,C_uniform.("Daily Avg R_LFP FS"))
    xlabel('Time')
    ylab = sprintf('Running FS (k=%d window size)',k);
    ylabel(ylab)

    yyaxis right
    right = [S.("Band of Interest (Hz)"){:}];
    right = right(1:2:end);
    plot(S.Date,right)
    ylabel('Band of Interest (Hz)')

    title('Daily Avg R_LFP FS')
    hold off

    subplot(2,1,2)
    hold on
    yyaxis left
    stimx = [S.("Date Time (HH:mm:ss)"){:}];
    stimy = [];
    for i = 1:length(S.("Stim Amp - R/CH1 (mA)"))
        stimy = [stimy;(S.("Stim Amp - R/CH1 (mA)"){i})];
    end
    scatter(stimx,stimy')
    ylabel('Stim Amp (mA)')

    yyaxis right
    cla
    right = [S.("Channel Names"){:}];
    right = right(1:2:end);
    cat = unique(right);
    y=zeros(1,length(right));
    %create categorical array
    for i = 1:length(cat)
        idx = find(right==cat(i));
        y(idx) = i;
    end
    c = categorical(y,1:length(cat),[cat]);
    right = [S.("Channel Names"){:}];
    right = right(1:2:end);
    plot(S.Date,c)
    ylabel('Lead Configuration')

    title('Stimulation and Channels')
    hold off
    if logical
        savefig("RightRunFS-Raw")
    else
        savefig("RightRunFS-MedicalLP")
    end

    %% Left Data
    figure
    subplot(2,1,1)
    hold on
    yyaxis left
    plot(C_uniform.Time,C_uniform.("Daily Avg L_LFP FS"))
    xlabel('Time')
    ylab = sprintf('Running FS (k=%d window size)',k);
    ylabel(ylab)

    yyaxis right
    left = [S.("Band of Interest (Hz)"){:}];
    left = left(1:2:end);
    plot(S.Date,left)
    ylabel('Band of Interest (Hz)')

    title('Daily Avg L_LFP FS')
    hold off

    subplot(2,1,2)
    hold on
    yyaxis left
    stimx = [S.("Date Time (HH:mm:ss)"){:}];
    stimy = [];
    for i = 1:length(S.("Stim Amp - L/CH2 (mA)"))
        stimy = [stimy;(S.("Stim Amp - L/CH2 (mA)"){i})];
    end
    scatter(stimx,stimy')
    ylabel('Stim Amp (mA)')

    yyaxis right
    cla
    left = [S.("Channel Names"){:}];
    left = left(1:2:end);
    cat = unique(left);
    y=zeros(1,length(left));
    %create categorical array
    for i = 1:length(cat)
        idx = find(left==cat(i));
        y(idx) = i;
    end
    c = categorical(y,1:length(cat),[cat]);
    left = [S.("Channel Names"){:}];
    left = left(1:2:end);
    plot(S.Date,c)
    ylabel('Lead Configuration')

    title('Stimulation and Channels')
    hold off
    if logical
        savefig("LeftRunFS-Raw")
    else
        savefig("LeftRunFS-MedicalLP")
    end

end