Observed Features and Unidentified Coordinates
A missing-data model can fit the observed data well and still leave important features of the complete-data distribution undetermined. The observed distribution fixes, for example, how often each response pattern occurs and the distribution of the outcomes actually recorded in each pattern. It does not, in general, fix the distribution of outcomes that were not recorded. A useful way to describe this gap is to choose sensitivity parameters: explicit coordinates for the unidentified features, with the observed-data features held fixed.
The word “coordinate” is deliberate. It means that a parameterization labels the possible complete-data laws consistent with the same observed-data law. It does not mean that the coordinate is estimated from the observed data, nor that every numerical value is permissible. Probability constraints, support restrictions, and any stated model assumptions restrict which values define a valid law.
This distinction builds on the earlier material on what the observed likelihood establishes. In population terms, the observed-data law determines the pattern probabilities and the distributions of \(Y_{\mathrm{obs}}\) within those patterns. The unobserved conditional distributions are additional ingredients required to specify a complete-data law. Any sensitivity analysis must state how those ingredients are parameterized and what values are considered plausible.
A Factorization That Separates What Is Fixed
For each pattern \(r\), let \(P_{R,Y_{\mathrm{obs}}}\) denote the observed-data probability measure for the pattern and its observed outcomes. Let \(g_r(dy_{\mathrm{mis}}\mid y_{\mathrm{obs}})\) be a conditional probability distribution for the missing outcomes given the observed outcomes in pattern \(r\). The full joint law is assembled by combining these pieces. This is a measure identity, so it applies whether observed outcomes are discrete, continuous, or a mixture of both.
Here \(A\) and \(B\) are measurable sets of possible observed and missing outcomes. The conditional distribution \(g_r\) is defined almost everywhere under the observed-data law for pattern \(r\); its values at observed-outcome values having probability zero do not affect the joint law. When densities exist, the same factorization can be written as a joint density \(f_{R,Y_{\mathrm{obs}}}(r,y_{\mathrm{obs}})\) multiplied by the conditional density of \(Y_{\mathrm{mis}}\) given \((R=r,Y_{\mathrm{obs}}=y_{\mathrm{obs}})\).
The observed-data law supplies the measure in this expression. A sensitivity model supplies \(g_r\), perhaps through one or more parameters. To be valid, the resulting conditional distributions must be normalized probability distributions and must respect the outcome’s support. If a pattern has probability zero, its conditional distribution is not learned from the data and its contribution to this observed-data factorization is zero; claims about that pattern require additional assumptions.
Worked Example: A Risk-Difference Coordinate
Worked Example: A Risk-Difference Coordinate
Suppose a binary outcome \(Y\) is recorded when \(R=1\), and is missing when \(R=0\). Imagine the observed-data law has \(P(R=1)=0.60\) and \(P(Y=1\mid R=1)=0.40\). These two quantities are treated as fixed. Define \(p_0=P(Y=1\mid R=0)\), and use the risk difference \(\delta=p_0-P(Y=1\mid R=1)\) as the sensitivity coordinate.
By definition, \(p_0=0.40+\delta\). Because a binary-outcome probability must lie between zero and one, the coordinate must satisfy
Thus the observed-data features do not select one value of \(\delta\), but probability validity rules out values outside this interval. The complete-data prevalence is obtained by averaging across the response patterns:
For a hypothetical sensitivity choice \(\delta=0.20\), the missing-group prevalence is \(0.40+0.20=0.60\), and the complete-data prevalence is \(0.40+0.40(0.20)=0.48\). The arithmetic is compatible with the fixed observed quantities: the response probability remains \(0.60\), and the respondent prevalence remains \(0.40\). Changing \(\delta\) changes only the assumed missing-group prevalence and the resulting complete-data prevalence.
This example also shows why a sensitivity parameter must come with a scale and interpretation. A risk difference is measured in probability points. Other coordinates may be more natural for different scientific questions, but they can have different constraints and different behavior near the boundaries of the parameter space.
Worked Example: An Odds-Ratio Coordinate
Worked Example: Re-expressing the Binary-Outcomes Assumption
Use the same fixed observed quantities, \(P(R=1)=0.60\) and \(p_1=P(Y=1\mid R=1)=0.40\). Instead of a risk difference, define the sensitivity coordinate \(\gamma\) as the odds of \(Y=1\) among those with missing outcomes divided by the odds among respondents:
For \(0<p_1<1\) and \(\gamma>0\), solve for \(p_0\). Multiplying through gives \(\gamma p_1(1-p_0)=p_0(1-p_1)\), so
If \(\gamma=2\), the missing group has twice the respondent odds of a positive outcome. Substituting the observed respondent prevalence gives
The implied complete-data prevalence is \(0.60(0.40)+0.40(0.5714)=0.4686\), or about 46.86%. The response rate and respondent prevalence have not changed. For interior \(p_1\), every positive finite \(\gamma\) gives a valid interior \(p_0\); \(\gamma=1\) gives \(p_0=p_1\). If the observed prevalence is exactly zero or one, the odds are on the boundary and this parameterization needs separate treatment rather than direct substitution into the odds formula.
The risk-difference and odds-ratio coordinates describe the same kind of unidentified contrast, but they are not interchangeable numerically. In the preceding example, \(\gamma=2\) corresponds to a risk difference \(0.5714-0.40=0.1714\). Choosing a coordinate therefore affects how experts express and vary assumptions; it does not make the underlying missing outcomes observable.
Worked Example: Coordinates Across Multiple Patterns
Worked Example: Pattern-Specific Missing-Outcome Coordinates
Consider two binary outcomes, \(Y_1\) and \(Y_2\), with three response patterns in a hypothetical population: \(R=(1,1)\), \(R=(1,0)\), and \(R=(0,1)\). Suppose their pattern probabilities are \(0.50\), \(0.30\), and \(0.20\), respectively. Within pattern \((1,1)\), the observed prevalence of \(Y_2=1\) is \(0.30\); within pattern \((0,1)\), it is \(0.70\). Within pattern \((1,0)\), the observed prevalence of \(Y_1=1\) is \(0.40\).
The missing \(Y_2\) in pattern \((1,0)\) is not specified by these observed quantities. Define two sensitivity coordinates: \(q=P(Y_2=1\mid Y_1=1,R=(1,0))\) and \(s=P(Y_2=1\mid Y_1=0,R=(1,0))\). Both must lie in \([0,1]\). The prevalence of \(Y_2=1\) in pattern \((1,0)\) is therefore \(0.40q+0.60s\), by conditioning on the observed \(Y_1\).
The overall prevalence of \(Y_2=1\) under this specification is
For the hypothetical choice \(q=0.20\) and \(s=0.60\), the pattern-\((1,0)\) prevalence is \(0.40(0.20)+0.60(0.60)=0.44\). The overall prevalence is \(0.15+0.30(0.44)+0.14=0.422\). For another valid choice, \(q=0.80\) and \(s=0.10\), the pattern-\((1,0)\) prevalence is \(0.40(0.80)+0.60(0.10)=0.38\), and the overall prevalence is \(0.15+0.30(0.38)+0.14=0.404\).
Both choices preserve every observed pattern probability and every observed outcome margin stated in the example. They differ only in the conditional distribution of the missing outcome in pattern \((1,0)\). This is the practical meaning of coordinates for unidentified information: they make explicit which aspects of the completed distribution change while the observed-data features remain fixed.
Derivation and Practical Checks
The factorization gives a direct justification for this way of working. Fix an observed-data law \(P_{R,Y_{\mathrm{obs}}}\). For each positive-probability pattern, choose a valid conditional distribution \(g_r\). Integrating \(g_r\) over all possible missing-outcome values gives one, so integrating the factorization over \(B\) equal to the whole missing-outcome space returns \(P(R=r,Y_{\mathrm{obs}}\in A)\). Thus the constructed law reproduces the fixed observed-data distribution. Conversely, any complete-data law has conditional distributions of the missing components given the pattern and observed components (up to almost-everywhere qualifications), and hence yields such a factorization.
This argument establishes a representation, not a claim that all conditional distributions are scientifically reasonable. A proposed sensitivity model can fail as a probability model or impose unintended implications. Useful checks include:
- Support and normalization: each conditional distribution must assign total probability one and respect the outcome’s possible values. For binary outcomes, every specified conditional risk must be in \([0,1]\).
- Observed-data preservation: changing sensitivity coordinates should not alter the fixed response-pattern probabilities or the distributions of recorded outcomes within patterns.
- Boundary behavior: odds, log odds, and other transformed coordinates may be undefined or infinite at probability boundaries. Inspect those cases instead of relying on an interior formula.
- Dependence structure: for several missing outcomes, marginal sensitivity parameters alone may not specify their joint distribution. The joint conditional law may require additional coordinates.
- Scientific interpretation: state what contrast a coordinate represents, its scale, and the values examined. A numerically convenient parameter is not automatically an interpretable or plausible one.
Comparisons and Limitations
A pattern-mixture specification uses the conditional distributions of outcomes within response patterns as its building blocks. A selection-model specification instead describes response probabilities conditional on outcomes. The earlier material on equivalent selection and pattern-mixture representations shows that these can be alternative factorizations of a joint law; neither factorization, by itself, makes an unidentified feature observable. Sensitivity coordinates can be used in either representation, but their meanings and convenient scales may differ.
A sensitivity parameter is also distinct from a parameter that the observed likelihood identifies. It is possible to estimate the observed-data features while the sensitivity coordinate remains unsupported by the observed likelihood. Assigning a prior distribution to such a coordinate, as in Bayesian analyses under nonidentification, can produce a posterior distribution, but that posterior’s information about the unidentified feature depends on the prior and model assumptions as well as the observed data.
The main limitation is fundamental: observed-data fit cannot determine whether a particular missing-outcome conditional distribution is correct when alternatives produce the same observed-data law. Sensitivity analysis makes the dependence of conclusions on such choices visible; it does not empirically verify the choices. Results should therefore report both the fixed observed features and the assumed sensitivity coordinates, with a scientifically justified range or set of scenarios. If a chosen parameterization rules out plausible distributions or fails to capture important dependence, conclusions drawn from it can be misleading.
Common Mistake
Wrong reasoning: “The observed records give us the respondent outcome distribution, so the missing group must have that distribution too.” That conclusion silently sets the missing-outcome conditional distribution equal to the respondent distribution. It is an assumption, not a consequence of the observed data.
Correction: Write down the observed quantities that are fixed, then identify which conditional distributions remain unspecified. If equality across response patterns is scientifically defensible, state it as a restriction and recognize that it is not tested by the observed-data distribution alone. Otherwise, define sensitivity coordinates and examine how conclusions change across valid values.
Check Your Understanding
Answer each question, then compare with the solutions below.
- A binary outcome has respondent prevalence \(0.25\), and the risk-difference coordinate is \(\delta=p_0-0.25\). What is the valid range for \(\delta\)?
- With \(P(R=1)=0.80\), respondent prevalence \(0.30\), and \(\delta=0.10\), calculate the missing-group prevalence and the complete-data prevalence.
- For respondent prevalence \(p_1=0.20\), calculate \(p_0\) when the missing-to-respondent odds ratio is \(\gamma=3\).
- Why does specifying a conditional distribution for missing outcomes not change the observed-data law in the factorization?
Solutions. (1) Since \(0\leq p_0=0.25+\delta\leq1\), \(-0.25\leq\delta\leq0.75\). (2) \(p_0=0.30+0.10=0.40\); the complete-data prevalence is \(0.80(0.30)+0.20(0.40)=0.24+0.08=0.32\). (3) \(p_0=\gamma p_1/(1-p_1+\gamma p_1)=3(0.20)/(0.80+0.60)=0.60/1.40=0.4286\). (4) The conditional distribution \(g_r\) integrates to one over the missing-outcome space. Marginalizing over missing outcomes therefore leaves the fixed measure for the pattern and its observed outcomes unchanged.