我进行了用树进行的田间试验，其中我们在随机完全区组设计（RCBD）后在一个小种植园中种植了不同的基因型。 现在我想在R中进行分析，但我对如何做到有一些疑问。 简而言之，我有3个区块和5个基因型，除了我们测量的几个变量之外，其中一个是HEIGHT。我用来进行ANOVA测试的代码是：
fit <- lm(HEIGHT~GENOTYPE+BLOCK,data=data)
anova(fit)
在某些网页中，我看到他们写道：
lm(HEIGHT~BLOCK+GENOTYPE,data=data)
我不知道哪个是完全不同的但是我已经尝试了两个线性模型（lm）并且结果不一样。问题很简单：为什么？当我写'＃34; Height~Genotype + Block＆＃34;时，我告诉R的确是什么？当我告诉＆＃34; Height~Block + Genotype＆＃34;？另一个问题是：我是否正确地进行阻断ANOVA？
非常感谢您提前!!
答案 0 :(得分：1)
在模型中的主效应之后应该考虑该块，因为您要确定主效应相对于阻塞因子描述了多少总变化。当首先考虑阻塞因子时，它实际上起主效应的作用，GENOTYPE效应成为阻塞因子。您也可以使用以下代码执行相同的分析。
fit <- aov(HEIGHT ~ GENOTYPE + BLOCK, data=data)
summary(fit)