[R] need 95% confidence interval bands on cubic extrapolation

[Liaw, Andy]

Look at the "interval" option in ?predict.lm.

Andy

From: James Salsman
> 
> Dear R experts:
> 
> I need to get this plot, but also with 95% confidence interval bands:
> 
>    hour <- c(1, 2, 3, 4, 5, 6)
>    millivolts <- c(3.5, 5, 7.5, 13, 40, 58)
> 
>    plot(hour, millivolts, xlim=c(1,10), ylim=c(0,1000))
> 
>    pm <- lm(millivolts ~ poly(hour, 3))
> 
>    curve(predict(pm, data.frame(hour=x)), add=TRUE)
> 
> How can the 95% confidence interval band curves be plotted too?
> 
> Sincerely,
> James Salsman
> 
> P.S.  I know I should be using data frames instead of parallel lists.
> This is just a simple example.
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! 
> http://www.R-project.org/posting-guide.html
> 
>