[R] predict problem

Alan T. Arnholt arnholt at math.appstate.edu
Mon Apr 9 16:08:37 CEST 2001


You need to change "new <- data.frame(x=80)" to read
"new <- data.frame(lot=80)"

Alan



On Mon, 09 Apr 2001 15:43:13 +0200 Troels Ring 
<tring at gvdnet.dk> 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)
> best wishes
> Troels 
> 
> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> 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
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

----------------------
Alan T. Arnholt
arnholt at math.appstate.edu

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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