Dane

Dane dotyczące zamachów samobójczych ze strony statystyka.policja.pl/st/wybrane-statystyki/zamachy-samobojcze/63803,Zamachy-samobojcze-od-2017-roku.html. Dane dotyczące liczby ludności z bazy Eurostatu (tabela demo_r_d2jan). Liczba ludności policzona dla roku 2021 (powinna być dla każdego roku oddzielnie, ale różnica nie jest duża a nie chciało mi się)

z <- read.csv("demo_r_d2jan_tabular.csv", sep = ';',  header=T, na.string="NA" )


nutsPL <- c ('PL21', 'PL22', 'PL41', 'PL42', 'PL43',
           'PL51', 'PL52', 'PL61', 'PL62', 'PL63', 'PL71',
           'PL72', 'PL81', 'PL82', 'PL84', 'PL91', 'PL92')


##levels(as.factor(z$tmp))
##levels(as.factor(z$age))

c <- read.csv("nuts-PL.csv", sep = ';',  header=T, na.string="NA" ) 

z1 <- z %>% mutate(age=recode(age, 'Y13'= "Y13.18", 'Y14'= "Y13.18", 'Y15'= "Y13.18",
                             'Y16'= "Y13.18", 'Y17'= "Y13.18", 'Y18'= "Y13.18" )) %>%
  filter (age == 'Y13.18' | age == 'TOTAL') %>%
  filter (year == 2021) %>%
  filter (geo %in% nutsPL) %>%
  ## Komendy są wg NUTS
  ## aggregate 91/92 together
  ##mutate(geo = recode (geo, 'PL91'='PL92')) %>%
  group_by(geo, sex, age) %>%
  summarise(n = sum(value)) %>%
  mutate (age = as.factor(age)) %>%
  pivot_wider(names_from = age, values_from = n) %>%
  mutate (p = Y13.18/TOTAL * 100) %>%
  left_join(c, by="geo") %>%
     ungroup()
## `summarise()` has grouped output by 'geo', 'sex'. You can override using the
## `.groups` argument.
## liczba ludności wg województw
z.total <- z1 %>% filter (sex == 'T')
z.total.PL <- z.total %>%
    group_by(sex) %>%
    summarise (T = sum(TOTAL), T13.18 = sum(Y13.18)) %>%
    ungroup()

##
z.total.with.total <- z.total %>%
  add_row(geo = 'Polska', TOTAL = z.total.PL$T, 
          Y13.18 = z.total.PL$T13.18,
          sex = 'T',
          name =  'Polska',
          p = Y13.18/TOTAL * 100)

Zamachy samobójcze (skuteczne/nieskuteczne)

v <- read.csv("zamachy_samobojcze_grupa_wiekowa_dzien_tygodnia_2017-2022.csv", sep = ';',  
              header=T, na.string="NA" )  %>%
  mutate (KWP=recode (KWP, 
            "KSPWarszawa"='PL91',
            "KWPBiałystok"='PL84',
            "KWPBydgoszcz"='PL61',
            "KWPGdańsk"='PL63',
            "KWPGorzów Wlkp."='PL43',
            "KWPKatowice" ='PL22',
            "KWPKielce"='PL72',
            "KWPKraków"='PL21',    
            "KWPLublin"='PL81',
            "KWPŁódź"='PL71',
            "KWPOlsztyn"='PL62',
            "KWPOpole"='PL52',
            "KWPPoznań"='PL41',
            "KWPRadom"='PL92',
            "KWPRzeszów"='PL82',
            "KWPSzczecin"='PL42',
            "KWPWrocław"='PL51'), KWP=as.factor(KWP))

v1 <- v %>%
  left_join(z.total.with.total, by=c('KWP'='geo')) %>%
  mutate (y_13.18p = y_13.18 / Y13.18 * 10000) %>%
  mutate (y_total = total / TOTAL * 10000)
p1 <- v1 %>%
  select (rok, name, y_13.18p) %>%
  ggplot(aes(x=rok, y=y_13.18p)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_13.18p, label= sprintf("%.1f",y_13.18p) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed') +
  ggtitle("Próby samobójcze (w grupie 13--18 lat)", 
          subtitle = "Na 10,000 osób w wieku 13--18 lat") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))
p1 

p2 <- v1 %>%
  select (rok, name, y_total) %>%
  ggplot(aes(x=rok, y=y_total)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_total, label= sprintf("%.1f",y_total) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed',
             #scales = 'free_y'
             ) +
  ggtitle("Próby samobójcze", 
          subtitle = "Na 10,000 osób") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))

p2 

ggsave(p1, file='proby-samobojcze-1318.png', height = 8, width=8)
ggsave(p2, file='proby-samobojcze-razem.png', height = 6, width=8)
v2 <- v1 %>% 
  group_by(name) %>%
  mutate (y_13.18p = y_13.18p/first(y_13.18p) * 100,
    y_total = y_total/first(y_total) * 100)


