2
我试图在R中运行anova()
并遇到一些困难。这是我迄今为止帮助解决我的问题所做的工作。按组拟合线性模型/方差分析
这里是我的数据到这一点的str()
。
str(mhw)
'data.frame': 500 obs. of 5 variables:
$ r : int 1 2 3 4 5 6 7 8 9 10 ...
$ c : int 1 1 1 1 1 1 1 1 1 1 ...
$ grain: num 3.63 4.07 4.51 3.9 3.63 3.16 3.18 3.42 3.97 3.4 ...
$ straw: num 6.37 6.24 7.05 6.91 5.93 5.59 5.32 5.52 6.03 5.66 ...
$ Quad : Factor w/ 4 levels "NE","NW","SE",..: 2 2 2 2 2 2 2 2 2 2 ...
列r是表示一个数值,其行中的字段中的各个描绘驻留 c栏是表示数值的列的单个情节驻留
柱四对应于地理位置在现场其中每个小区驻留
Quad <- ifelse(mhw$c > 13 & mhw$r < 11, "NE",ifelse(mhw$c < 13 & mhw$r < 11,"NW", ifelse(mhw$c < 13 & mhw$r >= 11, "SW","SE")))
mhw <- cbind(mhw, Quad)
我有适合lm()
如下
nov.model <-lm(mhw$grain ~ mhw$straw)
anova(nov.model)
这是整个田间的anova()
,它测试数据集中每个地块的谷物产量与秸秆产量。
我的麻烦是,我想运行一个单独的anova()
我的数据的四列来测试每个象限的粮食产量和秸秆产量。
也许一个with()
可能会解决这个问题。我以前从未使用过,现在正在学习R的过程中。任何帮助将不胜感激。
谢谢。这提供了数据的“系数”。如果我想产生完整的anova摘要来检验假设,“东部象限对西象限的谷物产量是否有差异”。我怎样才能用'fit'生产这个' – pc8807
Worked Perfect。非常感谢你的帮助。我对这些问题的简单性表示歉意。正如我所说的,我正在学习R的过程中,稍微回过头来看看我的语法。 – pc8807
该链接相当有帮助。这清除了我与R一起工作时出现的许多问题。 – pc8807