new;output file= c:\jae\univ\lottery.out reset;
cybar=518.958715596330;
csig=291.699727470969;
crho1=-0.120948622967393;
format /le 23,15;
load xx[]=c:\jae\univ\lottery.dat;
x0=reshape(xx,218,1);
y=x0;
dmean=y-meanc(y)';
/* rho1=autocor(dmean,0,2); */
yi=y[2:218];
ym1=y[1:217];
rho1=((yi-meanc(y))'(ym1-meanc(y)))/((y-meanc(y))'(y-meanc(y)));

"Univariate Analysis-Lottery";
{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   Lottery";
"Lottery" nadmu~nads~nadr ;
 "meanc2" nadmu2~nads2~nadr2 ;

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

end;


