This is promising, I very frequently end up writing one off functions to munge data on CIs between lme4 and ggplot2. However, what a prediction interval actual

This is promising, I very frequently end up writing one off functions to munge data on CIs between lme4 and ggplot2. However, what a prediction interval actually means is a slightly conentuous issue in the context of mixed models, and Douglas Bates had plenty of caution about constructing/interpreting them on the R mailing lists a few years back. This said, I’ve always found the tip-toeing around constructing them annoying; it just means that there’s little standard practice. Do you have any documentation on how you’re doing prediction intervals for mixed models, and what interpretations do and don’t make sense for your method?