sas >> How do you calculate Eta Squared using Proc Mixed?

by hodab » Thu, 17 Feb 2005 01:56:34 GMT

Hi all:

I have computed a repeated measures multilevel models analysis in SAS
using Proc mixed on some EMA data.
It was a pretty simple analysis looking at the degree of association
between two variables: mood and pain with time as the repeated measure.
Well, the editor of the journal I am sending the article to would like
me to report effect sizes for the analysis. I know how to get SAS to
print out a 95% confidence interval in Proc Mixed, but I cannot get it
to print out the Type III Sums of Squares that I would need to
calculate Eta squared. My problem, is that Proc Mixed uses MLE and
does not calculate Type III SS. My understanding is that you can get
it to calculate them if you put in a "Method = Type 3" statement, but,
because I have a repeated models analysis and am already using the
"REPEATED" statement, I cannot do both. Do you know any way for me to
get SAS to calculate Eta Squared? If not, do you think that reporting
the confidence interval along with the p value would be sufficient?

sas >> How do you calculate Eta Squared using Proc Mixed?

by stringplayer_2 » Thu, 17 Feb 2005 03:36:05 GMT

What covariance structure is fit by your REPEATED statement. If
you have a TYPE=CS covariance structure, then you can use an
equivalent RANDOM statement. The following two models produce
identical results:

proc mixed data=mydata;
class id time;
model pain = mood / s;
repeated time / subject=id type=cs;

proc mixed data=mydata;
class id;
model pain = mood / s;
random id;

You could specify method=type3 to the latter model to get the
variance partitioning that you need to construct eta.

Alternatively, you can fit the model with and without mood.
Each produces an estimate of the residual variance. If we have

V1=residual variance from model without mood
V2=residual variance from model with mood

then an estimate of eta would be 1 - (V2/V1).



Dale McLerran
Fred Hutchinson Cancer Research Center
Ph: (206) 667-2926
Fax: (206) 667-5977

Do you Yahoo!?
Yahoo! Mail - You care about security. So do we.

sas >> How do you calculate Eta Squared using Proc Mixed?

by hodab » Thu, 17 Feb 2005 06:23:37 GMT

Dale -

Thanks so much for your help! It worked! :-)

take care,


Hoda Badr, Ph.D.
The University of Texas M. D. Anderson Cancer Center

Similar Threads

1. how to calculate sum of squares for treatments when using proc mixed

3. how to calculate sum of squares for treatments when using

4. score chi-square & wald chi-square in logistic regression using sas

What is the difference between score chi-square & wald chi-square in
logistic regression?

5. r-squared in proc mixed

6. proc mixed / Latin squares+ repeated measures

Hi listers,

I would like to know if a Latin square design with repeated measures
within each cell can be set up with PROC MIXED. The Latin square has as
blocking effects, period (row) and animal (column) (both random
effects), the response of each animal at each treatment is measured four
times on time.
What term can be used as error for the between subjects design?
If period was fixed, then could I use the interaction treatment*animal
as error for the Latin square?
Thanks a lot !

7. R square in proc mixed?

8. Using proc freq : add column after doing freq

Hi Every one,

i am doing following

proc freq data= ****;
tables a*b*e/ list out= abc;

i have normal abnormal and not reported in B,
 were as i have only normal and abnormal  in E, i dont have not
here a is 1;

this is how my output  looks like

Normal  |Abnormal|  Total
9          0              9
                                             Abnormal       |
0        3               3
                                             Not Reported *
1         0              1
10        31            13
      in data set abc:

                                    obs  A
B                          E      COUNT    PERCENT

                                     1      1    Normal
Normal          99        75.5725
                                     2      1    Abnormal
Abnormal      31        23.6641
                                     3      1    Not Reported *
Normal         1          0.7634

How do i add not reported in E so my o/p looks like

Normal  |Abnormal|  not reported  Total
9          0                0               9
                                             Abnormal       |
0        3                0                3
                                             Not Reported *
1         0               0                  1
10        31               0            13