Courbe roc sous SAS
vis2
Messages postés
1
Statut
Membre
-
Poupou -
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
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:
- Courbe roc excel
- Liste déroulante excel - Guide
- Courbe excel - Guide
- Word et excel gratuit - Guide
- Si ou excel - Guide
- Déplacer colonne excel - Guide
4 réponses
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.
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.
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é.
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é.
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
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
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
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