[R] how to visualise the smoothers of a model averaged GAM ?

Youssef youssefoulhote at gmail.com
Wed Sep 5 18:57:40 CEST 2012


Hi,
If i understand what you mean, you can do 
Plot(M.gam.auto)
You Will get the centered predictions for smooths, you can also extract predictions for each term using (predict(M.gam.auto, type="terms")).

The results you got are thé coefficients for each of you smooths pieces (as k=4, you get 3), but they re not interpretable as they are.

Hope this could help

Youssef

Le 2012-09-05 à 11:50, fmoreira <fmoreira at isa.utl.pt> a écrit :

> Hi,
> 
> It is not clear to me what a model-averaged (sensu MuMIn package) GAM is...
> How is the smoother averaged and how can I see it ?
> 
> My model is:
> 
> M.gam.auto <- gam (Leksite ~ s(Altmean,k=4) + s(Class_unfav_prop, k=4) + 
>               + s(Visib2000, k=4) + s(Ddisturbance, k=4) + s(Steppe_prop,
> k=4)
>                  + s(autocov2, k=4), family=binomial, data= leks)
> 
> When i extract the confidence intervals for the average models (all models
> with delta AIC <2), I get:
> 
>> confint (avgmod)
>                             2.5 %       97.5 %
> (Intercept)           -0.705495919  0.090974395
> s(Altmean).1          -1.032000478  0.846255206
> s(Altmean).2          -0.131526964  3.181165030
> s(Altmean).3          -0.479038766  1.511470593
> s(autocov2).1         -0.594732400  5.030680753
> s(autocov2).2         -3.932830030 18.377288299
> s(autocov2).3         -5.859925892  3.269061133
> s(Ddisturbance).1     -0.324700903  0.609412974
> s(Ddisturbance).2     -0.016172651  5.880936668
> s(Ddisturbance).3     -0.881357264  0.862160976
> s(Steppe_prop).1      -0.002864215  0.002861629
> s(Steppe_prop).2      -0.022374707  0.022389396
> s(Steppe_prop).3       0.191842041  1.242861265
> s(Class_unfav_prop).1 -0.119628517  0.109562147
> s(Class_unfav_prop).2 -0.246237564  0.274136388
> s(Class_unfav_prop).3 -0.588583036  0.189353845
> s(Visib2000).1        -0.003592773  0.003593891
> s(Visib2000).2        -0.012680896  0.012686619
> s(Visib2000).3        -0.489773942  0.206553913
> 
> Not sure of this means nor how to visualize the smoothers...
> Help would be appreciated.
> 
> 
> Francisco Moreira
> Institute of Agronomy
> Technical University of Lisbon
> 
> 
> 
> --
> View this message in context: http://r.789695.n4.nabble.com/how-to-visualise-the-smoothers-of-a-model-averaged-GAM-tp4642302.html
> Sent from the R help mailing list archive at Nabble.com.
> 
> ______________________________________________
> R-help at r-project.org 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.




More information about the R-help mailing list