Tutorials › Regression Analysis (Undergraduate) › Missing Measurements in a Nutrition Cohort

Integrated Applied Case Studies · Tutorial 981 of 991

Missing Measurements in a Nutrition Cohort

Use a nutrition-cohort example to compare complete-case and imputed regression results, and learn what each analysis assumes about missing measurements.

Advanced 11 min read

What You'll Learn

  • Define complete-case analysis and calculate how much cohort data it retains.
  • Distinguish MCAR, MAR, and MNAR missingness assumptions.
  • Interpret a complete-case regression estimate and its uncertainty in context.
  • Pool estimates and uncertainty across multiple imputed data sets.
  • Compare analyses without treating imputation as proof that missing-data bias is absent.
  • Identify practical checks and sensitivity analyses for incomplete cohort measurements.

Why Missing Measurements Matter

A nutrition cohort may record participants’ usual diet, blood pressure, age, and body mass index, but some participants skip a diet questionnaire, miss a clinic visit, or leave a question blank. A regression fitted to the available records can still produce a coefficient, but the people contributing to that coefficient may differ systematically from those whose measurements are missing. The key issue is therefore not only how many values are absent, but also what process made them absent.

This tutorial considers a hypothetical cohort in which the research question is whether usual sodium intake is associated with systolic blood pressure after adjusting for age, sex, body mass index, and energy intake. The study is observational, so an estimated association does not by itself establish that changing sodium intake would change blood pressure.

Model: For participant \(i\), let \(Y_i\) be systolic blood pressure in mmHg and \(X_i\) be sodium intake in grams per day. The adjusted linear regression model is $$ Y_i=\beta_0+\beta_1X_i+\beta_2\mathrm{Age}_i+\beta_3\mathrm{Sex}_i+\beta_4\mathrm{BMI}_i+\beta_5\mathrm{Energy}_i+\epsilon_i. $$ Here \(\beta_1\) is the expected difference in systolic blood pressure, in mmHg, associated with a one-gram-per-day difference in sodium intake, among participants equal on the included adjustment variables. The model assumes an appropriately specified conditional mean, independent participant errors, and a suitable variance model for the reported standard errors; exact small-sample \(t\) inference additionally relies on normally distributed errors. These assumptions concern the regression, separately from assumptions about missingness.

What Complete-Case Analysis Does

A complete case is a participant with an observed value for the outcome and every variable included in the regression model. Complete-case analysis fits the model using only those participants. It is simple and transparent, but it can discard substantial data when missing values occur in different variables.

To reason about the consequences, let “observed” indicate that a value is recorded and “missing” indicate that it is not. Three standard descriptions of the missingness process are useful:

  • Missing completely at random (MCAR): the chance a measurement is missing is unrelated to the observed data and to the unobserved value. For example, a randomly lost batch of forms could produce MCAR missingness.
  • Missing at random (MAR): after conditioning on observed information, missingness does not additionally depend on the value that is missing. For example, if younger participants are less likely to complete a diet questionnaire, age may explain that pattern if age is observed.
  • Missing not at random (MNAR): even after accounting for observed information, the chance a value is missing still depends on that unobserved value. For example, participants with especially high sodium intake might be less likely to report it, even among people with the same recorded age and health measures.

These labels describe assumptions about the process, not patterns that can always be diagnosed from the observed data. In particular, observed associations between missingness and recorded variables can challenge MCAR, but cannot prove MAR: the missing values themselves are unavailable for a direct check.

Worked Example: Measuring How Much Complete-Case Analysis Retains

Suppose a hypothetical cohort enrolls 240 adults. Sodium intake is missing for 36, systolic blood pressure for 20, body mass index for 24, and energy intake for 30. A record audit shows that 168 participants have all variables required by the model. Calculate the sodium and outcome missingness percentages and the complete-case retention percentage.

The sodium missingness percentage is $$ \frac{36}{240}\times100\%=15.0\%. $$ The blood-pressure missingness percentage is $$ \frac{20}{240}\times100\%=8.33\%. $$ The proportion retained as complete cases is $$ \frac{168}{240}\times100\%=70.0\%. $$

Thus, 30.0% of the cohort is excluded from this particular regression. The missingness counts for separate variables cannot simply be added to find the number excluded, because the same participant may be missing more than one measurement. The audited joint count of 168 complete records is what determines the complete-case sample size.

