load moran_n1024_x0512.dat;
x = moran_n1024_x0512(:,1);
t = moran_n1024_x0512(:,2);
av=mean(t)
av=mean(x)
hist(t,20)
quantile(t,0.5)
tiling=0.05*(1:20)
xquant=quantile(t,tiling)
yquant=0.05*(1:20)
plot(xquant,yquant,'o')
yquant2=1-yquant
plot(xquant,yquant2,'o')
yquant3=log(yquant2+0.05)
plot(xquant,yquant3,'o')
x'*t
x'*t/sum(x)
(1-x')*t/(sum(1-x))