Pith. sign in

REVIEW 5 major objections 6 minor 40 references

Modeling large dimensional matrix time series with partially known and latent factors

T0 review · 5 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read Adding observable covariates to a matrix factor model need not slow the estimation of latent factors.

desk verdict A straightforward extension of matrix factor models to include known covariates, with plausible rates but no proofs and an untested exogeneity assumption that is likely violated in the real-data application. read the letter →

arxiv 2411.16192 v1 pith:XPACPMW2 submitted 2024-11-25 stat.ME

classification stat.ME MSC 62H2562M1062J05
keywords matrixtimeseriesfactormodelpartiallyknownfactorslatentleastsquaresestimationhigh-dimensionalconvergenceratestwo-step
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper proposes a matrix time series model that combines a regression term for observable covariates with a latent matrix factor model: $Y_t = A X_t + R F_t C' + E_t$. The central aim is to show that adding the known-factor term does not come at a statistical cost: after estimating $A$ by least squares and applying the standard matrix-factor estimation procedure to the residuals, the loading matrices converge at exactly the same rates as in the model without covariates. This matters because observable pricing factors are natural explanatory variables in finance, and a model that uses both known and latent factors is more interpretable without sacrificing accuracy for the latent structure. The paper also gives convergence rates for the coefficient matrix, the estimated signal, and eigenvalue-based dimension estimators, alongside simulations and a stock-return application.

What carries the argument

The workhorse is the lagged auto-covariance identity for the residual process. With $W_t = R F_t C' + E_t = Q_1 Z_t Q_2' + E_t$ and white-noise errors, the cross-covariance of columns of $W_t$ satisfies $\Omega_{w,ij}(h) = Q_1 \Omega_{zq,ij}(h) Q_1'$ for $h \ge 1$, so the statistic $M_1 = \sum_{h,i,j} \Omega_{w,ij}(h) \Omega_{w,ij}(h)'$ equals $Q_1$ times a sum of factor autocovariances times $Q_1'$. Consequently the eigenspace of $M_1$ coincides with the column space of the row loading matrix $Q_1$, and the analogous construction on transposed residuals gives $Q_2$. The QR decomposition first writes the loadings as semi-orthogonal matrices times nonsingular matrices, making the column spaces identifiable despite the rotational indeterminacy of $R F_t C'$.

What would settle it

Simulate from $Y_t = A X_t + R F_t C' + E_t$ with $X_t$ generated partly from $F_t$ (say $X_t = c F_t + \text{noise}$ with $c>0$), which violates the uncorrelatedness assumption; if the two-step estimator's loading error $\|\hat Q_i - Q_i\|_2$ stops shrinking at the claimed rate as $c$ grows, the assumption is doing the work.

Watch

Extended reading notes

Core claim

The central claim is that the two-step estimator, least squares for the coefficient matrix $A$ followed by eigen-analysis of a lagged cross-covariance statistic built from the residuals, recovers the latent factor loading spaces with error $O_p(p^{\delta_1} q^{\delta_2} T^{-1/2})$, identical to the rate proven for the matrix factor model with no known-factor term. The same rate holds for both the row and column loading spaces, and the eigenvalue ratios used to choose the number of latent factors remain valid. In other words, the regression term $A X_t$ is absorbed in the first step, and the residual process $W_t = R F_t C' + E_t$ behaves like the pure matrix factor model, so the presence of observable covariates does not degrade estimation of the latent factors.

Load-bearing premise

The entire argument depends on the observable covariates $X_t$ being uncorrelated with the latent factors $F_t$ and with the noise $E_t$, so that ordinary least squares isolates $A$ and leaves residuals that are exactly a matrix factor model plus white noise.

Editorial extensions

