stacksize(2e8)
getd ~/TANIT/SummerSchool15/PRACTICE
loadmatfile('~/TANIT/SummerSchool15/PRACTICE/Seizure1.mat');

plot(adat)

pplot(adat)

plot((0:20000)'/20,abs(fft(adat(10000:30000,11),-1)))
plot((0:20000)'/20,abs(fft(adat(90000:110000,11),-1)));

chn=11;
astart=1e5;
alength=1e4;
st=1;
orde=20;
ad=zeros(alength,orde);
for i=1:orde
ad(:,i)=adat(astart+(i-1)*st+1:astart+alength+(i-1)*st,chn);
end
at=adat(astart+(i+1)*st:astart+alength+(i+1)*st-1,chn);
w=ad\at;
//w=pinv(ad)*at;

ar=ad*w;
er=at-ar;0

plot(at)
plot(ar,'r')
plot(er,'g')
ser=sqrt(sum(er.^2)/alength)

plot(w($:-1:1))

chn=[11,32];
schn=size(chn,2);
astart=1e5;
alength=1e4;
st=1;
orde=5;
mser=zeros(2,2);
er2=zeros(alength,2);
for j=1:schn
ad=zeros(alength,orde);
for i=1:orde
ad(:,i)=adat(astart+(i-1)*st+1:astart+alength+(i-1)*st,chn(j));
end
at1=adat(astart+(i+1)*st:astart+alength+(i+1)*st-1,chn(j));
w=ad\at1;
//w=pinv(ad)*at;
ar1=ad*w;
er1=at1-ar1;
ser=stdev(er1);
mser(j,1)=ser;

ad=zeros(alength,2*orde);
for i=1:orde
ad(:,(i-1)*2+1:i*2)=adat(astart+(i-1)*st+1:astart+alength+(i-1)*st,chn);
end
at2=adat(astart+(i+1)*st:astart+alength+(i+1)*st-1,chn(j));
w=ad\at2;
//w=pinv(ad)*at;
ar2=ad*w;
er2(:,j)=at2-ar2;
er=at2-ar2;
ser=stdev(er);
mser(j,2)=ser;
end
cgr=log(mser(:,1))-log(mser(:,2));
co=er2(:,1)'*er2(:,2)/(alength-1);

ccgr=log(mser(1,1)*mser(1,2))-log(mser(1,2)^2-co^2)