Tutorials › Bioavailability and Bioequivalence › Sequence, Period, and Treatment Effects in Crossover ANOVA

Log-Scale Analysis and ANOVA · Tutorial 181 of 399

Sequence, Period, and Treatment Effects in Crossover ANOVA

Learn to distinguish the between-subject sequence comparison from the within-subject treatment contrast and the period comparison in a two-period crossover ANOVA.

Intermediate 10 min read

What You'll Learn

  • Interpret sequence, period, and treatment effects in a two-period crossover.
  • Derive the treatment contrast from sequence-specific within-subject differences.
  • Separate treatment and period contributions using balanced sequence means.
  • Translate a log-scale treatment contrast into an exposure ratio.
  • Recognize misleading comparisons and limits of the sequence effect.

Three Terms, Three Comparisons

In a two-period crossover, each participant receives both treatments, but the order differs by randomized sequence. For example, sequence TR means test treatment in period 1 and reference treatment in period 2; sequence RT means reference in period 1 and test in period 2. The ANOVA terms for sequence, period, and treatment therefore do not represent interchangeable versions of one treatment comparison. They refer to different features of the design.

This distinction matters especially when the endpoint is a pharmacokinetic measure analyzed on the log scale. A treatment contrast on that scale represents a difference in mean log exposure and, after exponentiation, a ratio of geometric means. A sequence contrast instead compares groups of participants assigned different orders. A period contrast compares measurements collected at different study times. Neither is, by itself, the treatment effect.

Definition: In a two-period, two-treatment crossover, sequence denotes the assigned order of treatments (TR or RT); period denotes when a measurement was collected (period 1 or period 2); and treatment denotes which treatment produced that measurement (T or R). The treatment effect is the contrast between T and R, adjusted for period and sequence under the crossover model.

Model and Meaning of the Terms

Let \(y_{ijk}\) be the log-transformed endpoint for participant \(j\) in sequence \(i\) and period \(k\). A conventional complete-data crossover ANOVA includes sequence, participant within sequence, period, treatment, and residual variation. To describe its fixed-effect means, write the expected log endpoint as

$$ E(y_{ijk})=\mu+\alpha_i+\pi_k+\tau_l $$

Here \(i\) identifies sequence, \(k\) identifies period, and \(l\) identifies the treatment received in that period. The terms \(\alpha_i\), \(\pi_k\), and \(\tau_l\) are the sequence, period, and treatment effects. A sum-to-zero convention, such as \(\alpha_{TR}+\alpha_{RT}=0\), \(\pi_1+\pi_2=0\), and \(\tau_T+\tau_R=0\), makes the parameters identifiable alongside the overall mean \(\mu\). The treatment contrast of interest is \(\delta=\tau_T-\tau_R\); the period contrast is \(\Delta_P=\pi_2-\pi_1\).

This mean structure is appropriate when treatment and period effects are additive: the treatment contrast is not assumed to change between periods. Its usual interpretation also relies on a properly randomized sequence assignment, complete measurements for the standard two-period analysis, independent participants, and no residual effect of a treatment from one period that changes the outcome in the next period. For statistical inference, the model’s error assumptions must also be reasonable: the relevant errors are centered around zero, have a suitable common variance, and are independent across participants. A log transformation can help make variability more stable, but it does not guarantee these assumptions.

The participant term is important conceptually: the two observations from one person are paired, not two independent observations. The sequence comparison is between groups of participants, whereas the treatment contrast uses the within-person treatment differences. The earlier material on Expected Mean Squares in the Two-Period Crossover ANOVA establishes why these comparisons belong to different error strata. Interpreting an effect estimate does not itself determine its test or error term.

Deriving the Treatment Contrast

The four sequence-by-period means make the distinction concrete. Denote by \(\bar y_{TR,1}\) and \(\bar y_{TR,2}\) the mean log outcomes in periods 1 and 2 for sequence TR, with corresponding notation for RT. Under the additive model, the mean difference oriented as T minus R is

$$ D_{TR}=\bar y_{TR,1}-\bar y_{TR,2}=\delta-\Delta_P $$

In sequence TR, T is measured in period 1 and R in period 2, so the treatment difference includes the period-1 versus period-2 difference. In sequence RT, orient the difference the same way—T minus R—even though T occurs second:

