Split-plot and other multi-stratum structures are widely used in factorial and response surface experiments and residual maximum likelihood (REML) and generalized least squares (GLS) estimation is seen as the state-of-the-art method of data analysis for nonorthogonal designs. We analyze data...