线性混合模型 - 矩阵排名不足

Question

我正在尝试将线性混合模型应用于我的数据集，并且我不断收到一条奇怪的消息：

"fixed-effect model matrix is rank deficient so dropping 1 column / coefficient"

现在我的型号代码是：

 m0 <- lmer(y ~ var1 + var2 + var3 + (1|var4))

其中y是数字响应变量，变量1-4是因子。

这可能意味着什么？

Answer 1

我的初始运行遇到了同样的问题所以我检查了我所包含的变量，其中一个（因子变量）只有1个级别。所以我删除它，它工作。还要确保变量中没有NA。

检查NA：

sapply(data[, c(variable_names), with=FALSE], function(x) sum(is.na(x)))

放弃NA：

data_no_NAs <- na.omit(data[, c(variable_names), with = FALSE])