inst/scripts/fpsim-votes-a.R

In fpsim: Compute Measures of Foreign Policy Similarity/Agreement

library(tidyverse)
library(fpsim)
library(peacesciencer)
library(foreach)
library(tictoc)
library(qs2)

half_cores  <- parallel::detectCores()/2
my.cluster <- parallel::makeCluster(
  half_cores,
  type = "PSOCK"
)

doParallel::registerDoParallel(cl = half_cores)
foreach::getDoParRegistered()

tic()
FPSIMVA <- foreach(
  #y = c(1946:1955) # debugging here...
  y = c(1946:1963, 1965:2022)
) %dopar% {

  print(paste("Starting", y, "on", Sys.time()))

  yearlyvotes <- qs_read(paste0("data-raw/un/y",y,".qs")) %>% select(-date)

  yearlyvotes %>%
    slice(1, .by=c(resid, ccode)) %>%
    group_split(resid) %>%
    map(~expand(., resid = resid, ccode1 = ccode, ccode2 = ccode)) %>%
    map(~left_join(., yearlyvotes,
                   by = c("ccode1" = "ccode",
                          "resid" = "resid"))) %>%
    map(~rename(., vote1 = vote)) %>%
    map(~left_join(., yearlyvotes,
                   by = c("ccode2" = "ccode",
                          "resid" = "resid",
                          "year" = "year"))) %>%
    map(~rename(., vote2 = vote)) %>%
    bind_rows() %>%
    filter(ccode1 < ccode2) -> UNDDY


  UNDDY %>%
    split(., paste(.$ccode1, .$ccode2, sep = "_")) -> yearsplits


  yearsplits %>%
    map(~bcai(.$vote1, .$vote2, distances = 'absolute', levels = 1:3)) %>%
    enframe(., name = "dyad", value = "avoteva") %>%
    unnest(avoteva) -> AVA

  yearsplits %>%
    map(~bcai(.$vote1, .$vote2, distances = 'squared', levels = 1:3)) %>%
    enframe(., name = "dyad", value = "avotevs") %>%
    unnest(avotevs) -> AVS

  full_join(AVA, AVS, by=c("dyad" = "dyad")) %>%
    separate(dyad, c("ccode1", "ccode2")) %>%
    mutate(ccode1 = as.numeric(ccode1),
           ccode2 = as.numeric(ccode2)) %>%
    mutate(year = y) %>%
    select(ccode1, ccode2, year, everything())  -> here_it_is


  print(paste("Ending", y, "on", Sys.time()))
  # ^ definitely don't end with this... Steve... okay...

  here_it_is

}

toc(log = TRUE) # and, time
parallel::stopCluster(cl = my.cluster) # close our clusters
rm(my.cluster)


#qs_save(FPSIMVPK, "docs/data/FPSIMVPK.qs")

FPSIMVA %>%
  bind_rows() %>%
  filter(ccode1 != ccode2) -> FPSIMVA

qs_save(FPSIMVA, "docs/data/fpsim-votes-a.qs")
saveRDS(FPSIMVA, "docs/data/fpsim-votes-a.rds")


sink(file = "inst/scripts/fpsim-votes-a.log")
timestamp()
tic.log()
sink()