上篇博文中的置信区间计算代码在使用过程中并不准确,本人没并没有搞清楚原因 - - 求大神解答:
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分布下的