Statistiques

tout sur les statistiques
 
AccueilAccueil  FAQFAQ  RechercherRechercher  S'enregistrerS'enregistrer  MembresMembres  GroupesGroupes  Connexion  

Partagez | 
 

 Régression logistique a effet aléatoire (proc nlmixed...)

Voir le sujet précédent Voir le sujet suivant Aller en bas 
AuteurMessage
Polo



Nombre de messages : 2
Date d'inscription : 03/06/2009

MessageSujet: Régression logistique a effet aléatoire (proc nlmixed...)   Mer 3 Juin - 16:56

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
Revenir en haut Aller en bas
Voir le profil de l'utilisateur
 
Régression logistique a effet aléatoire (proc nlmixed...)
Voir le sujet précédent Voir le sujet suivant Revenir en haut 
Page 1 sur 1
 Sujets similaires
-
» Régression logistique à effet aléatoire
» intéraction lors d'une régression logistique [RESOLU]
» Régression logistique
» Significativité coefficients en régression logistique
» régression logistique multiple

Permission de ce forum:Vous ne pouvez pas répondre aux sujets dans ce forum
Statistiques :: Statistiques-
Sauter vers: