求助！四因素三水平的正交实验编程

shiyiming

我这有一组四因素三水平的正交实验结果，要进行方差分析，需建立sas程序。我编了一个，但是结果不行。哪位大侠帮帮忙，指正指正！！

正交表省略了
sas自编程序：
data dx;
input x1 x2 x3 x4 y;
cards;
1 1 1 1 5.02
1 2 2 2 7.04
1 3 3 3 14.11
2 1 2 3 13.16
2 2 3 1 14.76
2 3 1 2 4.75
3 1 3 2 16.36
3 2 1 3 7.64
3 3 2 1 9.47
;
proc anova;
class x1 x2 x3 x4;
model y=x1 x2 x3 x4;
means x1 x2 x3 x4;
run;



运行结果：




                                  Analysis of Variance Procedure

Dependent Variable: Y

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

Model                    8                 5.79920000                0.72490000      .       .

Error                    0                  .                         .

Corrected Total          8                 5.79920000

                  R-Square                       C.V.                  Root MSE             Y Mean

                  1.000000                          0                                   5.83666667


Source                  DF                   Anova SS               Mean Square  F Value    Pr > F

X1                       2                 0.41086667                0.20543333      .       .
X2                       2                 0.10806667                0.05403333      .       .
X3                       2                 4.44740000                2.22370000      .       .
X4                       2                 0.83286667                0.41643333      .       .

shiyiming

no enough observations maybe.

shiyiming

这样说，四因素三水平的试验用l9（34）的表不行咯

shiyiming

Your design table is OK (9 is the number for 3^4), but you don't have adequate # of measurements for the proposed model.
Anyway, don't trust me, I'm still dangling at the bottom of the food chain!

shiyiming

So could you tell me how to Analysis it?
               and what to write

shiyiming

it's beyond my capacity. i don't know how to fit the model either.
my hunch is that you don't have enough data for the proposed analysis. the orthogonal design is a magic to reduce the # of runs, but several repetitions per run may still be desired to fit the model.

shiyiming

oh,Isee.You are wonderful! <!-- s:D --><img src="{SMILIES_PATH}/icon_biggrin.gif" alt=":D" title="Very Happy" /><!-- s:D -->  
   thank you very much!

shiyiming

四个因子三个水平，每个因子要有2个自由度，所以因子总共就要有8个自由度，而你总共只有9个观测值，即只有8个自由度可以分配，所以8-8＝0就没有自由度给 残差（Error）了。

你可以选择其他的正交表来满足有足够的自由度可以分配。

shiyiming

实验次数不够。