Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

  1. #1
    Membre à l'essai
    Sensibilité et spécificité avec proc logistic / option ctable
    Bonjour,

    J'aimerais pouvoir sortir les Sensibilités et spécificités d'un test par rapport à une maladie :
    ma maladie est "covid" en 1 (non malade) et 2 (malade)
    mon test est "toux" en 1(test positif) et 0 (test negatif)

    Je fais ce code sous SAS :

    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    proc logistic data  = covid.appel; where CAT_TEST = 1;
    title 'logistique sensibilite';
    class toux (ref = '0');
    model covid (ref = '1') = toux / ctable ;
    run;


    L'option ctable me permet d'obtenir les Se et Sp mais le problème c'est que les effectifs ne sont pas au bon endroit et il me donne 1-Se et 1-Sp en affichant qu'il s'agit de Se et de Sp, est-ce que vous auriez une idée du problème ?

    Merci beaucoup,

  2. #2
    Membre éprouvé
    proc logistic descending
    Bonjour,

    Vous pouvez essayer avec l'option descending :

    http://support.sas.com/resources/pap...2/087-2012.pdf

    Intuitively, we flag our data as zero for not having the event and one for the event of interest, however, PROC LOGISTIC will model on whatever value is the lowest. This means when the above basic syntax is ran, SAS will model the probability a patient does not have the disease. This is not usually what the researcher wishes to answer. We are most interested in the probability we have the disease and how the explanatory variables relate to it. So a simple solution is to either change the way you are flagging the event in the data set (e.g. “1” for disease and “2” for no disease) or include the DESCENDING option.
    Vous pouvez voir aussi cette page et le guide SAS The LOGISTIC procedure.

    Ce n'est pas propre à la proc logistic de SAS. En R par exemple, la fonction twoby2() du package Epi considère par défaut que quand une variable est codée 0/1, 0 signifie malade versus 1 non malade et 0 signifie avoir le facteur de risque versus 1 ne pas avoir le facteur de risque.

    Cordialement,