### Panel B: Improved Figure 4 - EOP Enlargement Rates The original Figure 4 showed counts, which is misleading when group sizes were controlled. We show proportions/rates instead. ```{r figure4_improved} # Calculate rates of enlarged EOP by age group and sex eop_rates = posture_df |> drop_na(eop_enlarged, age_group, sex) |> group_by(age_group, sex) |> summarize( n_total = n(), n_enlarged = sum(eop_enlarged == "Enlarged"), rate_enlarged = n_enlarged / n_total, .groups = "drop" ) fig4_improved = eop_rates |> ggplot(aes(x = age_group, y = rate_enlarged, fill = sex)) + geom_col(position = position_dodge(0.9), alpha = 0.8) + geom_text( aes(label = sprintf("%.1f%%", rate_enlarged * 100)), position = position_dodge(0.9), vjust = -0.5, size = 3 ) + scale_y_continuous( labels = scales::percent_format(), limits = c(0, 1), breaks = seq(0, 1, 0.2) ) + labs( title = "Panel B: Rate of Enlarged EOP by Age and Sex", subtitle = "Proportion of participants with enlarged EOP (categories 2-5)", x = "Age Group", y = "Proportion with Enlarged EOP", fill = "Sex" ) + theme( axis.text.x = element_text(angle = 45, hjust = 1), legend.position = "bottom" ) fig4_improved ```