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

3.6.7 La probabilité des estimates : comment produire directement sa page de rapport
En reprenant la table des probabilités Parameterestimates que je capture après sélection ainsi que celle
contenant les infos sur les variables et les modalités y compris celles mises en référence ClasslevelInfo
ods select Stat.Logistic.ParameterEstimates
Stat.Logistic.ClasslevelInfo;
ods output Stat.Logistic.ParameterEstimates=b.LogisticEst
Stat.Logistic.ClasslevelInfo=b.LogisticModRef ;
proc logistic data=b.d simple plots(only)=oddsratio;
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;
On la retravaille en calculant la formule de p (à valoir plus loin) et on crée un indicateur visuel vp de
significativité à partir de la probabilité du Chi 2.
data b.estimPR;
set b.logisticest;retain intercept; ordre=_n_-1;
If _n_=1 then do;p=1/(1+exp(-Estimate));
intercept=estimate;variable ="_Intercept";
end;
if _n_>1 then p=1/(1+exp(-Intercept-Estimate));
ppct=round(p*100,0.1);
If abs(probchisq*100) >10.00 then vp="- ";
else do;
If abs(probchisq*100)