根据CHOPIT模型计算校正后的自报健康水平,除了估计系数外,还需要进行预测。

完整的代码流程如下:

建立CHOPIT模型,估计系数

fo <- list(
  self = Q04 ~ 1,
  vign = cbind(Q01, Q02, Q03) ~ 1,
  tau = ~ 1,
  cpolr = ~ 1
)

out1 <- chopit(fo, data=mydata)

在原数据上进行预测,生成校正后的自报健康水平

mydata$corrected_Q04 <- predict(out1, mydata)

predict函数将利用模型参数,对原始自报健康水平进行校正,输出校正后的健康水平指标。

查看校正效果

summary(mydata$corrected_Q04)  

这样就可以直接获得CHOPIT模型校正后的自报健康水平。注意predict时需要传入原始的数据。

另外,也可以用fitted()函数获得拟合值。但predict()更适合生成校正后的健康水平指标。