Блог

Я использую пакет metafor для проведения мета-анализа в R, но у меня возникли проблемы с интерпретацией вывода функции fsn. В целом, «Наблюдаемый уровень значимости», по-видимому, не соответствует фактическому уровню значимости из предоставленной модели. Пример игрушки приведен ниже.

 #toy example originally from https://stats.stackexchange.com/questions/7766/meta-analysis-in- 
r-using-metafor-package
df <- structure(list(study = structure(c(1L, 5L, 3L, 4L, 2L), 
                                       .Label = c("Foo2000", "Pete2008",
                                                  "Pric2005", "Rota2008", 
                                                  "Sun2003"), 
                                       class = "factor"), 
                     mean1 = c(0.78, 0.74, 0.75, 0.62, 0.68), 
                     sd1 = c(0.05, 0.08, 0.12, 0.05, 0.03), 
                     n1 = c(20L, 30L, 20L, 24L, 10L), 
                     mean2 = c(0.82, 0.72, 0.74, 0.66, 0.68), 
                     sd2 = c(0.07, 0.05, 0.09, 0.03, 0.02), 
                     n2 = c(25L, 19L, 29L, 24L, 10L)), 
                .Names = c("study", "mean1", "sd1", "n1", "mean2", "sd2", "n2"),
                class = "data.frame", 
                row.names = c(NA, -5L)) 

library(metafor)
a <- rma(measure = "SMD", m1i = mean1, m2i = mean2, 
         sd1i = sd1, sd2i = sd2, n1i = n1, n2i = n2, 
         method = "REML", data = df)
a
fsn(a$yi, a$vi)
 

Почему наблюдаемый уровень значимости (.04) не соответствует фактическому уровню значимости (.3061) из модели? И почему я получаю отказоустойчивый N >0 (хотя это всего лишь 1), когда результаты модели уже не являются значимыми?

Спасибо за любую помощь. Я чувствую, что в этих случаях я получал «соответствующий» fsn (0), но некоторое время не использовал пакет.

fsn() (с type="Rosenthal" методом по умолчанию) использует «метод Стауффера для объединения значений p и описан в Rosenthal (1979)» (см. help(fsn) ). Это не то, что rma() нужно. Таким образом, значения p для этих двух методов отличаются. И поскольку значение p из метода Стауффера значительно (едва ли), N > 0.