Hi, I am trying to get the marginal difference in mean cost and 95% CI for an exposure with 5 groups using a GLM with log link and gamma distribution and have tried lsmeans, estimate, etc but all the results are off. I also tried %NLEstimate but it did not work for me. Any advice? Thank you! proc genmod data=t0 ; class expa(ref="1") agegroup(ref='1') female(ref='1') region(ref='1') flag_ad(ref='0') flag_ohter_comb(ref='0')/param=ref; model total_cost=expa agegroup female region flag_ad flag_ohter_comb/dist=gamma link=log type3; store p1; output out=c predicted=pred_c; run; expa is a categorical variable with 5 levels I am trying to get the exponential transformed cost and 95% CI for expa= 2 vs. 1, 3 vs. 1, 4 vs. 1, 5 vs. 1 /*Test incorrect results*/ proc genmod data=t0 ; class expa(ref="1") agegroup(ref='1') female(ref='1') region(ref='1') flag_ad(ref='0') flag_ohter_comb(ref='0')/param=ref; model total_cost=expa agegroup female region flag_ad flag_ohter_comb/dist=gamma link=log type3; estimate 'expa=1' intercept 1 expa 1 agegroup 1 female 1 region 1 flag_ad 0 flag_ohter_comb 0 / exp; estimate 'expa=2' intercept 1 expa 2 agegroup 1 female 1 region 1 flag_ad 0 flag_ohter_comb 0 / exp; estimate 'expa=3' intercept 1 expa 3 agegroup 1 female 1 region 1 flag_ad 0 flag_ohter_comb 0 / exp; estimate 'expa=4' intercept 1 expa 4 agegroup 1 female 1 region 1 flag_ad 0 flag_ohter_comb 0 / exp; estimate 'expa=5' intercept 1 expa 5 agegroup 1 female 1 region 1 flag_ad 0 flag_ohter_comb 0 / exp; lsmestimate expa "2 vs. 1" 2 1 / cl; lsmestimate expa "3 vs 1" 3 1/ cl; lsmestimate expa "4 vs 1" 4 1/ cl; lsmestimate expa "5 vs 1" 5 1/ cl; run;
... View more