Of course. If I find one I’ll post it here for posterity.
Edit: “method = "predict"
” accounts for the residual (observation-level) variance, the uncertainty in the fixed coefficients, and the uncertainty in the variance parameters for the grouping factors. “method = "fitted"
” accounts only for the uncertainty in the fixed coefficients and the uncertainty in the variance parameters for the grouping factors.