Question

我想在blme :: bglmer（）中使用一个固定效果，作为我的一个协变量显示完全分离的变量，遵循上一篇文章：Binomial glmm with a categorical variable with full successes的指导原则。但是，就我而言，还有一个协变量，没有显示出完全分离。

以下是我的数据集（RAMP_Subset）的概述。 ID标识观察到的每条鱼，拖曳标识每条采样的拖网（在渔船上，同一条拖网有多条鱼），死亡是我的响应变量，具有2个级别（渐近线处的死亡率，存活率0和死亡率1），R1 R2是我的协变量，每个变量有2个级别（每条鱼的反射率，存在反射时为0，没有反射时为1）：

ID   haul death R1 R2
2704    9     1  1  1
2705   14     1  1  1
2715   13     1  1  1
2718   13     1  1  1
2719    8     1  1  1
2722    9     1  1  1
...

让我们说R1是完全分离的协变量，我想为此添加一个固定的影响。因此，我的模型如下：

bglmer(death ~ R1 + (1|haul) + (1|ID), data=RAMP_Subset, family=binomial, fixef.prior = normal(cov = diag(9,2)))

现在让我们说我的第二个协变量R2没有显示出完全分离。 reference manual for the R package "blme"指出，固定效果的先验以与协方差先验相似的方式指定，因此我尝试了以下操作：

bglmer(death ~ R1 + R2 + (1|haul) + (1|ID), data=RAMP_Subset, family=binomial, fixef.prior = R1 ~ normal(cov = diag(9,2)))

但是，我的尝试收到一条错误消息：

valuateFixefPrior（fixefPrior，defnEnv，evalEnv）中的错误：无法识别的Fixef分布：“〜”

我在网上找不到任何其他示例，非常感谢您提供有关正确格式的建议。

谢谢，以斯帖

多个固定效果时，在blme :: bglmer（）中使用固定效果的格式优先

0 个答案: