radir (version 1.0.4)

ci.dose.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.

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

# NOT RUN {
### 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)
# }