Get risk total
get_risk_total.RdCalcul the total illness probability and total dalys that are the sum of every exposure event
Examples
library(dplyr)
library(purrr)
scenario_test <- create_scenario(system.file("input_1culture_2pop.xlsx", package = "ambre"))
scenario_conc <- inflow_concentration(scenario = scenario_test,
pathogenName = c("Campylobacter jejuni"))
#> Simulated pathogen: Campylobacter jejuni
#> Create 1000 random distribution(s): uniform (n: 60, min: 100.000000, max: 5000.000000)
#> Providing inflow events ... ok. (0.00 secs)
#> Providing inflow paras ... ok. (0.00 secs)
#> Simulated pathogen: Campylobacter jejuni
#> Create 1000 random distribution(s): uniform (n: 48, min: 100.000000, max: 5000.000000)
#> Providing inflow events ... ok. (0.00 secs)
#> Providing inflow paras ... ok. (0.00 secs)
scenario_volume <-
scenario_conc |>
mutate(
volume = map(
.x = config,
.f = ~ simulate_exposure(config = .x)
)
)
#> Simulated exposure: volume per event
#> Distribution set from 'triangle' to 'uniform' because 'min' equals 'max'
#> Create 1000 random distribution(s): uniform (n: 60, min: 0.001000, max: 0.001000)
#> Simulated exposure: volume per event
#> Distribution set from 'triangle' to 'uniform' because 'min' equals 'max'
#> Create 1000 random distribution(s): uniform (n: 48, min: 0.001000, max: 0.001000)
scenario_dose_ini <- initial_dose_calculation(scenario_volume)
scenario_scheme <- update_treatment_scheme(scenario_dose_ini)
scenario_with_logreduc_and_co <- scenario_scheme |>
mutate( log_reduction = map(config, simulate_treatment))
#> Create 1000 random distribution(s): uniform (n: 60, min: 1.000000, max: 2.000000)
#> Create 1000 random distribution(s): uniform (n: 60, min: 1.000000, max: 3.000000)
#> Create 1000 random distribution(s): uniform (n: 60, min: 5.000000, max: 7.000000)
#> Simulated treatment: Q.1 - Activated Sludge for Bacteria
#> Simulated treatment: Q.2 - Maturation Pond for Bacteria
#> Simulated treatment: Q.6 - Chlorination for Bacteria
#> Simulated treatment: Q.1 - Activated Sludge for Bacteria
#> Simulated treatment: Q.2 - Maturation Pond for Bacteria
#> Simulated treatment: Q.6 - Chlorination for Bacteria
#> Joining with `by = join_by(TreatmentID)`
#> Joining with `by = join_by(TreatmentID)`
#> Create 1000 random distribution(s): uniform (n: 48, min: 1.000000, max: 2.000000)
#> Create 1000 random distribution(s): uniform (n: 48, min: 1.000000, max: 3.000000)
#> Create 1000 random distribution(s): uniform (n: 48, min: 5.000000, max: 7.000000)
#> Simulated treatment: Q.1 - Activated Sludge for Bacteria
#> Simulated treatment: Q.2 - Maturation Pond for Bacteria
#> Simulated treatment: Q.6 - Chlorination for Bacteria
#> Simulated treatment: Q.1 - Activated Sludge for Bacteria
#> Simulated treatment: Q.2 - Maturation Pond for Bacteria
#> Simulated treatment: Q.6 - Chlorination for Bacteria
#> Joining with `by = join_by(TreatmentID)`
#> Joining with `by = join_by(TreatmentID)`
scenario_final_dose_and_co_test <- final_dose_calculation(scenario_with_logreduc_and_co)
scenario_inf_proba_and_co <- infection_probability_calculation(scenario_final_dose_and_co_test)
scenario_illness_proba_and_co <- illness_probability_calculation(scenario_inf_proba_and_co)
scenario_dalys_and_co <- dalys_calculation(scenario_illness_proba_and_co)
get_risk_total(scenario_dalys_and_co)
#> # A tibble: 2 × 32
#> CropName Area PopulationName nb_population PathName STEPtreatmentName
#> <chr> <dbl> <chr> <dbl> <chr> <chr>
#> 1 Tomato 10 Irrigation staff 1 Ingestion o… Q.1 - Activated …
#> 2 Corn seed 35 Maintenance staff 1 Ingestion f… Q.1 - Activated …
#> # ℹ 26 more variables: CollectiveTreatmentName <chr>, InitialProcessName <chr>,
#> # SupplementaryProcessName <chr>, nb_day_decay <dbl>, config <list>,
#> # CropID <dbl>, CropHeight <dbl>, PopulationID <dbl>, PathID <dbl>,
#> # MatrixID <dbl>, STEPtreatmentID <list>, CollectiveTreatmentID <list>,
#> # InitialProcessID <list>, SupplementaryProcessID <list>,
#> # InitialTrainName <list>, InitialTrainID <list>,
#> # SupplementaryTrainName <list>, SupplementaryTrainID <list>, …