code: def mean_cl_normal(series, confidence_interval=0.95):
      
      a = np.asarray(series)
      m = np.mean(a)
      se = scipy.stats.sem(a)
      h = se * scipy.stats.t._ppf((1+confidence_interval)/2, len(a)-1)
      return pd.DataFrame({'y': [m],
                           'ymin': m-h,
                           'ymax': m+h})
docs: Adapted from http://stackoverflow.com/a/15034143