在R中,我很难获得由poisson.test
函数计算的速率比的精确置信区间。此函数使用binom.test
中的方法计算两个速率之间的置信限,例如
poisson.test(x = c(10000,20000), T = c(15000,15000), r = 1, conf.level = 0.95)
当x
(事件)低于T
(观察)时,此方法正常。然而,在非常高的x
(事件)和相对较低的T
(观测)中,泊松分布的方差将接近无穷大。 poisson.test
函数不考虑这一点,因为binom.test
方法没有这样的限制。因此,随着事件发生率逐步提高,个别利率和这些利率的比率的置信区间逐渐变窄,而它们应逐渐扩大。
是否有人知道使用均匀最强大的方法测试两个非常高的比率的替代方法,并获得正确的置信限制?