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)