new;
output file=c:\jae\anova\simon6\simon6.out reset;
"Anova problems on NIST datasets";
cfstat= 2.00100000000000E+03;

load xx[]= "c:\\jae\\anova\\simon6\\simon6.dat";
x=reshape2(xx,9);
{withi,betw,f}=anov(x);
"F=" f;
"within= " withi;
"Between= " betw;

NAD= numacc(f,cfstat);
format /ld  14,1;
"NAD=" nad;
{withi,betw,f}=anov2(x);
"F=" f;
"within= " withi;
"Between= " betw;

NAD= numacc(f,cfstat);
format /ld  14,1;
"numerically improved NAD=" nad;
end;

