topic Effect size in GLIMMIX in Statistical Procedures

Effect size in GLIMMIX

gvandewalle — Tue, 29 May 2018 15:05:36 GMT

Dear Madam, dear Sir,

I would like to estimate effect size in PROC GLIMMIX (in SAS 9.4)

I am seeking for an association between a continuous dependent variable (normal distribution) and a series of independent predictors.

Predictors consist in 1 class variable (sex) and 6 continuous variables (age, education, etc.).

Some predictors are correlated with one another but not redundant.

How could I get something which ressemble an effect size estimation (r^squared, r, Cohen's, etc.) for each of the (significant) predictors (there are at least 2 significant predictors)?

We tried to standardize all variables prior to computing the model. We do get (beta) estimates that are in what seem to be the same scale. Can this be seen as some sort of effect size measure?

Thanks in advance for your help.

Regards,

Gilles

Re: Effect size in GLIMMIX

PaigeMiller — Tue, 29 May 2018 17:17:55 GMT

For categorical predictor variables, the effect size can be obtained by the LSMEANS statement.

For continuous predictor variables, the effect size (which is actually the slope) can be obtained by using the SOLUTION option in the MODEL statement.

You ask for R-squared and similar measures, but these are not "effect size" measures in any reasonable usage of the term "effect size"; furthermore GLIMMIX does not compute R-squared. If you really want R-squared, then PROC GLM could be used if it fits the requirements of the problem.

Re: Effect size in GLIMMIX

djw23 — Sun, 27 Jan 2019 21:14:07 GMT

Hi,

I'm trying to get effect size values for my mixed model using PROC GLIMMIX. I used SOLUTION in the MODEL statement as suggested but I'm not sure how to interpret the output of this. Could someone please help?

Thanks,

Dave

Re: Effect size in GLIMMIX

PaigeMiller — Sun, 27 Jan 2019 22:16:56 GMT

@djw23 wrote:

Hi,

I'm trying to get effect size values for my mixed model using PROC GLIMMIX. I used SOLUTION in the MODEL statement as suggested but I'm not sure how to interpret the output of this. Could someone please help?

I believe I explained this in my earlier post in this thread. If that doesn't help, please be more specific, show us the output, and then state exactly what you don't understand in the output.

Re: Effect size in GLIMMIX

djw23 — Tue, 29 Jan 2019 14:42:44 GMT

Hi there,

So some background to start.

I have an dataset with several fixed effects and I'm measuring all 2 and 3-way interactions. As I have a random effect too, I'm using proc glimmix to analyse my data in a generalised linear mixed model. I also have a gamma distribution with a log link.

What I wish to do is calculate effect size for each of my fixed effects and interactions. I've tried adding the oddsratio function into my code, but this is the message i get:

NOTE: Odds ratios are computed only for the logit, cumulative logit, or the generalized logit link function.

When I use the SOLUTION function that you suggested, what I wish to know is if the effect size corresponds to the estimate column and how I should interpret it? Also, why are there zeros next to certain levels of my fixed effects?

Thanks in advance

David

Re: Effect size in GLIMMIX

PaigeMiller — Tue, 29 Jan 2019 16:32:23 GMT

Also, why are there zeros next to certain levels of my fixed effects?

Let's suppose you have a variable that has two levels, male and female, and your analysis is of the heights (in inches) of the subjects. Then the results COULD BE that males have an effect of +5 and that females have an effect of –5. Or, the same information could be represented as males +10 and females 0; or it could be represented as males 0 and females –10. These are all the same representations of the results, there is no unique representation. So, SAS chooses one of these, using the rule that the level which is last alphabetically gets a zero.

This is why I recommend you really really really really really really really really really really want to look at the LSMEANS, which doesn't have this problem, and are much easier to interpret.

what I wish to know is if the effect size corresponds to the estimate column and how I should interpret it?

If the variable is categorical, don't interpret it. Use LSMEANS and interpret those. If the variable is continuous, then the effects are the slopes estimated in the model.

Re: Effect size in GLIMMIX

cedric5 — Thu, 19 Mar 2020 18:14:37 GMT

Hi,

What do you look at in the LSMEANS output to get the effect size?

Where do you find the slope for the continuous predictor in the Solution output?

Thank you for your help.

Cedric

Re: Effect size in GLIMMIX

PaigeMiller — Thu, 19 Mar 2020 18:20:02 GMT

@cedric5 wrote:

Hi,

What do you look at in the LSMEANS output to get the effect size?

Where do you find the slope for the continuous predictor in the Solution output?

Thank you for your help.

Cedric

Please start a new thread with your question, as it doesn't belong in this thread. In your new thread, please be sure to show us the output you are looking at.