new;output file=c:\jae\univ\lew.out reset;
cybar=-177.435;
csig=277.332168044316;
crho1=-0.307304800605679;
format /le 23,15;
load xx[]=c:\jae\univ\lew.dat;
x0=reshape(xx,200,1);
y=x0;
{nadmu,nads,nadr}=nistuniv(y,cybar,csig,crho1);
{nadmu2,nads2,nadr2}=nistuni2(y,cybar,csig,crho1);
"above uses numerically superior meanc2 stdc3";
"  ";
format /ld 8,1;
"Lew" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2;
output file=c:\jae\univ\tableI.out on;
format /ld 8,1;
"Lew" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2;
end;


