Programming the statistical procedures from SAS

distribution option in Proc Glimmix

Reply
Contributor
Posts: 66

distribution option in Proc Glimmix

[ Edited ]

Hi. I chose to run proc glimmix instead of proc mixed because when I test shapiro wilk's test in proc mixed, the normality assumption was violated even after tranformation.

However, when I choose distribution of gamma instead of default distribution (gaussian), the model does not converge.

Is there any solution for this?

Should I still seek a model to converge when the distribution is gamma not gaussian?

Or can I go with gaussian distribution even the normality assumption was violated?

Thanks for your advice in advance.

Super User
Posts: 10,787

Re: distribution option in Proc Glimmix

You could try different initial value. Check PARAM statement.

Esteemed Advisor
Posts: 5,540

Re: distribution option in Proc Glimmix

Given enough observations, the normality checking tests will always give significant results. First, make sure you don't have an outlier problem. Then, if your residuals distribution has a single mode and not too much skewness, go with the gaussian analysis.

 

If you really have to model your data as gamma, then post your code for more help.

PG
Contributor
Posts: 66

Re: distribution option in Proc Glimmix

So do you mean it is not a problem to do gaussian distribution even when the purpose of running proc glimmix is because I thought my data is non-normally distributed?

 

Thanks.

SAS Super FREQ
Posts: 4,245

Re: distribution option in Proc Glimmix

> So do you mean it is not a problem to do gaussian distribution even when the purpose of running proc glimmix is because I thought my data is non-normally distributed?

 

It is not a problem, but I think you are misquoting the purpose of glimmix.

The assumptions for general linear regression models are that the residuals (not the data) are normally (and independently) distributed with constant variance. Generalized linear models (eg, PROC GENMOD) support alternative residual structures such as Poisson regression where the variance is assumed to vary, as well as non-normal error distributions.  PROC GLIMMIX supports generalized linear models with random effects. 

Ask a Question
Discussion stats
  • 4 replies
  • 161 views
  • 0 likes
  • 4 in conversation