If this is right

  • Known factors can be added to a matrix factor model without slowing down latent-factor estimation: the loading error rate is the same whether or not the term $A X_t$ is present.
  • The least-squares estimator of the coefficient matrix is consistent when $p/T \to 0$, with rate $\|\hat A - A\|_F = O_p(p^{1/2} T^{-1/2})$.
  • The ratio-based eigenvalue estimator for the number of latent factors remains theoretically valid, though weak factors make correct dimension selection harder in finite samples, as the simulations show.
  • The signal part $A X_t + R F_t C'$ is consistently estimable as $p$ and $q$ grow, so the fitted model can be used for forecasting and out-of-sample evaluation.
  • In the stock-return application, the known-factor term accounts for a substantial share of explained variation, and the latent column factor has near-zero loading on the U.S. market once known factors are included.

Reading between the lines

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

  • Beyond the paper, the orthogonality assumption suggests a natural stress test: if $X_t$ is allowed to correlate with $F_t$ or $E_t$, the residual $W_t$ is no longer a valid matrix factor model, and one would expect the claimed rates to fail; extending the model to endogenous covariates would require instrumental-variable or control-function ideas.
  • The same two-step structure could be adapted to the two-sided regression form $A X_t B'$ with observable matrix covariates, at the cost of a bilinear rather than linear first step; the paper notes this form is more natural but does not pursue it.
  • The rate-matching result hints that latent-factor estimation is first-order unaffected by the regression term, so any efficiency loss from estimating $A$ is asymptotically negligible; a finite-sample comparison of standard errors could quantify this.
  • Because the real-data column loading on the U.S. is near zero after including known factors, one testable implication is that for markets with strong observable pricing factors, latent column structure may be superfluous; this could be checked by fitting the model with and without the known-factor term on other multi-market panels.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 6 minor

Summary. The paper proposes a matrix-variate time series factor model with both known covariates and latent factors, Y_t = A X_t + R F_t C' + E_t. The estimation is done in two steps: ordinary least squares for the coefficient matrix A, followed by the Wang et al. (2019) eigenanalysis method applied to the residuals to estimate the latent loading spaces and the number of factors. The authors claim convergence rates for A, the loadings, and the signal part that coincide with those of Wang et al. (2019), and they support the method with simulations and a real-data application using Fama-French factors as known covariates for daily stock returns.

Significance. If the theoretical claims are correct, the paper offers a natural and useful extension of high-dimensional matrix factor models by incorporating observable covariates, and the preserved convergence rates for the latent part would be an interesting result. The paper also provides numerical experiments and a real-data demonstration. However, the central theoretical results are stated without proofs, the key exogeneity assumption is not formally listed or tested, and there are internal inconsistencies in the conditions. These issues currently prevent verification of the paper's main claims, so the significance cannot be assessed from the manuscript as it stands.

