Naděje dožití cyklistů

Kromě absolutního počtu zemřelých a zraněných při dopravních nehodách je důležitým ukazatele tzv. DALY - Disability-adjusted life years. Tento ukazatel popisuje o kolik let zdravého života jedinec/společnost přichází. https://www.who.int/data/gho/indicator-metadata-registry/imr-details/158.

Pro zjednodušení jsem pracoval jen se smrtelnými nehodami cyklistů a počítal jsem počet let života, o který je jedinec připraven. DALY jsem vypočetl jako rozdíl mezi očekávanou dobou dožití (na základě věku oběti v daném roce a pohlaví) a věku v době nehody. Data o očekávané době dižití byla vzata od ČSÚ (https://www.czso.cz/csu/czso/obyut_cr)

Počet zemřelých cyklistů při dopravních nehodách má sice v rámci sledovaného období 2011 - 2021 sestupný trend.

Zajímavé je, že podíl mužů mezi smrtelnými nehodami se v čase zvyšuje a to tempem cca 0.62 % za rok.

Průměrný věk zemřelých mužů nižší než průměrný věk žen, ale ne tolik, aby byl rozdíl signifikantní.

## 
##  Welch Two Sample t-test
## 
## data:  vek by pohlavi
## t = -1.0835, df = 268.01, p-value = 0.2795
## alternative hypothesis: true difference in means between group m and group f is not equal to 0
## 95 percent confidence interval:
##  -4.312559  1.250825
## sample estimates:
## mean in group m mean in group f 
##        55.70443        57.23529

U mužů dochází v čase ke zvyšování průměrného věku. U žen je trend vyhodnocen jako statisticky nesignifikantní, ale to může být dáno vysokou variabilitou dat - roky 2013 a 2016 jsou výrazně odlišné od zbytku dat.

Korelace u mužů je statisticky signifiantní

## 
##  Pearson's product-moment correlation
## 
## data:  ageAvgs$rok[ageAvgs$pohlavi == "m"] and ageAvgs$avg_age[ageAvgs$pohlavi == "m"]
## t = 2.9623, df = 14, p-value = 0.01029
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.1805736 0.8537348
## sample estimates:
##       cor 
## 0.6207177

Korelace u žen není statisticky signifikantní

## 
##  Pearson's product-moment correlation
## 
## data:  ageAvgs$rok[ageAvgs$pohlavi == "f"] and ageAvgs$avg_age[ageAvgs$pohlavi == "f"]
## t = 1.5319, df = 14, p-value = 0.1478
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1438299  0.7363038
## sample estimates:
##       cor 
## 0.3788892

Průměrné zkrácení doby života u smrtelných nehod cyklistů je 24.39 roku.

Ale průměrné zkrácení doby života je signifikantně vyšší u žen o cca 2.5 roku, což je dáno vyšší průměrnou dobou dožití u žen.

## 
##  Welch Two Sample t-test
## 
## data:  lexp by pohlavi
## t = -2.0177, df = 249.88, p-value = 0.04469
## alternative hypothesis: true difference in means between group m and group f is not equal to 0
## 95 percent confidence interval:
##  -4.89058840 -0.05909419
## sample estimates:
## mean in group m mean in group f 
##        23.93786        26.41271

Dá se předpokládat, že s rostoucí oblibou elektrokol a se zvyšujícími se počty seniorů bude přibývat smrtelných nehod starších cyklistů, kteří mají kratší dobu dožití. Průměrný počet let, o který zemřelí cyklisté byli připraveni, by tedy v čase měl klesat. Pro každý rok bylo vypočteno průměrné zkrácení doby dožití.

## 
## Call:
## lm(formula = avg_lexp ~ rok, data = lextAvgs)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.7378 -0.8984 -0.5286  0.5657  3.7617 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)   
## (Intercept) 533.63463  169.41952   3.150  0.00709 **
## rok          -0.25307    0.08414  -3.008  0.00941 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.551 on 14 degrees of freedom
## Multiple R-squared:  0.3925, Adjusted R-squared:  0.3491 
## F-statistic: 9.046 on 1 and 14 DF,  p-value: 0.009407

Průměrné zkrácení doby života v čase klesá a je statisticky signifikantní. Každý rok dojde k poklesu v průměru o -0.25 roku.

Současně tedy dochází ke zvyšování průměrného věku zemřelých při dopravních nehodách a zároveň ke snižování zkrácení doby dožití. Lze to interpretovat tak, že dopad na společnost vyjádřený počtem ztracených let života u nehod cyklistů v čase pomalu klesá tím, jak stále více zemřelých tvoří senioři.

Výrazný posun v závažnosti následků nehod cyklistů ve vyšším věku ukazuje následující graf. Výsledný histogram je pro každý sloupec součtem histogramů pro čtyři úrovně následků. Každý z dílčích histogramů zabírá v grafu stejnou plochu - stejnou plochu zabírá histogram smrtelných nehod, které tvoří 1 % nehod cyklistů a stejnou plochu zabírá histogram nehod s lehkým zraněním, které tvoží 73 % nehod cyklistů. To umožňuje lepší vzájemné srovnání, než kdyby plocha odpovídala počtu nehod - méně četné kategorie nehod by byly ve výsledném histogramu nerozpoznatelné.

## Warning: Removed 14 rows containing non-finite values (`stat_bin()`).
## Warning: Removed 8 rows containing missing values (`geom_bar()`).