PROC Logistic : scénario

[kristofer]

Bonjour, comment on peut mesurer l'impact des scénarii sur des modèles logistiques? quelle procédure ou écriture faut-il utiliser ?

[RemiBousquet]

Je pense qu'un peu plus d'info sur ton besoin serait nécessaire pour que quelqu'un puisse t'aider, car là il manque des infos

[kristofer]

en fait, il s'agit d'un modèle explicatif où on a 2 variables explicatives.
on a utilisé une régression logistique.

on veut tester l'effet de variation de l'un des deux sur la variable à expliquer

[Manoutz]

Toujours pas clair, en tout vas pour moi. Qu'entends tu par

tester l'effet de variation de l'un des deux sur la variable à expliquer

. Quel type de variable explcative? Quel est ton objectif, plus globalement? pourquoi faire une régression logistiue?

[kristofer]

il s'agit de modéliser la demande de transport confrontant la voiture à l'autobus en fonction du cout et du temps des deux modes.
cad il faut estimer un modèle de choix binaire à l'aide de proc logistic .

cad

mode de transport {voiture} = a*cout voiture+b*cout bus+ c* temps voit+
{bus } d*temps bus

après l'estimation il faut calculer les élasticités directes et croisées (ou déja je suis bloqué) et puis l'effet de majoration de prix de carburant se traduisant par une augmentation de 30 % de cout d'utilisation de voiture et déduire l'effet de ce scénario sur la modede transport.

[MEGAMIND2]

Bonjour,
Essaies un test de Khi2 de chaque variable significative retenue avec la variable à modéliser (Y) et regarde lequel est le plus significatif.

Code :

1
2
3
PROC FREQ  DATA=TABLE_;
TABLES Y*Var/chisq;
run;

Cdt,

[MEGAMIND2]

Cette macro détermine la contribution de chaque variables endogènes retenues sur la variable exogène.
Quant à l'effet de variation, elle se mesure par les odds-ratio (l'exponentielle des bétas estimés).

Code :

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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
%********************************************************************;
/*base : table de construction du score                             */
/*beta : table contenant les coefficients associés à chaque variable*/
/* 		 option (outest=) de la proc logistic               */
/*quanti: liste des variables utilisées en continue ds le score     */
/*quali: liste des variables découpées ds le score                  */
%********************************************************************;
 
 
/*=====> les contributions sont restituées dans la table traitement*/
 
%macro r_contribution(base=,beta=,quanti=,quali=);
 
options nonotes;
 
proc sql noprint;/*calcul du nombre total de clients dans la base*/
SELECT count(*) INTO: nbtot FROM &base;
quit;
 
%LET i=1; %* calcul du nb de VARIABLES quanti;
%DO %while(%scan(&quanti,&i,%STR( )) ne %STR( ));
   %LET i=%EVAL(&i+1);
   %END;
%let nquanti=%eval(&i-1);
 
%LET i=1; %* calcul du nb de VARIABLES quali;
%DO %while(%scan(&quali,&i,%STR( )) ne %STR( ));
   %LET i=%EVAL(&i+1);
   %END;
%let nquali=%eval(&i-1);
 
%put ---------> &nquanti variable(s) quantitative(s); %* affichage du nb de variables quanti;
%put ---------> &nquali variable(s) qualitative(s); %* affichage du nb de variables quali;
%put ---------------------------------------------------------------;
%put ---------> les variables quantitatives:;
 
DATA traitement;
run;
 
%do i=1 %TO &nquanti;
%let var=%scan(&quanti,&i);  %* ième mot de la liste des VARIABLES quanti;
proc sql noprint;
SELECT &var INTO : valeur FROM β
SELECT css(&var) INTO : carre FROM &base;
quit;
%put ---------> la variable &i quantitative est : &var;
%put ---------> le coefficient associé est : &valeur;
 
DATA traitement_&i;
format facteur $10.;
facteur="&var";
variance=%sysevalf(&valeur*&valeur*&carre/&nbtot);
run;
 
DATA traitement;
SET traitement traitement_&i;
IF facteur="" then DELETE;
run;
 
proc datasets nolist;
DELETE traitement_&i;
quit;
 
%end;
%put ---------------------------------------------------------------;
%put ---------> les variables qualitatives:;
 
%do k=1 %TO &nquali;
%let var=%scan(&quali,&k);
 
DATA trav;
SET &base;
keep &var;
run;
 
proc sort DATA=trav;
BY &var;
run;
 
proc sql noprint;
SELECT DISTINCT &var INTO : modalite separated BY " " FROM trav;   %* insertion dans une macro-
variable des différentes modalités prises par la variable quali;
SELECT count(DISTINCT &var) INTO : nmodalite FROM trav;   %* nb de modalités de la variable quali;
quit;
 
%let nmod=&nmodalite;
%put ---------> la variable &k qualitative est : &var;
%let moyenne=0;
%let variance=0;
	%do j=1 %TO &nmod;
	%let nummod=%scan(&modalite,&j);
 
	proc sql noprint;
	SELECT &var&nummod INTO : valeur&nummod FROM β
	quit;
 
	%IF %datatyp(&nummod)=NUMERIC %then %do; /*verification du type de la modalite*/
	proc sql noprint;
	SELECT count(*) INTO : freq&nummod FROM trav WHERE &var=&nummod;
	quit;
	%end;
	%else %IF %datatyp(&nummod)=CHAR %then %do; 
	proc sql noprint;
	SELECT count(*) INTO : freq&nummod FROM trav WHERE &var="&nummod";
	quit;
	%end;
 
		%put ---------> le coefficient associé de la modalité &nummod de &var est : &&valeur&nummod;
	%let moyenne=%sysevalf(&moyenne+&&valeur&nummod*&&freq&nummod/&nbtot);
	%end;
 
	%do j=1 %TO &nmod;
	%let nummod=%scan(&modalite,&j);
	%let variance=%sysevalf(&variance+&&freq&nummod/&nbtot*(&&valeur&nummod-&moyenne)*(&&valeur&nummod-&moyenne));
	%end;
 
DATA traitement_&k;
format facteur $10.;
facteur="&var";
variance=&variance;
run;
 
DATA traitement;
SET traitement traitement_&k;
IF facteur="" then DELETE;
run;
 
proc datasets nolist;
DELETE traitement_&k trav;
quit;
 
%end;
 
/***calcul de la contribution de chacun***/
DATA traitement;
SET traitement;
std=sqrt(variance);
run;
 
proc sql noprint;
SELECT sum(std) INTO:ecart_tot FROM traitement;
quit;
 
DATA traitement;
SET traitement;
std_tot=&ecart_tot;
run;
 
DATA traitement;
SET traitement;
contribution=std/std_tot;
run;
 
options notes;
 
%MEND r_contribution;