189 zwangerschappen. In deze telling kregen de rokers vaker een baby met een laag geboortegewicht — de odds liggen ruim twee keer zo hoog. Maar is dat verband echt, of kan het toeval het namaken? Links de telling. Rechts: hoe ver komt het toeval alleen al?
| moeder \ baby | laag gewicht | normaal gewicht | % laag |
|---|---|---|---|
| rookt | 30 | 44 | 41% |
| rookt niet | 29 | 86 | 25% |
De muur is nog leeg. Schud de moeders los van de baby's en kijk waar de odds ratio landt.
Schudden koppelt roken los van geboortegewicht: elke baby houdt zijn eigen gewicht, maar krijgt een wildvreemde moeder toegewezen — de een rookt, de ander niet, puur op de gok. Zo bouw je werelden waarin roken en gewicht niets met elkaar te maken hebben. En let op: de niks-waarde van een odds ratio is 1, niet 0 — even vaak laag als niet-laag betekent OR = 1. Kijk hoe het bergje rond 1 danst, en waar die rode lijn op 2.02 staat.
odds_ratio(roken, sample(laag_geboortegewicht))
Eén woord verschil met de echte som: sample(). Dat husselt de kolom door elkaar.
Duizend werelden waarin roken en geboortegewicht niets met elkaar te maken hebben — de gemeten OR
danst rond 1, niet rond 0; de echte 2.02 valt er (net) buiten.
De regressieve ruggengraat · W12 — het schud-blok bij de odds ratio · speeltje bij de simulatie
Data: geboortegewicht (Hosmer & Lemeshow, birthwt) · 189 zwangerschappen · roken × laag geboortegewicht.