求教 不均衡数据SAS统计

shiyiming

试验为随机区组，
施氮量为600、900、1200 kg/ha，
灌水量为FC高量灌水和FR节水灌水，
水氮耦合处理为FCN1200、FCN900、FRN1200、FRN900和FRN600，
没有设计FCN600处理

1）原始SAS程序，如下：
input block N $ F$ yield@@;
CARDS;
1.00         N1200        FC        313.65
2.00         N1200        FC        326.87
3.00         N1200        FC        239.52
1.00         N900        FC        189.62
2.00         N900        FC        309.04
3.00         N900        FC        292.19
1.00         N1200        FR        246.19
2.00         N1200        FR        210.35
3.00         N1200        FR        270.48
1.00         N900        FR        229.70
2.00         N900        FR        181.44
3.00         N900        FR        188.64
1.00         N600        FR        157.04
2.00         N600        FR        177.46
3.00         N600        FR        154.14

;
proc mixed;
CLASS block N F;
MODEL yield = N F N*F block/ ddfm=kenwardroger ;
LSMEANS N F N*F/pdiff cl adjust = Tukey ;
RUN;

由于试验没有设计FCN600处理，严重不均衡导致部分结果无法进行统计分析。因此将该程序拆成两个运行，如下：


DATA YIELD;
input block N $ F$ yield@@;
CARDS;
1.00         N1200        FC        313.65
2.00         N1200        FC        326.87
3.00         N1200        FC        239.52
1.00         N900        FC        189.62
2.00         N900        FC        309.04
3.00         N900        FC        292.19
1.00         N1200        FR        246.19
2.00         N1200        FR        210.35
3.00         N1200        FR        270.48
1.00         N900        FR        229.70
2.00         N900        FR        181.44
3.00         N900        FR        188.64
1.00         N600        FR        157.04
2.00         N600        FR        177.46
3.00         N600        FR        154.14

;
proc mixed;
CLASS block N F;
MODEL yield = N*F/ ddfm=kenwardroger ;
LSMEANS N*F/pdiff cl adjust = Tukey ;
random block;
RUN;
proc mixed;
CLASS block N F;
MODEL yield = N F/ ddfm=kenwardroger ;
LSMEANS N F/pdiff cl adjust = Tukey ;
random block;
RUN;

请问该程序是不是能用？
因为我发现该程序显示水氮交互作用显著，但是数据上并未呈现出交互作用？

shiyiming

有什么特别的理由支持直接加入交叉项然而忽略main effect？
京剧

shiyiming

proc glm 可以吗 ？