[R-meta] rma.mv: why some var components change but others don't across 2 models

Stefanou Revesz @te|@noureve@z @end|ng |rom gm@||@com
Fri Oct 29 17:24:11 CEST 2021


Dear Wolfgang and Expert List Members,

Why `study` with 57 levels in model `res` gives `sigma^2.1 = 0.0200`
but `study` with 57 levels in model `res2` gives `sigma^2.1  =
0.0122`?
(SAME LEVELS BUT DIFFERENT RESULTS)

Why `outcome` with 4 levels in model `res` gives `sigma^2.2 = 0.0093`
but `outcome` with 68 levels in model `res2` gives `sigma^2.2  =
0.0093`?
(DIFFERENT LEVELS BUT SAME RESULTS)

For reproducibility, below are my data and code.

Many thanks to you all,
Stefanou

m <- read.csv("https://raw.githubusercontent.com/fpqq/w/main/c.csv")

res <- rma.mv(yi, vi, random = list(~ 1 | study, ~1|outcome, ~ 1 |
measure), data=m)
                    estim       sqrt  nlvls  fixed   factor
sigma^2.1  0.0200  0.1415     57     no    study
sigma^2.2  0.0093  0.0964      4     no  outcome
sigma^2.3  0.0506  0.2249      7     no  measure

res2 <- rma.mv(yi, vi, random = list(~ 1 | study/outcome, ~ 1 |
measure), data=m)
                    estim       sqrt  nlvls  fixed         factor
sigma^2.1  0.0122  0.1105     57     no          study
sigma^2.2  0.0093  0.0964     68     no  study/outcome
sigma^2.3  0.0363  0.1904      7     no        measure



More information about the R-sig-meta-analysis mailing list