MCAR is a sufficient, relatively strong condition under which complete cases are not selected according to the study measurements. It is not the only possible condition under which a particular regression coefficient may remain unbiased. For example, selection can sometimes preserve the conditional relationship being estimated. But MAR alone does not guarantee that complete-case regression is unbiased: exclusion related to recorded age or health can change the composition of the analyzed sample and distort the fitted relationship. The mechanism and the model both matter.

Reading a Complete-Case Result

Before interpreting a complete-case estimate, state which participants contributed and what assumptions support its uncertainty calculation. In the example below, suppose checks of the hypothetical complete-case data support the chosen linear mean form and there is no evident dependence between participants. Assume the conventional model-based standard error is appropriate; if the variance pattern makes it unsuitable, the uncertainty calculation would need an appropriate alternative. These checks do not establish that complete cases represent participants with missing measurements.

Worked Example: Interpreting the Complete-Case Slope

In the 168 complete cases, suppose the fitted sodium coefficient is \(1.60\) mmHg per gram per day, with standard error \(0.872\) mmHg per gram per day. There are six fitted parameters including the intercept, so \(df=168-6=162\). Test \(H_0:\beta_1=0\) against a two-sided alternative at \(\alpha=0.05\), and calculate a 95% confidence interval.

The test statistic is $$ t=\frac{1.60-0}{0.872}=1.835. $$ For \(df=162\), the two-sided 5% critical value is approximately \(1.975\). Since \(|1.835|<1.975\), the result does not reject \(H_0\) at the 5% level; the two-sided \(p\)-value is approximately \(0.068\).

The 95% interval is estimate \(\pm\) critical value times standard error: $$ 1.60\pm1.975(0.872)=1.60\pm1.7222. $$ Therefore the interval is approximately \(-0.122\) to \(3.32\) mmHg per gram per day. In this complete-case sample, the data are compatible with adjusted associations ranging from a small decrease to an increase of about 3.32 mmHg per one-gram-per-day difference in sodium intake. The interval includes zero, so the analysis does not establish a nonzero association at the 5% level. It does not show that the true association is zero.

Multiple Imputation as a Practical Alternative

Rather than deleting every incomplete record, multiple imputation creates several completed versions of the data. Each version fills in missing values with plausible draws that reflect uncertainty about those values. The regression is fitted separately to each completed data set, and the estimates and their uncertainty are combined. Unlike replacing every missing value once with a mean, proper multiple imputation carries information about imputation uncertainty into the final standard error.

For this regression, an imputation model should generally include the outcome, all variables in the analysis model, and useful observed auxiliary variables that help predict missingness or the missing measurements. It should respect variable types and preserve important features of the analysis, such as nonlinear terms or interactions if those are in the regression. Including systolic blood pressure when imputing missing sodium helps retain the sodium–blood pressure relationship; imputing a predictor without using the outcome can weaken that relationship.

A practical analysis commonly relies on MAR conditional on the information included in the imputation model, along with appropriate imputation models and adequate computation. Imputation does not make that assumption true. Implausible imputed values, poorly chosen predictors, or an imputation model inconsistent with the analysis can still produce misleading results.

Pooling rule: With \(m\) imputed data sets, let \(Q_j\) be the fitted coefficient and \(U_j=SE_j^2\) its estimated variance in data set \(j\). The pooled coefficient is \(\bar{Q}=\frac{1}{m}\sum_{j=1}^{m}Q_j\). The average within-imputation variance is \(\bar{U}=\frac{1}{m}\sum_{j=1}^{m}U_j\), and the between-imputation variance is \(B=\frac{1}{m-1}\sum_{j=1}^{m}(Q_j-\bar{Q})^2\). The total variance is \(T=\bar{U}+(1+\frac{1}{m})B\), so the pooled standard error is \(\sqrt{T}\). The added term accounts for uncertainty because different plausible completed data sets give different estimates.

Pooling an Imputation Analysis

The pooling calculation follows directly from separating two sources of uncertainty. There is sampling uncertainty within each completed data set, summarized by \(\bar{U}\), and there is variation among the completed-data estimates, summarized by \(B\). Because the average across \(m\) imputations estimates the missing-data uncertainty, its contribution is inflated by \(1+1/m\), giving \(T=\bar{U}+(1+1/m)B\). When a confidence interval is needed, a \(t\) critical value with an imputation-adjusted degrees of freedom is commonly used; with many imputations and ample data, it is often close to the normal critical value 1.96. The calculation below uses that close approximation for illustration.

