Comment imprimer the the Somers'D dans un fichier SAS?

Question

Comme le suggère le titre, je me demande s'il existe un moyen d'imprimer les statistiques de Somers'D et la valeur p du prédicteur x dans un ensemble de données.

Vous pouvez obtenir ces statistiques simplement en cours d'exécution:

ODS TRACE ON; 
PROC LOGISTIC DATA = BETTING.TRAINING_DUMMIES NOPRINT; 
    MODEL Z1 (EVENT = '1') = D_INT_LNGAP_1; 
     OPTIONS; 
RUN; 

ODS TRACE OFF; 

ODS OUTPUT FITSTATISTICS=FITDS; 
PROC LOGISTIC DATA = BETTING.TRAINING_DUMMIES NOPRINT; 
    MODEL Z1 (EVENT = '1') = D_INT_LNGAP_1; 
     OPTIONS; 
RUN; 

Si je lance un code similaire à celui proposé here, je reçois que l'AIC, le SIC et enfin la stat LR et dans le SAS log Je trouve:

10 ODS TRACE ON; 
11 
12 PROC LOGISTIC DATA = BETTING.TRAINING_DUMMIES NOPRINT; 
13  MODEL Z1 (EVENT = '1') = D_INT_LNGAP_1; 
14   OPTIONS; 
15 RUN; 

NOTE: PROC LOGISTIC is modeling the probability that z1=1. 
NOTE: Convergence criterion (GCONV=1E-8) satisfied. 
NOTE: There were 3968 observations read from the data set BETTING.TRAINING_DUMMIES. 
NOTE: PROCEDURE LOGISTIC used (Total process time): 
     real time   0.07 seconds 
     cpu time   0.04 seconds 


16 
17 ODS TRACE OFF; 

dans le premier morceau de code, alors que dans le second, je trouve ce qui suit:

18 ODS OUTPUT FITSTATISTICS=FITDS; 
NOTE: Writing HTML Body file: sashtml.htm 
19 PROC LOGISTIC DATA = BETTING.TRAINING_DUMMIES NOPRINT; 
20  MODEL Z1 (EVENT = '1') = D_INT_LNGAP_1; 
21   OPTIONS; 
22 RUN; 

NOTE: PROC LOGISTIC is modeling the probability that z1=1. 
NOTE: Convergence criterion (GCONV=1E-8) satisfied. 
NOTE: There were 3968 observations read from the data set BETTING.TRAINING_DUMMIES. 
NOTE: PROCEDURE LOGISTIC used (Total process time): 
     real time   0.04 seconds 
     cpu time   0.04 seconds 

WARNING: Output 'FITSTATISTICS' was not created. Make sure that the output object name, label, 
     or path is spelled correctly. Also, verify that the appropriate procedure options are 
     used to produce the requested output object. For example, verify that the NOPRINT 
     option is not used. 

Certains d'entre vous peuvent suggérer un moyen d'imprimer de telles statistiques dans un nouvel ensemble de données?

Toute aide sera appréciée.

Merci!

Answer 1

Je ne sais pas pourquoi vous n'obtenez pas de sortie ODS TRACE. Je voudrais redémarrer votre version SAS ou le signaler à SAS.

Les tables que vous souhaitez s'appellent Association et ParameterEstimates. Somer's D nécessite que l'instruction Odds Ratio soit créée.

ods trace on; 
ods output association=somers parameterestimates=pe; 
proc logistic data=sashelp.heart; 
model status=ageatstart; 
oddsratio ageatstart; 
run; 
ods trace off;