-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathfig2.m
141 lines (116 loc) · 3.4 KB
/
fig2.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
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
% Interactions between multiple sources of short term plasticity
% during evoked and spontaneous activity at the rat calyx of Held
% J Physiol, 2008
%
% Matthias H. Hennig, Michael Postlethwaite, Ian D. Forsythe, Bruce
% P. Graham
% MHH: [email protected]; BPG: [email protected]
%
% This code produces a set of plots as shown in the paper in
% Figs. 2 and 3.
clear
% initialise graphics
h1 = figure(1);
clf
fs = [8.5 14]*1.5;
set(h1, 'PaperOrientation','portrait');
set(h1, 'PaperType','a4');
set(h1,'PaperUnits','centimeters');
set(h1,'Units','centimeters');
set(h1,'PaperPosition',[0 0 fs]);
set(h1,'Position',[0 19 fs]);
ms = 8;
colours = repmat([0 0 0],4,1)+repmat([0.4 0.8 0.6 0],3,1)';
%%%%%%%%%%%%%%%%%%%%%%%%%%
% load experimental data
load 'mean_epscs.mat'
frec = 0;
for freq = [10 20 50 100 ],
frec = frec + 1;
time = 1:freq;
exptime = 1:freq;
isi = repmat(1/freq, freq+1,1);
% simulate...
[resps, pprel, n, pb, nr, pf, rdes, final, retrieved ] = releasef(isi);
% ...and plot the results
subplot(3,2,1)
p = plot(exptime./freq,epsc_exp{frec}(:,2),'x');
set(p,'Color',colours(frec,:),'LineWidth',1,'MarkerSize',ms);
hold on
plot(time./freq,resps/resps(1),'k-','LineWidth',2,'Color',colours(frec,:))
hold on
subplot(3,2,2)
loglog(exptime./freq,epsc_exp{frec}(:,2),'x','Color',colours(frec,:),'LineWidth',1,'MarkerSize',ms)
hold on
loglog(time./freq,resps/resps(1),'k-','LineWidth',2)
subplot(3,2,3)
plot(time./freq,pf(1:length(time)),'k.-','LineWidth',2,'Color',colours(frec,:))
hold on
subplot(3,2,4)
plot(time./freq,pprel(1:length(time)),'k.-','LineWidth',2,'Color',colours(frec,:))
hold on
subplot(3,2,5)
plot(time./freq,n(1:length(time)),'k.-','LineWidth',2,'Color',colours(frec,:))
hold on
subplot(3,2,6)
plot(time./freq,retrieved,'k.-','LineWidth',1,'Color',colours(frec,:))
hold on
end
% set appropriate axis labels etc.
subplot(3,2,1)
set(gca,'XLimMode','manual');
set(gca,'XLim', [0.0 1.0]);
set(gca,'YLimMode','manual');
set(gca,'YLim', [0 1]);
set(gca,'FontName','Helvetica-Narrow');
set(gca,'FontSize',8);
ylabel('Normalised EPSC')
xlabel('Time/s')
subplot(3,2,2)
set(gca,'XLimMode','manual');
set(gca,'XLim', [0.01 1]);
set(gca,'YLimMode','manual');
set(gca,'YLim', [0.1 1.0]);
set(gca,'FontName','Helvetica-Narrow');
set(gca,'FontSize',8);
set(gca,'XTick',[0.01 0.1 1])
set(gca,'XTickLabel',[0.01 0.1 1])
set(gca,'YTick',[0.1 1])
set(gca,'YTickLabel',[ 0.1 1])
xlabel('Time/s')
subplot(3,2,3)
set(gca,'XLimMode','manual');
set(gca,'XLim', [0.0 1.0]);
set(gca,'YLimMode','manual');
set(gca,'YLim', [0.8 1.2]);
set(gca,'FontName','Helvetica-Narrow');
set(gca,'FontSize',8);
ylabel('Calcium Transient')
xlabel('Time/s')
subplot(3,2,4)
set(gca,'XLimMode','manual');
set(gca,'XLim', [0.0 1.0]);
set(gca,'YLimMode','manual');
set(gca,'YLim', [0.1 0.4]);
set(gca,'FontName','Helvetica-Narrow');
set(gca,'FontSize',8);
ylabel('Release Probability')
xlabel('Time/s')
subplot(3,2,5)
set(gca,'XLimMode','manual');
set(gca,'XLim', [0.0 1]);
set(gca,'YLimMode','manual');
set(gca,'YLim', [0 1.2]);
set(gca,'FontName','Helvetica-Narrow');
set(gca,'FontSize',8);
ylabel('Release Pool Occupancy')
xlabel('Time/s')
subplot(3,2,6)
set(gca,'XLimMode','manual');
set(gca,'XLim', [0.0 1.0]);
%set(gca,'YLimMode','manual');
%set(gca,'YLim', [0.0 1.2]);
set(gca,'FontName','Helvetica-Narrow');
set(gca,'FontSize',8);
xlabel('Time/s')
ylabel('Vesicles Retrieved')