2 つのカテゴリ予測変数 (X1、X2) と 1 つの連続予測変数 (X3) を含む混合モデルがあります。
model <- lmer(z ~ x1 * x2 * x3 + (1|group), data = data)
連続予測子の高特性と低特性を比較することに興味があります。
私の計画は、X3 の単純な傾き (-1 SD、M、+1 SD) を対比することです。
私が理解している限り、これは次のように emmeans の emtrends() 関数を使用して実行できます。
sd1 <- mean(data$X3, na.rm = T) + sd(data$X3, na.rm = T)
mean <- mean(data$X3, na.rm = T)
sd2 <- mean(data$X3, na.rm = T) - sd(data$X3, na.rm = T)
mylist <- list(X3 = c(sd1, mean, sd2))
emtrends(model, ~ X1 * X2 | X3,
var = "X3",
at = mylist)
ただし、提供される係数は、X3 の 3 つの値に対して同じです。
誰でもこれに光を当てることができますか?