There are two random effects - η₁ and η₂ - and I would like to relabel the y-axis titles as η_CL and η_Vc. Can we do that in the following function?
empirical_bayes_vs_covariates(
fitted_model_inspection;
categorical = [:formulation, :racen, :gender],
)