new;

/*load X[9314,4]=C:\GAUSSwin\x\boe\gausstlk\FT30_C.PRN;*/
load X[9314,4]=FT30_C.PRN;
x[1,.]; x[314,.];

/*load xx[]= C:\GAUSSwin\x\boe\gausstlk\FT30_C.PRN;*/
load xx[]= FT30_C.PRN;
nn=rows(xx);
n = int(nn/4);
if nn - 4*n /= 0 ;
	"number " nn " of data points not divisible by 4";
	stop;
endif;

xx = reshape(x,n,4);

meanc(x)~meanc(xx);

xx4 = xx[.,4];
save xx4;
/*save path=c:\gausswin\x\boe\gausstlk xx4;*/
load xx4; meanc(xx4);