BookmarkSubscribeRSS Feed
sodonnell
Calcite | Level 5

I ran a model using proc logisitic and proc surveylogistic and upon comparing the results I found the OR estimates to be same (as one would expect) however, the parameter estimates were different.


Upon taking the EXP of the parameter estimates from the proc logistic output I was able to reproduce the ORs but I was not able to do the same with the parameter estimates from the proc survey logistic output.  As a result, I believe the parameter estimates from proc logistic to be true.

Does anyone know why the parameter estimates would differ?

4 REPLIES 4
Reeza
Super User

Code for each proc?

sodonnell
Calcite | Level 5

proc surveylogistic  data=Mergedset_rev /*VARMETHOD=BRR;*/

class sex age2 marital edu_r Inc_Quin immigt aborig employ_4 cc_alt area smoke alco;

model sub2 (event='2')=sex age2 marital edu_r Inc_Quin immigt aborig employ_4 cc_alt area smoke alco/EXPB clodds alpha=.05;

ods output CLOddsWald=OR ParameterEstimates=parm;

weight wts_s;

/*REPWEIGHTS bsw1-bsw500;*/

run;

proc logistic descending data=Mergedset_rev;

model sub2= female _15_29 _30_49 w_s_d single le_sec_r sec_r so_pos_r inc_q1 inc_q2 inc_q3 inc_q4 imm abor no_work students onetw_ch three_ch daily occfor alcoh rural;

estimate 'female' female 1/ exp;

estimate '_15_29' _15_29 1/exp;

estimate '_30_49' _30_49 1/exp;

estimate 'w_s_d' w_s_d 1/ exp;

estimate 'single' single 1/exp;

estimate 'le_sec_r' le_sec_r 1/ exp;

estimate 'sec_r' sec_r 1/ exp;

estimate 'so_pos_r' so_pos_r 1/ exp;

estimate 'inc_q1' inc_q1 1/ exp;

estimate 'inc_q2' inc_q2 1/ exp;

estimate 'inc_q3' inc_q3 1/ exp;

estimate 'inc_q4' inc_q4 1/ exp;

estimate 'imm' imm 1/ exp;

estimate 'abor' abor 1/ exp;

estimate 'no_work' no_work 1/ exp;

estimate 'students' students 1/ exp;

estimate 'onetw_ch' onetw_ch 1/ exp;

estimate 'three_ch' three_ch 1/ exp;

estimate 'daily' daily 1/ exp;

estimate 'occfor' occfor 1/ exp;

estimate 'alcoh' alcoh 1/ exp;

estimate 'rural' rural 1/ exp;

weight wts_s;

run;

ballardw
Super User

Since there is significant difference in the two model statements I wouldn't expect much to be the same...

Reeza
Super User

Assuming you've done everything else right, it might be the common mistake of not adding param=ref to the class statement in proc surveylogistic. This creates parameter estimates that via effect coding not referential coding which is probably what you're after.

Check the docs for further information. Also, if you're using weights you should be using surveylogistic, again check the docs for the reasons.

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 4 replies
  • 4680 views
  • 0 likes
  • 3 in conversation