1

我正在尝试运行一个由 2 个固定效应变量组成的混合效应模型，其中第一个固定效应有两个级别，而另一个代表连续数据。连续数据是 pH 值，我有 12 个与我的处理块的不同子图相关联的唯一值。因为我使用的是连续数据，所以使用 cld(emmeans()) 框架只能给我 12 个唯一值的平均 pH 值，所以我有兴趣使用 emtrends() 来确定我的响应之间关系的斜率变量和 pH 值与 0 显着不同。

我已经尝试了后面的代码和这个评论。当我运行我的模型时，我们发现 pH 值对响应变量有显着影响。但是，当我尝试确定回归线的斜率是否与零显着不同时，我收到以下错误：“emm_basis.merMod(model, attr(data, "terms"), RG@model.info$ xlev，：缺少参数“选项”，没有默认值”

数据框示例：

site   species     nitrogen_treat    pH        response_var
   A         B          nit_added    3.4               36.8
   A         B          nit_added    3.4               35.7
   B         A        no_nitrogen    5.6               32.1
   A         C        no_nitrogen    5.6               30.1
   B         D        no_nitrogen    4.3               30.2
   C         C          nit_added    3.1               37.2
   C         A          nit_added    7.2               10.2

所以我的 lmer 函数代码如下：

lm <- lmer(response_var ~ nitrogen_treat * pH + 
               (1 | site) + (1 | species))

然后我用它来尝试 emtrends() 函数，该函数不起作用：

emtrends(lm, ~1, var = "pH")

我还尝试将“~1”与“~nitrogen_treat”交换，但我不断收到我上面提到的错误消息。我将不胜感激任何有关我收到这些错误消息的原因的见解！

4

1 回答 1

2

这是我在 github 上的当前工作版本中发现并修复的一个新错误。但是，您可以通过添加调用options = list()来解决它。emtrends()

于 2019-10-09T03:22:50.317 回答