pr.dose.radir {radir} | R Documentation |
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 |
lod |
Lower dose considered. Its default value is 0. |
upd |
Upper dose considered. Its default value is the maximum dose in |
Value
The probability that the real dose is between lod
and upd
.
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 <david.morina@uab.cat>
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
### 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)