подавление некоторых графиков из вывода сюжета (allEffects (модель)) - PullRequest
Новая
       12

подавление некоторых графиков из вывода сюжета (allEffects (модель))

0 голосов
/ 24 мая 2018

Я не уверен, что это возможно, но мне интересно, можете ли вы удалить некоторые из не столь интересных графиков из вывода сюжета (allEffects (модель)).Я предполагаю, что это возможно, отредактировав некоторую фоновую механику, которую R использует для создания массива графиков.Я уже использовал многострочность, чтобы уменьшить количество графиков с 20, но меня интересуют только 2 из 10 моих оставшихся графиков, потому что они наиболее точно представляют один из моих эффектов взаимодействия из моей модели (обобщенная линейная смешанная модель).

Пример кода: 

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")

Надеемся, что есть способ отредактировать что-то, чтобы показать только сюжет Iменя интересует. Любые идеи будут оценены :)

1 Ответ

0 голосов
/ 25 сентября 2018

Используйте функцию 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)
...