# NOT RUN {
### copy data into 'dat'
dat <- dat.craft2003
dat
### construct dataset and var-cov matrix of the correlations
tmp <- rcalc(ri ~ var1 + var2 | study, ni=ni, data=dat)
V <- tmp$V
dat <- tmp$dat
### examine data for study 1
dat[dat$study == 1,]
V[dat$study == 1, dat$study == 1]
### examine data for study 6
dat[dat$study == 6,]
V[dat$study == 6, dat$study == 6]
### examine data for study 17
dat[dat$study == 17,]
V[dat$study == 17, dat$study == 17]
### multivariate random-effects model
res <- rma.mv(yi, V, mods = ~ var1.var2 - 1, random = ~ var1.var2 | study, struct="UN", data=dat)
res
# }
Run the code above in your browser using DataLab