manuscript/R/3-cow-leader-year-initiation.R

In svmiller/peacesciencer: Tools and Data for Quantitative Peace Science Research

library(tidyverse)
library(peacesciencer)

create_leaderyears(standardize = "cow", subset_years=c(1875:2010)) %>%
  add_gml_mids() %>%
  add_spells() %>%
  add_lead() %>%
  mutate(milservice_nocombat = ifelse(milservice == 1 & combat == 0, 1, 0)) %>%
  add_nmc() %>%
  add_democracy() -> Data

modLY<- list()
modLY[[1]] <- glm(gmlmidonset_init ~ gender + leaderage + yrinoffice + milservice_nocombat + combat + rebel +
                          warwin + warloss + rebelwin + rebelloss + xm_qudsest + cinc +
                          gmlmidinitspell + I(gmlmidinitspell^2) + I(gmlmidinitspell^3), data = Data,
                    family = binomial(link = "logit"))

saveRDS(modLY, "data/modLY.rds")