【问题标题】:Multiple imputation using mice package使用鼠标包进行多重插补
【发布时间】:2021-03-29 13:54:37
【问题描述】:

我不是 R 用户,但我仅使用 R 来使用 mouse 包进行多重插补。由于我将 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")```

【问题讨论】:

  • 样本数据丢失

标签: r r-mice survival


【解决方案1】:

语法是正确的,问题出在将估算数据从 R 调用到 STATA 时。需要记录和重新标记一些分类变量。

【讨论】:

    猜你喜欢
    • 2017-10-15
    • 1970-01-01
    • 2021-09-28
    • 2019-08-04
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多