$$ D_{RT}=\bar y_{RT,2}-\bar y_{RT,1}=\delta+\Delta_P $$

Adding these two equations cancels the period contribution. Subtracting them cancels the treatment contribution. Thus, in a balanced design,

$$ \hat{\delta}=\frac{D_{TR}+D_{RT}}{2}, \qquad \widehat{\Delta}_P=\frac{D_{RT}-D_{TR}}{2} $$

This is a useful justification for the crossover treatment contrast: each sequence supplies a within-person comparison, and averaging the oppositely ordered comparisons removes an additive period effect. It also explains why a treatment comparison based only on one sequence, or on one period, need not isolate treatment.

Worked Example: Separating the Three Effects

Suppose a hypothetical balanced study has four participants in each sequence and the following mean log exposure values. The values are in natural-log units; the table gives group means, not individual observations.

SequencePeriod 1Period 2
TR10.20 (T)10.00 (R)
RT9.90 (R)10.10 (T)

For TR, the T-minus-R difference is \(D_{TR}=10.20-10.00=0.20\). For RT, it is \(D_{RT}=10.10-9.90=0.20\). Therefore, the balanced treatment contrast is

$$ \hat{\delta}=\frac{0.20+0.20}{2}=0.20 $$

The period contrast is \(\widehat{\Delta}_P=(0.20-0.20)/2=0.00\) log units. This agrees with the pooled period means: period 1 is \((10.20+9.90)/2=10.05\), and period 2 is \((10.00+10.10)/2=10.05\).

The treatment means pooled across periods and sequences are T: \((10.20+10.10)/2=10.15\), and R: \((10.00+9.90)/2=9.95\). Their difference is \(10.15-9.95=0.20\) log units. Exponentiating gives \(\exp(0.20)=1.221\), so the estimated geometric-mean exposure ratio T/R is about 1.22. This is an estimate, not a significance conclusion: the table does not provide the participant-level variability needed for a standard error or confidence interval.

The sequence means, averaged across periods, are 10.10 for TR and 10.00 for RT, a difference of 0.10 log units. That is a comparison between the two assigned-order groups. It is not the T-versus-R contrast of 0.20.

Interpreting Each Effect Carefully

Treatment. The treatment term estimates the adjusted mean difference between log outcomes under T and R. In a balanced two-period design, it can be understood through the average of the two sequence-specific T-minus-R differences. When the endpoint is log exposure, exponentiating the estimated contrast gives the ratio of geometric means. For example, a contrast of 0.20 log units corresponds to a ratio of \(\exp(0.20)\), not an increase of 0.20 in the original exposure units.

Period. The period term describes a systematic difference between period 2 and period 1, averaging over the treatment allocation under the model. It may reflect study-time changes that affect outcomes independently of the treatment contrast. The term is not a treatment effect and does not, by itself, identify why the measurements differ by period.

Sequence. The sequence term compares participants assigned to TR with participants assigned to RT. It is a between-subject comparison, even though every participant receives both treatments. It can reflect differences between the randomized groups as well as design-related patterns represented by the model. It does not directly estimate the within-subject T-versus-R contrast. In particular, a sequence difference alone does not establish that carryover occurred, and a nonsignificant sequence result does not prove that carryover is absent.

Worked Example: A Period Difference Does Not Equal a Treatment Difference

Consider hypothetical balanced sequence means in log units: TR has 10.20 in period 1 under T and 10.00 in period 2 under R; RT has 9.80 in period 1 under R and 10.20 in period 2 under T. The TR treatment-oriented difference is \(D_{TR}=10.20-10.00=0.20\). The RT difference is \(D_{RT}=10.20-9.80=0.40\).

The average treatment contrast is \(\hat{\delta}=(0.20+0.40)/2=0.30\) log units. The period contrast is \(\widehat{\Delta}_P=(0.40-0.20)/2=0.10\) log units. Directly, the period 1 mean is \((10.20+9.80)/2=10.00\), and the period 2 mean is \((10.00+10.20)/2=10.10\), also a difference of 0.10.

Pooling by treatment gives T mean \((10.20+10.20)/2=10.20\) and R mean \((10.00+9.80)/2=9.90\), reproducing the 0.30-log-unit treatment contrast. Its exposure ratio is \(\exp(0.30)=1.350\). The period difference is real as a model contrast in this hypothetical table, but it has not replaced or changed the definition of the adjusted treatment contrast.

