Pith. Sign in

REVIEW

Functional linear regression analysis for longitudinal data

Statistics Theory · Submitted 6 Mar 2006 · 31 references

Pith's verdict: Solid, influential paper that correctly extends FPCA-based functional regression to sparse longitudinal data; the confidence-band theory is asserted rather than proved, but the consistency results hold up.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDmath/0603132 (version 1)
Pith ID4QE4GTWH, Pith's permanent record of this paper and its review
CitationFang Yao, Hans-Georg M\"uller, Jane-Ling Wang. "Functional linear regression analysis for longitudinal data." arXiv:math/0603132 (2006). https://pith.science/paper/4QE4GTWH
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Sparse data still support consistent functional regression.

The reading

The paper claims that functional linear regression—where a response curve is modeled as a linear functional of a predictor curve—can be made to work when each subject contributes only a few noisy measurements at irregular times. It proposes estimating the regression surface and predicting new response trajectories through functional principal component curves, recovering the sparse per-subject scores by conditional expectation rather than by integrating over each observed curve. The central result is consistency: as the number of subjects grows, the estimated regression function converges to the true one in integrated squared error, and predicted trajectories converge to the best predictor available from the sparse observations. This matters because longitudinal studies—clinic visits, liver disease follow-up, blood pressure and body mass index—routinely produce sparse, irregular data for which earlier functional regression methods required dense or complete trajectory sampling.

What carries the argument

The load-bearing device is the conditional-expectation estimate of a new subject's functional principal component scores, $\hat{\zeta}_m^* = \hat{\rho}_m \hat{\psi}_m^{*T} \hat{\Sigma}_{U^*}^{-1}(\tilde{U}^* - \hat{\mu}_X^*)$ from (11), which extends the best linear prediction formula (10) to sparse data. This step replaces the usual integral approximation of scores with a prediction that borrows strength from the whole sample, which is what makes the method work when only three to five measurements per trajectory are available. Around this step the paper assembles the basis representation $\beta(s,t)=\sum_{k,m}(\sigma_{km}/\rho_m)\psi_m(s)\phi_k(t)$, the smoothed cross-covariance surface $\hat{C}$, and the Gaussian variance identity $\Omega_M = D - H\Sigma_{U^*}^{-1}H^T$ that drives the confidence bands.

Extended reading notes

Core claim

The paper's central claim is that estimation and prediction in the functional linear model remain consistent under sparse, irregular longitudinal sampling. For the regression surface, Theorem 1 states that with $\hat{\beta}(s,t)$ defined in (8) from estimated eigenfunctions, eigenvalues, and a smoothed cross-covariance surface, $\int_T \int_S [\hat{\beta}(s,t)-\beta(s,t)]^2\,ds\,dt \to 0$ in probability as the number of subjects $n\to\infty$ under eigenvalue-separation and bandwidth conditions, with a uniform version under a stronger summability assumption. Theorem 2 states that the predicted trajectory $\hat{Y}^*_{KM}(t)$ in (12) converges to the target $\tilde{Y}^*(t)$, the conditional-expectation predictor built from the sparse measurements of a new subject's predictor trajectory. Under joint Gaussianity of the predictor's principal component scores and measurement errors, Theorem 3 gives asymptotic normality of the prediction error and justifies the pointwise confidence bands in (13). The paper also defines a functional coefficient of determination $R^2 = \sum_{k,m}\sigma_{km}^2/\rho_m / \sum_k \lambda_k$ and proves that its plug-in estimator is consistent.

Load-bearing premise

The load-bearing premise is that the predictor's principal component scores and measurement errors are jointly Gaussian, because this is what makes the best linear score predictor equal to the conditional expectation and supplies the variance formula $\Omega_M = D - H\Sigma_{U^*}^{-1}H^T$ behind the confidence bands.

Editorial extensions

