Pacote Meta e Metafor R -

Question

Atualmente estou realizando uma metaanálise em R usando o pacote "metafor". Fazendo minha pesquisa, encontrei um pacote diferente para metanálises em R, chamado "meta". Eu gosto mais do gráfico de floresta criado pelo último pacote (em termos de design), mas infelizmente alguns dos dados não são os mesmos do gráfico que criei com metafor.

Especificamente, os dados são diferentes apenas para I^2 e a estimativa agrupada.

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)

Alguém sabe por que essas diferenças surgem?

Answer 1

Portanto, consegui que o intervalo de previsão correspondesse às funções, mas não aos valores I ^ 2 (mesmo que a diferença seja de apenas 2%). Pode haver alguma correção estatística que um pacote está fazendo em comparação com o outro ou tem a ver com o tipo de abordagem de modelagem RE/FE.

De qualquer forma, espero que este código ajude a apontar na direção certa. Para obter os ICs correspondentes, você também deve usar o parâmetro method.tau.ciem 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)