# single year
library(dplyr)
m <- korea_gravity %>%
filter(year == 2020,
orig != dest) %>%
select(orig, dest, flow)
m
p <- korea_gravity %>%
filter(year == 2020) %>%
distinct(dest, dest_pop)
p
index_impact(m = m, p = p, pop_col = "dest_pop", reg_col = "dest")
# multiple years
library(tidyr)
library(purrr)
korea_gravity %>%
select(year, orig, dest, flow, dest_pop) %>%
group_nest(year) %>%
mutate(m = map(.x = data, .f = ~select(.x, orig, dest, flow)),
p = map(.x = data, .f = ~distinct(.x, dest, dest_pop)),
i = map2(.x = m, .y = p,
.f = ~index_impact(
m = .x, p = .y, pop_col = "dest_pop", reg_col = "dest", long = FALSE
))) %>%
select(-data, -m, -p) %>%
unnest(i)
Run the code above in your browser using DataLab