negative lsmeans glimmix proc

Je fais un GLMM (loi binomiale) sur SAS, mon expé est la suivante:
j'ai plusieurs groupes d'individus, dans chaque groupe j'ai des individus malades ou sains. A chaque groupe je présente 3 traitement successivement. Chaque expé dure 20 min, et je réalise des 'scans' toutes les min, où je note quels individus sont dans la zone du traitement. J'ai donc pour chaque individu pour chaque test 20données, '1' si l'individu est dans la zone, '0' s'il ne l'est pas.
Ma réponse est le nombre de présence dans la zone, et mes variables explicatives d'intérêt: statut (sain/malade), traitement, statut*traitement, effet aléatoire individu dans son groupe, et effet aléatoire test.

Mon problème est que j'obtiens des estimations de lsmean négative...
Quelqu'un pourrait-il m'expliquer pourquoi...?

Mille fois merci!

Bonjour,

utilise l'option ILINK avec le lsmeans pour transformer les estimations en probabilité.
Dis-moi si c'est ok.

Niaboc

Bonjour Niaboc,
et merci.
Je pourrai faire tourner le modèle lundi seulement, je te dirai ce que ca donne.
Mais en fait j'ai avancé dans ma réflexion... mes probabilités pour chaque individu d'être dans la zone du stimulus sont systématiquement inférieures à 0.5, donc en appliquant une transformations logit, toutes mes estimations sont de fait négatives...
Je ne suis pas bien sûre de comprendre dans quel cas utiliser et à quoi sert l'option ILINK?

L'option Ilink te permet simplement de transformer le logit de tes lsmeans en probabilité (que tu retrouverais donc ici en dessous de 0.5).