R混合模型中两个因素间的随机相互作用



我想知道如何写出两个相互作用因素的随机效应。例如,我有6个物种,它们被种植在48个地块上,并在两个区块中复制。6个物种总共有48个组合(所有可能的1、2、3、5和6个物种组合(。

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

我想把块和组合之间的相互作用作为随机效应。我不知道如何在R中的混合模型中正确写入。谢谢

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

这应该适用于R:nlmelme4brmsGLMMadaptive中的大多数混合模型包。。。

最新更新