Nested/SplitModel - RepeatedMeasures/MixedModel ANOVA:R 中的嵌套和脚本级别

机器算法验证 r 方差分析 重复测量 裂区 嵌套数据
2022-04-03 02:12:52

我的数据集有以下变量:

  • 治疗(4种-固定)
  • 位置(8 个位置 - 固定)
  • 位置位置(每个位置 3 个位置 - 固定)
  • 在每个位置取样(每个位置随机 3 个样本)
    • 时间(两个采样时间 - 固定)
    • 矿化率(作为样本分析的结果)

两个位置用于测试每个处理(即 4 个处理,每个处理 2 个位置,总共 8 个位置)。

我想在R使用上述变量时进行拆分图(/嵌套?)重复测量(/混合模型?)方差分析。

问题 1。这听起来合适吗?

我的目标是查看是否存在 1) 位置、2) 治疗、3) 时间和 4) 所有相互作用(即 pos*treat*、pos*time、treat*time、pos*treat*time)的影响矿化率。

Q 2. 位置是否嵌套在治疗中?样品是否嵌套到位?

Q 3. 什么是中间因素和内部因素?

Q 4. 主题/情节是什么?- 是位置还是位置或样本或比率?

Q 5. 如何在我的 R 公式中将时​​间作为重复度量?

Q 6. 你会使用 aov、lme 还是 ezANOVA?

Q 7. 如何将单独的自变量以及它们的相互作用编码到适当的 R 公式中?

几天来我一直试图弄清楚这一点,但我似乎找不到一个有意义的答案......

1个回答

棘手的问题!位置是固定的还是随机的?位置是固定的还是随机的?我假设样本是随机的。

  • 由于处理被分配到位置,因此位置是抽样单位​​。基本上,治疗之间的比较是在那个水平上进行的。n=8.
  • 测量单位是您在给定时间对“样本”进行的观察。
  • 位置不嵌套在处理中。治疗适用于该位置。
  • 位置嵌套在位置内。
  • 样本嵌套在内部位置。
  • 时间嵌套在 Sample 中。
  • 时间与治疗相交。

您有 3 级嵌套(样本内的时间、位置内的样本、位置内的位置)。

如果位置、位置和样本是随机的,我认为 R 公式将如下所示:

 Y ~ Treatment * Time +(1|location|position|sample)

每次对于每个样本观察,您的数据框中都有 1 行 - 为您的所有设计特征提供适当的编码。

将重复的测量值组合成一个分数,比如它们的平均值或它们的差值是否有效?这可以使模型更容易解释。