"场所" visreg中有多个预测因子的论证

时间:2017-05-06 18:01:50

标签: r plot data-visualization interaction mixed-models

我试图从使用包glmer生成的一组二项lme4个对象中绘制预测概率。它们的形式如下:

library(lme4)
model = glmer(
  LogFlag ~ FRI + Mound + woody + (1 | ID) + 
    (0 + FRI | ID) + (0 + Mound | ID) +
    (0 + woody | ID)  + Mound:FRI + Mound:woody, 
  family = "binomial", 
  data = BBmod
)

要制作我正在使用visreg调用的图表,如下所示:

visreg(
  model, xvar = "FRI", by = "Mound", scale = "response", 
  rug = 0, overlay = T, partial = F, gg = T, 
  breaks = c(2, 4, 6, 8, 10, 15, 20, 40, 80), 
  band = F, line = list(col = rainbow(9)), legend = F
)

这通常会产生非常漂亮的图,x轴上有FRI,y轴上的概率,以及breaks参数中指定的Mound值的9条单独的线。

我的问题是第三个预测因子woody的使用价值是多少?如果没有任何breaks参数,visreg默认使用未包含在调用中的预测变量的中位数,但我已为breaks提供了Mound。我如何指定或至少确认正在使用的木质的价值?

1 个答案:

答案 0 :(得分:1)

我混淆了breakscond个参数。文档显示了

cond = list(woody = Z)

我可以控制woody的值,其中Zwoody的值。也许这会对某人有所帮助。