Programming the statistical procedures from SAS

Mix procedures question

Reply
New Contributor
Posts: 2

Mix procedures question

Hi everyone,

I am having an issue how to account for the extra error term in our mixed statistical model using the SAS mixed procedure. I want to treat it as a split-plot design to come up with an ANOVA table, and one error term should not include the factor L (see equation below).

My guess is something like:

Model.png

Proc Mixed data=StrainGauges method=type3; class rep F R C L;

model S=F R C FR  FC RC FRC L FL RL CL FRL FCL RCL FRCL/ddfm=kr;

Random rep(A|B|C)

However, I am not sure how to write the random statement (e.g. "Random rep(A*B*C)" or "Random rep(A|B|C)"). Can someone help me to write this code?

Thanks in advanced!

Josue

Respected Advisor
Posts: 2,655

Re: Mix procedures question

This is a case where you have several random effects, so the rep(A*B*C) is obviously not what you need (and it is just as obvious that it is the residual error).

A much better approach would be:

random intercept F R C F*R F*C R*C F*R*C/subject=rep;

Also in your model statement, be sure to represent interactions with the '*' operator.  As it is now, this would treat FR, FC, RC and FRC as continuous (and undefined) variables.

Steve Denham

New Contributor
Posts: 2

Re: Mix procedures question

Thank you, Steve. Your comments will be very useful for our analysis!

Josue Bastos

Ask a Question
Discussion stats
  • 2 replies
  • 224 views
  • 1 like
  • 2 in conversation