Proc logistic : résultats différents avec param=ref et option noint , comparés à param=glm et noint

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
Pièce jointe 534878

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
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;