ma_dichotomous_fit: ma_dichotomous_fit - Fit a model averaged dichotomous BMD...

In ToxicR: Analyzing Toxicology Dose-Response Data

ma_dichotomous_fit - Fit a model averaged dichotomous BMD model.

Description

Fit a model averaged dichotomous BMD model.

Usage

ma_dichotomous_fit(
  D,
  Y,
  N,
  model_list = integer(0),
  fit_type = "laplace",
  BMD_TYPE = "extra",
  BMR = 0.1,
  point_p = 0.01,
  alpha = 0.05,
  samples = 21000,
  burnin = 1000
)

Arguments

D	doses matrix
Y	response matrix
N	number of replicates matrix
model_list	a list of configurations for the single models (priors, model type)
fit_type	the method used to fit (laplace, mle, or mcmc)
BMD_TYPE	BMD_TYPE specifies the type of benchmark dose analysis to be performed. For continuous models, there are four types of BMD definitions that are commonly used. - Standard deviation is the default option, but it can be explicitly specified with 'BMR_TYPE = "sd"' This definition defines the BMD as the dose associated with the mean/median changing a specified number of standard deviations from the mean at the control dose., i.e., it is the dose, BMD, that solves \mid f(dose)-f(0) \mid = BMR \times σ - Relative deviation can be specified with 'BMR_TYPE = "rel"'. This defines the BMD as the dose that changes the control mean/median a certain percentage from the background dose, i.e. it is the dose, BMD that solves \mid f(dose) - f(0) \mid = (1 \pm BMR) f(0) - Hybrid deviation can be specified with 'BMR_TYPE = "hybrid"'. This defines the BMD that changes the probability of an adverse event by a stated amount relitive to no exposure (i.e 0). That is, it is the dose, BMD, that solves \frac{Pr(X > x| dose) - Pr(X >x|0)}{Pr(X < x|0)} = BMR. For this definition, Pr(X < x|0) = 1 - Pr(X > X|0) = π_0, where 0 ≤q π_0 < 1 is defined by the user as "point_p," and it defaults to 0.01. Note: this discussion assumed increasing data. The fitter determines the direction of the data and inverts the probability statements for decreasing data. - Absolute deviation can be specified with 'BMR_TYPE="abs"'. This defines the BMD as an absolute change from the control dose of zero by a specified amount. That is the BMD is the dose that solves the equation \mid f(dose) - f(0) \mid = BMR
BMR	This option specifies the benchmark response BMR. The BMR is defined in relation to the BMD calculation requested (see BMD). By default, the "BMR = 0.1."
point_p	This option is only used for hybrid BMD calculations. It defines a probability that is the cutpoint for observations. It is the probability that observations have this probability, or less, of being observed at the background dose.
alpha	Alpha is the specified nominal coverage rate for computation of the lower bound on the BMDL and BMDU, i.e., one computes a 100\times(1-α)\% . For the interval (BMDL,BMDU) this is a 100\times(1-2α)\% . By default, it is set to 0.05.
samples	the number of samples to take (MCMC only)
burnin	the number of burnin samples to take (MCMC only)

Value

a model object containing a list of single models

• Individual_Model_X: Here X is a number 1≤q X ≤q n, where n is the number of models in the model average. For each X, this is an individual model fit identical to what is returned in 'single_continuous_fit.'

• ma_bmd: The CDF of the model averaged BMD distribution.

• posterior_probs: The posterior model probabilities used in the MA.

• bmd: The BMD and the 100\times(1-2α)\% confidence intervals.

Examples

mData <- matrix(c(0, 2,50,
                  1, 2,50,
                  3, 10, 50,
                  16, 18,50,
                  32, 18,50,
                  33, 17,50),nrow=6,ncol=3,byrow=TRUE)
D <- mData[,1]
Y <- mData[,2]
N <- mData[,3]
model = ma_dichotomous_fit(D,Y,N)

summary(model)

ToxicR documentation built on Dec. 28, 2022, 3:07 a.m.