Statistical Procedures

Programming the statistical procedures from SAS
BookmarkSubscribeRSS Feed
Dicarlis
Obsidian | Level 7

Good morning, I'm currently analyzing a much larger database than I'm used to.
For this analysis I need to use the proc mixed with repeated measures, however when I test the covariance matrices and use the command DDFM=KR or DDFM=kr2 I have two problems:
1- the script does not work due to lack of memory even using options memsize=; independent of how much memory you use (max obs used 20G)
2- The DEN F appear as infinite,
However, I was unable to understand the differences between DDFM= betwithin, contain, kenwardroger, kenwardroger2, residual or satterthwaite.
Here is the model used:
PROC MIXED;
CLASS id_Animal TRT day;
MODEL visit = TRT dday trt*ddai / ddfm= ?? ;
RANDOM id_Animal ;
REPEATED day /TYPE = ?? SUBJECT = id;
RUN;

2 REPLIES 2
Ksharp
Super User

You could use "ddfm= BETWITHIN " and "RANDOM int /subject=subjectid" to save memory .


PROC MIXED;
CLASS id_Animal TRT day;
MODEL visit = TRT dday trt*ddai / ddfm= BETWITHIN  ;
RANDOM int /subject=id_Animal ;
REPEATED day /TYPE =ar(1) SUBJECT = id_Animal;
RUN;
SteveDenham
Jade | Level 19

If you have 20G observations, and have DEN F coming out as Inf, think about what that means - you are out at the point where F tests are indistinguishable from chi squared tests, and the t distribution doesn't differ from the normal until you get out to several decimal places. So, rather than trying DDFM= basically anything, try using the CHISQ option, and use those results rather than the F tests.

 

If you ran the analysis in GLIMMIX, you could use ddfm=none.

 

SteveDenham

sas-innovate-white.png

Join us for our biggest event of the year!

Four days of inspiring keynotes, product reveals, hands-on learning opportunities, deep-dive demos, and peer-led breakouts. Don't miss out, May 6-9, in Orlando, Florida.

 

View the full agenda.

Register now!

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.

Discussion stats
  • 2 replies
  • 1324 views
  • 4 likes
  • 3 in conversation