ObsData$Pred.TE <- ObsData$Pred.initialQ0.Y1 - ObsData$Pred.initialQ0.Y0 mean(ObsData$Pred.TE)
## [1] 0.007043316
hist(ObsData$Pred.TE)