Singular-value decompositions of the nonparametric two-way regression kernel produce eigenfunction-based proxies that identify unit and time fixed effects: injectivity holds almost everywhere under mild nondegeneracy and holds uniformly (for a finite truncation) under a local full-rank/immersion condition, enabling identification of the regression function at the realised types.

Main Finding

The paper shows that the latent, two-way unobserved heterogeneity in the nonparametric panel regression Git = g(αi, γt) + εit can be identified (up to observational equivalence) by using the singular value decomposition (SVD) of the bivariate regression kernel g. The left/right singular functions (ur, vr) serve as injective, identified proxies for the cross-sectional and time heterogeneity respectively. Two identification regimes are established: (i) an almost-everywhere regime where a finite collection of leading eigenfunctions becomes injective on a set of probability arbitrarily close to 1, and (ii) a uniform regime where under a stronger local full-rank condition a finite number of leading eigenfunctions are globally injective on the whole support.

Key Points

• Model & goal

  • Observed transformed outcomes Git = G(Yit) with conditional mean g(αi, γt).
  • Aim: identify realised fixed effects (αi, γt) or injective proxies for them and identify g at those realised types under large N, large T.
  • Observational-equivalence simplification: distinguishability is only required up to values of α that induce distinct functions g(α, ·).

• SVD-based proxy construction

  • Treat g as a Hilbert–Schmidt kernel on Ωα × Ωγ and form the integral operator Tg.
  • SVD: g(α, γ) = Σr σr ur(α) vr(γ) with left singular functions ur ∈ Hα (RdG-valued) and right singular functions vr ∈ L2π(Ωγ).
  • Define eigenfunction maps UR(α) = (u1(α),...,uR(α))⊤ and VR(γ) similarly. These serve as proxies λi = UR(αi) and ft = VR(γt).

• Two identification regimes (Theorem 2)

  • Almost-everywhere regime (weaker assumption): If the gradient-information matrix M∞(α) := ∫ J(α,γ)⊤J(α,γ) dπγ is positive definite for almost every α, then for any δ>0 there exists finite R(δ) and large-probability compact sets on which UR and VR are injective for all R ≥ R(δ). Thus proxies can be made injective on sets of probability arbitrarily close to one, at the cost of increasing R.
  • Uniform regime (stronger assumption): If a finite “design” of time points yields a stacked Jacobian with full column rank uniformly in α (a local immersion that holds uniformly), then there exists finite R such that UR and VR are globally injective on the whole supports for every R ≥ R.

• Key technical ingredients

  • Use gradient expansion of J(α,·) = Dαg(α,·) in the vr basis and define MR(α) = Σ_{r≤R} σr^2 Dur(α)⊤Dur(α), with M∞ its limit.
  • Show M∞(α) equals ∫ J⊤J and relate positive-definiteness of M∞ to local injectivity of the infinite eigenfunction map U∞.
  • Reduce global injectivity of finite UR to a local first-order (immersion) condition plus topological compactness; avoids global Jacobian (univalence) conditions.

• Additional points

  • Theorem 1: If identified finite-dimensional proxies (λi, ft) are injective maps of (αi, γt), then g at realised types is identified and inherits regularity under smoothness of U,V,g.
  • The paper treats both α- and γ-sides symmetrically (γ-analogues required).
  • Additive structure in the kernel across coordinates of the types yields componentwise decomposition that can remove the curse of dimensionality (Section 5.8).

Data & Methods

• Mathematical / theoretical approach (no new empirical dataset)
  • Functional-analytic framework: Hilbert–Schmidt integral operators on L2 spaces, SVD of compact operators, orthonormal systems of singular functions.
  • Sobolev regularity: eigenfunction smoothness established under Hp assumptions (Appendix A); use of weighted/unweighted Sobolev equivalence given densities bounded away from 0 and ∞.
  • Differential calculations: Jacobian expansions, Bessel and dominated convergence arguments show convergence of differentiated SVD series and closed form of M∞.
  • Topological/compactness arguments: reduce global injectivity to local immersion plus compactness to get finite truncation injectivity.
