R包“ez”中的ezANOVA函数中的“wid”参数究竟是什么?

机器算法验证 r 方差分析
2022-03-19 21:52:26

我对这一切都很陌生,我正在测试不同的方法来对我的数据执行双向 III 型方差分析。

  • 在用;拟合线性回归后,我已经尝试过从anova()包中statslm()
  • 我已经Anova()从包中尝试过,使用相同的线性回归(这给了我与我使用时car相同的结果- 我认为默认使用 I SS 类型?)。anova()type="II"anova()
  • 我现在正试图ezANOVA()ez包中使用。

对于最后一个,我无法理解wid=.()参数是什么(即使阅读帮助),并且由于它不是可选的,我不能将其留空。我正在尝试使用如下,其结果是:

> attach(data)
> library("ez")
> ezANOVA(data=data, dv=.(AG.DW), wid=.(), within=.(Genotype, Treatment), type=3)
Warning: Converting "" to factor for ANOVA.
Error in sort.list(y) : 'x' must be atomic for 'sort.list'
Have you called 'sort' on a list?

这是正确的脚本吗?wid什么,我应该用什么填充它?

关于我的数据,列GenotypeTreatment是我的两个因素,我想看看在查看我的植物的地上干重(列AG.DW)时是否存在相互作用。我的数据是平衡的。

如果这里的信息丢失或不准确,我很抱歉:这是我在这里的第一个贡献,我现在只是发现统计数据(我看不到如何加入我的数据文件)。

1个回答

首先,你熟悉 R 的帮助系统吗?键入?ezANOVA将显示该函数的文档,您可以在其中读到wid参数是:

.() 对象指定 data 中包含指定 case/Ss 标识符的变量的列。

根据您迄今为止指定模型的方式,您的设计似乎是受试者内的设计,其中每个受试者(这是一个地块/植物/区域?)被分配到基因型和治疗的每个级别。然后,您的wid论点将是标识单个植物/地块的列,例如。wid=.(plant.number).

如果这些似乎都不符合您的实验设计,那么可能是您的基因型和治疗因素中的一个或两个实际上是受试者之间的变量,您应该将它们包含在between = .()参数中。