У меня есть следующие данные:
set.seed(3)
library(data.table)
library(lme4)
a <- rep(1:5, times = 20)
b <- rep(c(1,1,1,1,1,2,2,2,2,2), times = 50)
ppt <- rep(101:110, each = 10)
item <- rep(1:10, times = 10)
dv <- rnorm(n = 100)
contrasts(data$a) = contr.sum(4)
data <- data.table(cbind(ppt, item, a, b, dv))
data$ppt <- as.factor(data$ppt)
data$item <- as.factor(data$item)
data$a <- as.factor(data$a)
data$b <- as.factor(data$b)
Я хотел бы получить коэффициент для каждого уровня a. u / omsa_d00d и u / dead-Серьезные указали мне на идею запуска модели без перехвата.
Если я запускаю эту модель:
m1 <- lmer(dv ~ a + b -1 +(1|ppt) + (1|item), data = data)
Я получаю коэффициенты для каждого уровня a.
Однако, если я запускаю эту модель, в которой b стоит первым:
m2 <- lmer(dv ~ b + a -1 +(1|ppt) + (1|item), data = data)
Я получаю коэффициенты для каждого уровня b, но не a.
Что именно происходит в каждом случае?
Кроме того, выполняется ли m1 достаточно, чтобы получить эффект каждого уровня по сравнению со средним значением, и одновременно контролировать b ?
Имеет ли значение, если я имею в виду сначала центрировать мои предикторы?
Каковы различные значения коэффициента кодирования фиктивного и суммарного кодов a?