[python skill]利用python计算T分布下的置信区间
上篇博文中的置信区间计算代码在使用过程中并不准确,本人没并没有搞清楚原因 - - 求大神解答: import numpy as np from scipy import stats X1=np.array([14.65,14.95,8.49,9.51,10.23,2.75]) Xmean=X1.mean() Xstd=X1.std(ddof=1) interval=stats.t.interval(0.95,len(X1-1),Xmean,Xstd) print("置信区间为:",interval) 在计算一个文献中出现的数据时出现了问题: ↑为python代码输出结果 ↑为R输出结果 ↑正确结果。 并没有找到其他现成的包里面的命令用于直接计算T分布下的置信区间,通过查阅资料找到自定义函数实现置信区间的计算: def ci_t (data,confidence=0.95): sample_mean = np.mean(data) sample_std = np.std(data,ddof=1) sample_size = len(data) alpha = 1 - confidence...