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

