Bonjour à tous ,
J'ai un petit problème sur une méthode que je ne connais pas mais que je dois utiliser en stage...
C'est la régression logistique a effet aléatoire et je dois la réaliser à l'aide de SAS ...
Je cherche à étudier la relation entre la mortalité lié à une pathologie et des indicateurs de prise en charge de cette pathologie.
J'ai une base de données avec une colonne " établissement " et une colonne " patient ". Chaque patient appartient donc à un établissement. La base de données se compose aussi du sexe et de l'âge de chaque patient et enfin si celui-ci est décédé ou non.
Enfin, on trouve pour chaque établissement la valeur de l'indicateur.
Au niveau méthodologique, on a pensé à une régression logistique à effet aléatoire car les données sont dites " groupées " et il peut y avoir un effet " établissement ". La variable à expliquer serait le décés (oui ou non) et les variables explicatives seraient les indicateurs (même valeur pour chaque patient d'un même établissement) mais aussi le sexe et l'âge (propre à chaque patient).
Après quelques recherches, il me semble que la procédure NLMIXED sous SAS est adéquate ?
Cependant je ne comprends pas trop sa syntaxe et les exemples que je trouve sont surtout liés à des données répétées...
Voici ce que j'ai commencer à faire :
Proc nlmixed data = avc.individus ;
parms b0=0.1 b1=0.1 b2=0.1 b3=0.1 b4=0.1 b5=0.1 b6=0.1 b7=0.1 sc0=0.1 sce=0.1;
eta = b0 + b1*I1 + b2*I2 + b3*I3 + b4*I4 + b5*I5 + b6*I6 + b7*age +g0 ;
expeta = exp(eta);
p = expeta / (1+expeta);
model Deces ~ binary(p);
random g0 ~ normal (0,sc0) subject = id_etab;
run;
/* I1...I6 représentent mes indicateurs */
Je ne sais pas trop à quoi correspondent les paramètres en haut...
et SAS me met une erreur : "Quadrature accuracy of 0.0001000 could not be achieved in 31 points..."
Un grand merci à ceux qui pourront m'aider !
Polo
Accueil