major comments (5)
  1. [Section 4] Theorems 1-5 are stated without proofs, and no supplementary file or appendix is provided. This is a fundamental gap for a theory paper: the central claim that the convergence rates coincide with Wang et al. (2019) cannot be verified. The authors should provide complete proofs (or a detailed proof sketch with all steps) for these results.
  2. [Section 3] The exogeneity assumption Cov(vec(X_t), vec(F_t)) = 0 and Cov(vec(X_t), vec(E_t)) = 0 is used to justify the least squares estimator of A, but it is not listed among Conditions 1-7 in Section 4 and is not tested in the simulations or real data. If this assumption fails, \hat A is biased, the residual \hat W_t contains latent factor variation, and the rate in Theorem 2 is unsupported. In the real data application, using Fama-French factors as X_t for individual stock returns makes correlation with latent factors plausible. Please state this assumption formally and provide diagnostic evidence or a sensitivity analysis.
  3. [Section 4, Condition 4] Condition 4 states that E(X_t X_t') = P_t with (1/T)Σ P_t → P and λ_min(P) ≍ q. Since X_t is m × q with m fixed, X_t X_t' is m × m, so its eigenvalues cannot grow as q → ∞. This internal inconsistency also appears in Theorem 1, whose rate O_p(p^{1/2}T^{-1/2}) does not depend on q. Please clarify the intended normalization of the covariate matrix.
  4. [Section 5, Table 3] For δ1 = δ2 = 0.5, the relative frequency of correctly estimating the number of factors is 0.0 in nearly all configurations, even at T = 2pq. The theoretical condition p^{δ1} q^{δ2} T^{-1/2} = o(1) is not satisfied in these settings, which the paper does not mention. The authors should either present simulation settings where the theoretical condition holds and the method works, or explicitly discuss this limitation.
  5. [Section 4, Condition 3] Condition 3 states that the rank of Σ_f(h) is k* = max(k, r), but Σ_f(h) is the covariance matrix of vec(F_t), which has dimension kr × kr. The condition is thus unclear as written, and the subsequent conditions on individual row/column covariances are not connected to identifiability of the loading spaces. Please rewrite this condition precisely.
minor comments (6)
  1. [Abstract and Introduction] There are several grammatical errors, e.g., 'This article considers to model' should be 'This article considers modeling', and 'provides a extensive framework' should be 'provides an extensive framework'.
  2. [Section 3, equations (7)-(8)] The notation Ω_{zq,ij}(h) is used before Z_t has been properly defined, and the argument of Cov in (7) is difficult to parse. Please define Z_t and clarify the covariance arguments.
  3. [Section 6] The text says the estimated dimensions are k = 1 and r = 1 but then states 'we use k = 2 and r = 1 here for accessible illustration.' This is confusing and should be explained.
  4. [Section 6, Table 9] Table 9 is formatted ambiguously: the columns for R2_K, R2_U, R2_T, and the matrix factor model R2 are not clearly separated, and for the (2,1) row the values appear inconsistent (R2_U exceeds R2_T). Please fix the table layout.
  5. [References] Some references appear unrelated to the topic (e.g., Psychogios et al. 2012, Sahay 2005, Tian et al. 2015); these should be removed or replaced with relevant literature.
  6. [Section 6, equation (18)] The out-of-sample R2 notation is not fully defined: it is unclear how the test set is constructed, whether parameters are re-estimated on the training set, and how the fitted values are computed. Please clarify.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the estimation chain is a standard two-step least-squares plus PCA procedure benchmarked against an external result.

full rationale

The paper's central claim is that introducing the known-factor term A X_t does not alter the Wang et al. (2019) convergence rates for the estimated loading spaces. This is a comparison with an external benchmark (Wang et al. 2019), not a quantity derived from its own output. The estimation procedure is: OLS for A (Theorem 1), residual W-hat_t = Y_t - A-hat X_t, then eigenvalues of M-hat_1 built from autocovariances of W-hat_t, exactly the Wang et al. construction. The population identity Omega_{w,ij}(h) = Q_1 Omega_{zq,ij}(h) Q_1' follows from model (1) and the white-noise assumption on E_t; it is an identification argument, not a self-referential prediction. No parameter is fitted to the quantity later reported as a prediction; simulation data are generated from (16) independently and real data are separate. The cited identification uniqueness (Lam et al. 2011; Wang et al. 2019) is external and not by the present authors, so no self-citation chain is load-bearing. A genuine weakness exists: Theorem 2 is stated without proof, and the exogeneity assumptions Cov(vec(X_t), vec(F_t))=0 and Cov(vec(X_t), vec(E_t))=0 are asserted rather than tested, so the rate claim may fail in the Fama-French application; but this is a correctness/robustness risk, not circularity, because the assumptions are stated inputs rather than outputs of the derivation.

Assumptions & free parameters 2 free parameters · 7 assumptions · 0 invented entities

The model introduces no new latent objects beyond the standard matrix factor structure. The main input assumptions are the seven technical conditions and the exogeneity condition. The most fragile assumption is the zero-correlation between the observable covariates and the latent factors or errors, which is not tested in the paper.

free parameters (2)
  • h0 = 1
    Tuning parameter for the M1 and M2 statistics in Section 3; set to 1 in the simulations and real data analysis.
  • K = p/2 or p/3
    Upper range for the ratio-based factor number estimator in equation (14); chosen by the practitioner.
assumptions (7)
  • domain assumption The joint process (vec(Y_t), vec(X_t), vec(F_t)) is alpha-mixing with coefficients satisfying sum alpha(k)^{1-2/gamma} < infinity for some gamma > 2.
    Condition 1, Section 4: imposes temporal dependence structure needed for the concentration results.
  • domain assumption Bounded 2-gamma moments for the elements of F_t, X_t, and E_t.
    Condition 2, Section 4: tail conditions for the weak convergence arguments.
  • domain assumption The matrix Sigma_f(h) has rank max(k, r) and non-zero autocovariances for rows and columns of F_t.
    Condition 3, Section 4: ensures the latent factors are identifiable from the autocovariance structure.
  • domain assumption The limit of (1/T) sum E(X_t X_t') has minimum eigenvalue of order q.
    Condition 4, Section 4: requires the observable covariates to be uniformly strong across the q columns of X_t.
  • domain assumption The factor loadings satisfy ||R||_2^2 approximately p^{1-delta1} and ||C||_2^2 approximately q^{1-delta2}.
    Condition 6, Section 4: defines factor strength and is needed for the stated convergence rates.
  • domain assumption Cov(vec(X_t), vec(F_t)) = 0 and Cov(vec(X_t), vec(E_t)) = 0.
    Section 3, before the LS estimator: this exogeneity condition is required for the least squares estimator of A to be unbiased and for the residual to recover the latent factor structure.
  • domain assumption The matrices M1 and M2 have k and r distinct positive eigenvalues, respectively.
    Condition 7, Section 4: ensures the eigenvectors are unique and the eigen-space estimation is well-defined.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Modeling large dimensional matrix time series with partially known and latent factors." pith.science (2026). https://pith.science/paper/XPACPMW2

@misc{pith2026241116192,
  author       = {Pith},
  title        = {Pith review of: Modeling large dimensional matrix time series with partially known and latent factors},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XPACPMW2}},
  note         = {Machine review of arXiv:2411.16192}
}
read the original abstract

