GLMM vs GEE plot with 95% confidence intervals

Plotting 95% Confidence Bands in R

I am comparing estimates from subject-specific GLMMs and population-average GEE models as part of a publication I am working on. As part of this, I want to visualize predictions of each type of model including 95% confidence bands. First I had to make a new data set for prediction. I could have compared fitted values…