求助！生物统计三因素方差分析题！

我的答案
setwd("D:/RBOOK/RBOOK")
CL <- read.table(file = "li8-2(2).txt", header=TRUE)
CL$F时间 <- factor(CL$时间)
CL$F浓度 <- factor(CL$浓度)
CL$F区组 <- factor(CL$区组)
CL$F类型 <- factor(CL$类型)
aov.CL <- aov(CL$产量 ~ CL$F区组+ CL$F时间 + CL$F浓度 + CL$F时间 + CL$F浓度: CL$F类型 + CL$F时间: CL$F浓度+CL$F时间: CL$F类型+CL$F时间: CL$F浓度:CL$F类型, data = CL)
summary(aov.CL)
Df Sum Sq Mean Sq F value Pr(>F)
CL$F区组 1 0.7 0.7 1.291 0.265
CL$F时间 1 1198.0 1198.0 2195.892 < 2e-16 ***
CL$F浓度 4 904.6 226.2 414.558 < 2e-16 ***
CL$F浓度:CL$F类型 10 69.0 6.9 12.645 4.28e-08 ***
CL$F时间:CL$F浓度 4 404.6 101.2 185.433 < 2e-16 ***
CL$F时间:CL$F类型 2 1.3 0.7 1.195 0.317
CL$F时间:CL$F浓度:CL$F类型 8 29.8 3.7 6.832 4.92e-05 ***
Residuals 29 15.8 0.5
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

书上手算的答案
时间间 ( A ) F =1258.087** ，浓度间 ( B ) F =248.222** ，籽粒类型间 (C) F=7.980** ， A × B 的 F =98.204** ， A × C 的 F =0.424 ， B × C 的 F=5.017** ， A × B × C 的 F =2.433* 。

数据录入
区组时间浓度类型产量
1 1 1 1 7.0
1112 13.5
111310.7
112112.8
112213.2
112312.4
113122.0
113220.8
113322.6
114121.3
114219.0
114321.3
115124.2
115224.6
115324.5
1211 3.6
1212 4.7
1213 8.7
122110.7
1222 9.8
1223 9.6
1231 4.7
1232 2.7
1233 3.4
124112.4
124212.4
124313.0
125112.6
125214.0
125314.8
2111 6.5
211213.8
211310.3
212111.4
212214.2
212313.2
213121.8
213221.4
213321.8
214120.3
214219.6
214322.4
215123.2
215223.8
215324.2
2211 1.5
2212 4.9
2213 8.5
2221 8.8
222210.5
2223 9.7
2231 3.4
2232 4.2
2233 4.2
224110.5
224213.2
224312.7
225113.7
225214.2
225312.6

不知道哪里不对。。。

#创建样本测量值向量
obs<-c(7,6.5,13.5,13.8,10.7,10.3,3.6,1.5,4.7,4.9,8.7,3.5,12.8,11.4,13.2,14.2,12.4,13.2,10.7,8.8,9.8,10.5,9.6,9.7,22,21.8,20.8,21.4,22.6,21.8,4.7,3.4,2.7,4.2,3.4,4.2,21.3,20.3,19,19.6,21.3,22.4,12.4,10.5,12.4,13.2,13,12.7,24.2,23.2,24.6,23.8,24.5,24.2,13.6,13.7,14,14.2,14.8,12.6)
#创建A,B,C因子
A<-factor(gl(2,6,len=60))
B<-factor(gl(5,12))
C<-factor(gl(3,2,len=60))
#作方差分析
anova(aov(obs~A*B*C))
根据你写的函数，暂判断你对因子分析还未完全弄明白，很简单的问题，多看看例子，多想想就搞清楚了。