Courbe roc sous SAS

vis2 Messages postés 1 Statut Membre -  
 Poupou -
Bonjour,

Je cherche à faire des courbes ROC sous SAS pour déterminer le seuil opitimal de décision.
J'arrive à obtenir la courbe grâce au programme:

proc logistic data=ttt;
model acidose=dfg / outroc=roc1 roceps=0;
output out=outp p=p;
run;
proc print data=roc1;
run;
title "ROC plot for acidose = dfg";
%rocplot(outroc = roc1,
out = outp,
p = p,
id = alb tp totscore,
roffset = 8)

Comment trouver le seuil? Existe il une option?
Si vous connaissez des logiciels qui sont mieux, ca m'interesse aussi.

Merci beaucoup
A voir également:

4 réponses

Poupou
 
Bonjour,

Une solution est de calculer la distance entre chaque point de la courbe de coordonnées (1-spécificité ; sensibilité) au point de coordonnées (0 ; 1). Le seuil sera le point minimisant cette distance.
4
katy
 
Bonjour
Pour faire la courbe de ROC sous SAS, tu peux utliser la fonction gplot.
Avec la proc logistique tu récupères la sensibilité et la spécificité. Du coup tu peux tracer le graphe.

J'espère t'avoir aidé.
2
Bioinfo
 
Bonjour,

déjà je m'excuse car je ne connais pas la réponse à ta question principale, celle sur le seuil.
Je vois à peu près ce qu'est une roc curve, mais je ne suis pas sûr de savoir de quel seuil tu veux parler.

Par contre, pour les autres logiciels que tu pourrais utiliser, il y a R. Tu peux aussi allier R à Splus ou utiliser le langage S.
Mais vu ton niveau en SAS (qui me paraît plutôt bon), je suis sûr que tu connais aussi ces autres languages.
Je trouve que SAS possède plein de fonctionnalités intégrées qui permettent de faire des études statistiques poussées en quelques lignes, par contre il me paraît moins "modulable" et "programmable" que les autres langages que je t'ai cité. Tu trouveras, par exemple, une comparaison entre SAS et le langage S sur un pdf disponible à l'adresse suivante:
http://cran.r-project.org/doc/contrib/Alzola+Harrell-Hmisc-Design-Intro.pdf

J'aimerais te poser une question à mon tour: sais-tu si c'est possible de récupérer la likelihood value d'une proc logistic, ailleur que dans l'output de SAS? En fait, j'aimerais ne sortir que les résultats dont la likelihood value (ou pi-value) est inférieur à 0.05.

Merci d'avance
1
imen
 
bonjour,
actuellemnt, je suis entrain de préparer ma thèse en medécine et je cherche à construire des crourbe ROC, ^car c'est très utilisé dans les articles, le problème c'est que je ne comprends pas parfaitement la notion, et la fonçon d'utilisation de ce test statistique
merci de bien vouloir m'aider
0