Plan for adding 95% CI error bars to Kenny et al. original graphs

[yuexinz04]

On the foundation of adding a 95% CI band of the gam fit line to the original graphs, add a vertical line to each point of fitted value. The range of this vertical line is as follows:
{
B <- gam_fit_model()
Pred.b = Predict(B, se.fit=T) # predicted value and sd
Lowerbound= Pred.b$fit- 1.96pred.b$se.fit
Upperbound=Pred.b$fit+1.96pred.b$se.fit
(95 CI%)
}
Add-on to ggplot()
{
geom_errorbar(aes(ymin=fitted(m)-1.96se, ymax=fitted(m)+1.96se), width=.1)
}

@leahghazali @nami-jain @RuobinGong

[yuexinz04]

Complete

https://github.com/yuexinz04/2022-Summer-Statistical-Privacy-and-Differential-Privacy/tree/kennyetal_modification/Add_Range