Calcite | Level 5

Need help with PROC GENMOD /Gamma distribution with log link

Hi,

I am building a multivariable model using proc genmod (gamma and a log link), and I have a major explanatory variable (X1, which is 5-level category variable) and a few covariates.

When I got the output, I had estimates for intercept B0, B1(X1), and estimates for other covairables....but I could not find the overall mean value for the dependent variable in the whole model. Also, I want to calculate the group means by variable X1 in the model.

Are there any ways to calculate the overall mean of dependent variable and the group means by X1 in the gamma and a log link model?

Thanks a lot!
4 REPLIES 4
Pyrite | Level 9

Re: Need help with PROC GENMOD /Gamma distribution with log link

In order to produce estimates of the response for each level of X1, you need to specify X1 as a categorical variable and include an LSMEANS statement which names X1 as the effect of interest. Since you have specified a log link, you will want to include the ILINK option on the LSMEANS statement. I trust that you have a recent update to SAS so that the ILINK option is available to you.

Providing specific help would be much easier if you were to provide more details about your problem, including the version of SAS which you have available as well as code which you have submitted.
Calcite | Level 5

Re: Need help with PROC GENMOD /Gamma distribution with log link

Thank you very much for explanation. Now I understand how to get the group means!
Rhodochrosite | Level 12

Re: Need help with PROC GENMOD /Gamma distribution with log link

Here is one simple example, with trt for a factor (categorical variable) and x for a covariable (continuous).

proc genmod data=ra;
class trt;
model y = trt x / dist=gamma link=log ;
lsmeans trt / ilink at x=-1;
run;

With LSMEANS, you get estimated means (based on the model) for each covariable at its overall mean. At least you do with MIXED and GLIMMIX, and I assume this holds with GENMOD. You may want the estimated values at a different value for the covariable. Here I asked for the values when x=-1 (example).
Calcite | Level 5

Re: Need help with PROC GENMOD /Gamma distribution with log link

That works! Thank you so much!
Discussion stats
• 4 replies
• 2741 views
• 0 likes
• 3 in conversation