-
Notifications
You must be signed in to change notification settings - Fork 2
/
plotfigure9.m
77 lines (69 loc) · 1.87 KB
/
plotfigure9.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
getplotid
figure(9)
%% panel 2,4,6
if (VIR_ind ==4 && freq_ind ==4)
AMPK_AUC = AMPK_AUC./max(AMPK_AUC,[],'all');
mTORC1_AUC = mTORC1_AUC./max(mTORC1_AUC,[],'all');
mTORC2_AUC = mTORC2_AUC./max(mTORC2_AUC,[],'all');
subplot(3,2,2)
currentheatmap1 = heatmap(freq_vary,parscan,AMPK_AUC);
ylabel('\bf{V_{IR} Multiplier}')
xlabel('\bf{Stimulus Frequency [Hz]}')
title('Normalized AMPK AUC')
plotformatheatmap
currentheatmap1.FontSize = 12;
%currentheatmap1.ColorScaling ='log';
set(gca,'CellLabelColor','none')
caxis([0.3,1])
%% panel 2
subplot(3,2,4)
currentheatmap2 = heatmap(freq_vary,parscan,mTORC1_AUC);
ylabel('\bf{V_{IR} Multiplier}')
xlabel('\bf{Stimulus Frequency [Hz]}')
title('Normalized mTORC1 AUC')
plotformatheatmap
currentheatmap2.FontSize = 12;
%currentheatmap2.ColorScaling ='log';
set(gca,'CellLabelColor','none')
caxis([0.3,1])
%% panel 3
subplot(3,2,6)
currentheatmap3 = heatmap(freq_vary,parscan,mTORC2_AUC);
ylabel('\bf{V_{IR} Multiplier}')
xlabel('\bf{Stimulus Frequency [Hz]}')
title('Normalized mTORC2 AUC')
plotformatheatmap
caxis([0.3,1])
currentheatmap3.FontSize = 12;
%currentheatmap3.ColorScaling ='log';
set(gca,'CellLabelColor','none')
colormap(colorcet('L18'))
end
%% panel 1,3,5
subplot(3,2,1)
hold on
plot(t,y(:,pAMPK),'Color',colorvector(VIR_ind,:),'LineStyle',char(markervector(freq_ind)))
xlabel('Time [s]')
ylabel('Concentration [\muM]')
plotformat
xlim([-50,100])
xticklabels('auto')
subplot(3,2,3)
title('AMPK')
hold on
plot(t,y(:,pmTORC1),'Color',colorvector(VIR_ind,:),'LineStyle',char(markervector(freq_ind)))
xlabel('Time [s]')
ylabel('Concentration [\muM]')
plotformat
xlim([-50,100])
xticklabels('auto')
title('mTORC1')
subplot(3,2,5)
hold on
plot(t,y(:,pmTORC2),'Color',colorvector(VIR_ind,:),'LineStyle',char(markervector(freq_ind)))
xlabel('Time [s]')
ylabel('Concentration [\muM]')
plotformat
xlim([-50,100])
xticklabels('auto')
title('mTORC2')