Java Mailing List Archive

http://www.r-help.com/

Home » R Help for Statistical Computing »

[R] Prediction and confidence intervals from predict.drc

Brant Inman

2010-09-06

Replies: Find Java Web Hosting

Author LoginPost Reply
R-helpers,

I am using the package "drc" to fit a 4 parameter logistic model. When I
use the predict function to get prediction on a new dataset, I am not
getting the requested confidence or prediction intervals. Any idea what
is going on? Here is code to reproduce the problem:

---

library(drc)

# Fit model to existing dataset in package
spinach.model <- drm(SLOPE~DOSE, data = spinach, fct = LL.4())

#Generate new fake dataset
newdt <- data.frame(matrix(c(seq(0, 150, 0.1), rep(NA, 1501)), ncol=2,
byrow=F))
colnames(newdt) <- c('DOSE', 'SLOPE')

#Use predict function to get prediction and confidence intervals
pred <- predict(spinach.model, interval='prediction', newdata=newdt)
conf <- predict(spinach.model, interval='confidence', newdata=newdt)
head(pred); head(conf)

---

Examining the output shows the point estimates but not the intervals.
Would like the intervals.


I am using the most recent versions of R and drc on Windows XP.

Thanks,

Brant
 [[alternative HTML version deleted]]

______________________________________________
R-help@(protected)
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
©2008 r-help.com - Jax Systems, LLC, U.S.A.