If this is right

  • Regression-function estimation remains consistent when per-subject measurement counts are bounded and irregular, provided the number of subjects grows and the eigenstructure is well separated.
  • Response trajectories for new subjects can be predicted from sparse predictor measurements alone, with predictions converging to the best available conditional predictor.
  • Asymptotic pointwise confidence bands for predicted trajectories are available under joint Gaussianity of scores and errors; the paper's simulation evidence indicates the predictions remain useful when Gaussianity fails.
  • The functional coefficient of determination $R^2$ gives a consistent, interpretable measure of the variance of the response trajectories explained by the predictor, and the paper demonstrates it on two longitudinal datasets.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The conditional-expectation score predictor could be extended to non-Gaussian settings via an elliptical or rank-based version of the best linear predictor; consistency would likely carry over, but the exact normal bands would need a different derivation.
  • Because the method permits predictor and response measurements at different times, it is a natural template for regression between two longitudinally measured biomarkers when visit schedules and missingness differ between the two outcomes.
  • A practical danger point is the effective number of subjects with at least two measurements: the theory needs only a positive probability of multiple measurements, but users should verify that this subset is large enough to identify the covariance eigenfunctions.
  • A direct test of the Gaussian assumption's role would be to measure the coverage of the nominal bands (13) under skewed or heavy-tailed scores; the paper's robustness simulations address prediction error, not band coverage.
Original abstract

We propose nonparametric methods for functional linear regression which are designed for sparse longitudinal data, where both the predictor and response are functions of a covariate such as time. Predictor and response processes have smooth random trajectories, and the data consist of a small number of noisy repeated measurements made at irregular times for a sample of subjects. In longitudinal studies, the number of repeated measurements per subject is often small and may be modeled as a discrete random number and, accordingly, only a finite and asymptotically nonincreasing number of measurements are available for each subject or experimental unit. We propose a functional regression approach for this situation, using functional principal component analysis, where we estimate the functional principal component scores through conditional expectations. This allows the prediction of an unobserved response trajectory from sparse measurements of a predictor trajectory. The resulting technique is flexible and allows for different patterns regarding the timing of the measurements obtained for predictor and response trajectories. Asymptotic properties for a sample of $n$ subjects are investigated under mild conditions, as $n\to \infty$, and we obtain consistent estimation for the regression function. Besides convergence results for the components of functional linear regression, such as the regression parameter function, we construct asymptotic pointwise confidence bands for the predicted trajectories. A functional coefficient of determination as a measure of the variance explained by the functional regression model is introduced, extending the standard $R^2$ to the functional case. The proposed methods are illustrated with a simulation study, longitudinal primary biliary liver cirrhosis data and an analysis of the longitudinal relationship between blood pressure and body mass index.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 7,654 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.31 · human participants 0.90 · controlled comparison 0.02 · randomized assignment 0.01 · states quantitative result 0.17 · compares with prior results 0.08 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We propose nonparametric methods for functional linear regression which are designed for sparse longitudinal data, where both the predictor and response are functions of a covariate such as time. Predictor and response processes have smooth random trajectories, and the data consist of a small number of noisy repeated measurements made at irregular times for a sample of subjects. In longitudinal studies, the number of repeated measurements per subject is often small and may be modeled as a discrete random number and, accordingly, only a finite and asymptotically nonincreasing number of measurements are available for each subject or experimental unit. We propose a functional regression approach for this situation, using functional principal component analysis, where we estimate the functional principal component scores through conditional expectations. This allows the prediction of an unobserved response trajectory from sparse measurements of a predictor trajectory. The resulting technique is flexible and allows for different patterns regarding the timing of the measurements obtained for predictor and response trajectories. Asymptotic properties for a sample of $n$ subjects are investigated under mild conditions, as $n\to \infty$, and we obtain consistent estimation for the regression function. Besides convergence results for the components of functional linear regression, such as the regression parameter function, we construct asymptotic pointwise confidence bands for the predicted trajectories. A functional coefficient of determination as a measure of the variance explained by the functional regression model is introduced, extending the standard $R^2$ to the functional case. The proposed methods are illustrated with a simulation study, longitudinal primary biliary liver cirrhosis data and an analysis of the longitudinal relationship between blood pressure and body mass index.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We propose nonparametric methods for functional linear regression which are designed for sparse longitudinal data, where both the predictor and response are functions of a covariate such as time. contribution 1.00 0.91 0.14 asserted 0.99
s1 Predictor and response processes have smooth random trajectories, and the data consist of a small number of noisy repeated measurements made at irregular times for a sample of subjects. background 0.93 0.21 0.03 asserted 1.00
s2 In longitudinal studies, the number of repeated measurements per subject is often small and may be modeled as a discrete random number and, accordingly, only a finite and asymptotically nonincreasing number of measurements are available for each subject or experimental unit. background 1.00 0.10 0.03 hedged possibility 0.57
s3 We propose a functional regression approach for this situation, using functional principal component analysis, where we estimate the functional principal component scores through conditional expectations. method 0.90 0.93 0.09 asserted 1.00
s4 This allows the prediction of an unobserved response trajectory from sparse measurements of a predictor trajectory. interpretation 0.36 0.77 0.14 asserted 0.97
s5 The resulting technique is flexible and allows for different patterns regarding the timing of the measurements obtained for predictor and response trajectories. contribution 0.67 0.72 0.05 asserted 0.99
s6 Asymptotic properties for a sample of $n$ subjects are investigated under mild conditions, as $n\to \infty$, and we obtain consistent estimation for the regression function. result 0.85 0.94 0.26 asserted 1.00
s7 Besides convergence results for the components of functional linear regression, such as the regression parameter function, we construct asymptotic pointwise confidence bands for the predicted trajectories. contribution 0.72 0.96 0.13 asserted 1.00
s8 A functional coefficient of determination as a measure of the variance explained by the functional regression model is introduced, extending the standard $R^2$ to the functional case. contribution 0.97 0.94 0.07 asserted 1.00
s9 The proposed methods are illustrated with a simulation study, longitudinal primary biliary liver cirrhosis data and an analysis of the longitudinal relationship between blood pressure and body mass index. contribution 0.55 0.93 0.11 asserted 1.00

