Bonjour
Je chercher à évaluer l'effet d'une exposition au bruit sur la prise de médicaments. Mes sujets ont été recrutés aléatoirement dans 7 pays différents.
Je dois appliquer un effet aléatoire à la fois sur l'intercept et sur la pente de l'effet du bruit avec le pays comme subject.

Mon code est le suivant :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
proc glimmix data=data   method = quad; 
   class X1 X2 X3  ; 
   model medication = noise X1 X2 X3 / dist=binary link=logit  solution or ;
   random int / subject = country ; 
   random noise/ subject = country  solution cl ;
   covtest 'Need Random Intercept?' zeroG ;
   covtest "Need Random slope?" . 0 0 ;
run;
Est-ce correct selon vous ?
Comment obtenir les différents odds ratio pour l'effet du bruit sur médication dans chaque pays ?
Comment évaluer le fait de conserver ou non l'effet aléatoire sur l'intercept? sur la pente ?

Merci pour votre aide