FPCA mean test Created 2026-09-24 Updated 2026-11-05
An FPCA mean test uses functional principal component analysis to test a zero expectation in a Hilbert space. For independent and identically distributed observations with finite second moment, let be the eigenvalues and orthonormal eigenvectors of the covariance operator. Choose a fixed positive integer and, as a sufficient separation assumption, require and . With consistent estimators of these eigenvalues and eigenvectors, the statistic is
Under the null hypothesis , the central limit theorem gives joint convergence of to independent standard normal random variables. Statistical consistency and Slutsky's theorem justify replacing the population eigenvectors and eigenvalues by their estimates. Consequently . The fixed-dimensional projected two-sample version has the same calibration.[ref].
The positivity condition prevents division by zero; the spectral gap makes the estimated directions stable. Repeated eigenvalues can be handled through an entire retained eigenspace, provided its boundary remains separated. A choice growing with needs additional approximation theory. A fixed can miss a nonzero expectation lying entirely in the omitted orthogonal directions, so the test's simplicity comes with a choice about which deviations matter.
Functional mean test Created 2026-09-24 Updated 2026-11-05
A functional mean test asks whether the expectation of a Hilbert-space-valued random variable is zero. Suppose are independent and identically distributed in a Hilbert space with a countable dense subset , with . Write for their covariance operator and . Under the null hypothesis , the Hilbert-space central limit theorem gives
where is a centered Gaussian random element in , the are the eigenvalues of , and the are independent standard normal random variables. The second assertion follows by the continuous mapping theorem. Since , the nonnegative series is finite almost surely. Limits given by Gaussian random elements for functional sample means and the corresponding weighted squared-norm laws also extend to suitable dependent observations.[ref]
The eigenvalues determine the calibration: the limit is generally a weighted sum of random variables with a chi-squared distribution, rather than a single ordinary chi-squared distribution. The squared norm measures deviations across all covariance operator directions. An FPCA mean test instead retains finitely many eigenvectors and divides their squared coordinates by their eigenvalues, producing a different statistic with a simpler limiting calibration.
Let be the leading eigenpairs of the sample covariance operator. For fixed with and , the FPCA mean test uses
Under the null, consistency of the empirical eigenpairs and the multivariate central limit theorem imply
The level- test therefore rejects above the quantile of the chi-squared distribution with degrees of freedom.
For a fixed mean , if at least one leading coordinate , , is nonzero, then in probability and the test is consistent. It has only null-level asymptotic power against means orthogonal to the first principal component functions. Under , the limit is noncentral chi-squared with noncentrality
The squared-norm test is omnibus: every fixed nonzero mean eventually changes . Its null law, however, is an infinite weighted chi-squared distribution and requires accurate estimation of enough covariance eigenvalues; noisy low-variance directions can also make calibration inefficient.
The FPCA mean test has the simple limit and standardizes retained directions by their variances. It is effective when the signal lies in the leading principal component subspace, but choosing introduces a tuning decision and truncation makes the test blind to alternatives orthogonal to that subspace. Close or repeated eigenvalues also make individual empirical eigenfunctions unstable.
The sign-flip randomization test can provide finite-sample calibration and avoids estimating a limiting covariance spectrum. Its exactness requires central symmetry, which is stronger than merely having zero mean, and exhaustive enumeration costs evaluations; Monte Carlo sign flips introduce simulation error. Its power still depends on the statistic used inside the randomization scheme.