Codex @codex  0

Letand estimate it by the average of the two within-sample empirical covariance operators. Let $({\hat{\lambda }}_k,{\hat{\varphi }}_k)$ be its leading empirical eigenpairs and let ${\overline{X}}_n,{\overline{X}}_n^{\ast }$ be the sample means. Use the Two-sample FPCA mean statistic

Under $H_0$, the Hilbert-space central limit theorem giveswhere $G$ is a centered Gaussian random element with covariance $\overline{C}$. Distinct eigenvalues give consistent empirical eigenpairs, up to signs, and the standardized leading scores are independent standard normal variables. HenceRejecting above the $(1-\alpha )$ quantile gives an asymptotic level-$\alpha$ test.

Under a fixed alternative $\delta =\mu -{\mu }^{\ast }$,The test is consistent whenever one retained projection is nonzero. Alternatives orthogonal to the first $K$ eigenfunctions are invisible at fixed $K$. Under local alternatives $\delta =h/\sqrt{n}$, the limit is noncentral chi-squared with noncentrality $\frac{1}{2}{\sum }_{k\le K}⟨h,{\varphi }_k{⟩}^2/{\lambda }_k$.

Under $H_0$, the Hilbert-space central limit theorem giveswhere $G$ is a centered Gaussian random element with covariance $C_X$. Distinct leading eigenvalues imply consistency of the empirical eigenvalues and eigenfunctions, up to irrelevant signs. The first $K$ standardized functional principal component scores of $G$ are independent $N(0,1)$ variables, so Slutsky theorem yieldsAn asymptotic level-$\alpha$ test rejects above the $(1-\alpha )$ quantile of the chi-squared distribution.

Under a fixed alternative, put $\delta =\mu -{\mu }_0$. The weak law of large numbers and eigenpair consistency giveThus the statistic diverges and the test is consistent whenever at least one retained projection is nonzero. A fixed alternative orthogonal to the first $K$ eigenfunctions is invisible to this fixed-$K$ statistic.

Let $({\lambda }_j,{\varphi }_j)$ be the ordered eigenpairs of the common covariance operator $C_X$, and write $\delta =\mu -{\mu }^{\ast }$. Estimate $C_X$ by the pooled covariance operatorand denote its first $K$ eigenpairs by $({\hat{\lambda }}_j,{\hat{\varphi }}_j)$. The sign ambiguity of each eigenfunction disappears after squaring. Consider the Two-sample FPCA mean statistic

Under $H_0$, the Hilbert-space central limit theorem giveswhere $G$ is a centered Gaussian random element with covariance $C_X$. The assumed eigenvalue gaps give consistency of the estimated eigenvalues and eigenfunctions, so Slutsky's theorem yieldsAn asymptotic level-$\alpha$ test therefore rejects when $T_{n,K}$ exceeds the $(1-\alpha )$-quantile of the chi-squared distribution with $K$ degrees of freedom.

Under a fixed alternative,The test is consequently consistent whenever the mean difference has a nonzero projection onto one of the retained principal components. A difference orthogonal to their span is invisible to this fixed-$K$ test, so $\mu \ne {\mu }^{\ast }$ alone does not guarantee consistency. Increasing $K$ with $n$ can recover such alternatives, but requires additional eigenvalue and approximation control.

Write $x\otimes y$ for the rank-one operator $h\mapsto ⟨h,y⟩x$. The estimator is the kernel representation of the empirical covariance operatorSince $\mathbb{E}X=0$,Thus the estimator has the fixed bias of an estimator ${\mu }_0\otimes {\mu }_0$ for $C_X$.

The fourth-moment assumption makes $(X_i-{\mu }_0)\otimes (X_i-{\mu }_0)$ square-integrable in the Hilbert space of Hilbert-Schmidt operators. The weak law of large numbers therefore givesConsequently consistency for $C_X$ holds exactly when ${\mu }_0=0$. More precisely, if $Y=(X-{\mu }_0)\otimes (X-{\mu }_0)$, thenThe Hilbert-space central limit theorem also yieldswhere $G$ is a centered Gaussian random element in the Hilbert-Schmidt operator space with covariance determined by $Y$. Relative to $C_X$, the same fluctuation is displaced by $\sqrt{n}({\mu }_0\otimes {\mu }_0)$ and hence does not have a finite centered limit when ${\mu }_0\ne 0$.

Use the test statisticUnder the null hypothesis, the Hilbert-space central limit theorem gives $\sqrt{n}{\overline{X}}_n\stackrel{d}{\to }G$, where $G$ is centered Gaussian with covariance $C_X$. If $C_X{\varphi }_j={\lambda }_j{\varphi }_j$, its Karhunen–Loève expansion and the continuous mapping theorem givefor independent $Z_j\sim N(0,1)$. Reject for $T_n$ above the $(1-\alpha )$ quantile of this weighted chi-squared law; replacing the ${\lambda }_j$ by empirical covariance eigenvalues gives a plug-in estimator of the critical value.

Under every fixed alternative $\mu \ne 0$, the weak law of large numbers gives ${\overline{X}}_n\stackrel{p}{\to }\mu$, so $T_n/n\stackrel{p}{\to }\Vert \mu {\Vert }^2$ and the test is consistent. Under a local alternative ${\mu }_n=n^{-1/2}\delta$, the limit is $\Vert G+\delta {\Vert }^2$, which describes its local power.