Question

我正在尝试对治疗进行事后比较，但在运行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")))

Answer 1

您已经提到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

在这里，我们以_______ |A: 15| |A: 45| |D: 55| |A: 45| |D: 65| (...)为条件并刻画所有成对比较的特征，并使用Tukey方法校正多种假设检验。

对glm运行事后分析时出错

1 个答案: