data(regress_age)
# Homogeneous residual variance
metafit <- ipdlme(n,Y,S2)
summary(metafit)
fixef(metafit)
confint(metafit)
ranef(metafit)
labs <- paste("Center",c(1,10,11,2:9))
plot(metafit, y=labs)
#Heterogeneous residual variance
metafit <- ipdlme(n,Y,S2,equal=FALSE)
summary(metafit)
fixef(metafit)
confint(metafit)
ranef(metafit)
labs <- paste("Center",c(1,10,11,2:9))
plot(metafit, y=labs)
Run the code above in your browser using DataLab