Generalized reduced rank tests using the singular value decomposition

https://doi.org/10.1016/j.jeconom.2005.02.011Get rights and content

Abstract

We propose a novel statistic to test the rank of a matrix. The rank statistic overcomes deficiencies of existing rank statistics, like: a Kronecker covariance matrix for the canonical correlation rank statistic of Anderson [Annals of Mathematical Statistics (1951), 22, 327–351] sensitivity to the ordering of the variables for the LDU rank statistic of Cragg and Donald [Journal of the American Statistical Association (1996), 91, 1301–1309] and Gill and Lewbel [Journal of the American Statistical Association (1992), 87, 766–776] a limiting distribution that is not a standard chi-squared distribution for the rank statistic of Robin and Smith [Econometric Theory (2000), 16, 151–175] usage of numerical optimization for the objective function statistic of Cragg and Donald [Journal of Econometrics (1997), 76, 223–250] and ignoring the non-negativity restriction on the singular values in Ratsimalahelo [2002, Rank test based on matrix perturbation theory. Unpublished working paper, U.F.R. Science Economique, University de Franche-Comté]. In the non-stationary cointegration case, the limiting distribution of the new rank statistic is identical to that of the Johansen trace statistic.

Introduction

Tests of many economic hypotheses boil down to testing the rank of a matrix, see, for example, Lewbel (1991). Also the identification of parameters in econometric models often depends on the rank of a matrix. Therefore, since the early stages of the development of econometric methodology a literature focusing on rank tests has developed. An early contribution to this line of work is Anderson (1951). Anderson (1951) proposes a likelihood ratio rank statistic in a multivariate regression model based on canonical correlations. This canonical correlation rank statistic has a χ2 limiting distribution when the covariance matrix of the unrestricted matrix estimator has a Kronecker structure. A Kronecker covariance matrix is however not always valid. For example, it does not apply if one opts for the use of heteroscedasticity autocorrelation consistent covariance (HACC) matrix estimators, see, for example, Andrews (1991), Newey and West (1987), West (1997) and White (1980). Alternative approaches to test the rank of a matrix that overcome this deficit have therefore been proposed.

Gill and Lewbel (1992) and Cragg and Donald (1996) use a LDU decomposition (see, for example, Golub and van Loan, 1989) of the unrestricted matrix to construct a rank statistic. They use a (root-T) consistent estimator of the unrestricted matrix, which does not need to have a Kronecker covariance matrix. A disadvantage of their approach is that it can be sensitive to the ordering of the variables. To overcome this sensitivity, the estimator of the unrestricted matrix is permuted using Gaussian elimination. A high level assumption is however needed to ensure that the Gaussian elimination performs adequately.

Cragg and Donald (1997) propose a rank statistic based on a minimum χ2 criterion. They assume that the rank of the covariance matrix of the unrestricted matrix is maximal and minimize the objective function numerically. Numerical optimization is often difficult when the dimension of the optimization space is large. Hence, numerical optimization of the criterion function does not always perform satisfactorily when the number of columns or rows of the matrix of interest is large.

Recently, Robin and Smith (2000) construct a rank statistic on the basis of a (root-T) consistent estimator for the unrestricted matrix. The covariance matrix of this estimator does not need to have the Kronecker structure and its rank does not have to be maximal. The limiting distribution of the rank statistic is, however, a weighted average of χ2 distributions. The asymptotic critical values of the rank statistic are therefore not tabulated. The limiting distribution shows that this rank statistic is not always appropriately normalized so the range of values of the weights affects the power of the statistic.

To overcome the deficiencies of the rank statistics, we propose a novel rank statistic. The rank statistic involves a (root-T) consistent estimator of the unrestricted matrix, which does not need to have a Kronecker covariance matrix. We decompose the unrestricted matrix estimator using a singular value decomposition. If the limiting distribution of the unrestricted matrix estimator is normal, an orthogonal transformation of the smallest singular values has a normal limiting distribution. Our rank statistic is equal to the quadratic form of this orthogonal transformation with respect to the inverse of its covariance matrix, and hence it has a standard χ2 limiting distribution. The rank statistic simplifies to the canonical correlation rank statistic of Anderson (1951) in case of a Kronecker covariance matrix. The rank statistic can also be applied in case of non-stationary variables, for example, for cointegration testing. The asymptotic distribution is in that case a functional of Brownian motions and equal to the asymptotic distribution of the Johansen (1991) trace statistic.

The novel rank statistic is related to the rank statistic proposed by Ratsimalahelo, 2001, Ratsimalahelo, 2002. The limiting distribution of his rank statistic is however based on asymptotic normality of the smallest singular values. Since singular values are non-negative, they cannot have a mean zero normal limiting distribution. We therefore consider an orthogonal transformation of the smallest singular values which removes the non-negativity restriction on the singular values.