Worked Example: Pooling Five Imputed Regression Fits

Suppose five hypothetical completed data sets yield sodium coefficient estimates \(1.18,\ 1.34,\ 1.26,\ 1.42,\ 1.30\), each with standard error \(0.840\) mmHg per gram per day. Calculate the pooled estimate, standard error, and approximate 95% interval.

The pooled coefficient is $$ \bar{Q}=\frac{1.18+1.34+1.26+1.42+1.30}{5}=\frac{6.50}{5}=1.30. $$ Each within-imputation variance is \(0.840^2=0.7056\), so \(\bar{U}=0.7056\). The deviations from the pooled estimate are \(-0.12,\ 0.04,\ -0.04,\ 0.12,\ 0.00\). Thus, $$ B=\frac{(-0.12)^2+(0.04)^2+(-0.04)^2+(0.12)^2+0^2}{5-1} =\frac{0.0320}{4}=0.00800. $$ The total variance and pooled standard error are $$ T=0.7056+\left(1+\frac{1}{5}\right)(0.00800)=0.7152, \qquad SE=\sqrt{0.7152}=0.846. $$ Using 1.96 as the approximate critical value gives $$ 1.30\pm1.96(0.846)=1.30\pm1.658, $$ or an approximate 95% interval from \(-0.358\) to \(2.96\) mmHg per gram per day.

This imputation analysis estimates an adjusted association of 1.30 mmHg per gram per day, but its interval includes zero. It uses more participants’ partial information than the complete-case analysis, yet its conclusion remains uncertain. Five imputations are used here only to make the arithmetic visible; an actual analysis should use enough imputations for stable estimates and uncertainty, particularly when missing information is substantial.

Comparing Results and Checking Assumptions

The complete-case estimate, 1.60, and the pooled imputation estimate, 1.30, are not identical because they use information differently and need not target the same result under every missingness process. Their confidence intervals overlap substantially, and both include zero in this hypothetical example. That comparison is informative, but it is not a test proving that either missing-data assumption is correct. A report should state the sample sizes, missingness patterns, imputation variables and methods, and the assumptions on which interpretation depends.

Investigators can compare observed characteristics of participants with and without complete measurements, and tabulate missingness across relevant recorded groups. Such checks can identify patterns that make MCAR doubtful and guide the imputation model. They cannot determine whether missingness depends on the unobserved sodium value after conditioning on observed information. When MNAR is plausible, a sensitivity analysis can repeat the analysis under explicit alternative assumptions—for example, shifting imputed sodium values by a prespecified amount—and show whether the substantive conclusion changes. The shift and its rationale must be reported; sensitivity analysis measures dependence on assumptions rather than eliminating it.

Common Mistakes

  • Wrong: “Only 15% of sodium values are missing, so missingness cannot matter.” Correction: Even a moderate missing fraction can bias a coefficient if missingness relates to the outcome or exposure. Examine joint patterns and the likely reasons for missingness.
  • Wrong: “The imputed analysis uses everyone, so it is unbiased.” Correction: Imputation relies on a model and assumptions such as MAR conditional on included information; it is not a guarantee against bias.
  • Wrong: “Fill each gap with the observed mean, then run ordinary regression.” Correction: Single mean filling understates variability and can distort associations. Multiple imputation represents uncertainty across plausible completions and pools that uncertainty.
  • Wrong: “The two analyses differ, so one is proven correct.” Correction: Differences can signal sensitivity to missing-data handling, but observed data alone may not reveal the true missingness mechanism.
Key takeaway: Complete-case and imputation analyses answer the regression question using different assumptions about incomplete records. Report the missingness, model, method, and uncertainty, and interpret conclusions with the limits of those assumptions in view.

Check Your Understanding

Use the definitions and calculations in this tutorial to answer the following questions.

  1. A cohort enrolls 300 people, of whom 225 have every variable needed for a regression. What percentage are complete cases, and what additional information is needed to know how many were excluded because of each variable?
  2. A diet questionnaire is less often completed by younger participants, and age is observed for everyone. Which missingness mechanism might be plausible if age explains the pattern, and what can the observed data not establish?
  3. In multiple imputation, why does the total variance include both within-imputation and between-imputation variance?
  4. Why should systolic blood pressure generally be included when imputing missing sodium values for a regression of blood pressure on sodium?
  5. If a complete-case confidence interval and a multiple-imputation confidence interval both include zero, what can be concluded—and what cannot be concluded—about the adjusted association?