.
. Starting Step 3.2
anova diameter machine
.
delete stderr nim1 ucl lcl level ni levmean
skip 0
read dpst2f.dat junk1 level ni levmean
let n = size diameter
let stderr = ressd/sqrt(ni)
let ucl = levmean + stderr*tppf(0.975,resdf)
let lcl = levmean + stderr*tppf(0.025,resdf)
set write format f6.0,f8.0,f10.6,f10.6,f10.6,f10.6
print "Level  Number      Mean Std Error Lower 95% Upper 95%"
print level ni levmean stderr lcl ucl
set write format
.
. Finished Step 3.2: Page Back for Output