This article considers to model large-dimensional matrix time series by introducing a regression term to the matrix factor model. This is an extension of classic matrix factor model to incorporate the information of known factors or useful covariates. We establish the convergence rates of coefficient matrix, loading matrices and the signal part. The theoretical results coincide with the rates in Wang et al. (2019). We conduct numerical studies to verify the performance of our estimation procedure in finite samples. Finally, we demonstrate the superiority of our proposed model using the daily returns of stocks data.

Figures

Figures reproduced from arXiv: 2411.16192 by the authors.

Figure 1
Figure 1. Data structure of multi-region and multi-industry stock returns. [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
Figure 2
Figure 2. Time series of five regions of daily returns in Energy Industry. [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Time series of five regions of Fama-French three factors. [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: The eigenvalues and their ratios of Mc1 and Mc2 for row factors and column factors. 15 [PITH_FULL_IMAGE:figures/full_fig_p015_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references · 30 canonical work pages

  1. [1]

    , year 1997

    author Carhart, M.M. , year 1997 . title On persistence in mutual fund performance . journal The Journal of Finance volume 52 , pages 57--82

  2. [2]

    , author Guo, B

    author Chang, J. , author Guo, B. , author Yao, Q. , year 2015 . title High dimensional stochastic regression with latent factors, endogeneity and nonlinearity . journal Journal of Econometrics volume 189 , pages 297--312

  3. [3]

    , author He, J

    author Chang, J. , author He, J. , author Yang, L. , author Yao, Q. , year 2023 . title Modelling matrix time series via a tensor cp-decomposition . journal Journal of the Royal Statistical Society Series B: Statistical Methodology volume 85 , pages 127--148

  4. [4]

    , author Chen, R

    author Chen, E.Y. , author Chen, R. , year 2022 . title Modeling dynamic transport network with matrix factor models: an application to international trade flow . journal Journal of Data Science volume 21 , pages 490--507 . :10.6339/22-JDS1065

  5. [5]

    , author Fan, J

    author Chen, E.Y. , author Fan, J. , year 2023 . title Statistical inference for high-dimensional matrix-variate factor models . journal Journal of the American Statistical Association volume 118 , pages 1038--1055

  6. [6]

    , author Tsay, R.S

    author Chen, E.Y. , author Tsay, R.S. , author Chen, R. , year 2020 . title Constrained factor models for high-dimensional matrix-variate time series . journal Journal of the American Statistical Association volume 115 , pages 775--793

  7. [7]

    , author Hagmann, M

    author Connor, G. , author Hagmann, M. , author Linton, O. , year 2012 . title Efficient semiparametric estimation of the fama--french model and extensions . journal Econometrica volume 80 , pages 713--754

  8. [8]

    , author Linton, O

    author Connor, G. , author Linton, O. , year 2007 . title Semiparametric estimation of a characteristic-based factor model of common stock returns . journal Journal of Empirical Finance volume 14 , pages 694--717

Show all 40 references
  1. [9]

    , year 2002

    author Davidson, J. , year 2002 . title Establishing conditions for the functional central limit theorem in nonlinear and semiparametric time series processes . journal Journal of Econometrics volume 106 , pages 243--269

  2. [10]

    , author Granger, C.W

    author Engle, R.F. , author Granger, C.W. , author Rice, J. , author Weiss, A. , year 1986 . title Semiparametric estimates of the relation between weather and electricity sales . journal Journal of the American statistical Association volume 81 , pages 310--320

  3. [11]

    , author French, K.R

    author Fama, E.F. , author French, K.R. , year 1993 . title Common risk factors in the returns on stocks and bonds . journal Journal of Financial Economics volume 33 , pages 3--56

  4. [12]

    , author French, K.R

    author Fama, E.F. , author French, K.R. , year 2015 . title A five-factor asset pricing model . journal Journal of Financial Economics volume 116 , pages 1--22

  5. [13]

    , author Liao, Y

    author Fan, J. , author Liao, Y. , author Wang, W. , year 2016 . title Projected principal component analysis in factor models . journal Annals of statistics volume 44 , pages 219

  6. [14]

    , author Li, Q

    author Fan, Y. , author Li, Q. , year 1999 . title Root-n-consistent estimation of partially linear time series models . journal Journal of Nonparametric Statistics volume 11 , pages 251--269

  7. [15]

    , year 2007

    author Gao, J. , year 2007 . title Nonlinear time series: semiparametric and nonparametric methods . publisher Chapman and Hall/CRC

  8. [16]

    , author Lu, Z

    author Gao, J. , author Lu, Z. , author Tj stheim, D. , year 2006 . title Estimation in semiparametric spatial regression . journal The Annals of Statistics , pages 1395--1435

  9. [17]

    , author Tong, H

    author Gao, J. , author Tong, H. , year 2004 . title Semiparametric non-linear time series model selection . journal Journal of the Royal Statistical Society Series B: Statistical Methodology volume 66 , pages 321--336

  10. [18]

    , author Tsay, R.S

    author Gao, Z. , author Tsay, R.S. , year 2023 . title A two-way transformed factor model for matrix-variate time series . journal Econometrics and Statistics volume 27 , pages 83--101

  11. [19]

    , author Kong, X

    author He, Y. , author Kong, X. , author Trapani, L. , author Yu, L. , year 2023 . title One-way or two-way factor model for matrix sequences? journal Journal of Econometrics volume 235 , pages 1981--2004

  12. [20]

    , author Kong, X

    author He, Y. , author Kong, X. , author Yu, L. , author Zhang, X. , author Zhao, C. , year 2024 . title Matrix factor analysis: From least squares to iterative projection . journal Journal of Business & Economic Statistics volume 42 , pages 322--334

  13. [21]

    , author Wang, Y

    author He, Y. , author Wang, Y. , author Yu, L. , author Zhou, W. , author Zhou, W.X. , year 2022 . title Matrix kendall's tau in high-dimensions: A robust statistic for matrix factor model . journal arXiv preprint arXiv:2207.09633

  14. [22]

    , author Yao, Q

    author Lam, C. , author Yao, Q. , year 2012 . title Factor modeling for high-dimensional time series: Inference for the number of factors . journal The Annals of Statistics volume 40 , pages 694 -- 726

  15. [23]

    , author Yao, Q

    author Lam, C. , author Yao, Q. , author Bathia, N. , year 2011 . title Estimation of latent factors for high-dimensional time series . journal Biometrika volume 98 , pages 901--918

  16. [24]

    , author Linton, O

    author Li, D. , author Linton, O. , author Lu, Z. , year 2015 . title A flexible semiparametric forecasting model for time series . journal Journal of Econometrics volume 187 , pages 345--357

  17. [25]

    , author Chen, E.Y

    author Liu, X. , author Chen, E.Y. , year 2022 . title Identification and estimation of threshold matrix-variate factor models . journal Scandinavian Journal of Statistics volume 49 , pages 1383--1417

  18. [26]

    , author Chen, R

    author Liu, X. , author Chen, R. , year 2016 . title Regime-switching factor models for high-dimensional time series . journal Statistica Sinica , pages 1427--1451

  19. [27]

    , author Chen, S.X

    author Mao, X. , author Chen, S.X. , author Wong, R.K. , year 2019 . title Matrix completion with covariate information . journal Journal of the American Statistical Association volume 114 , pages 198--210

  20. [28]

    , author Hasnawati, S

    author Muhammad, H.N.S. , author Hasnawati, S. , author Fiska, H. , year 2020 . title Fama-french three factor model: A study on lq 45 companies in indonesia stock exchange . journal Journal of Economics and Finance volume 11 , pages 25--30

  21. [29]

    , author Atanasovski, J

    author Psychogios, A.G. , author Atanasovski, J. , author Tsironis, L.K. , year 2012 . title Lean six sigma in a service context: A multi-factor application approach in the telecommunications industry . journal International Journal of Quality & Reliability Management volume 2...

  22. [30]

    , year 1988

    author Robinson, P.M. , year 1988 . title Root-n-consistent semiparametric regression . journal Econometrica: Journal of the Econometric Society , pages 931--954

  23. [31]

    , year 2023

    author Rohuma, H.N. , year 2023 . title Fama and french (1993) three-factor model: Evidence from conventional and shariah-compliant portfolios in bursa malaysia . journal International Journal of Business and Management volume 17 , pages 1--66

  24. [32]

    , year 2013

    author Ross, S.A. , year 2013 . title The arbitrage theory of capital asset pricing , in: booktitle Handbook of the fundamentals of financial decision making: Part I . publisher World Scientific , pp. pages 11--30

  25. [33]

    , year 2005

    author Sahay, B. , year 2005 . title Multi-factor productivity measurement model for service organisation . journal International Journal of Productivity and Performance Management volume 54 , pages 7--22

  26. [34]

    , year 1964

    author Sharpe, W.F. , year 1964 . title Capital asset prices: A theory of market equilibrium under conditions of risk . journal The Journal of Finance volume 19 , pages 425--442

  27. [35]

    , author Lo, D

    author Tian, Y. , author Lo, D. , author Xia, X. , author Sun, C. , year 2015 . title Automated prediction of bug report priority using multi-factor analysis . journal Empirical Software Engineering volume 20 , pages 1354--1383

  28. [36]

    , year 2012

    author Viroli, C. , year 2012 . title On matrix-variate regression analysis . journal Journal of Multivariate Analysis volume 111 , pages 296--309

  29. [37]

    , author Liu, X

    author Wang, D. , author Liu, X. , author Chen, R. , year 2019 . title Factor models for matrix-valued high-dimensional time series . journal Journal of Econometrics volume 208 , pages 231--248

  30. [38]

    , author Tong, H

    author Xia, Y. , author Tong, H. , author Li, W.K. , year 1999 . title On extended partially linear single-index models . journal Biometrika volume 86 , pages 831--842

  31. [39]

    , author Li, D

    author Yu, C. , author Li, D. , author Jiang, F. , author Zhu, K. , year 2024 . title Matrix garch model: Inference and application . journal Journal of the American Statistical Association , pages 1--25

  32. [40]

    , author He, Y

    author Yu, L. , author He, Y. , author Kong, X. , author Zhang, X. , year 2022 . title Projected estimation for large-dimensional matrix factor models . journal Journal of Econometrics volume 229 , pages 201--217

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.