# NOT RUN {
# RUN THE MODEL
model <- mixmeta(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