# single year
library(dplyr)
m <- korea_gravity %>%
filter(year == 2020,
orig != dest)
m
p <- korea_gravity %>%
filter(year == 2020) %>%
distinct(dest, dest_pop)
p
index_intensity(mig_total = sum(m$flow), pop_total = sum(p$dest_pop*1e6), n = nrow(p))
# multiple years
library(tidyr)
library(purrr)
mm <- korea_gravity %>%
filter(orig != dest) %>%
group_by(year) %>%
summarise(m = sum(flow))
mm
pp <- korea_gravity %>%
group_by(year) %>%
distinct(dest, dest_pop) %>%
summarise(p = sum(dest_pop)*1e6,
n = n_distinct(dest))
pp
library(purrr)
library(tidyr)
mm %>%
left_join(pp) %>%
mutate(i = pmap(
.l = list(m, p, n),
.f = ~index_intensity(mig_total = ..1, pop_total = ..2,n = ..3, long = FALSE)
)) %>%
unnest(cols = i)
Run the code above in your browser using DataLab