| Title: | Quantitative Assessment of Zoonotic Disease Risk |
|---|---|
| Description: | Provides quantitative tools for assessing zoonotic disease risk across animal, human, environmental, and transmission interfaces. The package supports exposure and transmission risk estimation, spillover risk assessment, risk scoring, cross-species comparison, transmission-network analysis, Monte Carlo uncertainty simulation, and sensitivity analysis. The One Health framework underlying these assessments is described by World Health Organization, Food and Agriculture Organization of the United Nations, United Nations Environment Programme, and World Organisation for Animal Health (2022) <doi:10.4060/cc2289en>. |
| Authors: | Vinodhkumar Obli Rajendran [aut, cre], Keerthi Aaradhana [aut] |
| Maintainer: | Vinodhkumar Obli Rajendran <[email protected]> |
| License: | MIT + file LICENSE |
| Version: | 1.0.0 |
| Built: | 2026-09-24 15:26:30 UTC |
| Source: | https://github.com/cran/ZooRisk |
Compare risk across host species
cross_species(species, risk)cross_species(species, risk)
species |
Character species names. |
risk |
Numeric risk values between 0 and 1. |
A ranked data frame.
Exposure risk
exposure_risk(prevalence, exposure)exposure_risk(prevalence, exposure)
prevalence |
Probability or prevalence of pathogen occurrence. |
exposure |
Probability of exposure. |
Numeric exposure-associated risk.
Calculate host-specific risk
host_risk(prevalence, exposure)host_risk(prevalence, exposure)
prevalence |
Prevalence. |
exposure |
Exposure probability. |
Numeric host risk.
Summarizes the distribution of zoonotic risk obtained from
risk_simulate().
monte_carlo_risk(x, probs = c(0.025, 0.5, 0.975))monte_carlo_risk(x, probs = c(0.025, 0.5, 0.975))
x |
A data frame returned by |
probs |
Probability levels used to calculate the lower and upper quantiles. |
A data frame containing the mean, median, lower, and upper risk estimates.
sim <- risk_simulate( prevalence = c(0.1, 0.3), exposure = c(0.2, 0.6), transmission = c(0.3, 0.8), n_sim = 100 ) monte_carlo_risk(sim)sim <- risk_simulate( prevalence = c(0.1, 0.3), exposure = c(0.2, 0.6), transmission = c(0.3, 0.8), n_sim = 100 ) monte_carlo_risk(sim)
Calculates node-level centrality measures for a transmission network.
When the optional igraph package is available, degree,
betweenness, and closeness centrality are calculated. Otherwise,
degree is calculated from the network structure and the other
measures are returned as NA.
network_centrality(x)network_centrality(x)
x |
A transmission network created by |
A data frame containing node-level network centrality measures.
Calculates node-level connectivity in a transmission network using in-degree, out-degree, and total degree. Higher total degree indicates greater network connectivity and does not itself represent epidemiological risk or probability of transmission.
network_risk(x)network_risk(x)
x |
A transmission network created by |
A data frame containing the node identifier, out-degree, in-degree, and total degree, ordered by decreasing total degree.
Summarize One Health risk components
onehealth_summary(animal, human, environment)onehealth_summary(animal, human, environment)
animal, human, environment
|
Numeric risk components between 0 and 1. |
A data frame.
Plots a transmission network using igraph when available. Otherwise, a base R degree plot is produced.
plot_network(x, ...)plot_network(x, ...)
x |
A transmission network created by |
... |
Graphical arguments passed to the plotting method. |
Invisibly returns x.
Plot a zoonotic risk score
plot_risk(x, ...)plot_risk(x, ...)
x |
Numeric risk score or |
... |
Graphical arguments. |
Invisibly returns x.
Plot a two-dimensional risk matrix
plot_risk_matrix(x, y = NULL, ...)plot_risk_matrix(x, y = NULL, ...)
x |
Numeric risk values. |
y |
Optional second risk dimension. |
... |
Graphical arguments. |
Invisibly returns x.
Prints a quantitative zoonotic risk result, including the overall multiplicative risk estimate and its component probabilities.
## S3 method for class 'zoorisk' print(x, ...)## S3 method for class 'zoorisk' print(x, ...)
x |
A |
... |
Ignored. |
Invisibly returns the original zoorisk object of class
zoorisk. The object contains the overall risk estimate,
component probabilities, and the model specification.
Calculate component contributions
risk_contribution(probabilities, weights = NULL)risk_contribution(probabilities, weights = NULL)
probabilities |
Numeric component probabilities. |
weights |
Optional weights. |
Data frame of normalized contributions.
Classify risk into qualitative categories
risk_matrix(x, breaks = c(0, 0.25, 0.5, 0.75, 1))risk_matrix(x, breaks = c(0, 0.25, 0.5, 0.75, 1))
x |
Numeric risk scores between 0 and 1. |
breaks |
Five increasing cut points from 0 to 1. |
A factor of risk categories.
Rank zoonotic risks
risk_rank(x, decreasing = TRUE)risk_rank(x, decreasing = TRUE)
x |
Numeric risk scores. |
decreasing |
Sort from highest to lowest when TRUE. |
Integer ranks.
Calculate a weighted zoonotic risk score
risk_score(probabilities, weights = NULL)risk_score(probabilities, weights = NULL)
probabilities |
Numeric probabilities between 0 and 1. |
weights |
Optional non-negative weights. |
A numeric risk score between 0 and 1.
Generates random combinations of zoonotic risk components within user-specified bounds and calculates the resulting risk.
risk_simulate( prevalence, exposure, transmission, susceptibility = c(1, 1), n_sim = 10000, seed = NULL )risk_simulate( prevalence, exposure, transmission, susceptibility = c(1, 1), n_sim = 10000, seed = NULL )
prevalence |
Numeric vector of length two giving the lower and upper bounds for prevalence. |
exposure |
Numeric vector of length two giving the lower and upper bounds for exposure. |
transmission |
Numeric vector of length two giving the lower and upper bounds for transmission. |
susceptibility |
Numeric vector of length two giving the lower and upper bounds for susceptibility. |
n_sim |
Number of simulations. |
seed |
Optional random seed. |
A data frame containing simulated prevalence, exposure, transmission, susceptibility, and zoonotic risk values.
set.seed(123) sim <- risk_simulate( prevalence = c(0.1, 0.3), exposure = c(0.2, 0.6), transmission = c(0.3, 0.8), susceptibility = c(0.5, 1), n_sim = 100 ) head(sim)set.seed(123) sim <- risk_simulate( prevalence = c(0.1, 0.3), exposure = c(0.2, 0.6), transmission = c(0.3, 0.8), susceptibility = c(0.5, 1), n_sim = 100 ) head(sim)
Evaluates the effect of fractional changes to each risk component while holding the other components at their baseline values.
sensitivity_risk(base, changes = c(-0.2, 0.2))sensitivity_risk(base, changes = c(-0.2, 0.2))
base |
Named numeric vector of baseline risk components. Values must be between 0 and 1. |
changes |
Numeric vector of fractional changes applied to each component individually. |
A data frame containing the component and simulated risk
values for each specified change. The baseline risk is returned
as the baseline_risk attribute.
base <- c( prevalence = 0.2, exposure = 0.5, transmission = 0.6, susceptibility = 0.8 ) sensitivity_risk(base)base <- c( prevalence = 0.2, exposure = 0.5, transmission = 0.6, susceptibility = 0.8 ) sensitivity_risk(base)
Spillover risk
spillover_risk(prevalence, exposure, transmission, susceptibility = 1)spillover_risk(prevalence, exposure, transmission, susceptibility = 1)
prevalence |
Probability or prevalence in the reservoir. |
exposure |
Probability of relevant exposure. |
transmission |
Probability of transmission given exposure. |
susceptibility |
Probability of infection given transmission. |
Numeric spillover risk.
Summarize a ZooRisk result
## S3 method for class 'zoorisk' summary(object, ...)## S3 method for class 'zoorisk' summary(object, ...)
object |
A |
... |
Ignored. |
A data frame.
Constructs an edge-list representation of a transmission network.
transmission_network(from, to, directed = TRUE)transmission_network(from, to, directed = TRUE)
from |
Source node identifiers. |
to |
Destination node identifiers. |
directed |
Logical; whether edges are directed. |
An edge-list data frame of class zoorisk_network.
Transmission risk
transmission_risk(exposure, transmission)transmission_risk(exposure, transmission)
exposure |
Probability of exposure. |
transmission |
Probability of transmission given exposure. |
Numeric transmission-associated risk.
Calculates a transparent multiplicative risk estimate.
zoonotic_risk(prevalence, exposure, transmission, susceptibility)zoonotic_risk(prevalence, exposure, transmission, susceptibility)
prevalence |
Probability or prevalence of pathogen occurrence. |
exposure |
Probability of relevant exposure. |
transmission |
Probability of transmission given exposure. |
susceptibility |
Probability of infection given transmission. |
An object of class zoorisk.