四因素三水平无空列正交试验方差分析怎么弄啊？

shiyiming

四因素三水平无空列正交试验，我写了程序
data ortho;
input a b c d y;
cards;
1 1 1 1 1.1881
1 2 2 2 1.2184
1 3 3 3 1.2094
2 1 2 3 1.2543
2 2 3 1 1.1792
2 3 1 2 1.2250
3 1 3 2 1.2170
3 2 1 3 1.2522
3 3 2 1 1.1658
;
proc anova data=ortho;
class a b c d;
model y=a b c d;
means a b c d;
run;
但是这样的话，结果因为没有误差来源，没有F值：
Dependent Variable: y

                                                 Sum of
         Source                      DF         Squares     Mean Square    F Value    Pr > F

         Model                        8      0.00742832      0.00092854        .       .

         Error                        0      0.00000000       .

         Corrected Total              8      0.00742832


                          R-Square     Coeff Var      Root MSE        y Mean

                          1.000000           .               .      1.212156


         Source                      DF        Anova SS     Mean Square    F Value    Pr > F

         a                            2      0.00030354      0.00015177        .       .
         b                            2      0.00067300      0.00033650        .       .
         c                            2      0.00059608      0.00029804        .       .
         d                            2      0.00585571      0.00292785        .       .
我看到有人做出来的误差自由度为2，而且是有F值的，不知道这样的程序是怎么编写出来的？？
万分感激！！！

shiyiming

你把所有的自由度都控制死了， 还谈什么误差来源？它想犯错误都犯不了。
别的不知道，反正去掉一个因素，自由度就增加两个。

shiyiming

这个我是明白的，这样的话重复试验是可以有误差的，但我看到有人的试验安排也是这样的，但是有误差的，而且误差的自由度是2，就搞出清楚是怎么做出来的？