# NOT RUN {
### REPRODUCE THE RESULTS IN CHEN ET AL. (2012)
# INSPECT THE DATA
hsls
# FIXED-EFFECTS MODEL
S <- as.matrix(hsls[5:10])
model <- mvmeta(cbind(b1,b2,b3),S,data=hsls,method="fixed")
summary(model)
# MM MODEL
model <- mvmeta(cbind(b1,b2,b3),S,data=hsls,method="mm")
summary(model)
model$Psi
# }
Run the code above in your browser using DataLab