Interprétations GLM binomial et poisson
Bonjour,
j'ai utilisé R pour réaliser deux GLM. L'un "binomial" et l'autre "poisson". Seulement, si j'ai compris à quoi servent ces modèles je n'ai jamais utilisé ces fonctions en cours (pas de cours dessus du tout d'ailleurs). Du coup je suis incapable de comprendre les sorties.
Le but de l'expérience est d'étudier la migration d'insectes que l'on marque, lâche à un point précis et recapture (ou pas) grâce à des pièges disposés selon un schéma bien défini. Mon expé est la reconduite d'un protocole déjà mené l'an dernier. J'ai donc eu peu de choix quant aux tests à mener...
Le premier modèle a pour but de répondre à la question : le nombre d'individus recapturés (Nb.piege) est-il fonction de la distance au point de lâcher et de la direction de migration ?
Code:
1
2
3
4
5
6
7
| Call: glm(formula = Nb.piege ~ (Distance..m. + Direction), family = poisson(link = log))
Coefficients:
(Intercept) Distance..m. Direction
0.988594 -0.002222 -0.184426
Degrees of Freedom: 10 Total (i.e. Null); 8 Residual
Null Deviance: 4.558
Residual Deviance: 3.616 AIC: 33.5 |
Si je me trompe pas les coefficients correspondent à ça :
Nb.piege = 0.988594 - 0.002222 (Distance..m.) - 0.184426 (Direction) + epsilon
Comme je ne comprends pas le reste ça ne m'avance pas des masses... Les directions ont été codées avec des chiffres comme l'an dernier. L'expérimentateur avait réussi a déduire de ses résultats une direction préférentielle, perso ça ne me saute pas aux yeux. Quelqu'un a-t-il la moindre piste d'interprétation ?
Le second modèle devait tester l’influence de la couleur de marquage (c'est-à-dire de la date de lâcher) et du sexe sur la proportion d’individus recapturés :
Code:
1
2
3
4
5
6
7
8
| Call: glm(formula = cbind(Nb.mark.capt, Nb.mark - Nb.mark.capt) ~ (Date + Sx), family =
binomial(link = logit))
Coefficients:
(Intercept) // Date16/07/2009 // Date17/07/2009 // Date22/07/2009 // Date30/07/2009 // SxM
-4.3617 // -0.5541 // -0.1006 // 0.2508 // -0.4175 // 0.2128
Degrees of Freedom: 7 Total (i.e. Null); 2 Residual
Null Deviance: 2.395
Residual Deviance: 1.33 AIC: 32.32 |
Avec :
Nb.mark.capt = nombre d'insectes marqués recapturés, pour chaque couleur de marquage et chaque sexe
Nb.mark = nombre d'insectes marqués pour chaque couleur et chaque sexe
De même :
Nb.mark - Nb.mark.capt = - 4.3617 - 0.5541 (Date16/07/2009) - 0.1006 (Date17/07/2009) + 0.2508 (Date22/07/2009) -0.4175 (Date 30/07/2009) + 0.2128 (SxM)
Sachant qu'il me "manque" une date de recapture et que le SxF n'est pas représenté j'ai supposé qu'ils servaient de référence. Me trompe-je ?
Comme précédemment je ne comprend pas le reste des sorties. Toute explication pas trop théorique sera donc la bienvenue...
Je sais pas si mes explications sur mon but son très claires, mais j'espère sinon demandez j'essaierai d'éclaircir.
En tous cas, merci d'avance pour vos éventuelles réponses.
