使用mice软件包进行多重填充

Question

我不是R用户，但我只使用R来使用mice包进行多重补偿。由于我将R从版本3.6.3更新到版本4.0.4，因此当我从R调用到STATA时，我遇到了估算数据的问题。问题是，当我运行生存分析时，所有受试者都被视为已死(status=1)，而状态变量包括(alive和died)。如果你能查查我的R码，我将不胜感激。我已经将这个平台的变量命名为var1、var2、var3、var4。这是我的R脚本

library (MASS)
library (lattice)
library(splines)
library(survival)
library (Rcpp)
library (mice)

#reading data 
hla=read.dta("R:/mypath")

attach(dataname)
names(dataname)

#create data frame 
dataframe=data.frame(var1, var2, var3, var4, na) 

names(dataframe)

#generate automatic prediction matrix
predmatrix=quickpred(dataframe)
predmatrix

#ImputationSettings
maxIter=20
imputations=10

#Imputation settings
#set timer on
ptm=proc.time()      

imp10=mice(dataframe,m=imputations,maxit=maxIter,
method=c(
"",#var1
"pmm",#var2
"pmm",#var3
"logreg",#var4
""),#na

predictorMatrix=predmatrix,seed=100)

#extract original and imputed data in long form
ImputedData=complete(imp10,action="long",include=TRUE)
#export data to Stata file
write.dta(ImputedData,"R:/mypath/data.dta")```

Answer 1

语法是正确的，问题出在从R向STATA调用推算数据时。需要对一些分类变量进行记录和重新标记。