BookmarkSubscribeRSS Feed
Lyson
Quartz | Level 8

Hi friends,

I am trying to estimate a linear mixed model in SAS using prc mixed or proc glimmix. I have 3 categorical factors A, B and C where factor C is assumed random, and A,B are fixed. My data set looks like this:

Lyson_0-1721308307376.png

 

The  syntaxes that I am exploring are:

a) PROC GLIMMIX

proc glimmix data=FFRgroup outdesign=XZ;;
class factorA factorB factorC  groupA groupB groupC rep;
model yield= groupA | groupB / ddfm=satterthwaite;
random factorC factorA*factorC  factorB*factorC factorA*factorB*factorC /  V;
lsmeans groupA | groupB   / pdiff cl alpha=0.05 ;
Title "FFR Model";
run;

whch displays the following V matrix,

Lyson_2-1721307720866.png

 

b) PROC  MIXED:

proc mixed data=FFRgroup ;
class factorA factorB factorC  groupA groupB groupC rep;
model yield= groupA | groupB / ddfm=satterthwaite;
random factorC factorA*factorC  factorB*factorC factorA*factorB*factorC /  V;
lsmeans groupA | groupB   / pdiff cl alpha=0.05 ;
Title "FFR Model";
run;

displays the following V matrix,

Lyson_3-1721307827866.png

What could be the problem with these two approaches not producing the covariance structure that aligns with the model in question (block diagonal)? 

 

Thank you in advance for your assistance folks.

1 REPLY 1

sas-innovate-2026-white.png



April 27 – 30 | Gaylord Texan | Grapevine, Texas

Registration is open

Walk in ready to learn. Walk out ready to deliver. This is the data and AI conference you can't afford to miss.
Register now and lock in 2025 pricing—just $495!

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
  • 1 reply
  • 892 views
  • 0 likes
  • 2 in conversation