ci.dose.radir {radir}

Credible intervals for radiation dose

Description

The function allows the user to compute credible intervals for radiation doses objects of class radir.

Usage

ci.dose.radir(object, cr=0.95)

Arguments

object

the doses estimated by dose.distr function, an object of class radir.

cr

size of the credibility region. Its default value is 0.95.

Value

A vector with two elements containing the lower and upper bounds of the credible region.

Author(s)

David Moriña (Barcelona Graduate School of Mathematics), Manuel Higueras (Basque Center for Applied Mathematics) and Pedro Puig (Universitat Autònoma de Barcelona)

Mantainer: David Moriña Soler <[e-mail]>

References

Higueras M, Puig P, Ainsbury EA, Rothkamm K. A new inverse regression model applied to radiation biodosimetry. Proc R Soc A 2015;471, http://dx.doi.org/10.1098/rspa.2014.0588

See Also

radir-package, dose.distr, pr.dose.radir

Examples

### Example 3 (a)
f <- expression(b1*x+b2*x^2)
pars <- c("b1","b2")
beta <- c(3.126e-3, 2.537e-2)
cov  <- matrix(c(7.205e-06,-3.438e-06,-3.438e-06,2.718e-06),nrow=2)

### (a)
ex1.a <- dose.distr(f, pars, beta, cov, cells=1811, dics=102, 
m.prior="normal", d.prior="uniform", prior.param=c(0, Inf))

ci.dose.radir(ex1.a, 0.90)