On the degrees of freedom of reduced-rank estimators in multivariate regression

Biometrika. 2015;102(2):457-477. doi: 10.1093/biomet/asu067. Epub 2015 Feb 9.

Abstract

We study the effective degrees of freedom of a general class of reduced-rank estimators for multivariate regression in the framework of Stein's unbiased risk estimation. A finite-sample exact unbiased estimator is derived that admits a closed-form expression in terms of the thresholded singular values of the least-squares solution and hence is readily computable. The results continue to hold in the high-dimensional setting where both the predictor and the response dimensions may be larger than the sample size. The derived analytical form facilitates the investigation of theoretical properties and provides new insights into the empirical behaviour of the degrees of freedom. In particular, we examine the differences and connections between the proposed estimator and a commonly-used naive estimator. The use of the proposed estimator leads to efficient and accurate prediction risk estimation and model selection, as demonstrated by simulation studies and a data example.

Keywords: Adaptive nuclear norm; Degrees of freedom; Model selection; Multivariate regression; Reduced-rank regression; Singular value decomposition.