Calculate internal chemical dose — calc_internal

Estimate the internal dose from inhalation of a chemical given inhalation rate, time, and body weight

Usage

calc_internal_dose(C_ext, IR, time = 1, BW = 1, scaling = 1)

Arguments

C_ext: ambient chemical concentration in $\frac{mg}{m^3}$
IR: inhalation rate in $\frac{m^3}{day}$
time: total time in $days$
BW: body weight in $kg$
scaling: scaling factor encompassing any required unit adjustments

Value

internal chemical dose in $\frac{mg}{kg}$

Details

TODO Additional details... $$D_{int} = \frac{C_{ext} \,\times\, IR \,\times\, time}{BW}$$

Examples

n_chem <- 3
n_sample <- 5

# Single population
C_ext <- matrix(runif(n_sample * n_chem), ncol = n_chem)
IR <- runif(n_sample)
calc_internal_dose(C_ext, IR)
#>             [,1]       [,2]        [,3]
#> [1,] 0.015800366 0.09125914 0.171132967
#> [2,] 0.336685183 0.20087215 0.070595211
#> [3,] 0.038245313 0.01844701 0.002179853
#> [4,] 0.061107128 0.28487219 0.124534354
#> [5,] 0.007218509 0.75363169 0.392490442

# Multiple populations
C_ext <- list(
  "a" = matrix(runif(n_sample * n_chem), ncol = n_chem),
  "b" = matrix(runif(n_sample * n_chem), ncol = n_chem)
)
IR <- list(runif(n_sample), runif(n_sample))
calc_internal_dose(C_ext, IR)
#> $a
#>           [,1]       [,2]       [,3]
#> [1,] 0.1336944 0.34197969 0.01440300
#> [2,] 0.2138538 0.01621801 0.07110673
#> [3,] 0.1284426 0.09261534 0.05254789
#> [4,] 0.1041654 0.36988156 0.33832828
#> [5,] 0.4840307 0.33989673 0.23646425
#> 
#> $b
#>            [,1]       [,2]      [,3]
#> [1,] 0.33679479 0.53005683 0.5966388
#> [2,] 0.14423849 0.10185347 0.1571509
#> [3,] 0.67671195 0.45777230 0.7067715
#> [4,] 0.01176099 0.04306118 0.0632936
#> [5,] 0.07682741 0.03401241 0.1378512
#>