使用小鼠 package 進行多重插補

Question

我不是 R 用戶，但我使用 R 僅使用小鼠 package 朗姆酒多重插補。 由於我將 R 從版本 3.6.3 更新到版本 4.0.4，當我從 R 調用它到 STATA 時，我遇到了估算數據的問題。 問題是當我進行生存分析時，所有受試者都被視為死亡（狀態=1），而狀態變量包括（活着和死亡）。 請檢查我的 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 時。 需要記錄和重新標記一些分類變量。