Add confidence_intervals.R script