Pacchetto Meta e Metafor R -

Question

Attualmente sto conducendo una metaanalisi in R utilizzando il pacchetto "metafor". Facendo le mie ricerche mi sono imbattuto in un pacchetto diverso per le metaanalisi in R, vale a dire "meta". Mi piace di più la trama della foresta creata da quest'ultimo pacchetto (dal punto di vista del design) ma sfortunatamente alcuni dei dati non sono gli stessi della trama che ho creato con metafor.

Nello specifico, i dati sono diversi solo per I^2 e per la stima aggregata.

meta_1 <- rma(yi=yi, vi=vi, measure="SMD", method="ML", slab=Citation, data=dat)
forest(meta_1)

meta_2 <- metagen(yi,vi^.5,data = dat,studlab = paste(Citation), comb.fixed = FALSE, 
                 comb.random = TRUE, hakn = TRUE, method.tau = "ML", sm = "SMD")
forest(meta_2)

Qualcuno sa perché emergono queste differenze?

Answer 1

Quindi sono stato in grado di far corrispondere l'intervallo di previsione tra le funzioni ma non i valori I ^ 2 (anche se la differenza è ridotta solo del 2%). Potrebbe esserci qualche correzione statistica che un pacchetto sta facendo rispetto all'altro o ha a che fare con il tipo di approccio di modellazione RE/FE.

Comunque spero che questo codice ti aiuti a indirizzarti nella giusta direzione. Per far corrispondere gli elementi della configurazione devi anche utilizzare il parametro method.tau.ciin metagen().

library(meta)
library(metafor)
study<- c(1:10) 
yi<- c( -0.48965031,0.64970214, 0.11201680,0.07945655,-0.70874645 -0.54922759,0.66768916 , -0.45523574 ) 
vi <- c(0.10299697,0.14036855,0.05137812, 0.03255550, 0.34913525, 0.34971466, 0.07539957, 0.08428983) 
dat <- cbind(study, yi, vi) 
dat <- as.data.frame(dat)

meta_1 <- rma(yi=dat$yi, vi=dat$vi, measure="SMD", method="REML", slab=paste(study), data=dat) 
forest(meta_1) 

meta_2 <- meta::metagen(TE =dat$yi,seTE = dat$vi^.5, method.tau = 'REML',
                        method.tau.ci = 'BJ', comb.random = TRUE, comb.fixed = TRUE,
                        sm = 'SMD') 

forest(meta_2)