Słownik metodologiczny · Zależności i efekty

Paradoks Simpsona – co to jest i dlaczego związek odwraca się w podgrupach?

Paradoks Simpsona to odwrócenie kierunku związku po podziale próby na podgrupy. Pokazujemy, skąd się bierze, na przykładzie tabeli krzyżowej policzonej w R i wyjaśniamy, od czego zależy, której tabeli wierzyć.

Strona główna / Słownik / Paradoks Simpsona
W skrócie: Paradoks Simpsona to sytuacja, w której związek w całej próbie ma inny kierunek niż w każdej podgrupie – zwykle dlatego, że porównywane grupy mają różny skład pod względem trzeciej zmiennej. Którą wersję wyniku przyjąć, rozstrzyga wiedza o przyczynach, a nie sama statystyka.

Co to jest paradoks Simpsona?

Paradoks Simpsona (ang. Simpson's paradox) to sytuacja, w której związek dwóch zmiennych w całej próbie ma inny kierunek niż w każdej z jej podgrup. Kievit i in. (2013) pokazują, że w psychologii zdarza się częściej, niż się sądzi, zwłaszcza gdy wnioski przenosi się między poziomami: z populacji na podgrupy albo z grup na pojedyncze osoby. Najłatwiej go zobaczyć w tabeli krzyżowej z odsetkami.

Skąd się bierze

Odsetek w całej grupie to średnia ważona odsetków z podgrup:

odsetekrazem = w1 × odsetek1 + w2 × odsetek2

gdzie w1 i w2 to udziały podgrup w porównywanej grupie. Kierunek może się odwrócić, gdy trzecia zmienna silnie wiąże się z wynikiem, a porównywane grupy mają bardzo różny skład pod jej względem. Najczęściej jest nią zmienna zakłócająca.

Przykład

Dane 640 studentów z trzech roczników: udział w zajęciach wyrównawczych ze statystyki i zaliczenie egzaminu (dane przykładowe).

Odsetek zaliczeń egzaminu według udziału w zajęciach wyrównawczych i poziomu matury z matematyki
Matura z matematykiNa zajęciachBez zajęć
Podstawowa60.0%50.0%
Rozszerzona90.0%80.0%
Razem67.5% (n = 320)72.5% (n = 320)

W obu podgrupach uczestnicy zajęć zdawali częściej, a w całej próbie – rzadziej (iloraz szans OR = 0.79). Na zajęcia chodziły głównie osoby z maturą podstawową (75% uczestników wobec 25% pozostałych), które zdawały rzadziej niezależnie od zajęć: 0.75 × 60% + 0.25 × 90% = 67.5%, a 0.25 × 50% + 0.75 × 80% = 72.5%. Iloraz szans Mantela–Haenszela, liczony w podgrupach, wyniósł OR = 1.71, 95% CI [1.12, 2.62].

Jak zapisać w pracy

Podaj wynik łączny i wynik w podgrupach:

W całej próbie (N = 640) egzamin zaliczyło 67.5% uczestników zajęć wyrównawczych i 72.5% pozostałych osób, jednak po uwzględnieniu poziomu matury z matematyki odsetek zaliczeń był wyższy wśród uczestników zajęć (60.0% wobec 50.0% przy maturze podstawowej i 90.0% wobec 80.0% przy rozszerzonej), OR Mantela–Haenszela = 1.71, 95% CI [1.12, 2.62], χ²(1) = 5.65, p = .017.
Przykład reprodukowalny

Paradoks Simpsona w R

Kod liczy odsetki zaliczeń w tabeli łącznej i w podgrupach, skład grup oraz iloraz szans – surowy i Mantela–Haenszela. Liczby w przykładzie pochodzą z tego kodu.

Tabela 2 × 2 × 2 i test Mantela–Haenszela

Dane wejściowe: liczebności z przykładu (dane przykładowe, 640 studentów).

Co odczytać: odsetki zaliczeń razem i w podgrupach, liczebności i skład grup, surowy OR oraz wspólny OR Mantela–Haenszela z 95% CI.

R
# Zajęcia wyrównawcze a zaliczenie egzaminu: 640 studentów, warstwy według matury z matematyki (dane przykładowe)
tab <- array(c(144, 96, 40, 40,   72, 8, 192, 48), dim = c(2, 2, 2), dimnames = list(
  egzamin = c("zdany", "niezdany"), zajecia = c("tak", "nie"), matura = c("podstawowa", "rozszerzona")))
razem <- apply(tab, 1:2, sum)                                  # tabela łączna (bez podziału na warstwy)
round(100 * prop.table(razem, 2)["zdany", ], 1)               # odsetek zaliczeń: razem
round(100 * prop.table(tab, 2:3)["zdany", , ], 1)             # odsetek zaliczeń w warstwach
apply(tab, 2:3, sum)                                           # liczebności grup
round(100 * prop.table(apply(tab, 2:3, sum), 1), 1)           # skład grup według matury
round(c(OR_laczny = (razem[1, 1] * razem[2, 2]) / (razem[1, 2] * razem[2, 1])), 2)
mantelhaen.test(tab)                                           # wspólny iloraz szans w warstwach
Wynik w R
 tak  nie 
67.5 72.5 
       matura
zajecia podstawowa rozszerzona
    tak         60          90
    nie         50          80
       matura
zajecia podstawowa rozszerzona
    tak        240          80
    nie         80         240
       matura
zajecia podstawowa rozszerzona
    tak         75          25
    nie         25          75
OR_laczny 
     0.79 

	Mantel-Haenszel chi-squared test with continuity correction

data:  tab
Mantel-Haenszel X-squared = 5.6522, df = 1, p-value = 0.01743
alternative hypothesis: true common odds ratio is not equal to 1
95 percent confidence interval:
 1.121348 2.620752
sample estimates:
common odds ratio 
         1.714286 

Bibliografia oprogramowania

Wyniki na tej stronie pochodzą z uruchomienia pokazanego kodu w R 4.5.3 (pakiet stats). W pracy podaj wersje faktycznie użytych pakietów (sessionInfo()).

  1. R: R Core Team. (2026). R: A language and environment for statistical computing (Version 4.5.3) [Computer software]. R Foundation for Statistical Computing. https://www.R-project.org/

Na co uważać

  • Statystyka nie rozstrzyga, której tabeli wierzyć. Kievit i in. (2013) przypominają, że same dane obserwacyjne nie mówią, czy trzecia zmienna jest zakłócająca, czy jest mediatorem – a od tego zależy, czy ją kontrolować. Poziom matury powstał przed zajęciami, więc liczą się podgrupy; gdyby był skutkiem badanej przyczyny, właściwa byłaby tabela łączna.
  • Odwrócenie nie dowodzi skuteczności zajęć. To nadal dane z badania korelacyjnego; podobnie mogą działać niezmierzone różnice, a mocny wniosek daje dopiero randomizacja.
  • Klasyczny przykład. Bickel i in. (1975) pokazali, że w zbiorczych danych z Berkeley kobiety rzadziej przyjmowano na studia doktoranckie, choć w podziale na wydziały niewiele z nich istotnie odbiegało od oczekiwanych odsetków przyjęć kobiet, a mniej więcej tyle samo wydziałów sprzyjało kobietom, co mężczyznom. Kobiety częściej aplikowały na wydziały trudniejsze dla kandydatów obu płci.

To skrajna postać problemu opisanego w haśle korelacja pozorna: trzecia zmienna odwraca kierunek związku. Ten sam mechanizm w danych medycznych pokazuje poradnik statystyka dla lekarza.

Najczęstsze pytania

Skąd się bierze paradoks Simpsona?

Z różnego składu porównywanych grup. Gdy trzecia zmienna silnie wiąże się z wynikiem, a jedna grupa ma znacznie więcej osób z jej „gorszego” poziomu, odsetek łączny może odwrócić kierunek widoczny w każdej podgrupie.

Którym wynikom wierzyć przy paradoksie Simpsona?

To zależy od roli trzeciej zmiennej, a nie od statystyki. Jeśli jest wspólną przyczyną (zmienną zakłócającą), patrz na podgrupy; jeśli jest skutkiem badanej przyczyny, czyli mediatorem, właściwa jest tabela łączna.

Jak wykryć paradoks Simpsona w danych?

Policz związek w całej próbie i w podgrupach wyznaczonych przez zmienne, które mogą wiązać się z wynikiem, np. poziom wyjściowy. Dla tabel 2 × 2 w kilku podgrupach użyj testu Mantela–Haenszela: mantelhaen.test() w R.

Powiązane hasła

Poradniki, w których użyjesz tego pojęcia

Potrzebujesz wsparcia w analizie danych?

Wybierz ofertę dopasowaną do pracy dyplomowej albo profesjonalnego projektu naukowego.

Wyceń analizęOferta dla studentów

Bibliografia i dalsza literatura naukowa (APA 7)

  1. Bickel, P. J., Hammel, E. A., & O’Connell, J. W. (1975). Sex bias in graduate admissions: Data from Berkeley. Science, 187(4175), 398–404. https://doi.org/10.1126/science.187.4175.398
  2. Kievit, R. A., Frankenhuis, W. E., Waldorp, L. J., & Borsboom, D. (2013). Simpson's paradox in psychological science: A practical guide. Frontiers in Psychology, 4, Article 513. https://doi.org/10.3389/fpsyg.2013.00513 ● Otwarty dostęp

Cytowania BibTeX

@article{bickel1975sex, author={Bickel, P. J. and Hammel, E. A. and O'Connell, J. W.}, title={{Sex bias in graduate admissions: Data from Berkeley}}, journal={Science}, year={1975}, volume={187}, number={4175}, pages={398--404}, doi={10.1126/science.187.4175.398}}
@article{kievit2013simpson, author={Kievit, Rogier A. and Frankenhuis, Willem E. and Waldorp, Lourens J. and Borsboom, Denny}, title={{Simpson's paradox in psychological science: A practical guide}}, journal={Frontiers in Psychology}, year={2013}, volume={4}, eid={513}, doi={10.3389/fpsyg.2013.00513}}
R