def getConfidenceInterval(data, percent=0.95, distribution="t"):
n, min_max, mean, var, skew, kurt = stats.describe(data)
std = np.sqrt(var)
if distribution == "t":
R = stats.t.interval(percent, len(data) - 1, loc=mean,
scale=std / math.sqrt(len(data)))
else:
R = stats.norm.interval(
percent, loc=mean, scale=std / math.sqrt(len(data)))
error = (R[1] - R[0]) / 2
return mean, error
评论列表
文章目录