load('atlas22493_burst.mat')


%%
length(pre.burst.ind)/3
mean(pre.IBI)
std(pre.IBI)
quantile(pre.IBI,[0.25, 0.5, 0.75])

length(drug.burst.ind)/max(drug.time)*3600
mean(drug.IBI)
std(drug.IBI)
quantile(drug.IBI,[0.25, 0.5, 0.75])


%%
x1(1,:) = [length(pre.burst.ind)/3,mean(pre.IBI),std(pre.IBI),quantile(pre.IBI,[0.5, 0.25, 0.75])];
%x1(1,:) = [length(pre.burst.ind)/min(pre.time)*-3600,mean(pre.IBI),std(pre.IBI),quantile(pre.IBI,[0.5, 0.25, 0.75])];
x1(2,:) = [length(drug.burst.ind)/max(drug.time)*3600,mean(drug.IBI),std(drug.IBI),quantile(drug.IBI,[0.5, 0.25, 0.75])];
