Tutorials › Multi-Omics Data Integration › Summarizing Modality Availability Patterns

Missing Modalities and Partially Observed Subjects · Tutorial 63 of 597

Summarizing Modality Availability Patterns

Build and interpret modality-pattern frequency tables, with checks that the pattern counts and availability summaries agree.

Advanced 10 min read

What You'll Learn

  • Define a subject’s modality-availability pattern from its indicator vector.
  • Construct a frequency table that includes unobserved combinations.
  • Calculate modality-specific availability rates and complete-case proportions.
  • Check that pattern frequencies and marginal summaries are internally consistent.
  • Recognize what descriptive availability counts can—and cannot—establish.

Why Summarize Availability Patterns?

A multi-omics cohort can have many different combinations of observed modalities. One subject may have RNA and protein measurements, another may have only RNA, and a third may have every planned modality. A list of modality-specific sample sizes gives useful marginal information, but it does not show which measurements occur together. Pattern-frequency tables make that joint structure explicit.

The aim here is descriptive: count the observed combinations, compare their frequencies, and identify combinations that are common, rare, or absent in the cohort being summarized. These counts help assess how much of the cohort has each modality and how many subjects have all modalities. They do not, by themselves, explain why a modality is unavailable or establish that the observed subjects represent a wider population.

Definition: Suppose a cohort contains \(n\) subjects and \(K\) modalities. Let \(A_{ij}=1\) if modality \(j\) is available for subject \(i\), and let \(A_{ij}=0\) otherwise. Subject \(i\)'s availability pattern is the length-\(K\) vector \((A_{i1},\ldots,A_{iK})\). A pattern frequency is the number of subjects with a particular vector. A pattern is absent in the cohort if its frequency is zero; it is rare when its frequency is small relative to a stated descriptive threshold or comparison.

A pattern uses a fixed modality order. For example, if the order is RNA, protein, metabolite, then \(101\) means RNA and metabolite are available, while protein is not. The order must be stated and used consistently; changing the order changes the interpretation of every string.

Before counting, define the cohort and the modalities being summarized, and make sure each subject has one correctly aligned indicator for each modality. Availability is a binary summary at the subject-by-modality level: it records whether that modality is available for the subject under the chosen definition. It does not describe how many features within a modality are measured or the reason any measurement is unavailable.

From Indicator Matrix to Pattern Frequencies

With \(K\) modalities, there are \(2^K\) mathematically possible binary patterns. A complete table lists every one of them, including those with frequency zero. Listing only patterns seen in the data can hide the fact that other combinations were considered and found to be absent.

Let \(n_r\) be the number of subjects whose pattern is \(r\), where \(r\) ranges over all \(2^K\) possible patterns. The patterns divide the cohort into mutually exclusive groups: every subject belongs to exactly one group. Consequently, the frequencies must sum to the cohort size.

$$ \sum_{r\in\{0,1\}^K} n_r=n $$

The total number with modality \(j\) available is found by adding frequencies for patterns whose \(j\)th entry is 1. Dividing by \(n\) gives that modality’s availability proportion. The all-ones pattern gives the number and proportion of subjects with every listed modality available.

$$ n_j^{\mathrm{avail}}=\sum_{r:r_j=1}n_r,\qquad \hat p_j^{\mathrm{avail}}=\frac{n_j^{\mathrm{avail}}}{n},\qquad \hat p_{\mathrm{all}}=\frac{n_{(1,\ldots,1)}}{n} $$

These are descriptive proportions for the defined cohort. They can be reported as counts, proportions, or percentages, but the denominator should be explicit. A modality’s availability proportion and the complete-case proportion answer different questions: the first concerns one modality at a time; the second concerns the simultaneous availability of all \(K\) modalities.

Worked Example: Three-Modality Pattern Table

Consider a hypothetical cohort of 20 subjects with modalities ordered as RNA, protein, and metabolite. Suppose the observed patterns are \(111\) for 8 subjects, \(110\) for 4, \(101\) for 3, \(011\) for 2, and \(100\) for 3. No other pattern occurs. Fill in the full table and summarize availability.

