線性混合模型 - 矩陣的秩不足

我試圖線性混合模型應用到我的數據集，我不斷收到一個奇怪的消息：線性混合模型 - 矩陣的秩不足

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

現在我的模型代碼是：

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

其中y是數字響應變量，變量1-4是因子。

這可能意味着什麼？

您在您的數據NA值 –

這意味着'var1'，'var2'和'var3'（固定效應）是多重共線性的，即它們的某些線性組合是不變的。 –

http://stats.stackexchange.com/questions/35071/what-is-rank-deficiency-and-how-to-deal-with-it –

我與我的初始運行有同樣的問題，所以我檢查了我包含的變量，其中一個（因子變量）只有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])

2016-08-12 14:39:26

回答