# RUN THE MODEL AND SUMMARIZE THE RESULTS
model <- mvmeta(cbind(PD,AL)~pubyear,S=berkey98[5:7],data=berkey98,method="ml")
# MODEL FRAME
model$model
model.frame(model)
update(model,method="model.frame")
class(model.frame(model))
# MODEL MATRIX
model.matrix(model)
Run the code above in your browser using DataLab