Tracy-Widom statistic for the largest eigenvalue of autoscaled real matrices

Edoardo Saccenti*, Age K. Smilde, Johan A. Westerhuis, Margriet M.W.B. Hendriks

*Corresponding author for this work

Research output: Contribution to journalArticleAcademicpeer-review

12 Citations (Scopus)

Abstract

Eigenanalysis is common practice in biostatistics, and the largest eigenvalue of a data set contains valuable information about the data. However, to make inferences about the size of the largest eigenvalue, its distribution must be known. Johnstone's theorem states that the largest eigenvalues l1 of real random covariance matrices are distributed according to the Tracy-Widom distribution of order 1 when properly normalized to L1=l1-ηnpξnp, where ηnp and ξnp are functions of the data matrix dimensions n and p. Very often, data are expressed in terms of correlations (autoscaling) for which case Johnstone's theorem does not work because the normalizing parameters ηnp and ξnp are not theoretically known. In this paper we propose a semi-empirical method based on test-equating theory to numerically approximate the normalization parameters in the case of autoscaled matrices. This opens the way of making inferences regarding the largest eigenvalue of an autoscaled data set. The method is illustrated by means of application to two real-life data sets.

Original languageEnglish
Pages (from-to)644-652
JournalJournal of Chemometrics
Volume25
Issue number12
DOIs
Publication statusPublished - 1 Dec 2011

Keywords

  • Autoscaling
  • Covariance matrix
  • Eigenanalysis
  • Largest eigenvalue
  • Tracy-Widom distribution

Fingerprint

Dive into the research topics of 'Tracy-Widom statistic for the largest eigenvalue of autoscaled real matrices'. Together they form a unique fingerprint.

Cite this