最近给导师报告处理WISE数据的进展时,对于处理WISE光变导师希望我对每一个观测区间的星等数值给中值的置信度。
计算数据的中值95%的置信区间
假设有一组数据如下:
data = [-0.1, -2.4, -0.1, -0.7, -1.4, -0.9, -3.2, -0.2, -0.3, -0.6, -3.2, -5.5]
求中值的置信区间与求数据的均值置信区间的方法是类似的,在这里我不会讨论详细的数学原理,而是直接给出置信度为95%时,对应的数据的上限值与下限值。(有关python实现的重点在于代码)
下限:0.5n-0.98
上限:1+0.5n+0.98
上面公式给出的是理论值,具体应用到数据上时,要对得到的下限和上限值取整,下限值向上取整,上限值向下取整。
一般情况:
Lower lim =
-
*
Upper lim = 1 +
+
*
注&#