
library(meta)

N <- rep(100,7)
Mean <- c(1,2,4,8,16,32,64)
SD <- c(0.5,1,2,4,8,16,32)

res1 <- metamean(n = N,
                 mean = Mean,
                 sd = SD,
                 sm = "MRAW",
                 comb.fixed = FALSE,
                 comb.random = TRUE,
                 method.tau = "REML",
                 hakn = TRUE)
res1
forest(res1)

res2 <- metamean(n = N,
                 mean = Mean,
                 sd = SD,
                 sm = "MLN",
                 comb.fixed = FALSE,
                 comb.random = TRUE,
                 method.tau = "REML",
                 hakn = TRUE)
res2
forest(res2)
forest(update(res2, backtransf = FALSE))

res1$tau
res2$tau

res1$Q
res2$Q
