################################################# # # # exemple de régression logistique # # # # (gH) gilles.hunault@univ-angers.fr # # # ################################################# page 1 - Le Système SAS NOTE: Copyright (c) 2002-2003 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) 9.1 (TS1M3) Licensed to LICENCE GRATUITE DOMICILE ETUDIANTS PROFS SAS MENR, Site 0088785001. NOTE: La session est exécutée sur la plate-forme XP_PRO . NOTE: L'initialisation de SAS a utilisé : temps réel 15.36 secondes temps processeur 0.50 secondes 1 /* TimeStamp (dos) : 22 Mars 06 17:58 ; z:\rl.sas */ 2 3 OPTIONS LINESIZE=132 PAGESIZE=50 NOCENTER FORMDLIM=' ' FORMCHAR='-----------' ; 4 5 TITLE "Etude des données fictives PG (Petits et Grands)" ; 6 TITLE2 "================================================" ; 7 8 9 /* données */ 10 11 data CONNUES ; 12 infile datalines firstobs=3 ; 13 input ID $ TAILLE GROUPE $ ; 14 datalines ; NOTE: La table WORK.CONNUES a 30 observations et 3 variables. NOTE: L'étape DATA a utilisé (Durée totale du processus) : temps réel 0.01 secondes temps processeur 0.01 secondes 47 ; 48 49 data INCONNUES ; 50 infile datalines firstobs=3 ; 51 input ID $ TAILLE ; 52 datalines ; NOTE: La table WORK.INCONNUES a 3 observations et 2 variables. NOTE: L'étape DATA a utilisé (Durée totale du processus) : temps réel 0.00 secondes temps processeur 0.00 secondes 58 ; 59 60 run ; page 2 - Le Système SAS 61 62 /* régression logistique et scoring dans la foulée */ 63 64 proc logistic data=CONNUES descending ; 65 model GROUPE = TAILLE ; 66 output out=SORTIES1 predprobs=(i); 67 score data=INCONNUES out=SORTIES2; 68 69 run; NOTE: PROC LOGISTIC is modeling the probability that GROUPE='1'. NOTE: Convergence criterion (GCONV=1E-8) satisfied. NOTE: 30 observations copiées de la table WORK.CONNUES. NOTE: La table WORK.SORTIES1 a 30 observations et 7 variables. NOTE: La table WORK.SORTIES2 a 3 observations et 5 variables. NOTE: La procédure LOGISTIC a imprimé les pages 1-2. NOTE: La procédure LOGISTIC a utilisé (Durée totale du processus) : temps réel 0.07 secondes temps processeur 0.08 secondes 70 71 /* affichage des résultats */ 72 73 proc print data=SORTIES1 ; NOTE: 30 observations copiées de la table WORK.SORTIES1. NOTE: La procédure PRINT a imprimé la page 3. NOTE: La procédure PRINT a utilisé (Durée totale du processus) : temps réel 0.00 secondes temps processeur 0.01 secondes 74 proc print data=SORTIES2 ; 75 run ; NOTE: 3 observations copiées de la table WORK.SORTIES2. NOTE: La procédure PRINT a imprimé la page 4. NOTE: La procédure PRINT a utilisé (Durée totale du processus) : temps réel 0.00 secondes temps processeur 0.01 secondes 76 77 /* tracé non graphique des groupes */ 78 79 data GROUPEC ; 80 set INCONNUES ; page 3 - Le Système SAS 81 GROUPE = "?"; 82 NOTE: 3 observations copiées de la table WORK.INCONNUES. NOTE: La table WORK.GROUPEC a 3 observations et 3 variables. NOTE: L'étape DATA a utilisé (Durée totale du processus) : temps réel 0.00 secondes temps processeur 0.01 secondes 83 data ALL ; 84 set CONNUES GROUPEC ; 85 NOTE: 30 observations copiées de la table WORK.CONNUES. NOTE: 3 observations copiées de la table WORK.GROUPEC. NOTE: La table WORK.ALL a 33 observations et 3 variables. NOTE: L'étape DATA a utilisé (Durée totale du processus) : temps réel 0.00 secondes temps processeur 0.00 secondes 86 proc sort data=ALL out=ALLTRIE ; by TAILLE; 87 NOTE: 33 observations copiées de la table WORK.ALL. NOTE: La table WORK.ALLTRIE a 33 observations et 3 variables. NOTE: La procédure SORT a utilisé (Durée totale du processus) : temps réel 0.00 secondes temps processeur 0.01 secondes 88 proc plot data=ALLTRIE ; 89 plot TAILLE*ID=GROUPE ; 90 run ; 91 92 /* comparaison de logistic, probit, genmod et catmod */ 93 94 TITLE "COMPARAISON logistic, probit, genmod et catmod" ; 95 TITLE2 "==============================================" ; 96 NOTE: 33 observations copiées de la table WORK.ALLTRIE. NOTE: La procédure PLOT a imprimé la page 5. NOTE: La procédure PLOT a utilisé (Durée totale du processus) : temps réel 0.00 secondes temps processeur 0.00 secondes page 4 - Le Système SAS 97 proc logistic data=CONNUES descending ; 98 model GROUPE = TAILLE ; 99 NOTE: PROC LOGISTIC is modeling the probability that GROUPE='1'. NOTE: Convergence criterion (GCONV=1E-8) satisfied. NOTE: 30 observations copiées de la table WORK.CONNUES. NOTE: La procédure LOGISTIC a imprimé les pages 6-7. NOTE: La procédure LOGISTIC a utilisé (Durée totale du processus) : temps réel 0.02 secondes temps processeur 0.03 secondes 100 proc probit data=CONNUES ; 101 class GROUPE ; 102 model GROUPE = TAILLE / d = logistic ; 103 NOTE: PROC PROBIT is modeling the probabilities of levels of GROUPE having LOWER Ordered Values in the response profile table. NOTE: Algorithm converged. NOTE: La procédure PROBIT a imprimé les pages 8-9. NOTE: La procédure PROBIT a utilisé (Durée totale du processus) : temps réel 0.01 secondes temps processeur 0.02 secondes 104 proc genmod data=CONNUES ; 105 class GROUPE ; 106 model GROUPE = TAILLE / dist = bin 107 link = logit 108 lrci ; 109 NOTE: PROC GENMOD is modeling the probability that GROUPE='0'. One way to change this to model the probability that GROUPE='1' is to specify the DESCENDING option in the PROC statement. NOTE: Algorithm converged. NOTE: The scale parameter was held fixed. NOTE: La procédure GENMOD a imprimé les pages 10-11. NOTE: La procédure GENMOD a utilisé (Durée totale du processus) : temps réel 0.04 secondes temps processeur 0.03 secondes 110 proc catmod data=CONNUES ; 111 direct TAILLE ; 112 response logits ; page 5 - Le Système SAS 113 model GROUPE = TAILLE ; 114 115 run; NOTE: The default estimation method for this model is maximum-likelihood. NOTE: Maximum likelihood computations converged. NOTE: 30 observations copiées de la table WORK.CONNUES. NOTE: La procédure CATMOD a imprimé les pages 12-13. NOTE: La procédure CATMOD a utilisé (Durée totale du processus) : temps réel 0.01 secondes temps processeur 0.02 secondes NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414 NOTE: Le Système SAS a utilisé : temps réel 16.01 secondes temps processeur 0.78 secondes