kuiper_false_positive_probability#

astropy.stats.kuiper_false_positive_probability(D, N)[源代码]#

计算柯伊伯统计量的假阳性概率。

使用Paltani 2004中描述的四个公式;他们报告得出的函数从未低估假阳性概率,但在N=40..50的范围内可能有点高。(他们引用了1e-7级别的系数1.5。)

参数:
D : floatPython :浮点

柯伊伯测验分数。

N : floatPython :浮点

有效样本量。

返回:
fpp : floatPython :浮点

由零假设产生的分数如此之大的概率。

笔记

Paltani 2004的公式7似乎错误地引用了原始公式(Stephens 1965)。此函数实现原始公式,因为它产生的结果更接近蒙特卡罗模拟。

工具书类

[1]

Paltani,S.,“使用柯伊伯检验在X射线观察中寻找周期。应用于ROSAT PSPC档案”,天文学和天体物理学,v.240,p.789-790,2004。

[2]

Stephens,M.A.,“拟合优度统计VN:分布和显著性点”,Biometrika,v.52,p.3091965。