I have a cross-classified model as the following. I would like to plot the random effects of nonprofit on the industry level to a caterpillar plot. What are the steps? Thanks!
M <- lmer(salary ~ 1+education+male+nonprofit+(1+nonprofit|ind)+(1+nonprofit|occ), data, reml=TRUE)