coef_names <- coef_names[1:3] # Dropping intercept for plots plot_summs(fit, fit_b, fit_c, robust = "HC3", coefs = coef_names)
plot_summs( fit_c, robust = "HC3", coefs = coef_names, plot.distributions = TRUE )