Discussion :

[Hiliadan]

Bonjour,

J'ai utilisé la proc logistic et le statement "OUTPUT" avec l'option "PRED" en voulant calculer les probabilités prédites. Mais lorsque je compare les valeurs obtenues à ce que j'ai calculé par ailleurs en utilisant la formule suivante :


... je n'obtiens pas du tout la même chose.

Pourriez-vous me renseigner sur la manière dont sont calculés les valeurs fournies par PRED ?
J'ai cherché dans la doc SAS sans succès.


Voici pour information mon code ci-dessous avec "PRED" en gras.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
proc logistic data=&table descending 
class propri(ref='1') sexepr(ref='2') mena(ref='3') CSP(ref='2') age(ref='1') revenu(ref='1') taille_copro(ref='1') etages(ref='1')
part_loc(ref='1') date_cons(ref='4') ancien(ref='1') zoneclim(ref='3') pieces(ref='1') Q16_1(ref='0') Q16_2(ref='0') 
Q16_3(ref='0') Q16_4(ref='0') Q17_1(ref='0') Q17_2(ref='0') Q17_3(ref='0') Q17_4(ref='0') Q17_5(ref='0') chauffage(ref='2')
eauchaude(ref='2') clim(ref='2') /*DPE_priv(ref='2') DPE_coll(ref='2')*/ CPE(ref='2') /*Q42_1(ref='0') Q42_2(ref='0') 
Q42_3(ref='0') Q42_4(ref='0') Q42_5(ref='0') Q42_6(ref='0') Q42_7(ref='0') Q42_8(ref='0')*/ cons_synd(ref='2') 
implication(ref='2') AG(ref='2') type_syndic(ref='1') satis_syndic(ref='1') niveau_info(ref='1') Q51_1(ref='0') 
Q51_2(ref='0') Q51_3(ref='0') Q51_4(ref='0') Q51_5(ref='0') Q51_6(ref='0') Q51_7(ref='0') Q51_8(ref='0') 
Q51_9(ref='0') Q51_10(ref='0') calcul_chauff(ref='1') tvx_coll(ref='0') tvx_priv(ref='0') tvx_coll_env(ref='0')  
tvx_coll_ouv(ref='0')  tvx_coll_cha(ref='0')  tvx_coll_ecl(ref='0') 
tvx_coll_sys(ref='0') tvx_priv_env(ref='0') tvx_priv_ouv(ref='0') tvx_priv_cha(ref='0') tvx_priv_sys(ref='0') 
solution(ref='2') plan(ref='2') fond(ref='2') Q115_1(ref='0') Q115_2(ref='0') Q115_3(ref='0') 
Q115_4(ref='0') Q115_5(ref='0') Q115_6(ref='0') Q115_7(ref='0') Q115_8(ref='0') Q115_9(ref='0') 
Q116_1(ref='0') Q116_2(ref='0') Q116_3(ref='0') Q116_4(ref='0') Q116_5(ref='0') Q116_6(ref='0') Q116_7(ref='0') 
Q116_8(ref='0') Q116_9(ref='0')  ;
model &travaux = propri sexepr mena CSP age revenu taille_copro etages part_loc date_cons ancien zoneclim 
pieces Q16_1 Q16_2 Q16_3 Q16_4 Q17_1 Q17_2 Q17_3 Q17_4 Q17_5 chauffage eauchaude clim /*DPE_priv DPE_coll*/ 
CPE /*Q42_1 Q42_2 Q42_3 Q42_4 Q42_5 Q42_6 Q42_7 Q42_8*/ cons_synd implication AG type_syndic satis_syndic niveau_info 
Q51_1 Q51_2 Q51_3 Q51_4 Q51_5 Q51_6 Q51_7 Q51_8 Q51_9 Q51_10 calcul_chauff 
solution plan fond Q115_1 Q115_2 Q115_3 Q115_4 Q115_5 Q115_6 Q115_7 Q115_8 Q115_9 
Q116_1 Q116_2 Q116_3 Q116_4 Q116_5 Q116_6 Q116_7 Q116_8 Q116_9 / link=logit selection=stepwise alpha=0.10
OUTROC=&travaux._roc 
CTABLE PPROB=&seuil
;
OUTPUT OUT=&travaux._p
PRED=&travaux._p ;
weight pond ;
run ;

[Silkus]

Bonjour,

Voir si c'est ce que tu cherches :

http://support.sas.com/documentation...c.logisticppec

A+

[Hiliadan]

Malheureusement non.
Cette page correspond plus à ce que je cherche: http://support.sas.com/documentation...ic_sect018.htm

Mais comme vous pouvez le voir, il n'y a pas de détails sur ce que fait "PRED".

[mi893]

Mais lorsque je compare les valeurs obtenues à ce que j'ai calculé par ailleurs en utilisant la formule suivante

Avec quoi l'as-tu calculé? Est-ce que tu as utilisé les mêmes estimés que ceux de SAS?

[Hiliadan]

J'ai utilisé la formule en image dans mon premier post pour calculer les probabilités prédites. Et j'ai donc pris les coefficients (notés b0, b1, ... dans la formule) calculés par SAS pour calculer les probabilités prédites.

Est-ce bien ce que tu voulais dire par les "estimés" de SAS ?

[charlotteb]

Bonjour,

J'avais le même problème.
En fait, quand tu définis tes classes il y a différentes paramétrisations et tout est expliqué là : http://support.sas.com/documentation...od_sect033.htm

En gros, il faut qu'à la fin de ta ligne class, tu mettes "/param=ref" , dans ce cas, ça te donnera les mêmes coefficients que ceux auxquels tu t'attends dans la formule de la régression logistique, et en comparant la proba prédite et en la calculant toi même tu obtiendras la même chose.

Le fait de paramétriser les classes d'une manière ou d'une autre ne change pas les probabilités prédites (tu vas obtenir la même chose) mais change les coefficients donnés et donc leur signification.