[R] how to visualise the smoothers of a model averaged GAM ?
fmoreira
fmoreira at isa.utl.pt
Wed Sep 5 17:50:50 CEST 2012
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.
More information about the R-help
mailing list