Tutorials › Pattern-Mixture and Selection Models › Observed-Data Distributions and Identification Targets

Identifiability and Untestable Assumptions · Tutorial 53 of 499

Observed-Data Distributions and Identification Targets

Learn how missingness maps the joint law to observed data, and how to decide whether a target is identified by the observed-data distribution.

Advanced 12 min read

What You'll Learn

  • Define the observed-data record for a study with incomplete outcomes.
  • Derive the observed-data distribution by summing or integrating over unobserved outcomes.
  • Distinguish identifiable observed-data features from complete-data targets.
  • Interpret response rates and respondent outcome distributions without treating them as population outcome distributions.
  • Recognize when an assumption is needed to connect observed outcomes to missing outcomes.
  • Check identification claims by specifying the target and the model assumptions.

Why Separate the Observed Law From the Complete-Data Law?

A study with missing outcomes produces a distribution of what was actually recorded, not a direct sample from the distribution that would have been available if every outcome had been measured. Keeping those distributions distinct is essential. A response rate and an outcome distribution among respondents are features of the observed data. The overall outcome distribution, which includes outcomes that were not observed, is a feature of the complete data and may require assumptions beyond the recorded information.

This distinction is not a technicality about notation. It determines what the data alone can establish. The observed-data distribution describes the probabilities of observed records. A complete-data target, such as the population mean outcome, is identified only if its value is determined by that observed distribution under the model assumptions being used.

Defining the Observed-Data Distribution

Let \(X\) denote fully observed covariates, \(Y\) an outcome that may be missing, and \(R\) its response indicator, with \(R=1\) if \(Y\) is observed and \(R=0\) otherwise. A complete-data record is \((X,Y,R)\). The recorded outcome can be written \(Y_{\mathrm{obs}}\): it equals \(Y\) when \(R=1\) and is marked missing when \(R=0\). Thus an observed-data record is \(O=(X,R,Y_{\mathrm{obs}})\). This definition does not require that missingness be random or that any particular model for \(R\) hold.

Definition: The complete-data distribution is the joint distribution of all variables in the full-data record, such as \(P(X,Y,R)\). The observed-data distribution is the distribution of the record \(O\) that is actually observed. It is obtained from the complete-data distribution by retaining the response pattern and the outcomes observed under that pattern, and summing or integrating over outcomes that are missing.

For a discrete outcome, the observed-data probabilities follow directly from marginalization. For each covariate value \(x\), the probability of observing \(Y=y\) with \(R=1\) is the joint probability \(P(X=x,Y=y,R=1)\). The probability of \(X=x\) and \(R=0\) combines all possible values of the missing outcome:

$$ P(X=x,R=0)=\sum_y P(X=x,Y=y,R=0). $$

For a continuous outcome, replace the sum over missing \(y\) by an integral. More generally, if \(R\) indexes a pattern of observed and missing outcomes and \(Y_R\) denotes the components observed in pattern \(R\), then the observed-data probability or density is obtained by integrating the complete-data joint law over the missing components \(Y_{\bar R}\):

$$ p_{\mathrm{obs}}(x,r,y_r) = \int p_{\mathrm{full}}(x,y_r,y_{\bar r},r)\,dy_{\bar r}. $$

For discrete missing components, the integral is replaced by a sum. The resulting observed-data law includes the probabilities of response patterns as well as the distributions of outcomes observed within those patterns. It does not include the unrecorded values themselves.

What Does “Identified” Mean Here?

Identification is a statement about what follows from the observed-data distribution under a specified set of assumptions. It is not a statement about whether a particular sample is large or whether an estimator is precise. To make the statement precise, let \(\mathcal{M}\) be the set of complete-data distributions allowed by the model assumptions. Each law in \(\mathcal{M}\) induces an observed-data distribution by marginalization.

Definition: A target \(T(P_{\mathrm{full}})\) is identified from the observed data under \(\mathcal{M}\) if any two complete-data distributions in \(\mathcal{M}\) that induce the same observed-data distribution give the same value of the target. In symbols, whenever \(P_{\mathrm{obs}}^{(1)}=P_{\mathrm{obs}}^{(2)}\), the target must satisfy \(T(P_{\mathrm{full}}^{(1)})=T(P_{\mathrm{full}}^{(2)})\).

The observed-data distribution itself is, by definition, the distribution of the recorded record. Its features are therefore the natural starting point for identification. Examples include the probability of each response pattern and the distribution of \(Y\) among records with \(R=1\). A complete-data target such as \(P(Y=1)\) or \(E(Y)\) is different: it includes contributions from outcomes when \(R=0\). Whether that target is identified depends on whether the model assumptions determine those contributions from the observed law.

