从 plot(allEffects(model)) 输出中抑制一些图

Question

我不确定这是否可行，但我想知道您是否可以从 plot(allEffects(model)) 输出中删除一些不那么有趣的图。我认为这可以通过编辑 R 用来制作图形数组的一些背景机制来实现。我已经使用多线将图形数量从 20 个减少，但我只对剩余 10 个图中的 2 个感兴趣，因为它最准确地代表了我的模型（广义线性混合模型）的交互作用之一。

示例代码：

glmer1 <- glmer(dv ~ categoricalv1*dimensionalv1*categoricalv2*dimensionalv2 + (1 | subjectID), data = data1, family = binomial, contrasts=contrasts1, control = glmerControl(optimizer = "bobyqa"))

当我查看所有效果函数中的值时：

allEffects(model)
model: dv ~ categoricalv1 * dimensionalv2 * categoricalv2 * dimensionalv2

R 清楚地从输出中读取值以生成 10 个图表：

例如：

catv1*dimv1*catv2*dimv2 effect, , catv2 = 0, dimv2 = -4 ... 
           dimv1
catv1             -2      -0.7       0.5         2         3
dv level1        0.269329 0.4296272 0.5929832 0.7686877 0.8520351
dv level2        0.730671 0.5703728 0.4070168 0.2313123 0.1479649
... (x9 more times to show all the different combinations)

情节代码（到目前为止）：

 plot(allEffects(model), ylab="Proportion", xlab="Standardized Dimensional V1", main="",
 type="rescale", multiline=T, rescale.axis=F, ci.style="band")

希望有一种方法可以编辑某些内容以仅显示我感兴趣的情节。任何见解将不胜感激:)

Answer 1

使用函数selection=2,

plot(allEffects(model), ylab="Proportion", **selection=2**, xlab="Standardized Dimensional V1", main="", type="rescale", multiline=T, rescale.axis=F, ci.style="band")

for example;

model2<-plot(allEffects(model), selection=2)

model6<-plot(allEffects(model), selection=6)

如果你想合并不同的地块，你可以使用

library(gridExtra)

grid.arrange( model2, model6,  nrow=1,  ncol=2)

从 plot(allEffects(model)) 输出中抑制一些图

suppressing some of the plots from plot(allEffects(model)) output

plot

r

mixed-models