R: как указать доверительные интервалы в метапропе пакета мета

Question

library(meta)
event <- c(81, 15, 0, 1)
n <- c(263, 148, 20, 29)
#
m1 <- metaprop(event, n, sm="PLOGIT", method.ci="SA")

Мне интересно объединить пропорции в метаанализе.В приведенном выше примере у меня есть 4 исследования, каждое из которых сообщает пропорцию.И похоже, что metaprop рассчитывает CI для каждого из 4 исследований.Однако, поскольку в 4 исследованиях уже сообщается о КИ (в оригинальной статье), есть ли способ включить фактические, сообщенные КИ в расчеты мета-анализа, вместо того, чтобы metaprop вычислять их самостоятельно?

Я также открыт для изучения других пакетов, если они позволят мне указать CI.

Answer 1

Во-первых, имейте в виду, что вы можете совершить огромную ошибку.Я не вижу причин не брать CI из функции metaprop.Вы должны взять рассчитанные CI из пакета, особенно если вы используете модель случайных эффектов.Средние индивидуальные эффекты и CI будут рассчитываться иначе, чем фиксированный эффект.В любом случае, вот оно:

Нижние пределы доверия можно найти в:

> m1$lower
[1] 0.25219013 0.05273001 0.00000000 0.00000000

И верхние пределы доверия в:

> m1$upper
[1] 0.3637795 0.1499727 0.0000000 0.1008922

Вы также можете получить доступиндивидуальные средние эффекты в TE:

> m1$TE
[1] -0.8095575 -2.1822989 -3.7135721 -3.3322045

Затем вы можете вручную добавить свои собственные вычисленные значения в эти векторы, и они появятся на вашем лесном графике.

Допустим, вашсобственные нижние значения myLCL:

myLCL <- seq(-1, -.25, .25) # define lower confidence limits

m1$lower <- myLCL # assign them to the respective location
m1$lower
[1] -1.00 -0.75 -0.50 -0.25

А ваши верхние доверительные пределы myUCL:

myUCL <- seq(.25, 1, .25) # define upper confidence limits

m1$upper <- myUCL # assign them to the respective location
m1$upper
[1] 0.25 0.50 0.75 1.00