I’m trying to visualise my adjusted lmer model showing an interaction of education status using the predict function and ggplot, but I get a zig-zag output when I add in my covariates. I want to include 95% CI bands as well.