p3 <- v2 %>%
  select (rok, name, y_13.18p) %>%
  ggplot(aes(x=rok, y=y_13.18p)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_13.18p, label= sprintf("%.1f",y_13.18p) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed') +
  ggtitle("Próby samobójcze (w grupie 13--18 lat)", 
          subtitle = "2017 = 100%") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))
p3

p4 <- v2 %>%
  select (rok, name, y_total) %>%
  ggplot(aes(x=rok, y=y_total)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_total, label= sprintf("%.1f",y_total) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed',
             #scales = 'free_y'
  ) +
  ggtitle("Próby samobójcze", 
          subtitle = "2017 = 100%") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))

p4

ggsave(p3, file='proby-samobojcze-1318_p.png', height = 8, width=8)
ggsave(p4, file='proby-samobojcze-razem_p.png', height = 6, width=8)

Zamachy samobójcze zakończone zgonem

v <- read.csv("zamachy_samobojcze_zakonczone_ZGONEM_grupa_wiekowa_dzien_tygodnia_2017-2022.csv", sep = ';',  
              header=T, na.string="NA" )  %>%
  mutate (KWP=recode (KWP, 
            "KSPWarszawa"='PL91',
            "KWPBiałystok"='PL84',
            "KWPBydgoszcz"='PL61',
            "KWPGdańsk"='PL63',
            "KWPGorzów Wlkp."='PL43',
            "KWPKatowice" ='PL22',
            "KWPKielce"='PL72',
            "KWPKraków"='PL21',    
            "KWPLublin"='PL81',
            "KWPŁódź"='PL71',
            "KWPOlsztyn"='PL62',
            "KWPOpole"='PL52',
            "KWPPoznań"='PL41',
            "KWPRadom"='PL92',
            "KWPRzeszów"='PL82',
            "KWPSzczecin"='PL42',
            "KWPWrocław"='PL51'), KWP=as.factor(KWP))

v1 <- v %>%
  left_join(z.total.with.total, by=c('KWP'='geo')) %>%
  mutate (y_13.18p = y_13.18 / Y13.18 * 10000) %>%
  mutate (y_total = total / TOTAL * 10000)

source <- 'statystyka.policja.pl/st/wybrane-statystyki/zamachy-samobojcze/63803,Zamachy-samobojcze-od-2017-roku.html + Eurostat: demo_r_d2jan'
p1 <- v1 %>%
  select (rok, name, y_13.18p) %>%
  ggplot(aes(x=rok, y=y_13.18p)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_13.18p, label= sprintf("%.1f",y_13.18p) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed') +
  ggtitle("Próby samobójcze (w grupie 13--18 lat)", 
          subtitle = "Na 10,000 osób w wieku 13--18 lat") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))
p1 

p2 <- v1 %>%
  select (rok, name, y_total) %>%
  ggplot(aes(x=rok, y=y_total)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_total, label= sprintf("%.1f",y_total) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed',
             #scales = 'free_y'
             ) +
  ggtitle("Próby samobójcze", 
          subtitle = "Na 10,000 osób") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))

p2 

ggsave(p1, file='proby-samobojcze-Z-1318.png', height = 8, width=8)
ggsave(p2, file='proby-samobojcze-Z-razem.png', height = 6, width=8)
v2 <- v1 %>% 
  group_by(name) %>%
  mutate (y_13.18p = y_13.18p/first(y_13.18p) * 100,
    y_total = y_total/first(y_total) * 100)


p3 <- v2 %>%
  select (rok, name, y_13.18p) %>%
  ggplot(aes(x=rok, y=y_13.18p)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_13.18p, label= sprintf("%.1f",y_13.18p) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed') +
  ggtitle("Próby samobójcze (w grupie 13--18 lat)", 
          subtitle = "2017 = 100%") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))
p3
## Warning: Removed 1 rows containing missing values (geom_point).
## Warning: Removed 1 rows containing missing values (geom_text).

p4 <- v2 %>%
  select (rok, name, y_total) %>%
  ggplot(aes(x=rok, y=y_total)) +
  geom_line() +
  geom_point() +
  geom_text(aes(x=rok, y=y_total, label= sprintf("%.1f",y_total) ), 
            vjust=-.5,
            size=2.0, alpha=.5) +
  facet_wrap(~ name, scales = 'fixed',
             #scales = 'free_y'
  ) +
  ggtitle("Próby samobójcze", 
          subtitle = "2017 = 100%") +
  labs(caption=source) +
  ylab("liczba prób/10 tys") + xlab("")+
  expand_limits(y = 0) +
  theme(plot.caption=element_text(size=6), axis.text.x = element_text(size=6))

p4

ggsave(p3, file='proby-samobojcze-Z-1318_p.png', height = 8, width=8)
## Warning: Removed 1 rows containing missing values (geom_point).
## Warning: Removed 1 rows containing missing values (geom_text).
ggsave(p4, file='proby-samobojcze-Z-razem_p.png', height = 6, width=8)