# NOT RUN {
### REPRODUCE THE RESULTS OF EXAMPLE 3 IN JACKSON ET AL. (2011)
# INSPECT THE DATA
p53
# REML MODEL WITH INPUTTED CORRELATION EQUAL TO 0.95
model <- mixmeta(cbind(y1,y2), cbind(V1,V2), data=p53, control=list(Scor=0.95))
print(summary(model), digits=2)
# }
Run the code above in your browser using DataLab