Confidence Interval for the mean of a Normal distribution using scipy and Pythonn

confidence_interval.py

from scipy import stats
import numpy as np

def mean_confidence_interval(x, alpha=0.05):
    """Computes two-sided confidence interval for a Normal mean
    
    Assumes population variance is unknown.
    
    x is assumed to be a list or a 1-d Numpy array
    """
    n = len(x)
    xbar = np.mean(x)
    standard_error=stats.sem(x)
    return stats.t.interval(1 - alpha,
                            n - 1,
                            loc=xbar,
                            scale=standard_error)