给定这段代码作为另一个问题的答案:
def poisson_interval(k, alpha=0.05):
"""
uses chisquared info to get the poisson interval. Uses scipy.stats
(imports in function).
"""
from scipy.stats import chi2
a = alpha
low, high = (chi2.ppf(a/2, 2*k) / 2, chi2.ppf(1-a/2, 2*k + 2) / 2)
if k == 0:
low = 0.0
return low, high
这个片段返回了一个双侧置信区间,但是如果我希望它是单侧的,我该怎么做呢。这更为复杂,因为泊松分布是不对称的。如有任何帮助,我们将不胜感激。
我认为应该将a/2
更改为a
或0
,因为它指示间隔所在的位置:
def poisson_interval(k, alpha=0.05):
"""
uses chisquared info to get the poisson interval. Uses scipy.stats
(imports in function).
"""
from scipy.stats import chi2
a = alpha
low, high = (chi2.ppf(0, 2*k) / 2, chi2.ppf(1-a, 2*k + 2) / 2)
if k == 0:
low = 0.0
return low, high
告诉我它是否有效。