new;output file= c:\jae\univ\numacc4.out reset;
cybar=10000000.2;
csig=0.1;
crho1=-0.999;
format /ld 30,15;
load xx[]=c:\jae\univ\numacc4.dat;
x0=reshape(xx,1001,1);
y=x0;
{nadmu,nads,nadr}=nistuniv(y,cybar,csig,crho1);
" ";
"now numerically superior means and std dev";
{nadmu2,nads2,nadr2}=nistuni2(y,cybar,csig,crho1);
" above uses meanc2, stdc3 ";
format /ld 8,1;
"Numacc4" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2;
output file=c:\jae\univ\tableI.out on;
format /ld 8,1;
"Numacc4" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2;
end;
