Discussion :

[missipsa01]

Bonjour à toutes et à tous,
Je vous remercie de porter attention à ma présente problématique concernant l'utilisation de glmmADMB avec un modèle de Hurdle, car je n'ai jusqu'à présent pas encore trouvé de solution, d'autant plus que je suis nouveau dans R.
J'essaie d'analyser des données comportementales contenant un nombre d'essais dits ab. (design expérimental à mesures répétées. Effets fixe: Groupe et Délais; effet aléatoire l'id des sujets) (Données en pièce jointe).
Or il se trouve que différents sujets ont présenté plusieurs valeurs 0. Ce qui nécessite l'utilisation d'un Zero inflated ou d'un Hurdle model avec la fonction glmmADMB().
J'ai essayé de le faire suivant les indications de l'aide accompagnant le package de cette fonction, même si celle-ci me parait incomplète et non claire, mais je reçois des messages d'erreurs, en plus du fait que je ne sais que faire par la suite avec les deux modèles.
Voici les lignes de codes en question :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
datab$ID<-factor(datab$ID)
mod1<-glmmadmb(cbind(AB,Total-AB) ~ DELAY*GROUP+(1|ID), data = subset(datab, AB > 0), family = "truncnbinom1")
datab$nz <- as.numeric(datab$AB > 0)
mod2<-glmmadmb(nz ~ DELAY*GROUP+(1|ID), data = datab, family = "binomial")

Le message d'erreur :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
[1] AVIS: Warning in eval(expr, envir, enclos) :
sd.est not defined for this family
[2] ERREUR:
The function maximizer failed (couldn't find parameter file) Troubleshooting steps include (1) run with
'save.dir' set and inspect output files; (2) change run parameters

Qu'en pensez-vous ?
J'espère de tout cœur trouver de l'aide dans ce forum.
Je vous remercie d'avance.

Missipsa