For any real orthonormal basis , the argument of (b)(ii) gives almost sure convergence and convergence in L2 of . Finite linear combinations have normal distributions, and the Parseval identity for a Hilbertian basis gives
Thus and are centered Gaussian processes with the same covariance function. More explicitly, their finite-dimensional characteristic functions are both
This proves the basis invariance of an isonormal Gaussian process: equality of all finite-dimensional distributions, and hence equality in law of the indexed collections.
For the specified cosine orthonormal basis, put . Its coefficient against an interval indicator is
Consequently the Brownian half-integer sine expansion, a Karhunen–Loève expansion, is
For every fixed , the sum converges almost surely and in L2 space because the coefficient squares sum to . Its continuous modification is a standard Brownian motion by (b)(v). The expansion also holds in : the expectation of the squared tail is the integral of the coefficient-square tail, which tends to zero by the dominated convergence theorem, with bound .
To realize this expansion for the original , choose . These are Jointly Gaussian with covariance matrix equal to the identity, hence are independent standard normal random variables. The finite orthogonal projections of onto converge in , and the identity transfers that convergence to the displayed expansion for . With an independently chosen sequence the resulting Brownian motion has the same law; a pathwise identification uses the coupling just given.
Let . The Karhunen–Loève expansion and Gaussianity give
and
The conditional multivariate normal distribution formula yields
If is singular, the same formula uses its Moore-Penrose pseudoinverse.
The Karhunen–Loève expansion has the following Hilbert-space form. Let be a square-integrable Random element of a Hilbert space , with mean and covariance operator . There are nonnegative eigenvalues and orthonormal eigenvectors spanning the closure of the range of such that
in . The functional principal component scores
satisfy
If is a Gaussian random element, the scores are independent normal random variables.
For the proof, is positive, self-adjoint, and trace-class, with
It is therefore compact, so the spectral theorem for compact Hermitian operators supplies the eigenpairs. Their score covariance is
For the residual , Parseval identity and the trace formula give
The centered variable's projection onto has zero second moment and is therefore zero almost surely, which completes the mean-square expansion. When has a continuous covariance kernel, Mercer's theorem additionally expands that kernel as .
The stated kernel is the Brownian bridge covariance kernel. Its eigenvalue equation is
The right-hand side vanishes at and , and differentiating it twice gives
Thus the normalized eigenfunctions and eigenvalues are
The Karhunen–Loève expansion is consequently
with convergence in , where and . Covariance alone does not imply that the are independent or normal; they are independent standard normal variables when is Gaussian.
Use the test statistic
Under the null hypothesis, the Hilbert-space central limit theorem gives , where is centered Gaussian with covariance . If , its Karhunen–Loève expansion and the continuous mapping theorem give
for independent . Reject for above the quantile of this weighted chi-squared law; replacing the by empirical covariance eigenvalues gives a plug-in estimator of the critical value.
Under every fixed alternative , the weak law of large numbers gives , so and the test is consistent. Under a local alternative , the limit is , which describes its local power.
Let be the Hilbert-Schmidt operator with kernel , so the function-on-function linear model is . Independence and centering give . Let
and define the functional principal component scores
The cross-covariance operator identity gives, for every ,
Since is centered, the requested integrated variance is . Applying the Karhunen–Loève expansion to and the Parseval identity in the basis yields
Equivalently, if , the expression is .