使用 meta.summaries 来自 epi.2by2 的置信区间

Question

我正在寻找基于一系列值的平均值，每个值都有相关的错误。我发现 rmeta 包中的 meta.summaries 似乎工作得很好，但我对输入值的错误只是上限和下限，而不是标准错误测量，来自使用epi.2by2 个。 Meta.summaries 需要一个标准错误。

有没有办法让meta.summaries使用上限和下限，或者epi.2by2生成一个s.e？

attach(DATA)
m <- summary(meta.summaries(RR,SE), conf.level = 0.95)

Answer 1

也许你太依赖包裹了？从维基百科的Odd-Ratio page，我可以读到：

The distribution of the log odds ratio is approximately normal with:

L ~ N( log(OR), s^2 )

... 其中 s 可以近似为：

SE=sqrt(1 / n_11 + 1/n_01 + 1/n_00 + 1/n_10)

（其中 n_11 是接触药物的治愈人数，n_10 未接触药物的治愈人数，等等）因此您可以非常轻松地计算标准误差。假设您将所有观察结果打包为一个数据集：

n      <- 10                                       # number of odd-ratios
counts <- sample(20:100, replace = TRUE, size=4*n) # fake counts data
counts <- matrix(counts, cols=4)                   # one line = one experiment = one OR = 4

...然后您可以按如下方式计算标准误差：

se <- sqrt(rowSums(1/counts))

请注意：

This is an asymptotic approximation, and will not give a meaningful result if any of the cell counts are very small.

...更一般地说，平均奇比听起来很可疑，至少我也是这样。如果您有原始观察结果，那么您可以简单地将所有计数放在一起。如果你不这样做，那么你将在这个过程中超重小型研究。我没有详细调查，但似乎对此事有很多讨论：https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5352535/. You can also see what is said here 在 CrossValidated（用于统计的 Whosebug）上。

使用 meta.summaries 来自 epi.2by2 的置信区间

Using meta.summaries for confidence intervals from epi.2by2

meta

statistics

r