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