ProbitLogNormal-class {crmPack}R Documentation

Probit model with bivariate log normal prior

Description

This is probit regression model with a bivariate normal prior on the intercept and log slope. The covariate is the dose xx itself, potentially divided by a reference dose xx^{*}, or the logarithm of it:

Details

probit[p(x)]=α+βx/xprobit[p(x)] = \alpha + \beta \cdot x/x^{*}

or

probit[p(x)]=α+βlog(x/x)probit[p(x)] = \alpha + \beta \cdot \log(x/x^{*})

in case that the option useLogDose is TRUE. Here p(x)p(x) is the probability of observing a DLT for a given dose xx.

The prior is

(α,log(β))Normal(μ,Σ)(\alpha, \log(\beta)) \sim Normal(\mu, \Sigma)

The slots of this class contain the mean vector and the covariance matrix of the bivariate normal distribution, as well as the reference dose. Note that the parametrization inside the class uses alpha0 and alpha1.

This model is also used in the DualEndpoint classes, so this class can be used to check the prior assumptions on the dose-toxicity model - even when sampling from the prior distribution of the dual endpoint model is not possible.

Slots

mu

the prior mean vector μ\mu

Sigma

the prior covariance matrix Σ\Sigma

refDose

the reference dose xx^{*}

useLogDose

should the log of (standardized) dose be used?

Examples

model <- ProbitLogNormal(mu = c(-0.85, 1),
                           Sigma = matrix(c(1, -0.5, -0.5, 1), nrow = 2))

## we can also specify a reference dose, and use a log transformation of
## standardized dose in the model:
model <- ProbitLogNormal(mu = c(-0.85, 1),
                         Sigma = matrix(c(1, -0.5, -0.5, 1), nrow = 2),
                         refDose = 7.2,
                         useLogDose=TRUE)



[Package crmPack version 1.0.6 Index]