-
Notifications
You must be signed in to change notification settings - Fork 0
/
dipRMS_20150804.m
55 lines (45 loc) · 1.51 KB
/
dipRMS_20150804.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
%% set initial conditions
close all
mu = linspace(0,50,20);
inds = [1:20];
ws=1;
ws2=5000;
indstart=1;
indend=15e3;
for j = 1:length(Xd_pulse)
x = Xd_pulse{j} - smooth(Xd_pulse{j},ws2);
x = smooth(x,ws);
x = x(indstart:indend);
[dip(j),dipp(j)]=HartigansDipSignifTest(x,1);
RMSmag(j) = std(x);
end
subplot(2,1,1);plot(mu(inds),dip(inds));axis tight
set(gca,'XDir','reverse')
xlabel('Control parameter');ylabel('Dip statistic');
subplot(2,1,2);plot(mu(inds),RMSmag(inds));axis tight
set(gca,'XDir','reverse')
xlabel('Control parameter');ylabel('RMS mag.');
%% bootstrapping
Nboot = 1e3;
for j = 1:length(Xd_pulse)
x = Xd_pulse{j} - smooth(Xd_pulse{j},ws2);
x = smooth(x,ws);
x = x(indstart:indend);
[dip(j),dipp(j)]=HartigansDipSignifTest(x,Nboot);
RMSmag(j) = std(x);
RMSmagboot{j} = bootstrp(Nboot,@(x) std(x), x);
dipboot{j} = bootstrp(Nboot,@(x) HartigansDipSignifTest(x,1), x);
RMSmagsem(j) = std(RMSmagboot{j});
dipsem(j) = std(dipboot{j});
disp(num2str(j));
end
subplot(2,1,1);errorbar(mu(inds),dip(inds),dipsem(inds));axis tight
set(gca,'XDir','reverse')
xlabel('Control parameter');ylabel('Dip statistic');
subplot(2,1,2);errorbar(mu(inds),RMSmag(inds),RMSmagsem(inds));axis tight
set(gca,'XDir','reverse')
xlabel('Control parameter');ylabel('RMS mag.');
%% save
save('/Users/joshsalvi/Documents/Lab/Lab/Clamp Data/2015-08-10.01/Ear 1/Cell 1/dipRMSdata-Xd_pulse{1,3}.mat',...
'Xd_pulse','indstart','indend','ws','ws2','Fs','mu','dip','dipp','RMSmag','RMSmagboot','dipboot',...
'RMSmagsem','dipsem','Nboot','inds','-v7.3');