load data.dat;

data=data-median(data);


subplot(4,1,1);
plot(data);
ylabel("raw values - median")


subplot(4,1,2);
qqplot(data);
ylabel("CDF")


subplot(4,1,3);
[y,x]=hist(data,100,1);
z=normpdf(x,mean(data),std(data));
z2=normpdf(x,median(data),std(data));
plot(x,y,x,z,x,z2);
legend("samples","stdmean","stdmedi",2);
ylabel("stddev")


subplot(4,1,4);
y=y+1e-12;
z=z+1e-12;
z2=z2+1e-12;
plot(x,20*log10(y),x,20*log10(z),x,20*log10(z2));
ylabel("log of stddev")


print("out9.png","-color","-dpng");