GAM et interaction

[Jérémie Habasque]

Bonjour,

Je teste actuellement un GAM avec la formule suivante :

Code : Tout sélectionner

gam(W~s(X,by=log(Y)),data=dataset,family="binomial")

Est-ce que le fait d'introduire le paramètre "by" correspond d'une manière à une interaction X*log(Y) ?
Il y a peu de documentation sur ce point...

Merci d'avance,

[Jérémie Habasque]

Finalement j'ai opté pour la formule suivante :

Code : Tout sélectionner

gam(W ~ te(X, log(Y)),data=dataset,family="binomial")

En revanche, l'interprétation du graphique résultat du modèle est assez complexe...

http://www.casimages.com/img.php?i=14011406325417122511899120.png

Si quelqu'un a l'habitude de ce genre de GAM, je suis preneur !

Merci d'avance,

[Jérémie Habasque]

Bonjour,

Je relance le topic...

Merci,

[Jérémie Habasque]

Bonjour,

Petit update de ce post.
La solution trouvée pour l'interprétation du modèle a été de passer par une visualisation en 2D par l'utilisation de l'option 'plot.type="contour"' de la fonction vis.gam du package mgcv :

Code : Tout sélectionner

vis.gam(model,view=c("CHLA","SST"),xlim=c(0,15),ylim=c(10,30),xlab="CHL-a (mg.m-3)",ylab="SST (°C)",
        plot.type="contour",type="response",color="terrain",n.grid=200)

Maintenant je souhaiterais spécifier les couleurs utilisées avec la palette suivante :

Code : Tout sélectionner

# colors palette definition
jet.colors <-colorRampPalette(c("#00007F", "blue", "#007FFF", "cyan",
                                "#7FFF7F", "yellow", "#FF7F00", "red", "#7F0000"))
vis.gam(model,view=c("CHLA","SST"),xlim=c(0,15),ylim=c(10,30),xlab="CHL-a (mg.m-3)",ylab="SST (°C)",
        plot.type="contour",type="response",color=jet.colors(10),n.grid=200)

Mais le message d'erreur suivant apparaît :

Error in vis.gam(model, view = c("CHLA", "SST"), xlim = c(0, 15), ylim = c(10, :
color scheme not recognised

Si quelqu'un a une solution, je suis preneur !
Merci d'avance,