Tutorials › Causal Inference › Randomization Inference Under Interference

Interference and Spillover Effects · Tutorial 838 of 972

Randomization Inference Under Interference

Use the experimental assignment mechanism—not an assumed independent-treatment model—to construct valid randomization tests for direct and spillover effects.

Advanced 12 min read

What You'll Learn

  • Define sharp direct-effect, spillover-effect, and global null hypotheses using exposure-specific potential outcomes.
  • Explain why a component null may not permit unrestricted reassignment of treatments.
  • Construct conditional randomization distributions that hold the nuisance exposure dimension fixed.
  • Calculate assignment probabilities in a two-stage randomized saturation design.
  • Enumerate exact randomization distributions and interpret their p-values.
  • Identify when conditioning leaves too few assignments for a useful test.

Why Randomization Inference Requires Care Under Interference

In an ordinary randomized trial without interference, a sharp null hypothesis says that every unit’s outcome would be the same under every treatment assignment. That makes randomization inference direct: the observed outcomes can be treated as fixed while the assignments are redrawn according to the trial’s design. With interference, a person’s outcome may depend both on that person’s treatment and on other people’s treatments. A null hypothesis about only one of those dimensions may therefore leave outcomes unknown under some reassigned exposure conditions.

The key is to match the reassignment scheme to the null being tested. If testing a direct effect while allowing spillovers, reassignments must preserve the peer-exposure conditions relevant to each unit. If testing a spillover effect while allowing direct effects, they must preserve each unit’s own treatment. A global sharp null, by contrast, permits reassignment across both dimensions because it specifies no effect of either.

Definition: Let \(Z_i\) denote unit \(i\)’s own treatment and \(G_i\) its peer exposure under a stated exposure mapping. Write \(Y_i(z,g)\) for the potential outcome under own treatment \(z\) and peer exposure \(g\). A sharp null specifies enough potential outcomes to determine every outcome needed for the assignments included in the test.

