rdrm {drc} | R Documentation |
Simulation of a dose-response curve with user-specified dose values and error distribution.
rdrm(nosim, fct, mpar, xerror, xpar = 1, yerror = "rnorm", ypar = c(0, 1), onlyY = FALSE, seedVal = 20070527)
nosim |
numeric. The number of simulated curves to be returned. |
fct |
list. Any built-in function in drc or a list with similar components. |
mpar |
numeric. The model parameters to be supplied to fct . |
xerror |
numeric or character. The distribution for the dose values. |
xpar |
numeric vector supplying the parameter values defining the distribution for the dose values.
If xerror is a distribution then remember that the number of dose values also is part of this argument
(the first argument). |
yerror |
numeric or character. The error distribution for the response values. |
ypar |
numeric vector supplying the parameter values defining the error distribution for the response values. |
onlyY |
logical. If TRUE then only the response values are returned (useful in simulations). Otherwise both dose values and response values (and for binomial data also the weights) are returned. |
seedVal |
numeric. The seed to be used. |
The distribution for the dose values can either be a fixed set of dose values (a numeric vector) used repeatedly for creating all curves or be a distribution specified as a character string resulting in varying dose values from curve to curve.
The error distribution for the response values can be any continuous distribution
like rnorm
or rgamma
. Alternatively it can be the binomial distribution.
A list with up to 3 components (depending on the value of the onlyY
argument).
Christian Ritz
~put references to the literature/web site here ~
## Model fit m1 <- drm(rootl~conc, data = ryegrass, fct = LL.4()) ## 10 random dose-response curves based on the model fit sim10 <- rdrm(10, LL.4(), coef(m1), xerror = ryegrass$conc) sim10