Programming the statistical procedures from SAS

Proc mixed effect size: Does a solution exist?

Reply
Senior User
Posts: 1

Proc mixed effect size: Does a solution exist?

I have spent so much time trying to find out how to calculate an effect size for linear mixed models, but still success evades me. I am an avid SPSS user. Given the mixed model limitation in SPSS and lack of my success, I finally decided to transition to SAS in hopes of possibly overcoming this problem with the effect size. Is there anyone who has any clue how to calculate an effect size for lognitudinal linear mixed models with complex covariance structures like ARH1, ante(1), AR(1), etc?

 

Here is a sample code:

proc mixed data = sas.aLTP covtest scoring=5;
class dose Mouse Time;
model ln_aLTP = dose cTime cTime*cTime / solution ddfm=sat influence;
repeated Time /type=ante(1) subject=mouse;
random int cTime / subject = mouse type = arh(1) solution g v; 
lsmeans dose /adjust=dunnett pdiff=control('4') CL;
run;

I am aware of this thread for Cohen's d in proc mixed, which did not generate a solid resolution.

 

I'm also aware of this paper: A practical guide to calculating Cohen’s f2, a measure of local effect size, from PROC MIXED, but that covers only the simplest covariance structures. I have emailed various authors and  still have not found a practical solution. So I'm wondering if any SAS users have somehow overcome this issue, since the community is large.

 

Some more specific questions to guide the discussion:

 

  1. If I want to at least report Cohen's d for the difference between two means in lsmeans, does it sound resonable to run a t-test for those two means of interesy using the predicted values from proc mixed and then calculate Cohen's d?
  2. As far as converting SE to SD, are there actual sources that validate such an approach? The thread referenced above mentions SD=SE*SQRT(DF+1). Is this a justifyable approach?
  3. One of the authors who didn't have much time to explain his solution to me, said that he would compute the effect size as the ratio of the estimated fixed effect parameter to the sqrt of total marginal outcome (at a given timepoint if marginal variance isn't constant). It seems that SAS can produce the marginal variance output. Can anyone elaborate on this and perhaps give some tips about syntax? Do I need to add "marginal" or "var" somewhere in the syntax to get that output? I am not sure what would such an effect size be called in terms of the standardized variants, like Cohen's f sq, pseudo-r sq, or something else.

 

Frequent Contributor
Posts: 130

Re: Proc mixed effect size: Does a solution exist?

The reference to marginal variance likely refers to proc glimmix not proc mixed, but that's ok, both will estimate your models.

 

See

 

http://support.sas.com/resources/papers/proceedings16/11663-2016.pdf

 

for an example of specifying, exporting and calculating other statistics from proc glimmix covariance estimate output.

Frequent Contributor
Posts: 130

Re: Proc mixed effect size: Does a solution exist?

I think this posting is also relevant to your inquiry:

 

https://communities.sas.com/t5/SAS-Statistical-Procedures/Proc-Mixed/m-p/296000#U296000

Ask a Question
Discussion stats
  • 2 replies
  • 381 views
  • 0 likes
  • 2 in conversation