Pratique de Sas Windows ...9.3 Volume 2 - Ined

3.6.6 Rendre un tableau de résultats statistiques parlant visuellement : les rapports de cotes retravaillé
Nous avons marqué plus haut certaines cases dans le tableau des rapports de cotes (objet CLoddsWald) en
grisé, vert, rouge.
Essayons de le faire par programmation (sans macro), car nous avons suivi des règles précises pour déclarer
significatif ou non certaines lignes et opposé ce qui amoindrissait à ce qui augmentait la force de l’attachement.
* EXAMEN POUR METTRE EN COULEUR LE TABLEAU ODDS RATIO;
ods select Stat.Logistic.Cloddswald ;
ods output Stat.Logistic.Cloddswald=b.LogisticOR ;
proc logistic data=b.d simple plots=none;
class trajloc (ref= "Tjrs urbain") echstab_tot(ref="En Frce et étranger")
sexee (ref="Homme") sitfam(ref="Couple avec enfant")
LT1AMI1(ref="pas voisin ami") lienmig(ref="Autre") echdouvc(ref="Commune")
GRFAM(ref="Pas de famille") mproloc(ref="Autres") agepop(ref="18-35 ans")
/param=ref ;
weight poidsfc;
model glbio= sexee sitfam mproloc lienmig agepop echstab_tot trajloc
echdouvc LT1AMI1 grfam
/rsq link=logit alpha=0.05 clodds=wald ;
run;
ods select all;
ods output close;
A l’aide de fonctions caractère classiques, je distingue la variable (variable class) de ses modalités (variable
value) et j’isole la modalité correspondante après la chaîne 'vs' sous le nom de variable versus.
J’applique également les tests de validité pour construire une échelle soit la nouvelle variable vs.
data or;
set b.logisticor; order=_n_;
length class $ 11 value versus $ 20;
class=scan(effect,1,' ');
psmot=find(effect,scan(effect,2,' '));
pvs=find(effect,'vs');
value=substr(effect,psmot,pvs-(1+psmot));
versus=substr(effect,pvs+3);
If lowercl