Turn on suggestions

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

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Calculate standard deviations by proc mixed!!

Options

- 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 06-18-2014 02:20 PM
(12522 views)

Dear all,

Is there a way to obtain standard deviation (SD) by proc mixed, for each lsmean values instead (or plus) of standard error (SE)?

proc mixed data=have method=REML;

class id year province y;

model y = year province;

random id;

lsmeans year;

run;

Thanks in advance

zana

1 ACCEPTED SOLUTION

Accepted Solutions

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

Ah. Cohen's D.

Try looking at:

http://www.psy.mq.edu.au/psystat/documents/standardised_effect_size_in_mixed_ML_models.pdf

Also read Jake Westfall's blog on the subject:

It comes down to the notion that in mixed models, there is more than one source of variability, and a composite standard deviation is needed if you really want to calculate a standardized effect size.

Steve Denham

4 REPLIES 4

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

Sorry, i had a mistake in my model.

The true model equation is;

**proc mixed data=have method=REML;**

** class id year province;**

** model y = year province;**

** random id;**

** lsmeans year;**

** run;**

Regards

zana

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

The standard error reported is the standard deviation of the population of means, conditional on the random effects and marginal over differences in sample size. Additionally, under the model presented, you assume homogeneity of variances by year and province, so that standard errors of each year will be identical, as will the standard errors of each province (assuming a balanced design)..

I suppose it could be done by multiplying the standard error by sqrt(df + 1), but I am curious as to why you would do this. The value obtained has little to do with anything related to the sample or what might be considered the originating population. Why not just use the sample standard deviation reported by PROC MEANS?

In the end, do you want heterogeneous estimates of variability? That would involve use of the repeated statement--see Example 63.7 of The MIXED Procedure.

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

Hi Steve

Could you please explain this statement:

" The value obtained has little to do with anything related to the sample or what might be considered the originating population. "

Why wouldn't the SD have anything to do with the originating population. Wouldn't a person want to know the estimated mean and estimated SD to calculate estimated Cohen's Ds?

All the best.

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

Ah. Cohen's D.

Try looking at:

http://www.psy.mq.edu.au/psystat/documents/standardised_effect_size_in_mixed_ML_models.pdf

Also read Jake Westfall's blog on the subject:

It comes down to the notion that in mixed models, there is more than one source of variability, and a composite standard deviation is needed if you really want to calculate a standardized effect size.

Steve Denham

**Available on demand!**

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

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.