1 Le Système SAS 17:46 Wednesday, March 22, 2006 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.16 secondes temps processeur 0.43 secondes 1 /* TimeStamp (dos) : 22 Mars 06 12:03 ; z:\reglogi.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 "Z:\rlf.dar" firstobs=2 ; 13 input ID $ TAILLE GROUPE $ ; 14 NOTE: L'infile "Z:\rlf.dar" est : Nom de fichier=Z:\rlf.dar, RECFM=V,LRECL=256 NOTE: 30 enregistrements lus dans infile "Z:\rlf.dar". La longueur min. de l'enregistrement était 16. La longueur max. de l'enregistrement était 16. 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.02 secondes 15 data INCONNUES ; 16 infile "z:\pginc.dar" firstobs=2 ; 17 input ID $ TAILLE ; 18 19 run ; NOTE: L'infile "z:\pginc.dar" est : Nom de fichier=z:\pginc.dar, RECFM=V,LRECL=256 2 Le Système SAS 17:46 Wednesday, March 22, 2006 NOTE: 3 enregistrements lus dans infile "z:\pginc.dar". La longueur min. de l'enregistrement était 9. La longueur max. de l'enregistrement était 9. 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 20 21 /* régression logistique et scoring dans la foulée */ 22 23 proc logistic data=CONNUES ; 24 model GROUPE = TAILLE ; 25 output out=SORTIES1 predprobs=(i); 26 score data=INCONNUES out=SORTIES2; 27 28 run; NOTE: PROC LOGISTIC is modeling the probability that GROUPE='0'. One way to change this to model the probability that GROUPE='1' is to specify the response variable option EVENT='1'. AVERTISSEMENT: There is a complete separation of data points. The maximum likelihood estimate does not exist. AVERTISSEMENT: The LOGISTIC procedure continues in spite of the above warning. Results shown are based on the last maximum likelihood iteration. Validity of the model fit is questionable. 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.06 secondes temps processeur 0.07 secondes 29 30 /* affichage des résultats */ 31 32 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.00 secondes 33 proc print data=SORTIES2 ; 34 run ; 3 Le Système SAS 17:46 Wednesday, March 22, 2006 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 35 36 /* tracé non graphique des groupes */ 37 38 data GROUPEC ; 39 set INCONNUES ; 40 GROUPE = "?"; 41 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 42 data ALL ; 43 set CONNUES GROUPEC ; 44 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 45 proc sort data=ALL out=ALLTRIE ; by TAILLE; 46 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.00 secondes 47 proc plot data=ALLTRIE ; 48 plot TAILLE*ID=GROUPE ; 49 run ; 4 Le Système SAS 17:46 Wednesday, March 22, 2006 50 51 /* comparaison de logistic, probit et genmod */ 52 53 TITLE COMPARAISON logistic, probit et genmod ; 54 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 55 proc logistic data=CONNUES ; 56 model GROUPE = TAILLE ; 57 NOTE: PROC LOGISTIC is modeling the probability that GROUPE='0'. One way to change this to model the probability that GROUPE='1' is to specify the response variable option EVENT='1'. AVERTISSEMENT: There is a complete separation of data points. The maximum likelihood estimate does not exist. AVERTISSEMENT: The LOGISTIC procedure continues in spite of the above warning. Results shown are based on the last maximum likelihood iteration. Validity of the model fit is questionable. 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.01 secondes temps processeur 0.02 secondes 58 proc probit data=CONNUES ; 59 class GROUPE ; 60 model GROUPE = TAILLE / d = logistic ; 61 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.02 secondes temps processeur 0.03 secondes 62 proc genmod data=CONNUES ; 63 class GROUPE ; 64 model GROUPE = TAILLE / dist = bin 65 link = logit 66 lrci 5 Le Système SAS 17:46 Wednesday, March 22, 2006 67 ; 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. AVERTISSEMENT: Convergence not attained for at least one side of profile likelihood confidence interval for Prm1. Number of iterations = 50. AVERTISSEMENT: Convergence not attained for at least one side of profile likelihood confidence interval for Prm2. Number of iterations = 50. 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.05 secondes temps processeur 0.04 secondes 68 proc catmod data=CONNUES ; 69 direct TAILLE ; 70 response logits; 71 model GROUPE = TAILLE ; 72 73 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.02 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 15.73 secondes temps processeur 0.70 secondes