[R] predict.glm - how to?

Laurits Søgaard Nielsen laurits at lauritsnielsen.dk
Thu Mar 2 16:56:10 CET 2006


Hi

I have a little R problem. I have created a GLM model in R and now I want to
predict some values outside the values I have in the model (extrapolate).

I have this code:

fitted.model4 <- glm(Yval ~ time, family=gaussian, data=Fuel)

The question is - How do I predict a value of Yval ie with a value of time =
340 and also get confidence/prediction intervals for Yvar?

I have tried the predict.glm but cannot make it work - any suggestions how
to do this using predict.glm?

****

Yval are 312 fuel prices ranging from 60 to 140 and time is a sequence from
1 to 312.

The GLM summary output is as follows:

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-24.978   -4.033   -0.391    4.747   15.323  

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) 79.698229   0.871830   91.42   <2e-16 ***
seq(1:312)  0.180127   0.004828   37.31   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

(Dispersion parameter for gaussian family taken to be 59.00219)

Null deviance: 100408  on 311  degrees of freedom
Residual deviance:  18291  on 310  degrees of freedom
AIC: 2161.6

Number of Fisher Scoring iterations: 2




More information about the R-help mailing list