mc: Runs the Monte Carlo method on a baseline...

In eghuang/synergy: Tools for Synergy/Antagonism Modeling

Usage

mc(dose, LET, ratios, E, models, dE = NULL, n = 200, vcov = TRUE,
  interval_length = 0.95, seed = 100, check = TRUE, ...)

Arguments

dose	Numeric vector of all total dose values to be evaluated.
LET	Numeric vector of all LET values, must be length n.
ratios	Numeric vector of dose proportions applied on component DERs.
E	Vector of DER functions
models	List of nls models.
n	Numeric integer of the number of samples to be drawn.
vcov	Boolean for assessing inter-parameter correlation.
interval_length	Numeric double of the confidence interval width.
seed	Numeric value for pseudorandom generators.
...	Optional arguments to DER functions in E. Should be constant across all supplied E.

Details

Corresponding elements of ratios, LET should be associated with the same DER.

Value

Named list representing lower and upper bounds for a Monte Carlo confidence interval for a mixture DER over an interval of doses.

Author(s)

Yimin Lin, Edward Greg Huang <[e-mail]>

Examples

ion_data <- load_ion_data("one_ion.csv")
HZE_nte_der <- make_der(ion_data, TRUE, TRUE)
check_der(HZE_nte_der)
# Iron and silicon two-ion mixture.
mc(dose = 0:100, LET = c(193, 70), ratios = c(1/3, 2/3),
   E = rep(c(HZE_nte_der), 2),
   models = rep(list(HZE_nte_model), 2), n = 20, vcov = TRUE)

eghuang/synergy documentation built on Dec. 31, 2019, 8:55 a.m.