随手写程序-t检验计算置信区间
事实上,对于一个单一的样本和由它得来的单置信区间,我们无法得出总体均值落在这个区间内的概率是多大,下面的例子说明了这一事实。
如图所示:
深度讲解:下面使用t检验来计算置信区间。
有一个内置的函数scipy.stats可以计算置信区间,记住要指定自由度!
对于正态分布,还有一个内置的函数可以计算置信区间,这个函数不需要指定自由度。
下面是一些可视化图形上的置信区间的代码:
结果如下:
标准偏差,标准误差和置信区间的计算都依赖于一定的假设。如果违反这些假设,那么95%的置信区间的可信度将会降低。我们说,在这种情况下的置信区间是非校准的。下面是一个例子。
执行结果为:
事实证明,对于较大的样本量,我们应该看到样本均值渐近收敛到零。
结果如下:
从图中可以很明显的看出结果慢慢收敛,我们也可以查看一下样本均值的均值。