Я работаю над блеском, в котором переменная отклика является категориальной (результат: да, нет), а фиксированный коэффициент тоже категориален (sample.type: control, эксперимент).
здесь моймодель: # m1 <- glmer (результат ~ sample.type + (1 | dyad), семейство = бином, data = boxes_exp) </p>
#boxes_exp$sample.type:
[1] эксперимент эксперимент эксперимент эксперимент эксперимент эксперимент эксперимент контроль контроль контроль контроль....... (уровни: контроль, эксперимент)
#boxes_exp$outcome
[1] 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 1 0 1 0 0 0 0 0 0 1
1 1 1 1 0 1 1 1 1 0 1 1 1 1 0 1 0 0 0 0 0..... (levels: 0, 1)
Я пытался использовать эффект сюжета, висрег и мозаичный сюжет, мне было интересно, есть ли более подходящий способ представления таких данных
#plot(effect("sample.type", m1), x.var="sample.type",
main="outcome in relation to sample type",
multiline=TRUE,confint=TRUE,ci.style="none")
#Using package visreg:
visreg(m1, "sample.type", type = "conditional",main="Outcome in relation
to Sample type", xlab = "Sample type", ylab = "outcome")
#mosaicplot(table(boxes_exp$sample.type, boxes_exp$outcome))
Отображение в соответствии с результатами модели должно показывать, что результат 1 выше, чем результат 0 для эксперимента с типом образца, и наоборот для контроля типа образца.
Мне было интересноесли кто-то может помочь мне с более сложным способом представления результатов,
большое спасибо