Participation in activities of churches or other religious organisations by level of activity limitation, sex and age - hlth_ds040

Data - Eurostat

Author

age

Code

hlth_ds040 %>%
  left_join(age, by = "age") %>%
  group_by(age, Age) %>%
  summarise(Nobs = n()) %>%
  arrange(-Nobs) %>%
  {if (is_html_output()) print_table(.) else .}

age	Age	Nobs
Y16-29	From 16 to 29 years	348
Y25-34	From 25 to 34 years	348
Y35-44	From 35 to 44 years	348
Y_GE16	16 years or over	348
Y_GE65	65 years or over	348
Y16-24	From 16 to 24 years	347
Y16-64	From 16 to 64 years	347
Y45-54	From 45 to 54 years	347
Y55-64	From 55 to 64 years	347

unit

Code

hlth_ds040 %>%
  left_join(unit, by = "unit") %>%
  group_by(unit, Unit) %>%
  summarise(Nobs = n()) %>%
  arrange(-Nobs) %>%
  {if (is_html_output()) print_table(.) else .}

unit	Unit	Nobs
PC	Percentage	3128

geo

Code

hlth_ds040 %>%
  left_join(geo, by = "geo") %>%
  group_by(geo, Geo) %>%
  summarise(Nobs = n()) %>%
  arrange(-Nobs) %>%
  {if (is_html_output()) datatable(., filter = 'top', rownames = F) else .}

lev_limit

Code

hlth_ds040 %>%
  left_join(lev_limit, by = "lev_limit") %>%
  group_by(lev_limit, Lev_limit) %>%
  summarise(Nobs = n()) %>%
  arrange(-Nobs) %>%
  {if (is_html_output()) print_table(.) else .}

lev_limit	Lev_limit	Nobs
NONE	None	1566
SM_SEV	Some or severe	1562

yes_no

Code

hlth_ds040 %>%
  left_join(yes_no, by = "yes_no") %>%
  group_by(yes_no, Yes_no) %>%
  summarise(Nobs = n()) %>%
  arrange(-Nobs) %>%
  {if (is_html_output()) print_table(.) else .}

yes_no	Yes_no	Nobs
NO	No	1566
YES	Yes	1562

sex

Code

hlth_ds040 %>%
  left_join(sex, by = "sex") %>%
  group_by(sex, Sex) %>%
  summarise(Nobs = n()) %>%
  arrange(-Nobs) %>%
  {if (is_html_output()) print_table(.) else .}

sex	Sex	Nobs
F	Females	1043
T	Total	1043
M	Males	1042

Religious Share

Code

hlth_ds040 %>%
  filter(age == c("Y_GE16"),
         sex == "T",
         lev_limit == "NONE",
         yes_no == "YES") %>%
  left_join(geo, by = "geo") %>%
  right_join(europe_NUTS0, by = "geo") %>%
  filter(long >= -13.8, lat >= 33) %>%
  ggplot(., aes(x = long, y = lat, group = group, fill = values/100)) +
  geom_polygon() + coord_map() +
  scale_fill_viridis_c(na.value = "white",
                       labels = percent_format(accuracy = 1),
                       breaks = 0.01*seq(0, 90, 10),
                       values = c(0, 0.1, 0.2, 0.3, 0.4, 0.5, 1)) +
  theme_void() + theme(legend.position = c(0.15, 0.85)) + 
  labs(fill = "Religious Share")