Ошибка при запуске posthoc анализа для glm - PullRequest
Новая
       83

Ошибка при запуске posthoc анализа для glm

1 голос
/ 22 марта 2019

Я пытаюсь выполнить сравнение posthoc для моих обработок, но я продолжаю получать эту ошибку при запуске glht: «Ошибка в modelparm.default (model, ...): размеры коэффициентов и ковариационная матрица не совпадают»,

Есть ли лучший способ сделать несколько парных сравнений?Я также пытался использовать emmeans, но я не уверен, что это правильный метод.

Это подмножество моих данных:

mydata <- read.table(header=TRUE, text="
treatment    total.bites   hours   rep
                     A  10  3.1  a1
                     A  1   3.2  a2
                     A  1024   3.22 a3
                     B  0   3.13 a1
                     B  16  3.15 a2
                     B  1305  3.24 a3
                     C  0   3.13 a1
                     C  0  3.26 a2
                     C  0   3.11 a3
                     D  2  3.25 a1
                     D  0   3.17 a2 
                     D  3   3.21 a3
                     ")
mC4 <- glmmTMB(total.bites~treatment + offset(log(hours)) +(1|rep), ziformula=~0, family=nbinom1, data=mydata)
summary(mC4)
summary(glht(mC4, mcp(treatment = "Tukey")))

1 Ответ

0 голосов
/ 22 марта 2019

Как вы уже упомянули emmeans, вы могли бы сделать 

library(emmeans)
pairs(emmeans(mC4, "treatment"))
#contrast estimate       SE df t.ratio p.value
#A - B       0.323 8.94e-01  6  0.361  0.9824
#A - C      20.930 1.51e+04  6  0.001  1.0000
#A - D       0.892 9.05e-01  6  0.985  0.7631
#B - C      20.607 1.51e+04  6  0.001  1.0000
#B - D       0.569 9.92e-01  6  0.573  0.9365
#C - D     -20.038 1.51e+04  6 -0.001  1.0000
#
#Results are given on the log (not the response) scale.
#P value adjustment: tukey method for comparing a family of 4 estimates

Здесь мы включаем treatment и характеризуем все парные сравнения, исправляя тестирование множественных гипотез с использованием метода Тьюки.

...