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

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Re: generalized linear mixed model 'proc glimmix'

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

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 07-26-2021 03:39 PM
(921 views)

Hello,

I appreciate for any guidance that help me understand how to check the unconditional (marginal) distribution which is assumed to be known for maximum likelihood estimation in generalized linear mixed model (proc glimmix) for multilevel repeated measures data with random effect. If I'm right one assumption for proc glimmix is that the random effect variables should follow normal distribution, so does the variable that has random effect play the role of (marginal) distribution ? Is there any article or any source that explains how to check the normality of marginal distribution and assumptions in proc glimmix and what is the concept of marginal distribution in generalized linear mixed model ?

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

The first sentence of that article says that the discussion is about "how to check the assumptions of an **ordinary least squares** linear regression model." Some of the discussion is relevant to the **fixed** effects in generalized linear models.

*> The variables in a least squares regression model do not have to be normally distributed. Do the variables here mean independent variables and dependent variable ( response )?*

Yes, that is correct. Neither the explanatory variables nor the response variable is required to be normal in OLS. As discussed in the article, the residuals of the model are assumed to be normal

*> Does this rule apply to mixed model( proc mixed and proc glimmix )*

There is not a requirement that any of the **fixed-effect** variables (or response variable) in mixed models be normal. The estimation method does not matter. Generalized models do not require the normality of residuals.

The **random effects** in mixed models are assumed to be normally distributed with mean zero and a correctly specified covariance matrix. We usually estimate the variance-covariance components.

You can search the internet for details about the model assumptions for linear mixed models. Some of the course notes from professors (such as Haelwig's Linear Mixed-Effects Regression (umn.edu)) are good resources.

5 REPLIES 5

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Calling @lvm @StatDave @SteveDenham @Rick_SAS

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Marginal: population averaged

Conditional: subject 'averaged"

Marginal distribution: In GLIMMIX, this is specified as an R-side distribution of the residuals. It can be whatever you choose from the available distributions. The variance components are assumed to be normal.

Conditional distribution: In GLIMMIX, this is specified as a G-side distribution. Again, you may choose an appropriate distribution, but take care that it is appropriate at the subject level.

The Introduction to Mixed Models in the SAS/STAT documentation ought to help get you started.

SteveDenham

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

@SteveDenham can correct me if I'm wrong, but I don't see many people trying to "verify the assumptions" for random effects. The normality assumptions for random effects can't be checked as easily as the assumptions for fixed effects. One reason is that the random effects are often "groups" such as schools or hospitals, and there might only be a dozen groups. Another is that the model will be more heavily influenced by violations in the assumptions of the fixed effects, and those are rarely satisfied exactly (outside of simulation studies).

For a discussion of random effects and what happens if they do not satisfy the assumptions, see Schielzeth et all (2020). The paper concludes that "Our simulation analysis shows that the effect of violations of distributional assumptions of random effect variances and residuals is surprisingly small." They conclude that "mixed-effects models are largely robust even to quite severe violations of model assumptions."

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

according to this article " https://blogs.sas.com/content/iml/2018/08/27/on-the-assumptions-and-misconceptions-of-linear-regress... " ,**The variables in a least squares regression model do not have to be normally distributed**. Do the variables here mean independent variables and dependent variable ( response )? Dose this rule apply to mixed model( proc mixed and proc glimmix ) also, since mixed model is multilevel regression lines but estimation method can be maximum likelihood or (restricted) maximum likelihood?

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

The first sentence of that article says that the discussion is about "how to check the assumptions of an **ordinary least squares** linear regression model." Some of the discussion is relevant to the **fixed** effects in generalized linear models.

*> The variables in a least squares regression model do not have to be normally distributed. Do the variables here mean independent variables and dependent variable ( response )?*

Yes, that is correct. Neither the explanatory variables nor the response variable is required to be normal in OLS. As discussed in the article, the residuals of the model are assumed to be normal

*> Does this rule apply to mixed model( proc mixed and proc glimmix )*

There is not a requirement that any of the **fixed-effect** variables (or response variable) in mixed models be normal. The estimation method does not matter. Generalized models do not require the normality of residuals.

The **random effects** in mixed models are assumed to be normally distributed with mean zero and a correctly specified covariance matrix. We usually estimate the variance-covariance components.

You can search the internet for details about the model assumptions for linear mixed models. Some of the course notes from professors (such as Haelwig's Linear Mixed-Effects Regression (umn.edu)) are good resources.

Build your skills. Make connections. Enjoy creative freedom. Maybe change the world. **Registration is now open through August 30th**. Visit the SAS Hackathon homepage.

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.