This definition is relative to \(\mathcal{M}\). A target can be unidentified when only the observed-data law is specified, yet identified under a more restrictive model. The identifying restriction is then an assumption added to the observed information; it is not another observed-data feature. The distinction between what the records determine and what a model stipulates should remain explicit.

Worked Example: Response Rate and Respondent Success Rate

In a hypothetical sample of 100 people, the binary outcome \(Y\) is observed for 70 people. Among those 70, 42 have \(Y=1\). Suppose the sampled records are independent and are being used to describe a population with the same sampling scheme.

Observed response rate. The sample estimate is \(\hat{P}(R=1)=70/100=0.700\), or 70.0%. This is a feature of the observed response indicators.

Outcome distribution among respondents. The estimated conditional probability of success among respondents is \(\hat{P}(Y=1\mid R=1)=42/70=0.600\), or 60.0%. The denominator is the 70 people whose outcomes were observed, so this describes respondents, not automatically all 100 people.

Complete-data prevalence. The complete-data prevalence would count successes among both respondents and nonrespondents. The 30 missing outcomes do not tell us how many of those people had \(Y=1\). Consequently, the observed response rate and respondent success rate do not by themselves supply the complete-data prevalence. One would need an assumption or other information connecting the missing outcomes to the observed information.

The two calculated proportions are directly summaries of observed records. Treating 0.600 as the population success probability would silently replace the target \(P(Y=1)\) with the different target \(P(Y=1\mid R=1)\).

Deriving the Marginalization and Reading Its Consequences

The observed-data formula is a direct consequence of adding the probabilities of complete records that produce the same observed record. For example, when \(X=x\) and \(R=0\), the recorded information does not distinguish among the possible values of \(Y\). The disjoint events \(\{X=x,Y=y,R=0\}\), over all possible \(y\), make up the event \(\{X=x,R=0\}\). Additivity of probability gives:

$$ \begin{aligned} P(X=x,R=0) &=P\left(\bigcup_y\{X=x,Y=y,R=0\}\right)\\ &=\sum_y P(X=x,Y=y,R=0). \end{aligned} $$

The same argument applies to multiple missing outcomes: the observed record corresponds to all complete records that agree on the observed components, and their probabilities are combined. This derivation explains both what marginalization retains and what it loses. The observed law retains the total probability assigned to each observed record, but it does not reveal how that probability is divided among complete records that differ only in missing outcomes.

In an independent sample, the likelihood for the observed records is built from their probabilities under \(P_{\mathrm{obs}}\). A record with a missing outcome contributes the probability of its observed pattern and covariates, with the missing outcome marginalized out. It does not contribute a probability evaluated at an unobserved value. This is why an observed-data likelihood is not simply the complete-data likelihood with blank outcomes treated as known values.

Worked Example: An Observed Distribution for a Continuous Outcome

In a hypothetical sample of 10 people, \(Y\) is observed for four people, with values 2, 4, 6, and 8 units. The other six outcomes are missing. Assume the records are independent for interpreting sample proportions as estimates of population probabilities.

The observed response proportion is \(4/10=0.400\). The empirical distribution of \(Y\) among respondents assigns probability \(1/4\) to each of 2, 4, 6, and 8 units. Its mean is \[ \bar{x}_{R=1}=\frac{2+4+6+8}{4}=\frac{20}{4}=5.00 \] units.

The calculation establishes a respondent mean of 5.00 units. It does not establish that the mean among all 10 people is 5.00 units: that overall mean would include the six missing outcome values. Nothing in the four recorded values supplies those values by arithmetic alone. A procedure that replaces missing values with 5.00 would add a modeling or imputation rule; it would not be a consequence of the observed-data distribution itself.

Covariates and Pattern-Specific Features

When covariates are observed for everyone, the observed-data law can describe response and outcome distributions conditionally on \(X\). These conditional features can be important: response rates may differ across covariate groups, and the respondent outcome distribution may also differ across groups. Conditioning does not, by itself, turn a respondent distribution into the outcome distribution for everyone in that group.

Worked Example: Observed Features Within Covariate Groups

Consider a hypothetical sample of 80 people, with 40 in each of two fully observed covariate groups, \(X=0\) and \(X=1\). In group \(X=0\), outcomes are observed for 30 people, of whom 18 have \(Y=1\). In group \(X=1\), outcomes are observed for 20 people, of whom 8 have \(Y=1\).