• Main assumptions
  • Support assumptions: compact (and for α, convex and connected) supports; densities bounded away from 0 and ∞.
  • Smoothness: g continuously differentiable in α with joint continuity of Dαg on Ωα×Ωγ (so ur are C1).
  • Observational equivalence: the response map α ↦ g(α, ·) is injective after removing observationally equivalent points.
  • Assumption 4 (weaker): det M∞(α) = 0 only on a πα-null set (M∞ ≻ 0 a.e.).
  • Assumption 5 (stronger): existence of a finite list of evaluation points in Ωγ so that a stacked finite Jacobian has full column rank uniformly in α (gives uniform positive-definiteness of M∞).
• Results are identification proofs; estimation aspects (finite-sample error, rates) are discussed qualitatively (e.g., how truncation order R must grow) but detailed finite-sample theory/algorithms are not developed in this paper.

Implications for AI Economics

• What it enables
  • Rigorous identification of latent heterogeneous agent/time types in rich nonparametric panel models when heterogeneity affects conditional responses (i.e., when different types generate different conditional distributions of transformed outcomes).
  • A constructive route to build low-dimensional, injective proxies for latent types by spectral decomposition of the estimated regression kernel. Those proxies can be used as sufficient controls in downstream structural or causal analysis.
• Practical consequences for empirical work in AI economics
  • When studying heterogeneous effects of algorithmic policies, platform changes, or firm-level AI adoption over time, researchers can:
    • Estimate the conditional mean g nonparametrically (or its moments G(Yit)),
    • Compute an empirical SVD of the estimated kernel,
    • Use leading singular-function evaluations as proxies for cross-sectional/time types.
  • The paper gives conditions under which a finite number of leading components suffice to identify types (so using low-rank approximations is justified), and it clarifies the trade-off: under weak assumptions one can achieve injectivity on almost all types by letting the number of components R grow (possibly slowly), while under a stronger local full-rank condition a fixed finite R is enough globally.
• Guidance on choosing truncation R and interpretation
  • Under the almost-everywhere regime, R must grow with sample sizes to drive down the measure and diameter of the set of conflated types; in examples with bounded type density this can be logarithmic in NT, but no universal rate holds.
  • Under the uniform regime, a finite R* suffices; verifying the finite-design full-rank condition in applications amounts to checking that variation in g across some finite set of γ points spans the α-derivatives.
• Limitations and cautions for applied researchers
  • Identification is up to observational equivalence: types that produce exactly identical conditional responses cannot be distinguished by any procedure based on g; researchers should check whether types of interest can plausibly be distinguished via data transformations G(·).
  • The theory relies on smoothness, compact-support, and density bounds; violations may break eigenfunction regularity and injectivity conclusions.
  • The paper proves identification; it does not provide a full finite-sample estimation and inference recipe that quantifies sampling error from nonparametric estimation of g and of the operator SVD. Implementing the approach requires careful nonparametric estimation and numerical spectral analysis, with attention to discretisation and regularisation.
• Directions for empirical AI-economics research
  • Use eigenfunction proxies as controls in causal or structural estimators to account for two-way unobserved heterogeneity (e.g., heterogeneous response to algorithmic recommendations across user types and time regimes).
  • Exploit additive or separable kernel structure (when present) to mitigate curse-of-dimensionality concerns, per the componentwise decomposition in the paper.
  • Combine with asymptotic results for nonparametric estimators of g to develop finite-sample procedures (model selection for R, regularisation) that inherit the identification guarantees here.

Summary recommendation: the paper provides a clean, constructive identification strategy for two-way latent heterogeneity via kernel SVD and gives clear conditions (both weak and strong) under which a finite spectral truncation yields valid, interpretable proxies. Applied researchers in AI economics can leverage these proxies to control for rich unobserved heterogeneity, but must take care when estimating g nonparametrically and when verifying the assumptions required for the desired injectivity regime.