AES Consulting meeting on 06 Jun 2016
Split-plot-in-time
Considering a 2 block experiment with 3 experimental units per block each with a different treatment. Measurements are taken on each unit for each of 3 years, but the treatment is the same in each unit from year to year. Thus there are a total of (2x3x3=) 18 observations.
Full design
A full design would include these effects (degrees of freedom):
- (Intercept) (1)
- Treatment (2)
- Block (1)
- Year (2)
- Block x Treatment (2)
- Year x Treatment (4)
- Block x Year (2)
- Block x Treatment x Year (4)
This model has no error degrees of freedom left.
Design based on an experimental unit
Suppose we take averages for each experimental unit across the years. Then a design could be
- Treatment (2)
- Block (1)
- Year (2)
- Block x Treatment (2)
This still has no error degrees of freedom.
Split-plot-in-time design
Taking the design based on an experimental unit but expanding it to include year produces this design
- (Intercept) (1)
- Treatment (2)
- Block (1)
- Block x Treatment (2) (random effect for repeated measure)
- Year (2)
- Year x Treatment (4)
This design has 6 degrees of freedom which arise by eliminating
- Block x Year (2)
- Block x Treatment x Year (4)
Correlation structure
Treating the block x treatment interaction as a standard random effect will result in a compound symmetry structure for observations across years.
In SAS, specify this by using repeated / type = cs subject = Block x Treatment CS
.
This means that the correlation between every combination of years is the same which may not be reasonable since some of the years are more spread apart in time than others.
An alternative is to use an autoregressive structure, e.g. AR(1), which can be specified using type=ar(1)
.