Я использую пакет lsmeans/emmeans
в R для создания графика парных сравнений в ответе между уровнями TreatA (двоичная / факторная переменная). Я могу получить оценки разницы, используя lsmeans (контраст), но он предоставляет только SE для оценок, а не доверительные интервалы. Доверительные пределы предоставляются только для отдельных эффектов, а не контрастов. Кто-нибудь может помочь сгенерировать CL для средних различий (контрастов)?
model = glmer.nb(response ~ treatA * treatB + (1|random), data, family=nbinom1)
lsm <- lsmeans(model, pairwise ~ treatA*treatB)
contrast(lsm, method = "pairwise")
В конечном счете, я надеюсь создать такой сюжет, но с CL вместо SE.