PatternRNAProteinMetaboliteFrequency
111AvailableAvailableAvailable8
110AvailableAvailableUnavailable4
101AvailableUnavailableAvailable3
100AvailableUnavailableUnavailable3
011UnavailableAvailableAvailable2
000UnavailableUnavailableUnavailable0
001UnavailableUnavailableAvailable0
010UnavailableAvailableUnavailable0

The observed frequencies sum to \(8+4+3+3+2=20\), as required. The unlisted patterns are \(000,001,\) and \(010\). Each has frequency zero; no other pattern is absent.

RNA is available for patterns \(111,110,101,\) and \(100\), so its count is \(8+4+3+3=18\) and its availability proportion is \(18/20=0.900\), or 90.0%. Protein is available for \(111,110,\) and \(011\): \(8+4+2=14\), giving \(14/20=0.700\), or 70.0%. Metabolite is available for \(111,101,\) and \(011\): \(8+3+2=13\), giving \(13/20=0.650\), or 65.0%.

Only pattern \(111\) has all three modalities, so the complete-case count is 8 and the complete-case proportion is \(8/20=0.400\), or 40.0%. The marginal availability proportions are each larger than or equal to this proportion because having all three modalities necessarily means having any one of them.

A Short Justification for the Checks

The frequency and marginal summaries are linked by the definition of a pattern. To see why, define \(I_i(r)=1\) when subject \(i\) has pattern \(r\), and \(I_i(r)=0\) otherwise. Since every subject has exactly one pattern, the indicators over all possible patterns sum to 1 for each subject. Summing over subjects gives the total-frequency identity.

$$ \sum_{r\in\{0,1\}^K} n_r =\sum_{r\in\{0,1\}^K}\sum_{i=1}^n I_i(r) =\sum_{i=1}^n 1 =n $$

For modality \(j\), a subject contributes to the availability count exactly when their pattern has \(r_j=1\). Thus adding the frequencies of precisely those patterns counts each subject with modality \(j\) once and every other subject zero times. This also gives a practical audit: calculate marginal counts both from the indicator matrix and from the pattern table, then verify they match.

Key check: Pattern frequencies must sum to \(n\); each modality’s summed pattern count must match its column sum in the availability indicator matrix; and the all-ones frequency cannot exceed any individual modality availability count.

Worked Example: Four Modalities and a Rare Combination

In a hypothetical cohort of 12 subjects, use the order RNA, protein, metabolite, and lipid. Suppose the nonzero frequencies are \(1111: 4\), \(1110: 3\), \(1011: 2\), \(1000: 2\), and \(0111: 1\). The colon separates a pattern from its frequency.

First check the cohort total: \(4+3+2+2+1=12\). For this four-modality table, there are \(2^4=16\) possible patterns. The five listed patterns are observed; the other 11 patterns have frequency zero. The pattern \(0111\) is observed once, so it is rare relative to the cohort in the straightforward descriptive sense of one subject out of 12. Whether to label a frequency “rare” should be made explicit rather than treated as a universal statistical cutoff.

RNA is available in \(1111,1110,1011,\) and \(1000\), for a count of \(4+3+2+2=11\), or \(11/12=0.917\) (91.7%). Protein is available in \(1111,1110,\) and \(0111\), for \(4+3+1=8\), or \(8/12=0.667\) (66.7%). Metabolite is available in \(1111,1011,\) and \(0111\), for \(4+2+1=7\), or \(7/12=0.583\) (58.3%). Lipid is available in \(1111,1011,\) and \(0111\), also for \(4+2+1=7\), or 58.3%.

The complete-case count is the \(1111\) frequency, 4, so the complete-case proportion is \(4/12=0.333\), or 33.3%. These calculations show why a low-frequency combination can matter even when one or more marginal availability counts are relatively high: the margins do not, on their own, describe how modalities co-occur.

Comparing Groups and Choosing a Summary

When the cohort contains predefined groups, such as collection sites, summarize each group with its own denominator before reporting an overall table. This makes the group-specific pattern distribution visible and prevents a pooled count from obscuring differences in the composition of the groups. The example below is descriptive; it does not identify what caused the difference.

Worked Example: Availability by Collection Site

Suppose a hypothetical 16-subject cohort has two sites with 8 subjects each and modalities ordered RNA, protein, metabolite. At Site A, the frequencies are \(111:4\), \(110:2\), and \(100:2\). At Site B, they are \(111:1\), \(101:3\), \(011:2\), and \(001:2\).

