def CI(q, data, method = "mean",U = True, L = True):
    if q > 1: q = q/100
    norms = qnorm(q+(1-q)/2)*(std(data) / sqrt(len(data)) )
    if method == "mean": u = mean(data) + norms; l = mean(data) - norms
    if U == L == True: return (l,u)
    elif U == True: return u
    else: return l