8
$\begingroup$

Context:

My question concerns a typical design in my area – a researcher takes a group of subjects (say 10) and then applies three different conditions to them to measure the change in a response variable, e.g. vertical jump height performed after drinking a glucose drink, coloured plain water, and fruit juice (say). Every subject has every treatment, but in a random order with enough time between for effects to ‘wash out’.

Analysis:

Kuehl (2000) (Kuehl, R. O. (2009) Design of Experiments: Statistical principles of research design and analysis, Duxbury Press, CA, p497 2nd Ed.) states:

When each of the treatments is administered in a random order to each subject... then subjects are random blocks in a randomised complete block design”

and then shows the corresponding analysis.

In this case, the subject is a random effect, but a nuisance or blocking factor, and although our statistical model will test the significance of the block factor, we aren’t really interested in its significance. However, many researchers (and reviewers!) think that such a design should be analysed as a repeated measures design with a Mauchly test for the Huynh-Feldt condition (with the treatment as the repeated measure). However, this seems more appropriate when a time factor is being analysed – for example when observations are taken at 0 minutes, 10 minutes, 30 minutes and 60 minutes, for example. In this case the covariance between pairs of time points might reasonably be expected to change, particularly when unequal time intervals are used. [In fact, I use SAS to model different covariance structures in this case (e.g. autoregressive) and use the AIC to choose the best structure, though this is not an approach that is well received by many reviewers.]

I understood that when the subject is a block factor, and the different treatments are administered in a random order that is different for different subjects, this means that the correlation between observations is different for each subject so compound symmetry can be assumed.

Question:

  • How should repeated measures ANOVAs with 3 or more conditions presented in random order be analysed?
  • Is it reasonable to assume compound symmetry?
$\endgroup$
  • 5
    $\begingroup$ Could you specify what your question is? Is this about the relevance or validity of using repeated-measures ANOVA with appropriate test for the variance-covariance structure, alternative ways for analyzing RCBD, or so? $\endgroup$ – chl May 5 '11 at 12:32
  • $\begingroup$ I liked where you were going with this question. After 12 days with no response from you, I have updated your question, so that there is an actual question. Feel free to modify if I have misconstrued. $\endgroup$ – Jeromy Anglim May 17 '11 at 6:40
  • $\begingroup$ @chi @Jeromy Anglim - sorry I hadn't worked out how to use the site properly - I thought I'd get an email when someone replied. I suppose my question is, is this really a repeat measures design or is it a randomised block design. In a repeat measures design, if I measure a variable at 0 minutes, 10 minutes and 20 minutes, then the order of measurements and the time gap is same for all subjects. But if I randomise the order of presentation of three treatments, then Subject 1 may get T1 then T3 then T2, whereas Subject 2 may get T1 then T2 then T3, so the time gap between T1 and T3 is different $\endgroup$ – user5178 Jun 26 '11 at 16:35
  • $\begingroup$ @Sam I converted your reply as a comment, which is the way to do when you want to add information or reply to request for clarification (you can also update your question directly). (You'll get email notification if you ask for that; also, when you're logged, your SE message box should indicate whether you have new responses or comments.) $\endgroup$ – chl Jun 26 '11 at 20:53
  • $\begingroup$ I'm no expert in statistics, so a comment it is. Your reference is the first time I've come across the idea of not using an RM ANOVA for said design. Most textbooks I've read usually call it Randomized Blocks when you conduct a between-subjects study but match different subjects on a variable (e.g. athletic ability) to allow you to treat the values of multiple subjects as if they came from one subject. This leverages the power of an RM design when it is not appropriate to actually perform one (e.g. learning effects). But they still advocate analysing the results using an RM ANOVA in both cases $\endgroup$ – ThomasH Nov 28 '12 at 17:26
1
$\begingroup$

Repeated measures is kind of an overloaded term. To some people it refers to a particular statistical analysis method; to others it refers to the structure of the design.

This is a variant on a three period, three treatment crossover design.

It is a variant because usually in a crossover design you randomize subjects to sequences. In this case the sequence is determined randomly for each subject. Since there are six possible sequences, it might be that some sequences are not observed, especially with 10 subjects. Maybe this is formally the same as randomizing subjects to sequences, but I haven't looked at that yet.

The considerations for crossover designs are:

  • Carryover effects: Also known as residual effects, where prior treatment may affect response to current treatment. The goal of the washout periods is to remove this from consideration. You could also have (in theory) second-order residual effects, where the treatment given in the first period potentially affects the response to treatment given in the third period.

  • Period effects: Response to treatment(s) may change as the study goes on for a given subject.

  • Autocorrelation: Serial correlation in errors is usually an issue with more closely measured data. In simple balanced designs, having a random effect for subject is going to imply equal correlaation of errors from each subject.

  • Subject effects: Subjects may differ in mean response from each other regardless of treatments. You could conceive of a situation where measurement error was serially correlated separate from a random subject effect.

  • Sequence effect: In cases where you randomize subjects to sequences, subjects are considered nested in sequence.

A minimal analysis for this would be the suggested randomized complete block design. That is, a fixed effect for treatment and a random effect for subject. With a skimpy sample size that might be all you can really do.

I would argue for a bit more structure to the analysis, if possible. Assuming no carryover effects on scientific grounds, it seems like a good idea to have at fixed effects for treatment, period, and treatment $\times$ period interaction, and a random effect for subjects. For small data sets, if this model can't be fit, I would drop the treatment $\times$ period interaction first.

Period should be included because it represents a restriction on the randomization. You cannot "randomize" periods --- they always happen in the same order. Treatment $\times$ period interaction might be indicative of some sort of carryover effect.

With tons of data one could work up terms that would allow estimation of various specific carryover effects. My notes on this are gone, though I know I've seen it handled in some texts.

The strategy of additionally modelling the correlation structure on the R-side seems reasonable to me. That allows one to claim that one is handling the possible dependence structure induced by repeated measures on the same subject, which I would also probably claim about the random effect for subject if the analysis devolved to that level... It is also nice if various analysis strategies provide broadly or very similar results.

For implementation, I'd use PROC MIXED in SAS and likely nlme or lme4 in R.

I'll punt on the compound symmetry question, since that seems more like a holdover from the days where MANOVA was the only "correct" analysis for repeated measures.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.