我正在尝试通过使用 R 中的 miceadds 包从 SPSS 创建的多重插补数据库中汇集结果(我对 R 很陌生,如果术语有点偏离,我很抱歉)。当我尝试使用 as.mids2 代码将数据库更改为 mids 对象时,我收到以下错误消息:
Error in 1:(max(as.numeric(levels(data2[, .imp])))) :
result would be too long a vector
In addition: Warning messages:
1: In max(as.numeric(levels(data2[, .imp]))) :
no non-missing arguments to max; returning -Inf
2: In max(as.numeric(levels(data2[, .imp]))) :
no non-missing arguments to max; returning -Inf
有人知道这意味着什么吗?数据库是 28 个变量,有 2378 个 obs,所以不大。
这是我正在使用的脚本:
mydata <- read.csv("mydata.csv", na.strings = "999")
as.mids2 <- function(data2, .imp=1, .id=2){
ini <- mice(data2[data2[, .imp] == 0, - c(.imp, .id)], m = max(as.numeric(levels(data2[, .imp]))), maxit=0)
names <- names(ini$imp)
if (!is.null(.id)){
rownames(ini$data) <- data2[data2[, .imp] == 0, .id]
}
for (i in 1:length(names)){
for(m in 1:(max(as.numeric(levels(data2[, .imp]))))){
if(!is.null(ini$imp[[i]])){
indic <- data2[, .imp] == m & is.na(data2[data2[, .imp]==0, names[i]])
ini$imp[[names[i]]][m] <- data2[indic, names[i]]
}
}
}
return(ini)
}
mydata.mids <- as.mids2(mydata)
任何帮助将不胜感激。