turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Parameter estimates in proc glimmix?

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-21-2012 09:01 AM

My data are hospital cost data and I model them as a function of various cost drivers. I need to specify random effects and a lognormal or gamma distribution for my data, therefore I found proc glimmix convenient.My model converges with a lognormal distribution and link=id. I found it very difficult to interprete the parameter estimates. An OLS regression without any explanatory variables renders an intercept of around 77,000. My Glimmix model renders an intercept of about 10,500 (regardless of covariates and random effects). Clearly exponentiating brings these parameters close to the 'true' figure but firstly I can't find any literature telling me to exponentiate and secondly all parameters turn positive when exponentiating. I've spent some time on the internet trying to find solutions.

I hope you can help. Thanx

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to Marie

02-22-2012 08:30 AM

Do not exponentiate any single parameter estimate. The lognormal distribution refers to the distribution of the dependent variable, and so you would exponentiate the linear combination you get after plugging in appropriate values in the X matrix. The only time I could understand exponentiating the intercept is if all independent variables are zero. Given that, try centering any covariates and see what happens to the intercept. It would be good to think of the intercept then as estimating the median. See some recent posts in this community on GLIMMIX for more info.

Steve Denham

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to SteveDenham

02-22-2012 08:49 AM

Thanks Steve, it makes a lot of sense.

I gather from your reply that I basically can't use the individual parameter estimates, the way I normally do in proc reg and proc glm. Which is a pity since glimmix is otherwise quite approriate.

Thanks again

Marie

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to Marie

02-22-2012 08:57 AM

Well you can still use the parameter estimates, just remember that a shift of 1 due to an X variable, holding all other variables constant, leads to a shift of exp(beta) in the y variable--if you use a log link as is standard for a gamma distribution. I'll have to think about the dist=lognormal, as the link is the identity there. And by think, I mean go plug in some data and look at the results.

Hope this helps.

Steve Denham