Gilmore (2007) argues that the use of ANCOVA in Anstey et al. (2006) is invalid. Based on quotes in Winer (1971), he claims that covariates cannot be used to adjust interaction terms involving within-subject variables. Because the within-subject variable in Anstey et al. (2006) is time, he concludes that the results involving Time × Group interaction coefficient are invalid and consequently all analyses need to be re-done.
However, our approach was more cautious than that proposed by Winer (1971) and we therefore argue that our analysis published in this journal is acceptable. We agree that with constant covariates, the classical ANCOVA of Winer (1971) can only be used to adjust the between-subjects (BS) effects. To see this, in a simplest case with time (T) as within-subject factor and treatment group (G) as between-subject factor and there are only two levels for T and G, we can represent Winer's ANCOVA model in notational form as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170221052916903-0603:S135561770707049X:S135561770707049Xfrm001.gif?pub-status=live)
where yij refers to outcome for individual i measured at time Tj (0=pre-test,1=post-test), Gi is the treatment received by individual i (0= control, 1 = treatment), μ is the grand mean and xi is the covariate value of individual i and it is constant across time. We can see that since the G and x share the same subscript, the covariates only adjusts the group effects (BS effects) but not the other effects.
However, another class of ANCOVA can be created by having the covariates in the WS portion of the above model. This is achieved by adding interaction terms between the covariates and the WS variable (time). This is the default ANCOVA produced by SPSS (SPSS Resolution No. 22133). This model can now be written as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170221052916903-0603:S135561770707049X:S135561770707049Xfrm002.gif?pub-status=live)
It can be seen now that the new interaction term (x * T) does not adjust the main within subject effect T (time) but nonetheless it adjusts the Time × Group interaction effect. A significant λ can be interpreted as the change of scores between pre- and post-test for an individual is affected also by his/her covariates value x, apart from the treatment group he/she is allocated to. More importantly, this model is better than Winer's model because a significant λ will indicate that the homogeneity of slopes (HOS) assumption implicit in Winer's ANCOVA in Eq. 1 is violated (see also SPSS Resolution No. 22133). On the other hand, if λ is not significant, our approach and Winer's approach will yield virtually the same numerical solutions.
Hence, we conclude that our approach is different from Winer's ANCOVA and is a more cautious one and hence no re-analysis is needed.