TP SAS 6 : Régression logistique binomiale - Lsta
TP SAS 6 : Régression logistique binomiale - Lsta
TP SAS 6 : Régression logistique binomiale - Lsta
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Régression <strong>logistique</strong> de Satellite01 par le poids et la couleur<br />
13. Définir formellement le modèle de régression <strong>logistique</strong> qui explique la présence d’un satellite<br />
par les variables poids et dark.<br />
14. Utilisez le code ci-dessous pour ajuster ce modèle de régression <strong>logistique</strong> sur les données. Y-at-il<br />
un ≪ effet carapace foncée ≫ Comparez graphiquement les probabilités a posteriori pour les<br />
carapaces foncées et claires.<br />
proc logistic data = RegLog.Limule plots(only) = (roc oddsratio effect);<br />
title ’Régression <strong>logistique</strong> avec poids et dark’;<br />
class dark;<br />
model satellite01(event = ’1’ ) = poids dark / expb ctable pprob= 0.5;<br />
oddsratio dark;<br />
oddsratio poids;<br />
run;<br />
Régression <strong>logistique</strong> et sélection de variable<br />
15. Discutez les sorties fournies par la code ci-dessous :<br />
proc logistic data = RegLog.Limule plots(only) = (roc oddsratio effect);<br />
title ’Régression <strong>logistique</strong> avec les quatre variables’;<br />
class epine dark;<br />
model satellite01 = poids largeur epine dark / ctable pprob= 0.5;<br />
oddsratio dark;<br />
oddsratio poids;<br />
oddsratio largeur;<br />
oddsratio epine;<br />
effectplot slicefit(sliceby=epine plotby=dark);<br />
run;<br />
16. Effectuez une sélection de variable par un algorithme backward pour identifier un modèle de<br />
régression <strong>logistique</strong> parcimonieux en utilisant l’option selection = stepwise dans l’instruction<br />
model :<br />
proc logistic data = RegLog.Limule plots(only) = (roc oddsratio effect);<br />
title ’Régression <strong>logistique</strong> et sélection de variable’;<br />
class epine dark;<br />
model satellite01 = poids largeur epine dark / ctable pprob= 0.5<br />
selection = backward ;<br />
run;<br />
4