httkpop_mc {httk}

httk-pop: Correlated human physiological parameter Monte Carlo

Description

This is the core function for httk-pop correlated human physiological variability simulation as described by Ring et al. (2017) (doi:10.1016/j.envint.2017.06.004). This functions takes the data table of population biometrics (one individual per row) generated by httkpop_generate, and converts it to the corresponding table of HTTK model parameters for a specified HTTK model.

Usage

httkpop_mc(model, samples = 1000, httkpop.dt = NULL, ...)

Arguments

model

One of the HTTK models: "1compartment", "3compartmentss", "3compartment", or "pbtk".

samples

The number of Monte Carlo samples to use (can often think of these as separate individuals)

httkpop.dt

A data table generated by httkpop_generate. This defaults to NULL, in which case httkpop_generate is called to generate this table.

...

Additional arugments passed on to httkpop_generate.

Details

The Monte Carlo methods used here were recently updated and described by Breen et al. (submitted).

Value

A data.table with a row for each individual in the sample and a column for each parater in the model.

Author(s)

Caroline Ring and John Wambaugh

References

Ring CL, Pearce RG, Setzer RW, Wetmore BA, Wambaugh JF (2017). “Identifying populations sensitive to environmental chemicals by simulating toxicokinetic variability.” Environment International, 106, 105–118.

Breen M, Wambaugh JF, Bernstein A, Sfeir M, Ring CL (2022). “Simulating toxicokinetic variability to identify susceptible and highly exposed populations.” Journal of Exposure Science & Environmental Epidemiology, 32(6), 855–863.

Rowland M, Benet LZ, Graham GG (1973). “Clearance concepts in pharmacokinetics.” Journal of pharmacokinetics and biopharmaceutics, 1(2), 123–136.

Examples

set.seed(42)
indiv_examp <- httkpop_generate(method="d", nsamp=10)

httk_param <- httkpop_mc(httkpop.dt=indiv_examp, 
                        samples=10,
                        model="1compartment")