plot(gamm_intercept, select=2) # select=2 because the random effect appears as the second entry in the summary table.