0

我想知道如何编写两个相互作用因素的随机效应。例如，我有 6 个物种，它们被种植在 48 个地块中，并在两个街区中进行了复制。6 个物种共有 48 种组合（所有可能的 1、2、3、5 和 6 种组合）。

通过设计，组合实际上嵌套在多样性中：组合是每个多样性级别的正确复制；因此，多样性的影响应该针对多样性水平内组合之间的（随机）变化进行测试。由于多样性或物种丰富度是一个连续变量，我不能认为组合嵌套在多样性中。

我想将块和组合之间的交互作为随机效应。我不确定如何在 R 的混合模型中正确写入。谢谢

4

1 回答 1

1

在标准的 R 混合模型表示法（回到 Doug Bates 的nlme包）中，将块和组合之间的交互视为随机效应(1|block:combination)，即响应变量的基线（截距）值根据跨 {块的高斯分布而变化× 组合} 值。

这应该适用于 R: nlme, lme4, brms, GLMMadaptive...中的大多数混合模型包

于 2022-03-03T19:58:57.327 回答