Programming the statistical procedures from SAS

proc mixed split plot regression

Posts: 1

proc mixed split plot regression

Split plot design with a ordinal whole-plot factor.

I have an experimental design where the whole plot factor (call it A) has five unreplicated ordinal levels. Imagine this as five water tanks each maintained at one of five different temperatures.  Within each is a replicated treatment (call it B) with two levels and continuous dependent variable.  So instead of the typical split plot with each level of factor A [whole-plot] minimally replicated, I have five unreplicated ordinal levels (temperature) of Factor A each represented by a single tank.  So I would like to use a split plot regression approach using factors A and B (similar to that recommended by Jones and Nachstheim) but how do I specify AND interpret the output in Proc Mixed?   So this is split plot design with a ordinal whole-plot factor.  Any experience with this?

Respected Advisor
Posts: 2,655

Re: proc mixed split plot regression

No experience, so I'm just spitballing here.  If the temperatures differ from tank to tank but are a known continuous value, could you just include temperature as a covariate in the model?

Maybe something like:

proc mixed;

class B;

model response=B temperature B*temperature;


Comparisons of means for the levels of B can be controlled by using the AT= option in an LSMEANS statement. If you suspect heterogeneity of response due to temperature, then this may get more complicated, but I would think this would be appropriate as a start.

Steve Denham

Ask a Question
Discussion stats
  • 1 reply
  • 2 in conversation