EnviroPRA-package {EnviroPRA}R Documentation

Environmental Probabilistic Risk Assessment Tools

Description

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

Details

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 ~~

Author(s)

F. Barrio-Parra, with contributions from A. Dominguez-Castillo.

Maintainer: Fernando Barrio-Parra <fernando.barrio@upm.es>

Examples


#### 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))


[Package EnviroPRA version 1.0 Index]