python - python中的单边泊松置信区间

Question

鉴于此代码是在另一个问题中作为答案给出的:

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

此代码段返回一个双侧置信区间，但如果我想要单侧置信区间，我该怎么做。这更复杂，因为泊松分布是不对称的。任何帮助将不胜感激。

Answer 1

我认为您应该将 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

告诉我它是否有效。