mean_uncertainUCDP {uncertainUCDP} | R Documentation |
Mean, median, and quantiles of the parametric uncertainty distributions for UCDP events
Description
Mean, median, and quantiles of the parametric uncertainty distributions for UCDP events. The parametric uncertainty distributions are based on the reported-value inflation Gumbel mixture distribution. The median
and quantile
functions are shortcuts for the quncertainUCDP
function.
Usage
mean_uncertainUCDP(fatalities, tov = c("sb", "ns", "os", "any"))
median_uncertainUCDP(fatalities, tov = c("sb", "ns", "os", "any"))
quantiles_unceartainUCDP(probs, fatalities, tov = c("sb", "ns", "os", "any"))
Arguments
fatalities |
A vector of non-negative integers representing the number of fatalities of the UCDP events. Non-integer values are allowed but should be considered experimental. |
tov |
A character string representing the type of violence of the UCDP. Must be one of "sb", "ns", "os", or "any". The options are: * "sb" for state-based violence * "ns" for non-state violence * "os" for one-sided violence * "any" for parameters estimated across all type of violence. This is somewhat experimental and should be used with caution. This is possibly useful when the type of violence is unknown or when the user wants to combine all types of violence into a single category. |
probs |
A numeric vector of probabilities with values in [0,1]. The quantiles to calculate. |
Value
A numeric vector of the same length as the input vector of fatalities representing the means, medians, and quantiles of the parametric uncertainty distribution for each UCDP event.
Examples
data(ucdpged)
# Calculate the mean for an arbitrary UCDP event
mean_uncertainUCDP(fatalities = 100, tov = 'sb')
# Calculate the mean for the first event in the UCDP GED sample
mean_uncertainUCDP(ucdpged$best[1], tov = ucdpged$type_of_violence[1])
# Calculate the median for an arbitrary UCDP event
median_uncertainUCDP(fatalities = 100, tov = 'sb')
# Calculate the median for the first event in the UCDP GED sample
median_uncertainUCDP(ucdpged$best[1], tov = ucdpged$type_of_violence[1])
# Calculate the 90th percentile for an arbitrary UCDP event
quantiles_unceartainUCDP(probs = 0.9, fatalities = 100, tov = 'sb')
# Calculate the 90th percentile for the first event in the UCDP GED sample
quantiles_unceartainUCDP(ucdpged$best[1], 0.9, tov = ucdpged$type_of_violence[1])