Question

我不确定这是否可行，但我想知道你是否可以从绘图（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个图形：

E.g：

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））输出

1 个答案: