close all
clear all


 disp('Neuronfires QQPLOTS')
 lw = 2.5; 
 set(0, 'DefaultAxesFontSize', 16);
 fs = 15;
 msize = 10;

load 'C:\STAT\Estim\Estimdat\neuronfires.mat'
x = diff(neuronfires);

% We argued E are exponential
p =0.5/988:1/988:1  %
figure(1)
subplot(221)
plot( expinv(p, 1), sort(x),'*' )
subplot(222)
plot( chi2inv(p, 3), sort(x),'*' )
subplot(223)
plot( norminv(p), sort(x),'*' )
subplot(224)
plot( tinv(p, 3), sort(x),'*' )

figure(2)
%qqplot(expinv(p,1),x)
qqplot(gaminv(p, 1.4, 1), x)




