24.04.2014 Views

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

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

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

data b.logisticmref;<br />

set b.logisticmodref; length variable $ 11; retain variable ;<br />

ordre=_n_ ; if class ne "" then variable=class;<br />

drop control_var class;<br />

run;<br />

proc sort data=b.estimPR out=b.estimPRTRI; by variable classval0;run;<br />

proc sort data=b.logisticmref out=b.logistmreftri; by variable value;run;<br />

data b.PEassemble;<br />

merge b.logistmreftri(rename=(value=classval0)) b.estimPRTRI; by variable classval0;<br />

if estimate=. then do; estimate=0.00; vp="REF";<br />

end;<br />

run;<br />

Il nous reste une <strong>de</strong>rnière chose avant <strong>de</strong> construire notre table présentant le résultat classiquement présenté par<br />

les chercheurs dans leurs publications : habiller proprement les noms <strong>de</strong> variable, hérités du consortium Insee &<br />

<strong>Ined</strong> et trop abrégés pour être publiés :<br />

proc format;<br />

value $intitule<br />

"SexeE" = "Sexe"<br />

"sitfam" = "Situation familiale actuelle"<br />

"MPROLOC" = "Logement"<br />

"lienmig" = "Lien à la migration"<br />

"LT1AMI1" = "Relation avec le voisinage"<br />

"GRFAM" = "Présence <strong>de</strong> la famille dans la région"<br />

"echstab_tot" = "Trajectoire"<br />

"trajloc" = "Localisation <strong>de</strong> la trajectoire"<br />

"echdouvC" = "Echelle d'appartenance"<br />

"agepop" = "Age"<br />

"_Intercept" = " Constante"<br />

;<br />

run;<br />

Nous pouvons donc maintenant sortir notre résultat synthétique et l’habiller grâce aux outils habituels déjà vu en<br />

V1. Les 2 résultats (rtf et pdf) qui seront obtenus le sont sans fermer la <strong>de</strong>stination html <strong>de</strong> la work en 9.3 qui<br />

elle sera nettoyée à la fermeture <strong>de</strong> <strong>Sas</strong>.<br />

Un simple proc report fera l’affaire.<br />

Ayant à présenter le facteur Beta <strong>de</strong> l’équation logistique, j’utilise le supplément/correctif mentionné en V2p7 .<br />

Le startpage=no a été utilisé pour forcer le texte d’ods text= à rester sur la même page.<br />

Ces ods text= ont été enrichis <strong>de</strong> saut <strong>de</strong> ligne et ajustés grâce au choix d’une police et <strong>de</strong> sa taille combinés à<br />

un leftmargin <strong>de</strong> 2 cm. Tout cela est conforme à V1_p134.<br />

ods escapechar="^";<br />

options orientation=portrait nodate nonumber; options leftmargin=2cm;<br />

ods rtf file="d:\notes\logistic\resultat.rtf" style=journal startpage=no;<br />

ods pdf file="d:\notes\logistic\resultat.pdf" style=journal startpage=no;<br />

title1 "L'attachement au lieu <strong>de</strong> rési<strong>de</strong>nce <strong>de</strong> l’ensemble <strong>de</strong>s enquêtés" ;<br />

title2 "Paramètres estimés (*ESC*){unico<strong>de</strong> beta} pour l'attachement à son lieu <strong>de</strong> rési<strong>de</strong>nce vs un nonattachement<br />

^{newline} (modèle Logit)";<br />

proc report data=b.PEassemble nowd spanrows;<br />

column variable classval0 estimate vp ;<br />

<strong>de</strong>fine variable /or<strong>de</strong>r format= $intitule. "Variables et" ;<br />

<strong>de</strong>fine classval0 /display "modalités actives" ;<br />

<strong>de</strong>fine estimate /display format=12.2 "Val estimée" ;<br />

<strong>de</strong>fine vp /display "Significatif" format=$12.;<br />

run;<br />

ODS text= " ^{style[font_face='Times New Roman' font_size=1]Lecture : Les coefficients sont très<br />

significatifs *** à moins <strong>de</strong> 1 %, significatif ** entre 1% et 5% et peu significatif * entre 5% et<br />

10%.}<br />

^{newline} ^{style[font_face='Times New Roman' font_size=1]Au-<strong>de</strong>là <strong>de</strong> ce seuil, la valeur est non<br />

significative(-), selon le tableau <strong>de</strong> sortie Estimations par l'analyse du maximum <strong>de</strong> vraisemblance} ";<br />

ODS text= "^{style[font_face='Times New Roman' font_size=1]Champ : L'ensemble <strong>de</strong>s personnes qui ont<br />

déclaré être attaché au lieu <strong>de</strong> rési<strong>de</strong>nce soit 2890 individus. }";<br />

ODS text= "^{style[font_face='Times New Roman' font_size=1]Source : Enquête HdV, Insee-<strong>Ined</strong>, 2003.}";<br />

run;<br />

ods rtf close;<br />

ods pdf close;<br />

160 <strong>Pratique</strong> <strong>de</strong> <strong>Sas</strong> windows ….93 <strong>Volume</strong> 2 et Coll.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!