I have pooled results from the mitml package:
fit <- with(data = implist, gls(bsi_gsi ~ Mean_LPA + MEAN_MVPA + agl_ado_age + agl_ado_sex + season, method= "ML", na.action="na.omit"))
summary(pool(fit))
I'd like to plot the pooled results- any suggestions?