At Site A, the frequencies sum to \(4+2+2=8\). RNA is available for all 8 subjects; protein for \(4+2=6\); and metabolite for 4. The corresponding proportions are \(8/8=1.00\), \(6/8=0.750\), and \(4/8=0.500\). The complete-case proportion is \(4/8=0.500\).

At Site B, the frequencies sum to \(1+3+2+2=8\). RNA is available for \(111\) and \(101\), giving \(1+3=4\), or \(4/8=0.500\). Protein is available for \(111\) and \(011\), giving \(1+2=3\), or \(3/8=0.375\). Metabolite is available for \(111,101,011,\) and \(001\), giving \(1+3+2+2=8\), or \(1.00\). The complete-case proportion is \(1/8=0.125\).

Pooling the two sites gives \(111:5\), \(110:2\), \(100:2\), \(101:3\), \(011:2\), and \(001:2\), which sum to 16. The pooled availability counts are RNA \(=5+2+2+3=12\), protein \(=5+2+2=9\), and metabolite \(=5+3+2+2=12\). Their proportions are \(12/16=0.750\), \(9/16=0.5625\), and \(12/16=0.750\). The pooled complete-case proportion is \(5/16=0.3125\).

The pooled table is a valid description of these 16 subjects, but it hides the different site-specific distributions. Reporting the stratified counts alongside pooled summaries gives readers both views without asserting a reason for the difference.

Diagnostics, Alternatives, and Limitations

A pattern table is only as reliable as the subject alignment and availability definition behind it. Before interpreting the table, verify that subject identifiers refer to the same cohort across modalities, each subject appears once in the summary, and the order of modalities is documented. Then check the frequency total and independently reproduce each marginal count. If a pattern is not shown, determine whether it was deliberately omitted or has frequency zero; a complete table makes that distinction unambiguous.

A common mistake is to report only marginal sample sizes and assume they reveal the joint pattern. They do not: two cohorts can have identical counts for each modality but different combinations of modalities observed on the same subjects. Conversely, a long pattern table may be cumbersome when \(K\) is large, since the number of possible patterns grows as \(2^K\). In that setting, a concise report can show the most frequent patterns, the number of unobserved patterns, and key marginal and all-modality counts, while making clear what was omitted.

The main alternatives serve different purposes. Marginal availability counts are compact and useful when the question concerns each modality separately, but conceal co-occurrence. The complete-case count is compact when the question is specifically how many subjects have every modality, but discards information about partial combinations. Pattern frequencies retain that joint descriptive information, at the cost of a potentially large table. None of these summaries establishes why a modality is absent, whether availability is related to subject characteristics, or what an analysis should do with partially observed data. Those are separate questions.

Common pitfall: “RNA is available for 90% of subjects and protein for 70%, so 70% have both.” This does not follow from the marginal proportions. The overlap must be counted from the joint patterns; in the three-modality example, RNA and protein are both available for \(8+4=12\) of 20 subjects, or 60%.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. For three modalities in a fixed order, how many possible availability patterns are there? What does pattern \(010\) mean?
  2. A cohort has 30 subjects. The frequencies of all listed patterns sum to 28. Give two checks to make before interpreting the table.
  3. In a cohort of 20, patterns \(111:6\), \(110:5\), \(101:3\), and \(100:6\) are the only observed patterns. Find the RNA availability count and proportion, and the complete-case proportion.
  4. Why can modality-specific availability proportions not determine the frequency of a particular pair of modalities being jointly available?

Solutions. (1) There are \(2^3=8\) patterns. With the stated order, \(010\) means only the second modality is available. (2) Check whether subjects were omitted or duplicated, and whether the table includes all possible patterns or has an arithmetic/transcription error; by definition the frequencies must sum to 30. (3) RNA is available in \(111,110,101,\) and \(100\), so the count is \(6+5+3+6=20\), and the proportion is \(20/20=1.00\). The complete-case proportion is \(6/20=0.300\). (4) Marginal proportions count availability one modality at a time and do not show which subjects share both. The joint frequency must be obtained by summing patterns in which both specified entries are 1.

Clinical Biostats