[R] question about predict.gam
Simon Wood
s.wood at bath.ac.uk
Fri Oct 5 18:03:39 CEST 2007
On Friday 05 October 2007 16:25, Michela Leone wrote:
> I'm fitting a Poisson gam model, say
>
> model<-gam(a65tm~as.factor(day.week
> )+as.factor(week)+offset(log(pop65))+s(time,k=10,bs="cr",fx=FALSE,by=NA,m=1
>),sp=c( 0.001),data=dati1,family=poisson)
>
>
>
> Currently I've difficulties in obtaining right predictions by using
> gam.predict function with MGCV package in R version 2.2.1 (see below my
> syntax).
--- What mgcv version is it?
>
> previ<-predict.gam(model,dati2,type="response")
--- is this supposed to be `dati2' or `dati1'? It's hard to understand the
following comment if the data frame is different in fitting and prediction...
>
> I expect to get the predicted values in the same scale of the response
> variable but they are far what I expect. (i.e 1.3 when my response variable
> is 40), while when I use predict.gam with gam package I have right
> prediction (in the response scale).
>
> Does anybody know if I missed anything important? Is there anyway to make
> it work correctly or reliably?
>
--- It works on my test examples and no one has reported a problem until now.
A quick fit of your model structure to simulated data also produces nothing
unusual. Could you perhaps try again with the current mgcv version? (Not that
there is anything in the changeLog likely to make a difference.) If there is
still a problem would you be able to send me data and model fitting code
showing the problem, and I'll look into it (data confidential, of course)?
> Moreover, if it's possible, could anybody explain me why using gam function
> with GAM package I obtain different predictions from those obtained using
> gam function with MGCV package?
>
--- They'll always be somewhat different as the representation of the models
is slightly different in gam::gam and mgcv::gam. See ?gam
best,
Simon
--
> Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK
> +44 1225 386603 www.maths.bath.ac.uk/~sw283
More information about the R-help
mailing list