# NOT RUN {
# RUN THE MODEL
model <- mvmeta(cbind(PD,AL)~pubyear,S=berkey98[5:7],data=berkey98)
# COEFFICIENTS
model$coef
coef(model)
coef(model,format="matrix")
summary(model)$coef
# (CO)VARIANCE MATRIX
vcov(model)
# }
Run the code above in your browser using DataLab