new;output file= c:\jae\univ\michelso.out reset;
cybar=299.8524;
csig=0.0790105478190518;
crho1=0.535199668621283;
format /le 23,15;
load xx[]=c:\jae\univ\michelso.dat;
x0=reshape(xx,100,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;
"Michelso" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2;

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

end;


