0

我正在尝试通过使用 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)

任何帮助将不胜感激。

4

0 回答 0