Section 10 Box 8: Update the initial outcome model prediction based on targeted adjustment of the initial predictions using the PS model

ObsData$Pred.updateQ1.Y1 <- plogis(qlogis(ObsData$Pred.initialQ0.Y1) +  
                                   epsilon["H.A1L"]*ObsData$H.A1L)
ObsData$Pred.updateQ1.Y0 <- plogis(qlogis(ObsData$Pred.initialQ0.Y0) + 
                                   epsilon["H.A0L"]*ObsData$H.A0L)
summary(ObsData$Pred.updateQ1.Y1)

##        V1         
##  Min.   :0.00100  
##  1st Qu.:0.04354  
##  Median :0.05579  
##  Mean   :0.05441  
##  3rd Qu.:0.06610  
##  Max.   :0.11170

boxplot(ObsData$Pred.updateQ1.Y1)

summary(ObsData$Pred.updateQ1.Y0)  

##        V1         
##  Min.   :0.00100  
##  1st Qu.:0.03691  
##  Median :0.04874  
##  Mean   :0.04744  
##  3rd Qu.:0.05853  
##  Max.   :0.10932

boxplot(ObsData$Pred.updateQ1.Y0)