/*  TimeStamp  (dos) : 22 Mars 06 12:03  ;  z:\reglogi.sas */

OPTIONS LINESIZE=132 PAGESIZE=50 NOCENTER FORMDLIM=' ' FORMCHAR='-----------' ;

TITLE  "Etude des données fictives PG (Petits et Grands)" ;
TITLE2 "================================================" ;


/* données */

data CONNUES ;
     infile "Z:\rlf.dar" firstobs=2 ;
     input ID $  TAILLE  GROUPE $ ;

data INCONNUES ;
     infile "z:\pginc.dar" firstobs=2 ;
     input ID $  TAILLE  ;

run ;

/* régression logistique et scoring dans la foulée */

proc logistic data=CONNUES ;
  model GROUPE = TAILLE ;
  output out=SORTIES1 predprobs=(i);
  score data=INCONNUES out=SORTIES2;

run;

/* affichage des résultats  */

proc print data=SORTIES1 ;
proc print data=SORTIES2 ;
run ;

/* tracé non graphique des groupes */

data GROUPEC ;
     set INCONNUES ;
     GROUPE = "?";

data ALL ;
     set CONNUES GROUPEC ;

proc sort data=ALL out=ALLTRIE ; by TAILLE;

proc plot data=ALLTRIE ;
     plot TAILLE*ID=GROUPE ;
run ;

/* comparaison de logistic, probit et genmod */

TITLE COMPARAISON logistic, probit et genmod ;

proc logistic data=CONNUES ;
  model GROUPE = TAILLE ;

proc probit data=CONNUES ;
  class GROUPE ;
  model GROUPE = TAILLE / d = logistic ;

proc genmod data=CONNUES ;
  class GROUPE ;
  model GROUPE = TAILLE / dist = bin
                          link = logit
                          lrci
                          ;
proc catmod data=CONNUES ;
  direct TAILLE ;
  response logits;
  model GROUPE = TAILLE ;

run;