The outline of the paper is as follows. In Section 2 we discuss the singular value decomposition of a matrix. We show how this decomposition can be used to obtain a parameter that reflects rank reduction. In Section 3 we construct a rank statistic based on this parameter. We derive its limiting distribution and show that the statistic is a generalization of the canonical correlation rank statistic. We limit ourselves in this section to stationary variables and compare the rank statistic with existing rank statistics. The non-stationary cointegration case is discussed in Section 4. In this case, the limit distribution of our rank statistic is identical to that of the Johansen trace statistic. In Section 5, we apply the rank statistic to test the rank of a matrix that governs the identification of the parameters in the stochastic discount factor model of Jagannathan and Wang (1996). Furthermore, we use a factor model which stems from this application to illustrate the validity of the limiting distribution and to conduct a power comparison of different rank statistics. We end with some concluding remarks and topics for further research in Section 6.

Throughout the paper we use the notation: a=vec(A) for the column vectorization of the k×m matrix A such that for A=(A1Am), vec(A)=(A1Am), Im is the m×m identity matrix and tr(A) stands for the trace of a square matrix A. Furthermore, “p” stands for convergence in probability and “d” for convergence in distribution.

Section snippets

Singular value decomposition and rank reduction

To construct our rank test, we use the fact that a k×m real matrix C can be decomposed asC=AqBq+Aq,ΛqBq,,with Aq a k×q matrix, Bq a q×m matrix, Λq a (k-q)×(m-q) matrix, Aq, a k×(k-q) matrix, Bq, a (m-q)×m matrix and where AqAq,0, Bq,Bq0, Aq,Aq,Ik-q and Bq,Bq,Im-q with q<min(k,m), see Kleibergen and van Dijk (1998) and Kleibergen and Paap (2002). If Λq=0, the rank of the matrix C is determined by the rank of AqBq. If both Aq and Bq have full rank, the rank of C is then equal to

Rank test: the standard (stationary) case

We are concerned with testing the rank of the k×m matrix of parameters Π. We assume that the matrix is non-symmetric in case k=m and hence Π is not allowed to be a covariance matrix. Under our null hypothesis, the rank of the matrix Π is equal to q with q<min(k,m), that is, H0:rank(Π)=q. To derive the limiting distributions of our test statistics and estimators, we make an assumption about the limiting behavior of an estimator of the unrestricted value of Π denoted by Π^, where we use the

Rank test: the non-stationary cointegration case

In vector autoregressive (VAR) models, cointegration implies a reduced rank value of the long-run multiplier, see, for example, Engle and Granger (1987) and Johansen, 1991, Johansen, 1995. A convenient specification of a VAR(1) model to reflect cointegration is the so-called error correction specificationΔyi=Πyi-1+εifori=1,,T,where yi is a k×1 vector that contains time-series observations of the variable y at time i, Δyi=yi-yi-1, εi is a k×1 vector of disturbances which we assume to be

Testing rank in stochastic discount factor models

Jagannathan and Wang (1996) base the stochastic discount factor (SDF) in a SDF model on a conditional capital asset pricing model (CAPM). A SDF model for asset returns is represented by, see, for example, Campbell et al. (1997),E[rji+1si+1|Ii]=1,j=1,,k,where rji is the return at time i on asset j, j=1,,k, Ii is the information set at time i and si is the SDF at time i. Based on the conditional CAPM, Jagannathan and Wang (1996) specify the SDF si+1 as a linear function of a set of additional

Concluding remarks

To overcome some of the deficiencies of statistics that test the rank of a matrix, we propose a novel one. The rank statistic requires a consistent estimator of the unrestricted matrix, which does not need to have a Kronecker covariance matrix. This allows one, for example, to use HACC estimators for the covariance matrix. The unrestricted matrix is decomposed using a singular value decomposition. The rank statistic is a quadratic form of an orthogonal transformation of the smallest singular

Acknowledgements

We thank three anonymous referees and participants of ESEM 2003 in Stockholm and the EC2 meeting in Bologna for helpful comments and suggestions. The first authors’ research documented in this article has been partially funded by the NWO Vernieuwingsimpuls research grant “Empirical Comparison of Economic Models”.

References (31)

  • R.F. Engle et al.

    Co-integration and error correction: representation, estimation and testing

    Econometrica

    (1987)
  • L. Gill et al.

    Testing the rank and definiteness of estimated matrices with applications to factor, state-space and ARMA Models

    Journal of the American Statistical Association

    (1992)
  • G.H. Golub et al.

    Matrix Computations

    (1989)
  • L.P. Hansen

    Large sample properties of generalized method moments estimators

    Econometrica

    (1982)
  • R.A. Horn et al.

    Topics in Matrix Analysis

    (1991)
  • Cited by (1668)

    • Indirect effects of trading restrictions

      2024, Journal of Corporate Finance
    View all citing articles on Scopus
    View full text