/*---------------------------- | Example: die -----------------------------*/ data die; input face $ count; datalines; 1 3 2 7 3 5 4 10 5 2 6 3 ; run; proc freq; weight count; tables face/ chisq; /*tables face /nocum all chisq testp=(16.67 16.67 16.67 16.67 16.67 16.67);*/ run; /********computation of residuales, deviance residuals, X2 and G2 ****/ data cal; set die; pi=1/6; ecount=30*pi; res=(count-ecount)/sqrt(ecount); devres=sqrt(abs(2*count*log(count/ecount)))*sign(count-ecount); run; proc print;run; proc sql; select sum((count-ecount)**2/ecount) as X2, 1-probchi(calculated X2,5) as pval1, 2*sum(count*log(count/ecount)) as G2, 1-probchi(calculated G2,5) as pval2 from cal; quit;