我正在尝试使用lmer来安装模型,并且我不断收到错误:
Error : (p <- ncol(X)) == ncol(Y) is not TRUE
该模型如下:
model.covars.im.adapt <- lmer (data$im.adpt ~ data$pnr + data$female_c + data$age_c + data$age_c2 + (1 | data$nomem_encr) + (1 | data$nohouse_encr))
这个问题可能是由于NAs,因为以下模型,其中因变量没有NA,运行良好:
model.covars.rw <- lmer (data$rw ~ data$pnr + data$female_c + data$age_c + data$age_c2 + (1 | data$nomem_encr) + (1 | data$nohouse_encr))
然而,当我编码排除NA时,我不断收到同样的错误:
data <- data[-is.na(data$im.adpt), ]
model.covars.con <- lmer (data$im.adpt ~ data$pnr + data$female_c + data$age_c + data$age_c2 + (1 | data$nomem_encr) + (1 | data$nohouse_encr))
以下是数据摘要$ rw:
V1
Min. :-2.431773
1st Qu.:-0.574774
Median :-0.110525
Mean : 0.000063
3rd Qu.: 0.817975
Max. : 2.210723
数据汇总$ im.adpt:
V1
Min. :-2.75092
1st Qu.:-0.69883
Median : 0.32721
Mean :-0.00002
3rd Qu.: 0.32721
Max. : 1.35326
NA's :316
你推荐什么?我有能力提供更多相关信息。