R4CSS/V06_Association.Rmd at master · p4css/R4CSS · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
```{r include=FALSE}
knitr::opts_chunk$set(cache = T, warning = F, message = F,
                      class.output = "output", out.width='100%',
                      fig.asp = 0.5, fig.align = 'center')
library(tidyverse)
```

# ASSOCIATION

## 等比例座標軸

### UNICEF-Optimistic (WGOITH)

<https://www.nytimes.com/2021/11/17/upshot/global-survey-optimism.html> <https://changingchildhood.unicef.org/about>

```{r}
plot.opt <- read_csv("data/unicef-changing-childhood-data.csv") %>%
    select(country = WP5, age = WP22140, bw = WP22092) %>%
    mutate(country = ordered(country,
                             levels=c(1, 3, 4, 10, 11, 12,
                                      13, 14, 17, 29, 31,
                                      33, 35, 36, 60, 61,
                                      77, 79, 81, 87, 165),
                             labels=c("USA", "Morocco", "Lebanon",
                                      "Indonesia", "Bangladesh",
                                      "UK", "France", "Germany",
                                      "Spain", "Japan", "India",
                                      "Brazil", "Nigeria", "Kenya",
                                      "Ethiopia", "Mali", "Ukraine",
                                      "Cameroon", "Zimbabwe",
                                      "Argentina", "Peru"))) %>%
    count(country, age, bw) %>%
    group_by(country, age) %>%
    mutate(perc = n/sum(n)) %>%
    ungroup() %>%
    filter(bw == 1) %>%
    select(country, age, perc) %>%
    spread(age, perc) %>%
    rename(`15-24y` = `1`, `40+y` = `2`)

plot.opt %>% head(10) %>% knitr::kable()
```

```{r}
plot.opt %>%
    ggplot() + aes(`40+y`, `15-24y`, label = country) +
    geom_point(color = "skyblue", size = 2) +
    xlim(0, 1) + ylim(0,1) +
    geom_text(hjust = -0.1, vjust = -0.5) +
    geom_abline(intercept = 0, slop = 1,
                color="lightgrey", alpha=0.5, linetype="dashed") +
    theme_minimal() +
    theme(aspect.ratio=1)

```