R/add_gml_mids.R
In svmiller/peacesciencer: Tools and Data for Quantitative Peace Science Research
Defines functions
add_gml_mids
Documented in
add_gml_mids
#' Add Gibler-Miller-Little (GML) Militarized Interstate Dispute (MID) data to a
#' data frame
#'
#' @description
#'
#' `r lifecycle::badge("deprecated")`
#'
#' \code{add_gml_mids()} merged in GML's MID data to a (dyad-year, leader-year,
#' leader-dyad-year, state-year) data frame. The current version of the GML MID
#' data is 2.2.1. However, this function is deprecated at the request of one of
#' the data set's maintainers, insisting that the use of the data be done
#' manually with design decisions made 100% at the discretion of the researcher.
#' The function now returns a stop communicating this development.
#'
#' @return
#'
#' \code{add_gml_mids()} now returns a stop communicating one of the maintainer's
#' requests to reject all software that facilitates the use of the data in this
#' fashion. Before, it took a (dyad-year, leader-year, leader-dyad-year,
#' state-year) data frame and adds dispute information from the GML MID data.
#' If the data are dyad-year, the return is a laundry list of information about
#' onsets, ongoing conflicts, and assorted participant- and dispute-level
#' summaries. If the data are leader-dyad-year, these are carefully matched to
#' leaders as well. If the data are state-year or leader-year, the function
#' returns information about ongoing disputes (and onsets) and whether there
#' were any ongoing disputes (and onsets) the state (or leader) initiated.
#'
#' @details
#'
#' Dyads are capable of having multiple disputes in a given year, which can
#' create a problem for merging into a complete dyad-year data frame.
#' Consider the case of France and Italy in 1860, which
#' had three separate dispute onsets that year (MID#0112, MID#0113,
#' MID#0306), as illustrative of the problem. This merging process employs
#' several rules to whittle down these duplicate dyad-years for merging
#' into a dyad-year data frame.
#'
#' The function will also return a message to the user about the
#' case-exclusion rules that went into this process. Users who are
#' interested in implementing their own case-exclusion rules should
#' look up the "whittle" class of functions also provided in this
#' package.
#'
#' Determining "initiation" for state-year summaries of inter-state disputes
#' is possible since there is an implied directionality of "initiation."
#' In about half of all cases, this is straightforward. You can use the
#' participant summaries and determine that if the dispute was bilateral
#' and the dispute did not escalate beyond an attack, the state on Side A
#' initiated the dispute. For multilateral MIDs, these conditions still hold
#' at least for originators. However, there is *considerable* difficulty for
#' cases where 1) participant-level summaries suggested actions at the level of
#' clash or higher, 2) the participant was a joiner and not an originator.
#' The effort required to flesh this out is enormous, and perhaps
#' forthcoming in a future update.
#'
#' \code{add_gml_mids()} allows you to make one of three judgment
#' calls here (see the arguments section of the documentation).
#' If it were my call to make, I would say you should probably use the option
#' \code{"sidea-all-joiners"}. My review of the MID data with Doug Gibler
#' suggests most states that join a dispute are not roped into a conflict
#' (i.e. targeted by some other state) after the first incident. They
#' routinely initiate their entry into the conflict, which is
#' what this concept of "initiation" is supposed to capture in the
#' literature. There are no doubt cases where some third state is brought into
#' the dispute by the actions of some other state even as the original MID
#' coding rules place a high barrier on coding that type of dispute entry.
#' However, the time required to individually assess whether a state initiated
#' their entry into a MID under something other than the simplest of cases
#' (e.g. bilateral cases where the highest participant action fell short of
#' a clash) would be too time-consuming. It would require an audit of almost
#' half of all participant-level summaries in the data. In a forthcoming
#' publication, Gibler and Miller offer excellent coverage here with a
#' new data set on militarized events. However, this would include only
#' confrontations after World War II.
#'
#' @author Steven V. Miller
#'
#' @param data a data frame with appropriate \pkg{peacesciencer} attributes
#' @param keep an optional parameter, specified as a character vector,
#' applicable to just the dyad-year data, and passed to the function in a
#' \code{select(one_of(.))} wrapper. This allows the user to discard unwanted
#' columns from the directed dispute data so that the output does not consume
#' too much space in memory. Note: the Correlates of War system codes
#' (\code{ccode1}, \code{ccode2}), the observation year (\code{year}), the
#' presence or absence of an ongoing MID (\code{gmlmidongoing}), and the
#' presence or absence of a unique MID onset (\code{gmlmidonset}) are *always*
#' returned. It would be foolish and self-defeating to eliminate those
#' observations. The user is free to keep or discard anything else they see fit.
#'
#' If \code{keep} is not specified in the function, the ensuing output returns
#' everything.
#'
#' @param init how should initiators be coded? Applicable only to state-year,
#' leader-dyad-year, and leader-year data. This parameter accepts one of
#' three possible values (\code{"sidea-orig"}, \code{"sidea-with-joiners"},
#' \code{"sidea-all-joiners"}). \code{"sidea-orig"} = a state initiates a MID
#' (which appears as a summary return in the output) if the state was on Side A
#' at the onset of the dispute. \code{"sidea-with-joiners"} = a state initiates
#' a MID (which appears as a summary return in the output) if the state was on
#' Side A at the onset of the dispute or if the state joined the MID on Side A.
#' \code{"sidea-all-joiners"} = a state initiates a MID (which appears as a
#' summary return in the output) if the state was on Side A at the onset of
#' the dispute or if it joined at any point thereafter. See details section for
#' more discussion. The default is \code{"sidea-all-joiners"}.
#'
#' @references
#'
#' Gibler, Douglas M., Steven V. Miller, and Erin K. Little. 2016. “An Analysis of the Militarized
#' Interstate Dispute (MID) Dataset, 1816-2001.” International Studies Quarterly 60(4): 719-730.
#'
#' @examples
#'
#' \dontrun{
#' # just call `library(tidyverse)` at the top of the your script
#' library(magrittr)
#' cow_ddy %>% add_gml_mids()
#'
#' # keep just the dispute number and Side A/B identifiers
#' cow_ddy %>% add_gml_mids(keep=c("dispnum","sidea1", "sidea2"))
#' }
#'
add_gml_mids <- function(data, keep, init = "sidea-all-joiners") {
if (init %nin% c("sidea-all-joiners", "sidea-with-joiners", "sidea-orig")) {
stop("init must be one of 'sidea-all-joiners', 'sidea-with-joiners', or 'sidea-orig'. You may have mispelled something. This argument applies only to state-year or leader-year analyses.")
}
if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "dyad_year") {
if (!all(i <- c("ccode1", "ccode2") %in% colnames(data))) {
stop("add_gml_mids() merges on two Correlates of War codes (ccode1, ccode2), which your data don't have right now. Make sure to run create_dyadyears() at the top of the pipe. You'll want the default option, which returns Correlates of War codes.")
} else {
if (nrow(data %>% filter(.data$ccode1 > .data$ccode2)) == 0) {
message("Dyadic data are non-directed and initiation variables make no sense in this context.")
if (missing(keep)) {
gml_mid_ddydisps %>% select(everything()) -> dirdisp
} else {
gml_mid_ddydisps %>% select(one_of("ccode1", "ccode2", "year", "gmlmidonset", "gmlmidongoing", keep)) -> dirdisp
}
dirdisp %>%
left_join(data, .) %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing"), ~ifelse(is.na(.) & between(.data$year, 1816, 2010), 0, .)) -> data
} else { # Assuming data are directed...
if (init == "sidea-orig") {
gml_mid_ddydisps %>%
mutate(init1 = ifelse(.data$sidea1 == 1 & .data$orig1 == 1, 1, 0),
init2 = ifelse(.data$sidea2 == 1 & .data$orig2 == 1, 1, 0)) -> hold_this
} else if (init == "sidea-with-joiners") {
gml_mid_ddydisps %>%
mutate(init1 = ifelse(.data$sidea1 == 1 | (.data$orig1 == 0 & .data$sidea1 == 1), 1, 0),
init2 = ifelse(.data$sidea2 == 1 | (.data$orig2 == 0 & .data$sidea2 == 1), 1, 0)) -> hold_this
} else if (init == "sidea-all-joiners") {
gml_mid_ddydisps %>%
mutate(init1 = ifelse(.data$sidea1 == 1 | (.data$orig1 == 0), 1, 0),
init2 = ifelse(.data$sidea2 == 1 | (.data$orig2 == 0), 1, 0) ) -> hold_this
}
if (missing(keep)) {
hold_this %>% select(everything()) -> dirdisp
} else {
hold_this %>% select(one_of("ccode1", "ccode2", "year",
"gmlmidonset", "gmlmidongoing",
"init1", "init2",
"sidea1", "sidea2", "orig1", "orig2",
keep)) -> dirdisp
}
dirdisp %>%
left_join(data, .) %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing"), ~ifelse(is.na(.) & between(.data$year, 1816, 2010), 0, .)) -> data
}
message("add_gml_mids() IMPORTANT MESSAGE: By default, this function whittles dispute-year data into dyad-year data by first selecting on unique onsets. Thereafter, where duplicates remain, it whittles dispute-year data into dyad-year data in the following order: 1) retaining highest `fatality`, 2) retaining highest `hostlev`, 3) retaining highest estimated `mindur`, 4) retaining highest estimated `maxdur`, 5) retaining reciprocated over non-reciprocated observations, 6) retaining the observation with the lowest start month, and, where duplicates still remained (and they don't), 7) forcibly dropping all duplicates for observations that are otherwise very similar.\nSee: http://svmiller.com/peacesciencer/articles/coerce-dispute-year-dyad-year.html")
return(data)
}
} else if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "state_year") {
if (!all(i <- c("ccode") %in% colnames(data))) {
stop("add_gml_mids() for leader-year data merges on Correlates of War state codes (ccode), which you don't have.")
}
gml_part %>%
filter(.data$allmiss_leader_start == 0 & .data$allmiss_leader_end == 0) %>%
rowwise() %>%
mutate(year = list(seq(.data$styear, .data$endyear)),
gmlmidonset = list(ifelse(.data$year == min(.data$year), 1, 0))) %>%
unnest(c(.data$year, .data$gmlmidonset)) %>%
mutate(gmlmidongoing = 1) -> hold_this
if (init == "sidea-orig") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & (.data$sidea == 1 & .data$orig == 1), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) %>%
select(.data$dispnum:.data$ccode, .data$year, .data$gmlmidonset:ncol(.)) -> hold_this
} else if (init == "sidea-with-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & (.data$sidea == 1 | ( .data$orig == 0 & .data$sidea == 1)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) %>%
select(.data$dispnum:.data$ccode, .data$year, .data$gmlmidonset:ncol(.)) -> hold_this
} else if (init == "sidea-all-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & (.data$sidea == 1 | ( .data$orig == 0)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) %>%
select(.data$dispnum:.data$ccode, .data$year, .data$gmlmidonset:ncol(.)) -> hold_this
}
# ,
# gmlmidongoing_dispnums = paste0(.data$dispnum, collapse = "; ")
hold_this %>%
group_by(.data$ccode, .data$year) %>%
summarize(gmlmidongoing = sum(.data$gmlmidongoing),
gmlmidonset = sum(.data$gmlmidonset),
gmlmidongoing_init = sum(.data$gmlmidongoing_init),
gmlmidonset_init = sum(.data$gmlmidonset_init)) %>%
mutate_at(vars(c(.data$gmlmidongoing, .data$gmlmidonset, .data$gmlmidongoing_init, .data$gmlmidonset_init)), ~ifelse(. >= 1, 1, 0)) %>%
ungroup() -> hold_this
data %>%
left_join(., hold_this) -> data
data %>%
mutate_at(vars(c(.data$gmlmidongoing, .data$gmlmidonset, .data$gmlmidongoing_init, .data$gmlmidonset_init)), ~ifelse(is.na(.), 0, .)) -> data
} else if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "leader_year") {
if (!all(i <- c("ccode") %in% colnames(data))) {
stop("add_gml_mids() for leader-year data merges on Correlates of War state codes (ccode), which you don't have.")
}
leaderdays <- create_leaderdays(standardize = "cow")
gml_part %>%
filter(.data$allmiss_leader_start == 0 & .data$allmiss_leader_end == 0) %>%
mutate(stdate = as.Date(paste0(.data$styear,"/", .data$stmon,"/", .data$dummy_stday)),
enddate = as.Date(paste0(.data$endyear,"/", .data$endmon,"/", .data$dummy_endday))) %>%
rowwise() %>%
mutate(date = list(seq(.data$stdate, .data$enddate, by = "1 day")),
gmlmidonset = list(ifelse(date == min(.data$date), 1, 0))) %>%
unnest(c(.data$date, .data$gmlmidonset)) %>%
mutate(gmlmidongoing = 1) %>%
select(.data$dispnum:.data$ccode, .data$obsid_start, .data$sidea, .data$orig, .data$date, .data$gmlmidongoing, .data$gmlmidonset) %>%
left_join(leaderdays, .) -> hold_this
if (init == "sidea-orig") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start == .data$obsid & (.data$sidea == 1 & .data$orig == 1), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-with-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start == .data$obsid & (.data$sidea == 1 | ( .data$orig == 0 & .data$sidea == 1 )), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-all-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start == .data$obsid & (.data$sidea == 1 | (.data$orig == 0)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
}
hold_this %>%
# left_join(leaderdays, .) %>%
mutate(year = .pshf_year(.data$date)) %>%
group_by(.data$obsid, .data$year) %>%
summarize(gmlmidongoing = sum(.data$gmlmidongoing, na.rm=T),
gmlmidonset = sum(.data$gmlmidonset, na.rm=T),
gmlmidongoing_init = sum(.data$gmlmidongoing_init, na.rm=T),
gmlmidonset_init = sum(.data$gmlmidonset_init, na.rm=T)) %>%
ungroup() -> hold_this
hold_this %>%
mutate_at(vars("gmlmidongoing", "gmlmidonset", "gmlmidongoing_init", "gmlmidonset_init"), ~ifelse(.data$year >= 2011, NA, .)) %>%
mutate_at(vars("gmlmidongoing", "gmlmidonset", "gmlmidongoing_init", "gmlmidonset_init"), ~ifelse(. >= 1, 1, 0)) -> hold_this
data %>%
left_join(., hold_this) -> data
} else if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "leader_dyad_year") {
if (!all(i <- c("ccode1", "ccode2") %in% colnames(data))) {
stop("add_gml_mids() merges on two Correlates of War codes (ccode1, ccode2), which your data don't have right now. Make sure to run create_leaderdyadyears(system='cow') at the top of the pipe.")
}
if (nrow(data %>% filter(.data$ccode1 > .data$ccode2)) == 0) {
message("The leader-dyadic data are non-directed and initiation variables make no sense in this context.")
data %>%
left_join(., gml_mid_ddlydisps) -> hold_this
hold_this %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing"), ~ifelse(is.na(.) & between(.data$year, 1816, 2010), 0, .)) -> data
} else { # Assuming the data are directed....
data %>%
left_join(., gml_mid_ddlydisps) -> hold_this
if (init == "sidea-orig") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start1 == .data$obsid1 & (.data$sidea1 == 1 & .data$orig1 == 1), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-with-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start1 == .data$obsid1 & (.data$sidea1 == 1 | ( .data$orig1 == 0 & .data$sidea1 == 1 )), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-all-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start1 == .data$obsid1 & (.data$sidea1 == 1 | (.data$orig1 == 0)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
}
hold_this %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing",
"gmlmidonset_init", "gmlmidongoing_init"), ~ifelse(is.na(.) & !(is.na(obsid1) & is.na(obsid2)) & between(.data$year, 1816, 2010), 0, .)) -> data
}
} else {
stop("add_gml_mids() requires a data/tibble with attributes$ps_data_type of state_year, leader_year, or dyad_year. Try running create_dyadyears(), create_leaderyears(), or create_stateyears() at the start of the pipe.")
}
return(data)
}
svmiller/peacesciencer documentation built on April 3, 2025, 1:15 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions add_gml_mids
Documented in add_gml_mids
#' Add Gibler-Miller-Little (GML) Militarized Interstate Dispute (MID) data to a
#' data frame
#'
#' @description
#'
#' `r lifecycle::badge("deprecated")`
#'
#' \code{add_gml_mids()} merged in GML's MID data to a (dyad-year, leader-year,
#' leader-dyad-year, state-year) data frame. The current version of the GML MID
#' data is 2.2.1. However, this function is deprecated at the request of one of
#' the data set's maintainers, insisting that the use of the data be done
#' manually with design decisions made 100% at the discretion of the researcher.
#' The function now returns a stop communicating this development.
#'
#' @return
#'
#' \code{add_gml_mids()} now returns a stop communicating one of the maintainer's
#' requests to reject all software that facilitates the use of the data in this
#' fashion. Before, it took a (dyad-year, leader-year, leader-dyad-year,
#' state-year) data frame and adds dispute information from the GML MID data.
#' If the data are dyad-year, the return is a laundry list of information about
#' onsets, ongoing conflicts, and assorted participant- and dispute-level
#' summaries. If the data are leader-dyad-year, these are carefully matched to
#' leaders as well. If the data are state-year or leader-year, the function
#' returns information about ongoing disputes (and onsets) and whether there
#' were any ongoing disputes (and onsets) the state (or leader) initiated.
#'
#' @details
#'
#' Dyads are capable of having multiple disputes in a given year, which can
#' create a problem for merging into a complete dyad-year data frame.
#' Consider the case of France and Italy in 1860, which
#' had three separate dispute onsets that year (MID#0112, MID#0113,
#' MID#0306), as illustrative of the problem. This merging process employs
#' several rules to whittle down these duplicate dyad-years for merging
#' into a dyad-year data frame.
#'
#' The function will also return a message to the user about the
#' case-exclusion rules that went into this process. Users who are
#' interested in implementing their own case-exclusion rules should
#' look up the "whittle" class of functions also provided in this
#' package.
#'
#' Determining "initiation" for state-year summaries of inter-state disputes
#' is possible since there is an implied directionality of "initiation."
#' In about half of all cases, this is straightforward. You can use the
#' participant summaries and determine that if the dispute was bilateral
#' and the dispute did not escalate beyond an attack, the state on Side A
#' initiated the dispute. For multilateral MIDs, these conditions still hold
#' at least for originators. However, there is *considerable* difficulty for
#' cases where 1) participant-level summaries suggested actions at the level of
#' clash or higher, 2) the participant was a joiner and not an originator.
#' The effort required to flesh this out is enormous, and perhaps
#' forthcoming in a future update.
#'
#' \code{add_gml_mids()} allows you to make one of three judgment
#' calls here (see the arguments section of the documentation).
#' If it were my call to make, I would say you should probably use the option
#' \code{"sidea-all-joiners"}. My review of the MID data with Doug Gibler
#' suggests most states that join a dispute are not roped into a conflict
#' (i.e. targeted by some other state) after the first incident. They
#' routinely initiate their entry into the conflict, which is
#' what this concept of "initiation" is supposed to capture in the
#' literature. There are no doubt cases where some third state is brought into
#' the dispute by the actions of some other state even as the original MID
#' coding rules place a high barrier on coding that type of dispute entry.
#' However, the time required to individually assess whether a state initiated
#' their entry into a MID under something other than the simplest of cases
#' (e.g. bilateral cases where the highest participant action fell short of
#' a clash) would be too time-consuming. It would require an audit of almost
#' half of all participant-level summaries in the data. In a forthcoming
#' publication, Gibler and Miller offer excellent coverage here with a
#' new data set on militarized events. However, this would include only
#' confrontations after World War II.
#'
#' @author Steven V. Miller
#'
#' @param data a data frame with appropriate \pkg{peacesciencer} attributes
#' @param keep an optional parameter, specified as a character vector,
#' applicable to just the dyad-year data, and passed to the function in a
#' \code{select(one_of(.))} wrapper. This allows the user to discard unwanted
#' columns from the directed dispute data so that the output does not consume
#' too much space in memory. Note: the Correlates of War system codes
#' (\code{ccode1}, \code{ccode2}), the observation year (\code{year}), the
#' presence or absence of an ongoing MID (\code{gmlmidongoing}), and the
#' presence or absence of a unique MID onset (\code{gmlmidonset}) are *always*
#' returned. It would be foolish and self-defeating to eliminate those
#' observations. The user is free to keep or discard anything else they see fit.
#'
#' If \code{keep} is not specified in the function, the ensuing output returns
#' everything.
#'
#' @param init how should initiators be coded? Applicable only to state-year,
#' leader-dyad-year, and leader-year data. This parameter accepts one of
#' three possible values (\code{"sidea-orig"}, \code{"sidea-with-joiners"},
#' \code{"sidea-all-joiners"}). \code{"sidea-orig"} = a state initiates a MID
#' (which appears as a summary return in the output) if the state was on Side A
#' at the onset of the dispute. \code{"sidea-with-joiners"} = a state initiates
#' a MID (which appears as a summary return in the output) if the state was on
#' Side A at the onset of the dispute or if the state joined the MID on Side A.
#' \code{"sidea-all-joiners"} = a state initiates a MID (which appears as a
#' summary return in the output) if the state was on Side A at the onset of
#' the dispute or if it joined at any point thereafter. See details section for
#' more discussion. The default is \code{"sidea-all-joiners"}.
#'
#' @references
#'
#' Gibler, Douglas M., Steven V. Miller, and Erin K. Little. 2016. “An Analysis of the Militarized
#' Interstate Dispute (MID) Dataset, 1816-2001.” International Studies Quarterly 60(4): 719-730.
#'
#' @examples
#'
#' \dontrun{
#' # just call `library(tidyverse)` at the top of the your script
#' library(magrittr)
#' cow_ddy %>% add_gml_mids()
#'
#' # keep just the dispute number and Side A/B identifiers
#' cow_ddy %>% add_gml_mids(keep=c("dispnum","sidea1", "sidea2"))
#' }
#'
add_gml_mids <- function(data, keep, init = "sidea-all-joiners") {
if (init %nin% c("sidea-all-joiners", "sidea-with-joiners", "sidea-orig")) {
stop("init must be one of 'sidea-all-joiners', 'sidea-with-joiners', or 'sidea-orig'. You may have mispelled something. This argument applies only to state-year or leader-year analyses.")
}
if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "dyad_year") {
if (!all(i <- c("ccode1", "ccode2") %in% colnames(data))) {
stop("add_gml_mids() merges on two Correlates of War codes (ccode1, ccode2), which your data don't have right now. Make sure to run create_dyadyears() at the top of the pipe. You'll want the default option, which returns Correlates of War codes.")
} else {
if (nrow(data %>% filter(.data$ccode1 > .data$ccode2)) == 0) {
message("Dyadic data are non-directed and initiation variables make no sense in this context.")
if (missing(keep)) {
gml_mid_ddydisps %>% select(everything()) -> dirdisp
} else {
gml_mid_ddydisps %>% select(one_of("ccode1", "ccode2", "year", "gmlmidonset", "gmlmidongoing", keep)) -> dirdisp
}
dirdisp %>%
left_join(data, .) %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing"), ~ifelse(is.na(.) & between(.data$year, 1816, 2010), 0, .)) -> data
} else { # Assuming data are directed...
if (init == "sidea-orig") {
gml_mid_ddydisps %>%
mutate(init1 = ifelse(.data$sidea1 == 1 & .data$orig1 == 1, 1, 0),
init2 = ifelse(.data$sidea2 == 1 & .data$orig2 == 1, 1, 0)) -> hold_this
} else if (init == "sidea-with-joiners") {
gml_mid_ddydisps %>%
mutate(init1 = ifelse(.data$sidea1 == 1 | (.data$orig1 == 0 & .data$sidea1 == 1), 1, 0),
init2 = ifelse(.data$sidea2 == 1 | (.data$orig2 == 0 & .data$sidea2 == 1), 1, 0)) -> hold_this
} else if (init == "sidea-all-joiners") {
gml_mid_ddydisps %>%
mutate(init1 = ifelse(.data$sidea1 == 1 | (.data$orig1 == 0), 1, 0),
init2 = ifelse(.data$sidea2 == 1 | (.data$orig2 == 0), 1, 0) ) -> hold_this
}
if (missing(keep)) {
hold_this %>% select(everything()) -> dirdisp
} else {
hold_this %>% select(one_of("ccode1", "ccode2", "year",
"gmlmidonset", "gmlmidongoing",
"init1", "init2",
"sidea1", "sidea2", "orig1", "orig2",
keep)) -> dirdisp
}
dirdisp %>%
left_join(data, .) %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing"), ~ifelse(is.na(.) & between(.data$year, 1816, 2010), 0, .)) -> data
}
message("add_gml_mids() IMPORTANT MESSAGE: By default, this function whittles dispute-year data into dyad-year data by first selecting on unique onsets. Thereafter, where duplicates remain, it whittles dispute-year data into dyad-year data in the following order: 1) retaining highest `fatality`, 2) retaining highest `hostlev`, 3) retaining highest estimated `mindur`, 4) retaining highest estimated `maxdur`, 5) retaining reciprocated over non-reciprocated observations, 6) retaining the observation with the lowest start month, and, where duplicates still remained (and they don't), 7) forcibly dropping all duplicates for observations that are otherwise very similar.\nSee: http://svmiller.com/peacesciencer/articles/coerce-dispute-year-dyad-year.html")
return(data)
}
} else if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "state_year") {
if (!all(i <- c("ccode") %in% colnames(data))) {
stop("add_gml_mids() for leader-year data merges on Correlates of War state codes (ccode), which you don't have.")
}
gml_part %>%
filter(.data$allmiss_leader_start == 0 & .data$allmiss_leader_end == 0) %>%
rowwise() %>%
mutate(year = list(seq(.data$styear, .data$endyear)),
gmlmidonset = list(ifelse(.data$year == min(.data$year), 1, 0))) %>%
unnest(c(.data$year, .data$gmlmidonset)) %>%
mutate(gmlmidongoing = 1) -> hold_this
if (init == "sidea-orig") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & (.data$sidea == 1 & .data$orig == 1), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) %>%
select(.data$dispnum:.data$ccode, .data$year, .data$gmlmidonset:ncol(.)) -> hold_this
} else if (init == "sidea-with-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & (.data$sidea == 1 | ( .data$orig == 0 & .data$sidea == 1)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) %>%
select(.data$dispnum:.data$ccode, .data$year, .data$gmlmidonset:ncol(.)) -> hold_this
} else if (init == "sidea-all-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & (.data$sidea == 1 | ( .data$orig == 0)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) %>%
select(.data$dispnum:.data$ccode, .data$year, .data$gmlmidonset:ncol(.)) -> hold_this
}
# ,
# gmlmidongoing_dispnums = paste0(.data$dispnum, collapse = "; ")
hold_this %>%
group_by(.data$ccode, .data$year) %>%
summarize(gmlmidongoing = sum(.data$gmlmidongoing),
gmlmidonset = sum(.data$gmlmidonset),
gmlmidongoing_init = sum(.data$gmlmidongoing_init),
gmlmidonset_init = sum(.data$gmlmidonset_init)) %>%
mutate_at(vars(c(.data$gmlmidongoing, .data$gmlmidonset, .data$gmlmidongoing_init, .data$gmlmidonset_init)), ~ifelse(. >= 1, 1, 0)) %>%
ungroup() -> hold_this
data %>%
left_join(., hold_this) -> data
data %>%
mutate_at(vars(c(.data$gmlmidongoing, .data$gmlmidonset, .data$gmlmidongoing_init, .data$gmlmidonset_init)), ~ifelse(is.na(.), 0, .)) -> data
} else if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "leader_year") {
if (!all(i <- c("ccode") %in% colnames(data))) {
stop("add_gml_mids() for leader-year data merges on Correlates of War state codes (ccode), which you don't have.")
}
leaderdays <- create_leaderdays(standardize = "cow")
gml_part %>%
filter(.data$allmiss_leader_start == 0 & .data$allmiss_leader_end == 0) %>%
mutate(stdate = as.Date(paste0(.data$styear,"/", .data$stmon,"/", .data$dummy_stday)),
enddate = as.Date(paste0(.data$endyear,"/", .data$endmon,"/", .data$dummy_endday))) %>%
rowwise() %>%
mutate(date = list(seq(.data$stdate, .data$enddate, by = "1 day")),
gmlmidonset = list(ifelse(date == min(.data$date), 1, 0))) %>%
unnest(c(.data$date, .data$gmlmidonset)) %>%
mutate(gmlmidongoing = 1) %>%
select(.data$dispnum:.data$ccode, .data$obsid_start, .data$sidea, .data$orig, .data$date, .data$gmlmidongoing, .data$gmlmidonset) %>%
left_join(leaderdays, .) -> hold_this
if (init == "sidea-orig") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start == .data$obsid & (.data$sidea == 1 & .data$orig == 1), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-with-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start == .data$obsid & (.data$sidea == 1 | ( .data$orig == 0 & .data$sidea == 1 )), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-all-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start == .data$obsid & (.data$sidea == 1 | (.data$orig == 0)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
}
hold_this %>%
# left_join(leaderdays, .) %>%
mutate(year = .pshf_year(.data$date)) %>%
group_by(.data$obsid, .data$year) %>%
summarize(gmlmidongoing = sum(.data$gmlmidongoing, na.rm=T),
gmlmidonset = sum(.data$gmlmidonset, na.rm=T),
gmlmidongoing_init = sum(.data$gmlmidongoing_init, na.rm=T),
gmlmidonset_init = sum(.data$gmlmidonset_init, na.rm=T)) %>%
ungroup() -> hold_this
hold_this %>%
mutate_at(vars("gmlmidongoing", "gmlmidonset", "gmlmidongoing_init", "gmlmidonset_init"), ~ifelse(.data$year >= 2011, NA, .)) %>%
mutate_at(vars("gmlmidongoing", "gmlmidonset", "gmlmidongoing_init", "gmlmidonset_init"), ~ifelse(. >= 1, 1, 0)) -> hold_this
data %>%
left_join(., hold_this) -> data
} else if (length(attributes(data)$ps_data_type) > 0 && attributes(data)$ps_data_type == "leader_dyad_year") {
if (!all(i <- c("ccode1", "ccode2") %in% colnames(data))) {
stop("add_gml_mids() merges on two Correlates of War codes (ccode1, ccode2), which your data don't have right now. Make sure to run create_leaderdyadyears(system='cow') at the top of the pipe.")
}
if (nrow(data %>% filter(.data$ccode1 > .data$ccode2)) == 0) {
message("The leader-dyadic data are non-directed and initiation variables make no sense in this context.")
data %>%
left_join(., gml_mid_ddlydisps) -> hold_this
hold_this %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing"), ~ifelse(is.na(.) & between(.data$year, 1816, 2010), 0, .)) -> data
} else { # Assuming the data are directed....
data %>%
left_join(., gml_mid_ddlydisps) -> hold_this
if (init == "sidea-orig") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start1 == .data$obsid1 & (.data$sidea1 == 1 & .data$orig1 == 1), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-with-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start1 == .data$obsid1 & (.data$sidea1 == 1 | ( .data$orig1 == 0 & .data$sidea1 == 1 )), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
} else if (init == "sidea-all-joiners") {
hold_this %>%
mutate(gmlmidongoing_init = ifelse(.data$gmlmidongoing == 1 & .data$obsid_start1 == .data$obsid1 & (.data$sidea1 == 1 | (.data$orig1 == 0)), 1, 0),
gmlmidonset_init = ifelse(.data$gmlmidonset == 1 & .data$gmlmidongoing_init == 1, 1, 0)) -> hold_this
}
hold_this %>%
mutate_at(vars("gmlmidonset", "gmlmidongoing",
"gmlmidonset_init", "gmlmidongoing_init"), ~ifelse(is.na(.) & !(is.na(obsid1) & is.na(obsid2)) & between(.data$year, 1816, 2010), 0, .)) -> data
}
} else {
stop("add_gml_mids() requires a data/tibble with attributes$ps_data_type of state_year, leader_year, or dyad_year. Try running create_dyadyears(), create_leaderyears(), or create_stateyears() at the start of the pipe.")
}
return(data)
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.