312final/data_exploration.R at master · susantmtran/312final · 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
library(tidyverse)

uni_data <- read_excel("uni_data.xlsx") %>% clean_names()

ggplot(uni_data, aes(end_endow)) +
  geom_histogram(aes(y = ..density..), bins = 50) +
  stat_function(fun = dnorm,
                args = list(mean = mean(uni_data$end_endow),
                            sd = sd(uni_data$end_endow)),
                color = "pink",
                size = 1) +
  theme_light() +
  labs(x = "endowment at the end of fy19 ($)",
       y = "") +
  theme(axis.text.y = element_blank(),
        axis.ticks.y = element_blank())

ggplot(uni_data, aes(m_index)) +
  geom_histogram(aes(y = ..density..)) +
  stat_function(fun = dnorm,
                args = list(mean = mean(uni_data$m_index),
                            sd = sd(uni_data$m_index)),
                color = "pink",
                size = 1) +
  theme_light() +
  labs(x = "mutual information index",
       y = "") +
  theme(axis.text.y = element_blank(),
        axis.ticks.y = element_blank())

ggplot(uni_data, aes(h_index)) +
  geom_histogram(aes(y = ..density..)) +
  stat_function(fun = dnorm,
                args = list(mean = mean(uni_data$h_index),
                            sd = sd(uni_data$h_index)),
                color = "pink",
                size = 1) +
  theme_light() +
  labs(x = "theil's information index",
       y = "") +
  theme(axis.text.y = element_blank(),
        axis.ticks.y = element_blank())

ggplot(uni_data, aes(end_endow, h_index)) +
  geom_jitter() +
  geom_smooth(method = "lm") +
  theme_light() +
  labs(x = "endowment at the end of fy19 ($)",
       y = "theil's information index")