clc
clear all
aa=load('../Données/LOD.txt');
t=aa(:,1);
lod=aa(:,2);
bb=load('../Données/soi.dat');
for i=1:19
   ts((i-1)*12+1:i*12)=bb(i,1)+(0:11)/12;
   soi((i-1)*12+1:i*12)=bb(i,2:end);
end
soi=soi(ts>min(t));
ts=ts(ts>min(t));
soi=soi(ts<max(t));
ts=ts(ts<max(t));
lod=interp1(t,lod,ts);
lodp=lod(soi>0);
lodm=lod(soi<0);
vp=var(lodp);
vm=var(lodm);
np=length(lodp);
nm=length(lodm);
y=testvar(vp,vm,np,nm);