Bonjour à tous,
Mise en contexte : Je m'intéresse au nombre d'oeufs par pontes (variable explicative) en fonction du volume de repas sanguin (Hématine) chez les moustiques selon qu'ils soient infectés par un parasite ou traité par un antibiotique.
Je dispose donc d'un jeu de données contenant plusieurs variables => Hématine, Jour de ponte (Jponte), Traitement, Infection,Nombre d'oeufs par ponte. Je dispose également d'une variable aléatoire = Effet cage
Comme ma variable explicative est une variable de comptage, j'ai pensé à effectuer un glmer suivant une transformation logarithmique de poisson. Cependant mon modèle ne fonctionne pas et malgré toutes les informations que j'ai pu trouvées sur les forums, ça ne résout pas mon problème.
Voici mon code :
oeufmod0<-glmer(log(Ponte)~Hématine+Infe ction+Traitement+Jponte+Hémati ne*Infection+Hématine*Traiteme nt+Traitement*Infection+Jponte *Traitement+Jponte*Infection+H ématine*Jponte+Jponte*Infectio n*Traitement*Hématine+(1|Cage) ,data=Hématine2,family=poisson ) => Il s'agit d'un modèle maximal avec tous les effets d'interaction entre les variables.
Voici le summary de mon modèle :
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
Family: poisson ( log )
Formula: log(Ponte) ~ Hématine + Infection + Traitement + Jponte + Hématine *
Infection + Hématine * Traitement + Traitement * Infection +
Jponte * Traitement + Jponte * Infection + Jponte * Infection *
Traitement * Hématine + (1 | Cage)
Data: Hématine2
AIC BIC logLik deviance df.resid
Inf Inf -Inf Inf 337
Random effects:
Groups Name Std.Dev.
Cage (Intercept) 1
Number of obs: 354, groups: Cage, 6
Fixed Effects:
(Intercept) Hématine
1.5599295 0.0013092
InfectionNI TraitementContrôle
-0.0189709 0.0677743
Jponte Hématine:InfectionNI
0.0263973 0.0010510
Hématine:TraitementContrôle InfectionNI:TraitementContrôle
-0.0020560 -0.0832584
TraitementContrôle:Jponte InfectionNI:Jponte
-0.0385340 -0.0197733
Hématine:Jponte InfectionNI:TraitementContrôle :Jponte
-0.0005961 0.0681721
Hématine:InfectionNI:Jponte Hématine:TraitementContrôle:Jp onte
0.0004158 0.0011747
Hématine:InfectionNI:Traitemen tContrôle Hématine:InfectionNI:Traitemen tContrôle:Jponte
0.0026395 -0.0021242
convergence code 0; 1 optimizer warnings; 84252 lme4 warnings
>
J'obtiens de nombreux warnings et je ne sais pas du tout comment les interpréter. Merci d'avance pour vos réponses, cela me sera d'une bien grande aide
-----