0

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 つの値に対して同じです。

誰でもこれに光を当てることができますか?

4

1 に答える 1