Hi Henrique, Thank you for your suggestion.
Actually, I have already tried it, but I am confused because the plot I get is not the same as the output of plot(model) or plot.gam(model). The yaxis is different On the other hand, if I build a more complex model, as for example: model<- gam(x ~ s(lat,long) + s(temperature)) I would like to extract the information to build the effects for each explanatory factor (one graph for s(lat,long) and another for s(temperature)), as R does when you use 'plot(model)' and you press return for subsequent pages. My final aim is to plot the influence of s(lat,long) as a contourplot superposed on a geographical map. Maybe there is an easier way to do it... Thank you very much lucía _____________________________________________ Lucia Zarauz AZTI - Tecnalia / Unidad de Investigación Marina Herrera kaia portualdea z/g 20110 Pasaia (Gipuzkoa) Tel: 943 004 800 - Fax: 943 004 801 e-mail: [EMAIL PROTECTED] www.azti.es ; www.tecnalia.info _____________________________________________ **************** LEGE OHARRA **************** AVISO LEGAL **************** DISCLAIMER **************** Mezu hau pertsonala eta isilpekoa da eta baimenik gabeko erabilera debekatua dago legalki. Jasotzailea ez bazara ezabatu mezua, bidali eta kontserbatu gabe. Este mensaje es personal y confidencial y su uso no autorizado está prohibido legalmente. Si usted no es el destinatario, proceda a borrarlo, sin reenviarlo ni conservarlo. This message is personal and confidential, unauthorised use is legally prohibited. If you are not the intended recipient, delete it without resending or backing it. ________________________________________ De: Henrique Dallazuanna [mailto:[EMAIL PROTECTED] Enviado el: martes, 24 de julio de 2007 13:36 Para: Lucia Zarauz CC: r-help@stat.math.ethz.ch Asunto: Re: [R] plotting gam models see ?predict.gam -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O On 24/07/07, Lucia Zarauz <[EMAIL PROTECTED]> wrote: Hi everybody, I am working with gams and I have found some questions when plotting gams models. I am using mgcv, and my model looks something like this: model<- gam(x ~ s(lat,long)) I can plot the output of the model using plot(model) or plot.gam(model) and I get a surface plot. That is ok, but what I want to do now is to extract the data used to perform the surface plot. Like that I would be able to superpose them to a geographical map, and to plot these data using other programs. Thank you very much in advance Lucía _____________________________________________ Lucia Zarauz AZTI - Tecnalia / Unidad de Investigación Marina Herrera kaia portualdea z/g 20110 Pasaia (Gipuzkoa) Tel: 943 004 800 - Fax: 943 004 801 e-mail: [EMAIL PROTECTED] www.azti.es ; www.tecnalia.info _____________________________________________ **************** LEGE OHARRA **************** AVISO LEGAL **************** DISCLAIMER **************** Mezu hau pertsonala eta isilpekoa da eta baimenik gabeko erabilera debekatua dago legalki. Jasotzailea ez bazara ezabatu mezua, bidali eta kontserbatu gabe. Este mensaje es personal y confidencial y su uso no autorizado está prohibido legalmente. Si usted no es el destinatario, proceda a borrarlo, sin reenviarlo ni conservarlo. This message is personal and confidential, unauthorised use is legally prohibited. If you are not the intended recipient, delete it without resending or backing it. ______________________________________________ R-help@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 and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@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 and provide commented, minimal, self-contained, reproducible code.