如何在R中计算poisson.test的标准误差？

Question

我有这样的数据集，

ID   Freq.x Freq.y
 1      1    8
 2      5    3
        ...

我使用R编程语言poisson.test函数计算了Freq.x和Freq.y两个速率参数之间的比率，但我想计算标准误差。我该怎么办？

Answer 1

您的问题中没有可重复的数据，所以让我们做一些：

set.seed(69)

x <- rpois(100, lambda = 7)
y <- rpois(100, lambda = 8)

您可以像这样获得这两个变量的标准误差：

se_x <- sqrt(mean(x)) / length(x)
se_y <- sqrt(mean(y)) / length(y)

se_x
#> [1] 0.02638181

se_y
#> [1] 0.02840775

，您可以将两者进行比较，以确定基本利率是否存在显着差异，如下所示：

poisson.test(c(sum(x), sum(y)))
#> 
#>  Comparison of Poisson rates
#> 
#> data:  c(sum(x), sum(y)) time base: 1
#> count1 = 696, expected count1 = 751.5, p-value = 0.004533
#> alternative hypothesis: true rate ratio is not equal to 1
#> 95 percent confidence interval:
#>  0.7781748 0.9556714
#> sample estimates:
#> rate ratio 
#>  0.8624535

尽管不清楚poisson.test的标准错误是什么意思。