Prediction
Predicted values for models of class 'drc' or class 'mrdrc'.
## S3 method for class 'drc'
predict(object, newdata, se.fit = FALSE,
interval = c("none", "confidence", "prediction"),
level = 0.95, na.action = na.pass, od = FALSE, vcov. = vcov, ...)
## S3 method for class 'mrdrc'
predict(object, newdata, se.fit = FALSE,
interval = c("none", "confidence", "prediction"),
level = 0.95, pava = FALSE, ...)object |
an object of class 'drc'. |
newdata |
An optional data frame in which to look for variables with which to predict. If omitted, the fitted values are used. |
se.fit |
logical. If TRUE standard errors are required. |
interval |
character string. Type of interval calculation: "none", "confidence" or "prediction". |
level |
Tolerance/confidence level. |
na.action |
function determining what should be done with missing values in 'newdata'. The default is to predict 'NA'. |
od |
logical. If TRUE adjustment for over-dispersion is used. |
vcov. |
function providing the variance-covariance matrix. |
pava |
logical. If TRUE the fit is monotoniosed using pool adjacent violators algorithm. |
... |
further arguments passed to or from other methods. |
For the built-in log-logistics and Weibull-type models standard errors and confidence/prediction intervals can be calculated. At the moment it only works for the situations where all observations are assumed to have a common variance.
A matrix with as many rows as there are dose values provided in 'newdata' or in the original dataset (in case 'newdata' is not specified) and columns with fitted, standard errors, lower and upper limits of confidence intervals.
Christian Ritz
For details are found in the help page for predict.lm.
## Fitting a model
spinach.model1 <- drm(SLOPE~DOSE, CURVE, data = spinach, fct = LL.4())
## Predicting values a dose=2 (with standard errors)
predict(spinach.model1, data.frame(dose=2, CURVE=c("1", "2", "3")), se.fit = TRUE)
## Getting confidence intervals
predict(spinach.model1, data.frame(dose=2, CURVE=c("1", "2", "3")),
interval = "confidence")
## Getting prediction intervals
predict(spinach.model1, data.frame(dose=2, CURVE=c("1", "2", "3")),
interval = "prediction")Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.