Assumptions, Diagnostics, and Alternatives

The additive mean structure is a central condition for the simple decomposition above. If treatment behaves differently by period, the two sequence-specific treatment differences need not combine into a single common treatment contrast in the intended way. The complete, balanced example also relies on comparable numbers of participants and observations in the two sequences. With missing periods or unequal sequence sizes, the simple averages shown here may not equal the model-based estimates; the chosen sums-of-squares or model approach then matters.

Useful checks include reviewing the sequence-by-period cell means, plotting participant-level paired outcomes or differences, and inspecting residuals from the fitted model. Large discrepancies between the two sequence-specific treatment contrasts can signal a possible treatment-by-period pattern, influential observations, or chance variation. Marked changes in spread across groups, skewed residuals, or a few extreme log outcomes may make the usual error model questionable. These checks flag concerns; they do not by themselves prove which assumption failed.

The main alternative summaries answer different questions. Comparing the two treatment means after pooling all observations is straightforward in a balanced design, but it should be understood as the model-adjusted treatment contrast only when the design and model justify that pooling. Comparing T and R within one period uses different sequence groups and can mix treatment with sequence differences. Comparing sequence averages answers a between-group question, not a within-person treatment question. A parallel-group comparison, in turn, has no within-person treatment contrast because each participant receives only one treatment.

The classical two-period ANOVA is also limited by its assumptions and design. It does not make an additive model true, repair measurements affected by prior treatment, or supply information about participants who lack a required period. A sequence result should not be used as a shortcut for diagnosing carryover or deciding that the treatment comparison is valid. Conclusions about treatment should be based on the treatment contrast, its appropriate uncertainty, and the credibility of the model and study conduct.

Worked Example: A Misleading Single-Period Comparison

Return to the first hypothetical table. Comparing T in period 1 for sequence TR with R in period 1 for sequence RT gives \(10.20-9.90=0.30\) log units. It may be tempting to call 0.30 the treatment effect because the values are from the same period.

That comparison is between different people in different sequences. The sequence means averaged over periods were 10.10 for TR and 10.00 for RT, so the groups differ by 0.10 log units in this table. The single-period difference of 0.30 therefore does not isolate the within-subject treatment contrast, which the balanced calculation gave as 0.20. The appropriate correction is to orient both sequence-specific comparisons as T minus R and average them, or to use the corresponding fitted model contrast.

Common Pitfalls and Practical Use

  • Calling every ANOVA term a treatment effect. Sequence compares assigned-order groups, period compares study times, and treatment compares T with R.
  • Comparing treatments in just one period. This uses different sequence groups and can mix treatment and between-group differences.
  • Reversing the contrast orientation. In RT, T is in period 2; calculate T minus R as period 2 minus period 1 so both sequence contrasts have the same sign convention.
  • Interpreting a sequence p-value as a carryover test. A sequence comparison does not by itself identify carryover, and failure to detect a sequence difference is not evidence that carryover is impossible.
  • Reporting a log contrast as a raw-unit difference. A log contrast is in log units; exponentiate it to report a geometric-mean ratio on the original exposure scale.

In practice, first state which term is being interpreted and what units its contrast uses. Then identify whether its comparison is within participants or between sequences. For the treatment term, report the log-scale estimate and, when relevant, its exponentiated ratio with uncertainty. Keep the period and sequence terms as distinct model comparisons rather than treating either as a substitute for the treatment result.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. In sequence RT, which subtraction gives the treatment-oriented difference T minus R?
  2. If \(D_{TR}=0.10\) and \(D_{RT}=0.30\) log units, calculate the treatment and period contrasts.
  3. What comparison does the sequence term represent?
  4. Why does a nonsignificant sequence result not prove that carryover is absent?

Solutions. (1) Use period 2 minus period 1, because T is given in period 2 and R in period 1. (2) The treatment contrast is \((0.10+0.30)/2=0.20\) log units; the period-2-minus-period-1 contrast is \((0.30-0.10)/2=0.10\) log units. (3) It compares participants assigned TR with participants assigned RT, so it is a between-sequence comparison. (4) A nonsignificant result means the data did not provide sufficient evidence against the tested sequence contrast at the chosen significance level. It does not show that the contrast is exactly zero or that carryover cannot be present.