Bonjour,
J'ai comparé les résultats de différentes paramétrisation d'une logistique, en particulier la paramétrisation ref ou glm (je n'ai pas testé 'effect'), et l'ajout de l'option noint :
a) Param=ref ref=first
b) Param=glm ref= first
c) Param= ref ref=first avec l'option /noint dans le model statement
d) Param= glm ref=first avec l'option /noint dans le model statement
Le modéles A, B et D sont consistants. A et B montrent presque la même chose ; dans A) la catégorie de référence est implicite, alors que dans B) elle est visible.
Dans le modèle D (glm et pas d'intercept) , l'intercept s'ajoute aux valeurs de beta , et finalement on obtient les mêmes Odds Ratios et contrastes.
Par contre, je ne m'explique pas pourquoi le modèle C (param=ref , option noint) donne des résultats différents. Ces résultats sont en contradiction avec les autres, et je suppose que ma paramétrisation doit être incorrecte ? ou faut-il éviter d'utiliser param=ref avec noint?
Ci-dessous le code et les résultats dans le lien . Merci beaucoup pour vos suggestions, Françoise
Summary_examples_logistic.htm
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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51 title 'A. param=ref; model /expb'; PROC SURVEYLOGISTIC DATA=test_allyr; WEIGHT WFIN; CLUSTER NUM_MENA;STRATA PROVW; where hc04=2; class year age8/param=ref ref=first; MODEL &ind (event='1') = year age8 /EXPB ; contrast '2001' intercept 1 year 0 0 0 age8 &agecoef/estimate=exp e; contrast '2004' intercept 1 year 1 0 0 age8 &agecoef/estimate=exp e; contrast '2008' intercept 1 year 0 1 0 age8 &agecoef/estimate=exp e; contrast '2013' intercept 1 year 0 0 1 age8 &agecoef/estimate=exp e; ods output parameterestimates=param_a oddsratios=or_a surveylogistic.contrastestimate=contrastyear_A; run ;title; title 'B. param=glm and ref=first; model /expb '; PROC SURVEYLOGISTIC DATA=test_allyr; WEIGHT WFIN; CLUSTER NUM_MENA;STRATA PROVW; where hc04=2; class year age8/param=glm ref=first; MODEL &ind (event='1') = year age8 /EXPB ; contrast '2001' intercept 1 year 0 0 0 1 age8 &agecoef/estimate=exp ; contrast '2004' intercept 1 year 1 0 0 0 age8 &agecoef/estimate=exp ; contrast '2008' intercept 1 year 0 1 0 0 age8 &agecoef/estimate=exp ; contrast '2013' intercept 1 year 0 0 1 0 age8 &agecoef/estimate=exp ; ods output parameterestimates=param_b oddsratios=or_b surveylogistic.contrastestimate=contrastyear_b; run ;title; title 'C. param=ref ref=first ; model /noint expb';/*NO*/ PROC SURVEYLOGISTIC DATA=test_allyr; WEIGHT WFIN; CLUSTER NUM_MENA;STRATA PROVW; where hc04=2; class year age8/param=ref ref=first; MODEL &ind (event='1') = year age8 /noint EXPB ; contrast '2001' year 0 0 0 age8 &agecoef/estimate=exp e; contrast '2004' year 1 0 0 age8 &agecoef/estimate=exp e; contrast '2008' year 0 1 0 age8 &agecoef/estimate=exp e; contrast '2013' year 0 0 1 age8 &agecoef/estimate=exp e; ods output parameterestimates=param_c oddsratios=or_c surveylogistic.contrastestimate=contrastyear_c; run ;title; title 'D. param=glm ref firs; model /noint expb '; PROC SURVEYLOGISTIC DATA=test_allyr; WEIGHT WFIN; CLUSTER NUM_MENA;STRATA PROVW; where hc04=2; class year age8/param=glm ref=first; MODEL &ind (event='1') = year age8 /noint EXPB ; contrast '2001' year 0 0 0 1 age8 &agecoef/estimate=exp e; contrast '2004' year 1 0 0 0 age8 &agecoef/estimate=exp e; contrast '2008' year 0 1 0 0 age8 &agecoef/estimate=exp e; contrast '2013' year 0 0 1 0 age8 &agecoef/estimate=exp e; ods output parameterestimates=param_d oddsratios=or_d surveylogistic.contrastestimate=contrastyear_d; run ;title;
Partager