r - glm の事後分析の実行中にエラーが発生しました

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