EnviroPRA-package {EnviroPRA} | R Documentation |
A collection of functions employed in environmental risk assessment to model exposure to a toxicant and predicting health effects, allowing to characterize variability and uncertainty in risk estimations
The DESCRIPTION file:
Package: | EnviroPRA |
Type: | Package |
Title: | Environmental Probabilistic Risk Assessment Tools |
Version: | 1.0 |
Date: | 2017-02-22 |
Author: | F. Barrio-Parra, with contributions from A. Dominguez-Castillo. |
Maintainer: | Fernando Barrio-Parra <fernando.barrio@upm.es> |
Description: | Methods to perform a Probabilistic Environmental Risk assessment from exposure to toxic substances - i.e. USEPA (1997) <https://www.epa.gov/risk/guiding-principles-monte-carlo-analysis> -. |
License: | GPL |
Imports: | MASS, kSamples, stats, fitdistrplus, truncdist |
Index of help topics:
AD Dermal conctact with chemicals in soil ADboot Dermal conctact with chemicals in soil by bootstrap AIR Inhalation of airborne chemicals AIRboot Inhalation of airborne chemicals by bootstrap DWIR Chemical intake by Drinking Water DWIRboot Chemical intake by Drinking Water by bootstrap EnviroPRA-package Environmental Probabilistic Risk Assessment Tools Fit_dist_parameter Returns adjusted distribution parameters HI Hazard Index HIdermal Hazard Index for dermal contact HIinhal Hazard Index for inhalation of vapors INH Inhalation of resuspended soil particles RISK Risk RISKInhal Risk for inhalation of vapors RISKdermal Risk for dermal contact SIR Chemical intake by accidental soil ingestion SIRboot Chemical intake by accidental soil ingestion by bootstrap VI Chemical intake by ingestion of vegetables VIboot Chemical intake by ingestion of vegetables by bootstrap condition p-value significance checking function extr_par Extracts the fitted distribution parameters to be introduced in other function fit_dist_test Summary of Godness-of-fit tests plot_fit_dist Graphical representation of data fitting to a distribution random_number_generator Random number generator sampler Execute sampling with replacement sig Significance level cheking function
~~ An overview of how to use the package, including the most important functions ~~
F. Barrio-Parra, with contributions from A. Dominguez-Castillo.
Maintainer: Fernando Barrio-Parra <fernando.barrio@upm.es>
#### Performs Deterministic Environmental Risk Assessment #####
# Example of dermal contact with a chemical in swiming water
# Estimate the dermal absorbed dose during swiming in waters with a carcinogenic chemical
# (water concentration of 250 mg/m^3)
DWIR ( CW = 250)
# For a systemic effect:
DWIR ( CW= 250, AT=24*365)
# Specifying all the parameters for the carcinogenic case
I = DWIR ( CW=250, IR=1.5, EF = 300, ED = 24, BW = 85)
# Chemical Slope factor
SFAs = 1.5
# Dermal Absorption Factor
ABSAs = 3e-02
# Gastrointestinal Absorption Factor
GIAs = 1
# Risk Estimation
RISKdermal (AD = I, SF = SFAs, GI = GIAs)
#### Perform a test to assess the fitness of a theorical distribution to empirical data ####
set.seed(123)
a <- rnorm(n=100, mean =1.5, sd = 0.25)
b <- rnorm(n = 15, mean = 300, sd = 15)
fit_dist_test(a)
fit_dist_test(b)
# Graphical representation of data fitting to a distribution
plot_fit_dist(a, "norm")
plot_fit_dist(b, "norm")
#### Perform a Probabilistic Environmental Risk Assessment ####
Fita <- Fit_dist_parameter(a)
Fitb <- Fit_dist_parameter(b)
IRr <-random_number_generator(n = 10000, Fited = Fita,
dist = "norm", a =0.8, b = 2.1)
EFr <-random_number_generator(n = 10000, Fited = Fitb,
dist = "norm", a =250, b = 330)
I = DWIR ( CW=250, IR=IRr, EF = EFr, ED = 24, BW = 85)
# Risk Estimation
Risk <- RISKdermal (AD = I, SF = SFAs, GI = GIAs)
hist (Risk)
quantile (Risk, c (0.05, 0.25, 0.5, 0.75, 0.95))