ANOVA with a random effect from clusters

[fengyuwuzu]

For clustered data (multiple test sites from the same subject), if I want to do one way ANOVA to compare the baseline values in each group, is the followint the right way to add the random effect from clustering in the ANOVA model ?

proc glm 
  data=mydata;
  class treatment subjid;
  model change = treatment;
  random subjid;
run;

[lvm · Accepted solution]

You really (REALLY) don't want to use GLM for any situation with random effects. This was OK before procedures like MIXED and GLIMMIX were written, but GLM should be avoided now. For your example, use:

proc mixed data=mydata;
  class treatment subjid;
  model change = treatment;  
  random subjid;
lsmeans treatment / diff;
run;