What kind of paper the classifier reports this is: method or tool 1.00 · evidence from proof or derivation 0.42.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Left panel: Observed individual trajectories (solid ) and the smooth estimate of the mean function for albumin (thick solid ). Right panel: Corresponding observed individual trajectories (solid ) and the smooth estimate of the mean function for prothrombin time (thick solid ), for the primary biliary cirrhosis (PBC) data
Figure 2. Left panel: Smooth estimates of the first (solid ) and second (dashed ) eigenfunc￾tions for albumin, accounting for 87% and 8% of total variation. Right panel: Smooth estimates of the first (solid ) and second (dashed ) eigenfunctions for prothrombin time, accounting for 54% and 33% of total variation, for the PBC data
Figure 3. Estimated regression function (8), where the predictor (albumin) time is s (in days), and the response ( prothrombin) time is t (in days), for the PBC data
Figure 4. Observed values (circles) for prothrombin times (not used for prediction), pre￾dicted curves (solid ) and 95% pointwise bands (dashed ), for four randomly selected pa￾tients, where bands and predicted curves are based on one-curve-leave-out analysis and sparse noisy albumin measurements (not shown) of the predictor function
Figure 5. Estimated pointwise coefficient of determination Rb 2 (t) (19) as a function of time t for the PBC data (left panel, time is duration of study in days) and for the BLSA data (right panel, time is age in years). Research Center bi-annually. Time corresponds to age of each subject and is measured in years. On each visit, systolic blood pressure (in mm Hg) and body mass index (BMI = weight in kg/height in m2 ) were assessed. Since both measurements are highly variable, the data are noisy, and as many par￾ticipants missed scheduled visits, or were seen at other than the scheduled times,
Figure 7. Left panel: Smooth estimates of the first (solid ), second (dashed ) and third (dash-dot) eigenfunctions for BMI, accounting for 90%, 6% and 3% of total variation. Right panel: Smooth estimates of the first (solid ), second (dashed ), third (dash-dot) and fourth (dotted ) eigenfunctions for SPB, accounting for 76%, 15%, 4% and 2% of total variation, for the BLSA data
Figure 8. Estimated regression function (8), where the predictor (BMI) time is s (in years), and the response (SBP) time is t (in years), for the BLSA data. then, as n tends to infinity, remainder terms such as θn and ϑn in (41) will disappear. Appropriately modified versions of the following theoretical results still hold. Since in this case the true K and M are finite, most of the technical assumptions that are needed to handle the infinite case would not be needed, such as (A1)–(A3), (A6)–(A7) and (B5), and K(n),M(n) would be assumed to converge to the true K and M instead of infinity, as…
Figure 9. Observed values (circles) of systolic blood pressure (SBP) that are not used for the prediction, predicted curves (solid ) and 95% pointwise bands (dashed ), for four randomly selected patients, where bands and predicted curves are based solely on sparse and noisy measurements of the predictor function (BMI, not shown). (B5) in Appendix A.1. For simplicity, we suppress the dependency of M and K on n in the notation. The consistency of βˆ (8) is obtained as follows. Theorem 1. Under (A1) and the assumptions of Lemma A.1 and (B5) (see Appendix A.1), limn→∞ Z T Z S [βˆ(s,t) − β(s,t)]

Showing 8 of the paper's 9 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:math/0603132 by the authors.

References
How the literature uses it

1 arXiv paper cite this in 1 sentence

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it

Cites it as background 1 paper

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.