Is there an R package/function for performing two-stage residual inclusion or two-stage predictor substitution?

Question

I have to ultimately estimate patient health index (continuous variable bound between 0 and 1) using hospital record data. Usually patients that have poor health are also the on