我有三个级别的固定效果,我使用以下对比在我的数据集中创建了 2 列
data = data %>%
mutate(contrast_B = ifelse(group == "A", -0.5,
ifelse(type == "B", 0.5, 0))) %>%
mutate(contrast_C = ifelse(group == "A", -0.5,
ifelse(type == "C", 0.5, 0)))
我将它们与其他固定效果一起包含在我的模型中,这些效果都是二进制的,具有如下两个级别。
model <-lmer(log(dependent variable) ~ contrast_B * contrast_C *
fixed effect * fixed effect * fixed effect +
(1|subject)+ (1|item), data = data....
但我收到了一条消息
固定效应模型矩阵秩不足,因此删除 12 列/系数"
有人知道如何解决吗?