Discussion :

[indox]

Bonjour, j'ai obtenue les coefficients des variables à l'aide de la proc logistique. Je souhaiterais écrire dans SEG 4 à la main la formule permettant d'avoir la probabilité Y. Comment je dois écrire cette formule? Le problème est par exemple pour la variable sexe (homme/femme) pour laquelle je n'ai que le coeff de la modalité femme.

Merci pour votre aide.

[suistrop]

Bonjour,

Par défaut SAS n'affiche pas la modalité de référence dans les sorties. Tu peux utiliser l'option GLM pour avoir toutes tes valeurs en sortie:

Code :

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

class V1 V2 V3 V4 /  param = glm;

Cordialement,

[AllXS7117]

A partir de la sortie de ta proc logistic tu peux faire :

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
proc transpose data=sortie1 out=sortie2;
run;
data _null_;
set sortie1;
call symput("Intercept",Intercept);
call symput("proba0",100*1/(1+exp(-Intercept)) );
 
/*macro %proba0 : estimer proba de la situation de référence, et la proba de s'en écarter*/
run;
data sortie2;
set sortie2;
if _n_=1 then probe=1/(1+exp(-&Intercept.))*100;
         else do;
probe=1/(1+exp(-&Intercept.-var_dependante ))*100;
d=probe-&proba0.;
end;
 variable_class=compress(_label_);
run;
 
 
/*déterminer les odds ratio et les modalités significatives au seuil de 5%*/
data odds1;
  set odds;
  odds = exp(estimate);
  if ProbChiSq > 0.05 then Test="NS";else Test="-";
  variable_class=compress(variable||" "||ClassVal0);
run;