Draw parameters fixed effects with biomarker-quantity-dependent boosting

Adds biomarker quantity ceiling effects to the draw_parameters_fixed_fx function. Here, an individual’s realized biomarker boost is dependent on their biomarker quantity at the time of the exposure event.

Usage

draw_parameters_fixed_fx_biomarker_dep(
  i,
  t,
  x,
  b,
  demography,
  biomarker_states,
  model_pars,
  ...
)

Arguments

i

individual

t

time

x

exposure

b

biomarker

demography

demography information

biomarker_states

an array of true biomarker quantities for all individuals across all time steps and biomarkers

model_pars

tibble of biomarker (antibody) kinetics parameters with variables: 1) exposure_id: numeric exposure ID; 2) biomarker_id: numeric biomarker ID; 3) name: the character name of the parameter; 4) mean: numeric mean of this parameter distribution; 5) sd: the numeric standard deviation of the parameter distribution; 6) distribution: character description of the parameter distribution type (e.g., log-normal, normal)

...

Additional arguments

Value

A tibble with the simulated parameters for this exposure event

Examples

model_pars <- reformat_biomarker_map(example_model_pars_biphasic)
draw_parameters_fixed_fx_biomarker_dep(2,100,1,1,example_demography, 
example_biomarker_states_wide, model_pars)
#> # A tibble: 8 × 7
#>       i     t     x     b name                        value realized_value
#>   <dbl> <dbl> <dbl> <dbl> <chr>                       <dbl>          <dbl>
#> 1     2   100     1     1 boost_long                  0.5           0.0500
#> 2     2   100     1     1 boost_short                 0.25          0.0250
#> 3     2   100     1     1 wane_long                   0.01          0.01  
#> 4     2   100     1     1 wane_short                  0.033         0.033 
#> 5     2   100     1     1 biomarker_ceiling_threshold 1.7           1.7   
#> 6     2   100     1     1 biomarker_ceiling_gradient  0.529         0.529 
#> 7     2   100     1     1 biomarker_prot_midpoint     0.075         0.075 
#> 8     2   100     1     1 biomarker_prot_width        0.1           0.1