radir (version 1.0.4)

pr.dose.radir: Probability between doses

Description

This function allows the user to compute the probability between two radiation doses.

Usage

pr.dose.radir(object, lod = 0, upd = object[[2]][length(object[[2]])])

Arguments

object

An object of class radir containing the estimated doses.

lod

Lower dose considered. Its default value is 0.

upd

Upper dose considered. Its default value is the maximum dose in object.

Value

The probability that the real dose is between lod and upd.

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, ci.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))

pr.dose.radir(ex1.a, 1, 1.4)
# }