R中的MICE包:被动归罪

Question

我的目标是用多重填充来处理缺失值，然后用混合线性模型进行分析。

我被"BMI“(身体质量指数)和"BMI类别”的被动指责所堆积。根据身高和体重计算"BMI“，并将其归入"BMI类别”。

如何计算“BMI类别”？

数据库如下所示：

 sub_eu_surf[1:5, 3:12]

    age gender smoking exercise education sbp dbp height weight      bmi
    1  41      1       1        2        18 120  80    185    107 31.26370
    2  46      1       3        2        18 130  70    182    102 30.79338
    3  46      1       3        2        18 130  70    182    102 30.79338
    4  47      1       1        2        14 130  80    178     78 24.61810
    5  47      1       1        1        14 150  80    175     85 27.75510

由于'bmi类别‘不能预测我的估算，所以我决定在估算之后创建它。详细信息如下：

1. To define

predictor

ini<-mice(sub_eu_surf, maxit=0)
meth<-ini$meth
meth["bmi"]<-"~I(weight/(height/100)^2)"

pred <- ini$predictorMatrix
pred[c("pm25_global", "pm25_eu", "pm10_eu", "no2_eu"),  ]<-0 
pred[,c("bmi", "hba1c", "pm25_eu", "pm10_eu")]<-0
pred[,"tc"]<-0
pred[c("smoking", "exercise", "hdl", "glucose"), "tc"]<-1
pred[c("smoking", "exercise", "hdl", "glucose"), "ldl"]<-0
vis <- ini$vis
imp_eu<-mice(sub_eu_surf, meth=meth, pred=pred, vis=vis, seed=200, print=F, m=5, maxit=5)
long_eu<- complete(imp_eu, "long", include=TRUE)
long_eu$bmi_category<-cut(as.numeric(long_eu$bmi), breaks=c(0, 18.5, 25, 30, 72))
complete_eu<-as.mids(long_eu)

但是我在分析我的数据时收到了一个错误：

    test1<-with(imp_eu, lme(sbp~pm25_global+gender+age+education+bmi_category, random=~1|centre))
    Error in eval(expr, envir, enclos) : object 'bmi_category' not found

这是怎么发生的？

Answer 1

您是在原始mids object imp_eu上运行分析，而不是在修改后的complete_eu上运行分析。尝试：

test1<-with(complete_eu, lme(sbp~pm25_global+gender+age+education+bmi_category, random=~1|centre))