Field note · October 22, 2023

bmi by sex: a 1% spread

Real group means for one column across 2 segments, and what the pooled average conceals.

1 min read ·Analytics practice ·statistics practice

Breaking bmi down by sex on clinical-trial, because the headline average is 27.13 and no segment is actually there.

  • M — mean 27.32, median 27.50 (402 rows)
  • F — mean 26.97, median 26.90 (498 rows)

Top to bottom that is 27.32 against 26.97, a spread of 1.3%. The pooled average is 27.13.

sql
select sex,
  count(*)                                         as rows,
  round(avg(bmi)::numeric, 2)                      as mean,
  percentile_cont(0.5) within group (order by bmi) as median
from clinical_trial
group by 1
order by mean desc;

The groups are close enough that the pooled average is a fair summary. That is worth confirming rather than assuming: the check costs one query and the failure mode is invisible.

Notice the mean and median columns disagree only slightly here. Always compute both in the group-by. The comparison between them per segment is free and tells you whether you are looking at a level difference or a tail difference.

This is the setup for Simpson's paradox — the case where every segment moves one way and the total moves the other.