[R] predict problem

Uwe Ligges ligges at statistik.uni-dortmund.de
Mon Apr 9 16:05:21 CEST 2001



Troels Ring wrote:
> 
> Windows 98
> R : Copyright 2001, The R Development Core Team
> Version 1.2.1 (2001-01-15)
> Dear friends.
> How comes this works and produce a single prediction:
> x <- rnorm(15)
> y <- x + rnorm(15)
> predict(lm(y ~ x))
> new <- data.frame(x = seq(-3, 3, 0.5))
> predict(lm(y ~ x), new, se.fit = TRUE)
> pred.w.plim <- predict(lm(y ~ x), new, interval="confidence")
> new1 <- data.frame(x=3)
> predict(lm(y ~ x), new1, interval="confidence")
> 
> while this refuses to take the "new" and predict ?
> lot <- c(30,20,60,80,40,50,60,30,70,60)
> hours <- c(73,50,128,170,87,108,135,69,148,132)
> z1 <- lm(hours~lot)
> new <- data.frame(x=80)
> predict(z1,new,interval="confidence",level=90)


What you want is level=0.9 ???

Uwe Ligges
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list