r - 事後カスタム対比 - poly() を使用した glmmTMB

Question

glmmTMB多項式予測子を使用してモデルの一連のカスタム コントラストを実行する方法を理解する必要があります。モデルを操作できるようにするためにここで与えられた答えに従いましたが、まだ失敗しています.2つの段階にあるようです. 私の実際のモデルは非常に複雑で、間違いなくアプローチを台無しにするので、list() アプローチを使用する必要があります。glhtglmmTMBmatrixglmmTMB

library(contrast)
library(glmmTMB)

set.seed(1)
df <- structure(list(x = 1:20, y = c(4L, 7L, 11L, 12L, 23L, 21L, 42L, 
56L, 70L, 80L, 95L, 120L, 152L, 187L, 224L, 280L, 326L, 374L, 
438L, 500L), z = structure(c(2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 
1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L), .Label = c("a", 
"b"), class = "factor"), group = structure(c(1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L
), class = "factor", .Label = "1")), class = "data.frame", row.names = c(NA, 
-20L))

#just a toy example with a linear model and no poly()
m <- lm(y ~ x * z, data = df)
contrast(m, 
        a = list(x = 1, z = "a"),
        b = list(x = 10, z = "b"))

# a (big) simplification of what I'm after:
m1 <- glmmTMB(y ~ poly(x, 3) * z + (1|group), data = df, family = poisson)
contrast(m1, 
        a = list(x = 1, z = "a"),
        b = list(x = 10, z = "b"))

編集 1つの追加 - 私の実際のモデルはゼロインフレであるため、このアプローチはそれをどのように処理しますか? 「応答」値 (つまり、ゼロインフレとカウントの組み合わせ) をテストすることにのみ関心があります。