new;output file= c:\jae\univ\mavro.out reset;
errv3=zeros(3,1);
cybar=2.001856;
csig=0.000429123454003053;
crho1=0.937989183438248;
format /le 23,15;
load xx[]=c:\jae\univ\mavro.dat;
x0=reshape(xx,50,1);
y=x0;
{nadmu,nads,nadr}=nistuniv(y,cybar,csig,crho1);
" "; 
{nadmu2,nads2,nadr2}=nistuni2(y,cybar,csig,crho1);
"numericallyimproved mean and std dev above";
" ";
format /ld 8,1;
"Mavro" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2;

output file=c:\jae\univ\tableI.out on;
format /ld 8,1;
"Mavro" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2;

end;


