Dane

Dane dotyczące zgonów pobrano ze strony https://dane.gov.pl/pl/dataset/1953,liczba-zgonow-zarejestrowanych-w-rejestrze-stanu-cywilnego/resource/48464/table?page=1&per_page=20&q=&sort= Nazwa zbioru: Liczba zgonów zarejestrowanych w Rejestrze Stanu Cywilnego w okresie od 1 września 2015 r. - dane tygodniowe Rodzaj licencji: CC0 1.0 Link do licencji: https://creativecommons.org/publicdomain/zero/1.0/legalcode.pl

Biblioteki

install.packages("tidyverse")
## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.3'
## (as 'lib' is unspecified)
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.2     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ ggplot2   3.4.2     ✔ tibble    3.2.1
## ✔ lubridate 1.9.2     ✔ tidyr     1.3.0
## ✔ purrr     1.0.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

Utworzenie ramki danych

zgony <- read.csv("/cloud/project/Liczba_zgonów_w_Rejestrze_Stanu_Cywilnego_w_latach_2015-2023_-_dane_tygodniowe_ogółem.csv")

Analiza ramki danych

head(zgony)
##   Nr.tygodnia X2015 X2016 X2017 X2018 X2019 X2020 X2021 X2022 X2023
## 1           1    NA  8246  9321  8406  8706  8473 11496 12036 10442
## 2           2    NA  7942 10270  8107  8613  8275 10996 10808  9560
## 3           3    NA  7702  9902  8663  8740  8091 11038 10020  8641
## 4           4    NA  7999  9966  8811  8833  8042 10432 10184  8316
## 5           5    NA  8000 10102  9009  9381  8578  9903 10587  8281
## 6           6    NA  8160  9464  9118  9303  8363  9950 10718  8278
colnames(zgony)
##  [1] "Nr.tygodnia" "X2015"       "X2016"       "X2017"       "X2018"      
##  [6] "X2019"       "X2020"       "X2021"       "X2022"       "X2023"

Analiza

Przygotowanie zbioru do analizy

zgony1 <- rename(zgony, Rok_2015=X2015, Rok_2016=X2016, Rok_2017=X2017,Rok_2018=X2018,Rok_2019=X2019, Rok_2020=X2020,Rok_2021=X2021,Rok_2022=X2022,Rok_2023=X2023)
zgony_tyg <-c(rowSums(zgony1, na.rm=TRUE))

Utworzenie ramki danych

nr_tyg <- c(1:53)
zgony_tyg_lata <- data.frame(Nr_tyg=nr_tyg, zgony_tyg=zgony_tyg)

Wykres

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

# Zgony w ujęciu rocznym

Przygotowanie danych

zgony1_bez_tyg <- zgony1 %>% select(Rok_2015, Rok_2016, Rok_2017, Rok_2018, Rok_2019, Rok_2020, Rok_2021, Rok_2022, Rok_2023)
zgony_suma <- colSums(zgony1_bez_tyg, na.rm=TRUE)
rok <- c(2015,2016,2017,2018,2019,2020,2021,2022,2023)
liczba_zgon_rok <- c(133057,388426,404057,412859,409726,486430,520166,448712,186491)

Ramka danych

zgony_lata <- data.frame(Rok=rok, liczba=liczba_zgon_rok)

Wykres liczba zgonów w ujęciu rocznym

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.