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

3.6.9 Produire un forest plot avec l’odds ratios par une proc Sgplot
Pour clore ce parcours et retomber sur nos pattes (statistiques), je présente ici une autre ressource Sas :
la FAQ stat dont j’extrais l’article Odds Ratios, qui montre … qu’il y en a partout !
Odds Ratios
Model-based methods: By default, Sas/STAT PROC LOGISTIC produces odds ratio point estimates and confidence intervals
(CLODDS= to select Wald or profile likelihood intervals) for predictors not involved in interactions. Beginning in Sas 9.2, the
ODDSRATIO statement provides estimates and confidence intervals even when predictors are involved in interactions. The
ODDSRATIO option, along with the DIFF= or SLICEDIFF= option, in the LSMEANS statement of PROC GLIMMIX
provides point estimates (add CL option for confidence intervals). Custom odds ratios can be computed in the Sas/STAT
procedures CATMOD, LOGISTIC, GENMOD, GLIMMIX, or SURVEYLOGISTIC by determining the contrast of logistic
model parameters that represents the desired difference in log odds for two groups and then using the CONTRAST or
ESTIMATE statement with the ESTIMATE=EXP or EXP option. Beginning with Sas/STAT 9.22 in Sas 9.2 TS2M3 when
GLM parameterization of CLASS variables is used in the LOGISTIC, GENMOD, GLIMMIX, or SURVEYLOGISTIC
procedures, you can also use the EXP option in the LSMESTIMATE statement, or the DIFF and EXP options in the
LSMEANS and SLICE statements.
Nonmodeling methods: The RELRISK or MEASURES option in Base Sas PROC FREQ for 2x2 tables and the CMH option
for stratified 2x2 tables provide point estimates and confidence intervals. The OR and COMOR options in the EXACT
statement provide exact confidence intervals and tests. Also, OR option in Sas/Genetics PROC CASECONTROL.
http://support.sas.com/kb/30/333.html
Il y a même du flou autour du terme hors Sas http://www.amstat.org/publications/jse/v20n3/fulton.pdf
mentionné par Rob in http://blogs.sas.com/content/Sastraining/2013/02/14/what-are-the-odds-of-winning-the-lottery/
Cependant, comme pour le graphique du lifetest, le citron odds-ratio, n’aura pas été pressé à fond, si nous ne
montront pas comme essayer de les représenter par une proc SG.
L’odds-ratio en representation graphique par une proc SG semble possible, mais en s’inspirant du code lié à la
figure 9.1.8 du «petit livre rouge » page 242.
/*--mise en macro-variables d'info basées sur la table
et utilisées pour les décalages Offset sur l'axe Y--*/
data _null_;
pct=0.70/nobs;
height=14;
dpi=100;
call symputx("pct", pct);
call symputx("pct2", 2*pct);
call symputx("dpi", dpi);
call symputx("height", height);
call symputx("heightin", height || "cm");
call symputx("thickness", floor(height*dpi*pct));
set b.or nobs=nobs;
run;
/*--Enrichissement de la table déjà récupérée des Odds-ratio--*/
data b.oram;
set b.or;
OR='OR'; LCL='LCL'; UCL='UCL'; Sig="Signif";
format oddsratioest lowercl uppercl 6.3;
run;
/*--Création de notre premier template de style graphique modifié--*/
proc template;
define style styles.ForestColor93;
parent = Styles.htmlBlue;
style GraphFonts from GraphFonts /
'GraphDataFont' = (", ",7pt)
'GraphValueFont' = (", ",7pt);
style GraphData1 from GraphData1 /
contrastcolor = GraphColors('gcdata2')
color = GraphColors('gdata2');
style GraphData2 from GraphData2 /
contrastcolor = GraphColors('gcdata1')
color = GraphColors('gdata1');
end;
run;
170 Pratique de Sas windows ….93 Volume 2 et Coll.