For \(X=0\), the estimated response rate is \(30/40=0.750\). Among respondents in that group, the estimated success probability is \(18/30=0.600\). As proportions of the entire sample, the observed cells \((X=0,R=1,Y=1)\) and \((X=0,R=1,Y=0)\) have estimates \(18/80=0.225\) and \(12/80=0.150\), respectively.

For \(X=1\), the estimated response rate is \(20/40=0.500\), and the respondent success probability is \(8/20=0.400\). The observed sample proportions for \((X=1,R=1,Y=1)\) and \((X=1,R=1,Y=0)\) are \(8/80=0.100\) and \(12/80=0.150\).

These calculations describe response and outcome features in the observed records, including their variation by \(X\). They do not calculate \(P(Y=1\mid X=0)\) or \(P(Y=1\mid X=1)\) for everyone in the respective groups: the outcomes of nonrespondents remain unobserved. Stratification makes the observed description more detailed, but does not by itself resolve that missing-outcome information.

Assumptions, Diagnostics, and Common Mistakes

The mapping from the complete-data law to the observed-data law is definitional. It requires no missing-at-random assumption. Assumptions become necessary when the goal is to relate the observed law to features of the complete-data law that include unobserved outcomes. Before interpreting such a target, state the target and the restrictions being used. The earlier material on factorizing joint distributions helps express a joint law in a convenient order; changing the factorization does not change which records were observed.

Useful checks begin with the observed record itself. Verify that the response indicator and missingness codes are defined consistently, tabulate response patterns, and check that each reported outcome summary uses the correct observed subgroup as its denominator. With covariates, inspect pattern frequencies within covariate groups; an empty observed pattern in a group can prevent estimation of features that require that pattern. These checks can reveal coding problems or sparse data. They cannot verify the values of outcomes that were never recorded.

A particularly common mistake is to reason, “The mean among respondents is 5, so the population mean is 5.” The correction is to name the conditional target: the calculation estimates \(E(Y\mid R=1)\), not necessarily \(E(Y)\). Equating these quantities needs a justified assumption that connects response and outcome, or other identifying information. A second mistake is to treat an assumption as if it were a feature visible in the observed-data distribution. Assumptions can restrict the model and thereby identify a target, but their validity is not established merely because the observed records fit the model.

QuantityDefined fromWhat it describes
\(P(R=1)\)Observed response indicatorsProbability an outcome is recorded
\(P(Y=y\mid R=1)\)Observed outcomes among respondentsOutcome distribution among respondents
\(P(Y=y)\)Complete-data distributionOutcome distribution for the full population
\(E(Y)\)Complete-data distributionPopulation mean including missing outcomes

The main alternatives differ in what they target and how they bridge missing information. A complete-case summary uses only records with observed outcomes and directly describes that selected subset. An observed-data likelihood uses the full distribution of recorded patterns and outcomes, with missing components marginalized. A model-based complete-data analysis may additionally impose assumptions about the relation between response and the unobserved outcome. These approaches should not be confused: a more elaborate model does not change the observed-data law, though it may make additional targets identifiable conditional on its assumptions.

Key takeaway: First write down the distribution of the records actually observed. Then state the complete-data target separately and ask whether the allowed assumptions determine it from that observed distribution. A respondent summary is not automatically a full-population summary.

Check Your Understanding

Answer each question by identifying the relevant observed or complete-data quantity. Solutions follow the questions.

  1. In a sample of 50 people, outcomes are observed for 35. What is the sample response proportion, and what population feature does it describe?
  2. Among those 35 respondents, 21 have a binary outcome equal to 1. Calculate the respondent success proportion. Is it automatically the success proportion for all 50 people?
  3. For discrete \(Y\), express \(P(X=x,R=0)\) in terms of the complete-data joint probabilities.
  4. State the condition that makes a complete-data target identified under a specified model class.

Solutions. (1) The sample response proportion is \(35/50=0.700\), or 70.0%; it estimates the response probability under the sampling scheme. (2) The respondent success proportion is \(21/35=0.600\), or 60.0%. It describes respondents and is not automatically the proportion for all 50, because the 15 missing outcomes are not included. (3) Add over all possible outcome values: \(P(X=x,R=0)=\sum_y P(X=x,Y=y,R=0)\). (4) The target is identified if any two complete-data distributions allowed by the model that induce the same observed-data distribution give the same target value.

Clinical Biostats