new;
output file=c:\jae\anova\simon2\simon2.out reset;
"Anova problems on NIST datasets";
#include "c:\jae\anova\simon2\simon2.dat";
{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;

