Why Stratum Sizes Matter
When randomization is conducted separately within strata, each stratum receives a specified number of treatment assignments. If stratum sizes differ, the treatment count and the chance of treatment for an individual can differ from one stratum to another. The total number assigned to treatment may be fixed even though the overall treatment proportion is not exactly the same in every stratum.
The key is to distinguish the stratum's size from its treatment count. The randomization rule chooses a fixed-size subset of each stratum for treatment; it does not choose one fixed-size subset from the trial as a whole. We will calculate the resulting allocation probabilities and show how unequal stratum sizes affect them.
The design is conditional on the stratum rosters and the counts \(m_h\). It assumes that each participant belongs to exactly one stratum, the specified count is feasible in every stratum, and every subset of the required size within a stratum has the same chance of selection. Independence across strata means that a choice in one stratum does not alter the randomization probabilities in another. These are assumptions about the assignment mechanism; no outcome model is needed to calculate its probabilities.
Counting Allocations and Their Probabilities
There are \(\binom{n_h}{m_h}\) possible treatment subsets in stratum \(h\). Since the choices are independent across strata, the total number of complete allocations consistent with the stratum counts is the product of these numbers. Each such allocation has the same probability.
An allocation that violates even one required stratum count has probability zero. The total number assigned to treatment, \(M\), is fixed, but the stratum-specific counts are also fixed, so treatment counts cannot shift from one stratum to another under this design.
Worked Example: Counting Allocations in Three Unequal Strata
Suppose a hypothetical trial has three strata with sizes \(n_1=8\), \(n_2=5\), and \(n_3=3\). The design assigns \(m_1=4\), \(m_2=2\), and \(m_3=1\) participants to treatment, respectively. Find the total treatment count and the probability of one specified allowed allocation.
The total sample size is \(N=8+5+3=16\), and the total treatment count is \(M=4+2+1=7\). The numbers of treatment subsets in the strata are \(\binom{8}{4}=70\), \(\binom{5}{2}=10\), and \(\binom{3}{1}=3\). Therefore, the number of allowed complete allocations is
Because the randomization is uniform within each stratum and independent across strata, any one particular allocation with exactly 4, 2, and 1 treated participants in the respective strata has probability
The treatment proportions within the three strata are \(4/8=0.5\), \(2/5=0.4\), and \(1/3\approx0.3333\). The overall treatment proportion is \(7/16=0.4375\). These are different quantities: the overall proportion is a size-weighted result of the stratum-specific counts, not a treatment probability that applies identically to every participant.
Individual and Joint Assignment Probabilities
Let \(A_{hi}=1\) indicate that participant \(i\) in stratum \(h\) is assigned to treatment, and let \(A_{hi}=0\) indicate control. Within a stratum, uniform selection of exactly \(m_h\) participants gives each participant the same marginal treatment probability, \(m_h/n_h\). To see this by counting, when \(1\leq m_h\leq n_h\), the number of allowed treatment subsets containing a specified participant is \(\binom{n_h-1}{m_h-1}\). Dividing by the total number of subsets gives the probability. The boundary case \(m_h=0\) must be handled separately: treatment is impossible, so the probability is \(0=m_h/n_h\).
Thus, the probability is determined by the treatment fraction in that participant's own stratum. If two strata have unequal treatment fractions, participants in those strata have unequal treatment probabilities. The same expression also covers \(m_h=n_h\), where every participant is assigned to treatment and the probability is 1.
For two distinct participants in the same stratum, both can be treated only if \(m_h\geq2\). When \(2\leq m_h\leq n_h\), counting subsets that contain both participants gives the joint probability below. If \(m_h=0\) or \(m_h=1\), the event is impossible and has probability zero; the expression on the right is also zero. This separate treatment of the boundary cases avoids using a binomial coefficient with a negative lower argument.
For \(n_h\geq2\) and \(m_h\in\{0,1\}\), the joint probability is \(0\), which equals \(m_h(m_h-1)/[n_h(n_h-1)]\). The formula on the right can therefore be used for any \(m_h\) when \(n_h\geq2\), with the counting derivation applying only when \(m_h\geq2\). There are no two distinct participants in a stratum of size one, so a same-stratum pair probability is not defined there.
For participants in different strata, independence of the stratum selections gives a product of their individual probabilities. This differs from the same-stratum case, where the fixed treatment count makes assignments dependent.
Worked Example: Two Participants in the Same Stratum
In a hypothetical stratum of \(n_h=6\) participants, exactly \(m_h=2\) are assigned to treatment. Calculate the probability that a specified participant is treated and the probability that two specified distinct participants are both treated.
There are \(\binom{6}{2}=15\) possible treatment subsets. For the specified participant to be treated, the other treated participant can be any one of the remaining five. Therefore,
For both specified participants to be treated, the remaining four participants cannot be selected; there is exactly one treatment subset containing the specified pair. Equivalently, the pair formula gives
As a check, once the first specified participant is known to be treated, one of the five remaining participants must receive the second treatment assignment. Thus the conditional probability that the second specified participant is treated is \(1/5\), and \((1/3)(1/5)=1/15\), matching the joint calculation.
Worked Example: Joint Probability Across Two Strata
Suppose a hypothetical trial has two strata. In the first, \(n_1=4\) and \(m_1=2\); in the second, \(n_2=3\) and \(m_2=1\). Find the probability that a specified participant in each stratum is assigned to treatment. Also find the probability of one particular complete allocation.
The individual treatment probabilities are \(2/4=1/2\) in stratum 1 and \(1/3\) in stratum 2. Since the selections in the two strata are independent, the probability that both specified participants are treated is
There are \(\binom{4}{2}=6\) possible treatment subsets in the first stratum and \(\binom{3}{1}=3\) in the second. Hence there are \(6\times3=18\) allowed complete allocations, each with probability \(1/18\approx0.05556\). The joint event just calculated includes more than one complete allocation, which is why its probability is larger than \(1/18\).
Justification and Dependence Within a Stratum
The individual probability follows from a simple counting identity. For \(1\leq m_h\leq n_h\), the number of subsets containing a specified participant is the number of ways to choose the other \(m_h-1\) treated participants from the remaining \(n_h-1\). The ratio of this count to the total number of subsets simplifies as follows:
The boundary \(m_h=0\) is not covered by that counting ratio, because it would involve \(\binom{n_h-1}{-1}\). Instead, no participant can be treated when the fixed treatment count is zero, so the probability is zero, as stated above. The same care is needed for the pair calculation: its counting ratio is valid for \(2\leq m_h\leq n_h\), while the cases \(m_h=0\) and \(m_h=1\) have probability zero directly.
Within a stratum, assignments are generally not independent. If one participant is treated, fewer treatment assignments remain for the others. For example, in the six-person stratum above, the probability that a second specified participant is treated given that the first is treated is \(1/5\), less than the unconditional probability \(1/3\). Across strata, by contrast, the design specifies independent selections, so the corresponding probabilities multiply.
Comparing the Design With Complete Randomization Overall
Under complete randomization over the entire trial, a fixed total of \(M\) participants would be selected uniformly from all \(\binom{N}{M}\) possible subsets. Stratified complete randomization restricts the allowed subsets further by fixing each stratum's count at \(m_h\). The number of allowed subsets is \(\prod_h\binom{n_h}{m_h}\), which is no larger than \(\binom{N}{M}\). The restricted design fixes the treatment count within every stratum; complete randomization overall fixes only the total.
In the two-stratum example, stratified randomization fixes counts at 2 of 4 and 1 of 3, with 18 allowed allocations. Complete randomization over all 7 participants with 3 assigned to treatment would allow \(\binom{7}{3}=35\) allocations. Under that alternative, the treatment count in either stratum could vary, and each participant's treatment probability would be \(3/7\), rather than the stratum-specific values \(1/2\) and \(1/3\).
The choice between these designs depends on the intended allocation rule. Stratification can ensure specified counts in important groups, while complete randomization overall allows those group counts to vary by chance. Neither comparison changes the arithmetic: the probabilities must be calculated from the actual set of allocations the design permits.
Diagnostics, Limitations, and Common Mistakes
Before interpreting assignment probabilities, check that the design is internally consistent. Each stratum count must satisfy \(0\leq m_h\leq n_h\); the stratum sizes must sum to the stated total \(N\); and the treatment counts must sum to the stated total \(M\). If a proposed allocation has a different count in any stratum, it is not an allowed allocation under this design. If a schedule is presented, count its assignments separately within each stratum rather than checking only the overall total.
A common mistake is to use \(M/N\) as the treatment probability for every participant. That is generally wrong: the probability for someone in stratum \(h\) is \(m_h/n_h\). The overall fraction \(M/N\) is a weighted average of those stratum-specific fractions. For example, the first worked example has \(M/N=7/16\), while the stratum-specific probabilities are \(1/2\), \(2/5\), and \(1/3\).
Another mistake is to multiply probabilities for two participants in the same stratum as if their assignments were independent. The fixed-count constraint creates dependence; use the pair formula or condition on the first assignment. Multiplication is justified for assignments in different strata because independence across strata is part of the design definition.
This design fixes stratum counts, but it does not guarantee that all participant characteristics are balanced within a stratum, nor does it make treatment and control groups identical on every feature. Its probabilities also rely on the stated strata and treatment counts: changing either changes the randomization distribution. In practice, unequal stratum sizes can make the stratum-specific treatment fractions differ when counts must be whole numbers, so the planned counts should be checked against the intended allocation proportions.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- A hypothetical stratum has 10 participants and 4 treatment assignments. How many treatment subsets are possible, and what is one participant's treatment probability?
- In a stratum of 5 participants, exactly 1 is treated. What is the probability that two specified distinct participants are both treated? Explain without using a binomial coefficient with a negative lower argument.
- Two independent strata have sizes 6 and 4, with treatment counts 3 and 1. What is the probability that one specified participant in each stratum is treated?
- Why is \(M/N\) not necessarily the treatment probability for each participant?
Solutions. (1) There are \(\binom{10}{4}=210\) subsets. Each participant has probability \(4/10=0.4\) of treatment. (2) The event is impossible: only one participant in the stratum is treated, so the joint probability is 0. (3) Independence across strata gives \((3/6)(1/4)=1/8=0.125\). (4) \(M/N\) is the overall treatment fraction, while a participant in stratum \(h\) has treatment probability \(m_h/n_h\). Those fractions can differ across strata when treatment counts or stratum sizes differ.