docs: def mean_interval(data, alpha=_alpha):
      """
      Interval assuming gaussian posterior.
      """
      mean =np.mean(data)
      sigma = np.std(data)
      scale = scipy.stats.norm.ppf(1-alpha/2.)
      return interval(mean,mean-scale*sigma,mean+scale*sigma)