Máte dva datasety - pochází ze sčítání lidu v roce 1930 (poslední převálečné) a v roce 1980 (poslední socíkové). Oba bydlí v knihovně RCzechia.
Prvorepublikový dataset je v členění soudních okresů (n = 330), socíkový dataset je v členění okresů z dob ČSSR (n = 76; pro fanoušky: neexistoval okres Jeseník, Harrachov byl v Semilech).
Z vizualizace se pocitově zdá (a lze matematicky prokázat, což bude předmětem budoucí výuky), že podíl důchodců (definováno jako podíl obyvatel ve věku 65 let a starších ku celkovému počtu obyvatel přítomných) nebyl v roce 1980 rozdělen po republice rovnoměrně.
Vaším úkolem je vybrat ze dvou alternativních hypotéz tu, která věrněji vysvětluje tento nepoměr:
hypotéza Modrých zón: tato hypotéza tvrdí, že existují oblasti vyznačující se mimořádnou dlouhověkostí. Příčinou dlouhověkosti je tradiční způsob života, konzumace luštěnin a ořechů a funkční rodinné a společenské vazby. Podle hypotézy modrých zón očekáváme, že podíl důchodců v roce 1980 v regionu bude souviset s podílem důchodců v témže regionu v roce 1930.
hypotéza Postupim: tato hypotéza tvrdí, že po druhé válce proběhly v rámci Československa velké přesuny obyvatelstva. Pohraničí bylo opětovně osidlováno převážně mladými lidmi z vnitrozemí. Podle hypotézy Postupim očekáváme, že podíl důchodců v roce 1980 v regionu bude souviset s podílem občanů německé národnosti v témže regionu v roce 1930.
Oba datasety jsou pro zjednodušení transformovány do souřadnicového systému ing. Křováka.
Oba datasety obsahují absolutní počty obyvatel, ne podíly. Počty jsou tzv. extenzivní veličina – když okres rozdělíte, rozdělí se i jeho obyvatelé. Mějte to na paměti při počítání podílů.
Kontrola: přenášená metrika (obyvatelé, důchodci a / nebo Němci) by se přenosem z jedněch okresů na druhé neměla (výrazně) změnit. Tuto kontrolu vykonejte před vlastní regresí.
library(RCzechia)library(ggplot2)library(dplyr)# census 1980okresy_1980 <- RCzechia::historie("okresy_1980") %>%mutate(duchodci_80 =`obyvatelstvo celkem 65+`,obyvatele_80 =`počet obyvatel přítomných`) %>%select(duchodci_80, obyvatele_80) %>%st_transform(5514)# náhled na strukturuglimpse(okresy_1980)