Three null hypotheses illustrate the distinction. A sharp null of no direct effect states \(Y_i(1,g)=Y_i(0,g)\) for every unit \(i\) and relevant exposure \(g\), while allowing outcomes to vary across peer exposures. A sharp null of no spillover effect states \(Y_i(z,g)=Y_i(z,g')\) for every unit \(i\), own treatment \(z\), and relevant peer exposures \(g,g'\), while allowing outcomes to vary with own treatment. The global sharp null states \(Y_i(z,g)=Y_i(z',g')\) for all relevant combinations.

These are finite-population hypotheses about the potential outcomes of the units in the experiment. A randomization test evaluates how unusual the observed statistic would be if the specified null held, under the assignments the design could actually have produced. It does not require a sampling model for the outcomes.

Assignment Probabilities and Exact Tests

Let \(\Omega\) be the set of assignments allowed by the design, with assignment probability \(\Pr(A=a)\) for \(a\in\Omega\). For a statistic \(T(a)\), an upper-tail randomization p-value is the total design probability of assignments whose statistic is at least as large as the observed statistic. A lower-tail test reverses the inequality; a two-sided test needs a stated measure of extremeness, such as \(|T(a)|\) when the null value is zero.

$$ p_{\mathrm{upper}} = \sum_{a\in\Omega} \Pr(A=a)\, \mathbf{1}\{T(a)\geq T(a_{\mathrm{obs}})\}. $$

For a conditional test, restrict \(\Omega\) to assignments satisfying the chosen conditioning event \(C\), and renormalize their original design probabilities. The statistic and its outcomes must be imputable under the null for every assignment in this restricted set. Using equal weights is correct only when the assignments in the set are equally likely under the original design.

Theorem: Under a sharp null, if the observed assignment is drawn from the stated randomization distribution and the test statistic is computed for every assignment using outcomes imputed under that null, the resulting randomization test controls its type I error at the nominal level (with discreteness and ties making the usual inclusive-tail test potentially conservative). The same statement holds conditionally on an event of positive design probability, using the assignment distribution conditional on that event.

The justification is a rank argument. Under the null, the statistic that would be observed for each allowed assignment is known. The realized assignment is distributed according to the design probabilities, so its statistic has the very distribution used to calculate the tail probability. The probability of landing in the most extreme fraction of that distribution is no greater than that fraction; ties can only make the inclusive-tail rejection set larger in steps, so a nonrandomized test may reject less often than its nominal level.

Two-Stage Randomized Saturation Designs

In a two-stage saturation design, clusters are first assigned to saturation levels, and then units within each cluster are assigned treatment according to the cluster’s level. For example, the first stage may assign some clusters to a low probability \(q_L\) and others to a high probability \(q_H\). The second stage might assign each unit independently with its cluster’s probability, or assign a fixed number of treated units. The test must reproduce whichever second-stage rule the experiment used.

For a design that selects a fixed number of clusters for each saturation level and then treats each unit independently with probability \(q_{S_c}\) in cluster \(c\), the probability of a complete assignment is the first-stage probability multiplied by the conditional second-stage probabilities. If the first-stage allocation is uniform over permitted allocations, its probability is the reciprocal of their number.

$$ \Pr(A=a) = \Pr(S=s)\prod_c q_{s_c}^{\,\sum_{i\in c} z_i} (1-q_{s_c})^{\,n_c-\sum_{i\in c} z_i}, $$

Here \(S\) is the vector of cluster saturation assignments, \(s_c\) is cluster \(c\)’s assigned level, \(n_c\) is its size, and \(z_i\) is unit \(i\)’s treatment. For fixed-count assignment within clusters, replace the Bernoulli factor by the probability of the permitted treatment subset—for example, \(1/{n_c\choose m_c}\) when exactly \(m_c\) units are selected uniformly.

An exposure mapping must also be specified. In the examples below, \(G_i\) is the saturation level assigned to unit \(i\)’s cluster. This is a design-defined measure of peer exposure, not necessarily the realized fraction of that unit’s peers treated. If the analysis instead defines peer exposure from realized peer assignments, it must calculate that exposure anew for every reassignment.

Testing a Direct-Effect Null by Conditioning on Peer Exposure

Suppose the null is no direct effect at a fixed peer exposure, but spillovers are not ruled out. Under an unrestricted reassignment, a unit might move from one peer-exposure condition to another. Its observed outcome would then not reveal its outcome under that new peer exposure. A conditional test avoids this problem by retaining only assignments with the same \(G_i\) for every unit. Within that restricted set, the direct-effect null lets us impute each unit’s outcome under either own-treatment value at its fixed exposure.

Worked Example: A Conditional Test of a Direct Effect

Consider one hypothetical cluster of four units assigned saturation \(G_i=0.5\). The design treats exactly two units, chosen uniformly. The observed outcomes are 2, 4, 6, and 8 outcome units; units with outcomes 2 and 8 received treatment. Test the sharp null of no direct effect at this fixed saturation using the treated-minus-untreated mean as the statistic.

Under the null, each unit’s outcome at \(G_i=0.5\) is unchanged by its own treatment, so the four observed outcomes can be used for every permitted treatment pair. There are \({4\choose2}=6\) equally likely assignments. For the observed pair, the statistic is \(T_{\mathrm{obs}}=(2+8)/2-(4+6)/2=0\) outcome units.

The six possible treated pairs have statistics \(-4,-2,0,0,2,4\) outcome units. For a two-sided test using \(|T|\), every assignment is at least as extreme as the observed statistic, since \(|T_{\mathrm{obs}}|=0\). Thus \(p=6/6=1.000\). These data provide no evidence against this sharp direct-effect null. The test is conditional on the single saturation level; it does not test whether outcomes vary across saturation levels.

This example also shows why “conditioning” is more than a computational choice. It defines the assignments over which the direct-effect null supplies imputable outcomes. If fixing every unit’s peer exposure leaves no treatment variation, the conditional test cannot assess the direct effect for that design and exposure condition.

Testing a Spillover-Effect Null by Conditioning on Own Treatment

For a sharp null of no spillover effect, outcomes are known across peer-exposure changes only when own treatment is held fixed. A conditional spillover test therefore retains assignments with the observed \(Z_i\) for every unit and varies the assigned peer exposures according to the design. The design probabilities of those assignments may not be equal; they must be inherited from and normalized under the original two-stage assignment mechanism.

Worked Example: A Conditional Test of a Spillover Effect

Suppose four hypothetical clusters each contain two units. Exactly two clusters receive high saturation \(q_H=0.75\), and the other two receive low saturation \(q_L=0.25\). Within clusters, treatment is assigned independently with the stated probability. In the observed trial, all eight units are untreated. The four cluster mean outcomes are 1, 2, 8, and 9 outcome units. Test the sharp spillover null among these untreated units using the high-saturation minus low-saturation difference in cluster means.

Condition on all eight own-treatment values being zero. There are \({4\choose2}=6\) ways to choose the two high-saturation clusters. Each such assignment has the same conditional probability: two clusters have probability \((1-0.75)^2\) of both units being untreated and two have probability \((1-0.25)^2\). Their product is the same whichever clusters receive the high level. Therefore each of the six assignments has conditional probability \(1/6\).

Under the sharp spillover null, because every unit remains untreated, each cluster’s observed mean outcome stays fixed across these six assignments. The possible sums of the two high-saturation cluster means are 3, 9, 10, 10, 11, and 17. The total of all four means is 20, so for a high-cluster sum \(s\), the difference in means is \(s/2-(20-s)/2=s-10\). The six statistics are consequently \(-7,-1,0,0,1,7\) outcome units.

The observed high-saturation clusters have means 1 and 2, giving \(T_{\mathrm{obs}}=(1+2)/2-(8+9)/2=-7\) outcome units. The two-sided p-value based on \(|T|\) is \(2/6=0.3333\), since only the assignments with statistics \(-7\) and \(7\) are at least as extreme. This is an exact conditional test for the specified spillover null among units whose own treatment is zero, not a test of spillovers at other own-treatment values.

Using the Full Two-Stage Design Under a Global Sharp Null

When the null states that outcomes do not depend on own treatment or peer exposure, the observed outcome for each unit is known under every assignment. The test can then use the full two-stage assignment distribution without conditioning on \(G_i\) or \(Z_i\). The next example calculates an exact tail probability using the stage-one saturation lottery and stage-two Bernoulli assignments.

Worked Example: Enumerating a Full-Design Randomization Distribution

Consider two hypothetical clusters, each with two units. One cluster is assigned \(q_H=0.75\) and the other \(q_L=0.25\), with the two possible allocations equally likely. Within each cluster, the two units are independently treated with its assigned probability. Under the global sharp null, the four fixed outcomes are 1 and 2 in cluster A and 4 and 8 in cluster B. Use \(T=\sum_i Z_iY_i\), the sum of outcomes among treated units, for an upper-tail test.

Suppose A is assigned high saturation and both A units are treated, while B is assigned low saturation and neither B unit is treated. This observed assignment has probability \((1/2)(0.75^2)(0.75^2)=0.1582\), and its statistic is \(T_{\mathrm{obs}}=1+2=3\) outcome units. The observed assignment probability is not itself the p-value; the p-value sums probabilities over all assignments with \(T\geq3\).

First consider A high and B low. In A, the treated-outcome sum is 0, 1, 2, or 3 with probabilities 0.0625, 0.1875, 0.1875, and 0.5625. In B, the corresponding sum is 0, 4, 8, or 12 with probabilities 0.5625, 0.1875, 0.1875, and 0.0625. The only way for the total to be below 3 is for B’s sum to be 0 and A’s sum to be at most 2. Thus \(P(T<3\mid A\text{ high})=0.5625(0.0625+0.1875+0.1875)=0.2461\), and \(P(T\geq3\mid A\text{ high})=0.7539\).

If B is high and A low, the only way for \(T<3\) is for B’s sum to be zero and A’s sum to be at most 2. These probabilities give \(P(T<3\mid B\text{ high})=0.0625(0.5625+0.1875+0.1875)=0.05859\), so \(P(T\geq3\mid B\text{ high})=0.9414\). Averaging over the equally likely first-stage assignments gives \(p=(0.7539+0.9414)/2=0.8477\). Under the global sharp null, the observed statistic is not unusual in the upper tail of the full design distribution.

What Makes a Randomization Test Valid—and What Can Go Wrong

The essential validity conditions are concrete. The assignment mechanism must be known and represented correctly, including fixed totals, cluster restrictions, and unequal probabilities. The exposure mapping and potential-outcome null must be stated. Finally, every outcome needed to compute the statistic under every included assignment must be determined by the null and observed data. Under a component null, that last condition is often met only after conditioning on the other exposure dimension.

Common mistake: “The direct-effect null says treatment has no effect, so I can freely reassign treatment and keep all observed outcomes fixed.” That reasoning silently assumes outcomes are also invariant to peer exposure. Under interference, a unit whose peer exposure changes may have a different outcome even if its own-treatment effect is zero. The correction is to condition on the relevant peer exposures, or to test a broader sharp null that also rules out the exposure changes’ effects. For a spillover test that allows direct effects, the corresponding safeguard is to condition on own treatment.

Conditioning can reduce the number of assignments substantially. The resulting test can be exact but uninformative if few assignments remain, or impossible if only the observed assignment satisfies the condition. A small p-value is evidence against the particular sharp null and test statistic, not proof of a specific alternative mechanism. Likewise, failure to reject does not establish that direct or spillover effects are absent.

Key takeaway: Randomization inference under interference is design-based, but it is not automatically valid for every effect-specific null. Use the full assignment distribution only when the null imputes outcomes across every reassignment; otherwise, condition on the exposure dimension that the null does not eliminate and use the design’s conditional assignment probabilities.

Check Your Understanding

For each question, identify the null and assignment set carefully before calculating a statistic.

  1. Why does a sharp null of no direct effect alone not generally justify reassigning peer exposure while holding observed outcomes fixed?
  2. For a spillover test that allows direct effects, what should be held fixed across the assignments in the conditional randomization distribution?
  3. In a two-stage design, why might assignments within a conditional set require unequal weights?
  4. If conditioning on every unit’s peer exposure leaves only one permitted assignment, what does that imply about the usefulness of the conditional direct-effect test?
  5. What additional assumption allows a global sharp-null test to use assignments that change both own treatment and peer exposure?
Clinical Biostats