new;
output file=c:\jae\univ\pidigits.out reset;
cybar=4.53480000000000;
csig=2.86733906028871;
crho1=-0.00355099287237972;
format /le 23,15;  x0=0;
load xx[]=c:\jae\univ\pidigits.asc;
x0=reshape(xx,5000,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 ";
"  ";

"pidigits" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2 ;

output file=c:\jae\univ\pdigits.out off;
output file=c:\jae\univ\TableI.out reset;
format /ld 8,1;

"pidigits" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2 ;

end;


