RV coefficient

In statistics, the RV coefficient{{Cite journal

| author1 = Robert, P.

| author2 = Escoufier, Y.

| title = A Unifying Tool for Linear Multivariate Statistical Methods: The RV-Coefficient

| journal = Applied Statistics

| year = 1976

| volume = 25

| issue = 3

| pages = 257–265

| doi = 10.2307/2347233

| jstor=2347233

}}

is a multivariate generalization of the squared Pearson correlation coefficient (because the RV coefficient takes values between 0 and 1).{{Cite book

| last = Abdi

| first = Hervé

| title = RV coefficient and congruence coefficient

| publisher = Thousand Oaks

| year = 2007

|editor1-first= Neil J

|editor1-last= Salkind

| isbn = 978-1-4129-1611-0

}}

It measures the closeness of two set of points that may each be represented in a matrix.

The major approaches within statistical multivariate data analysis can all be brought into a common framework in which the RV coefficient is maximised subject to relevant constraints. Specifically, these statistical methodologies include:

:*principal component analysis

:*canonical correlation analysis

:*multivariate regression

:*statistical classification (linear discrimination).

One application of the RV coefficient is in functional neuroimaging where it can measure

the similarity between two subjects' series of brain scans{{Cite journal

| title = Group analysis in functional neuroimaging: selecting subjects using similarity measures

| author1 = Ferath Kherif

| author2 = Jean-Baptiste Poline

| author3 = Sébastien Mériaux

| author4 = Habib Banali

| author5 = Guillaume Plandin

| author6 = Matthew Brett

| journal = NeuroImage

| volume = 20

| issue = 4

| year = 2003

| pages = 2197–2208

| doi = 10.1016/j.neuroimage.2003.08.018

| pmid=14683722

| url = https://hal-cea.archives-ouvertes.fr/cea-00371054/file/Kherifetal_NeuroImage.pdf

}}

or between different scans of a same subject.{{Cite journal

| title = How to compute reliability estimates and display confidence and tolerance intervals for pattern classiffers using the Bootstrap and 3-way multidimensional scaling (DISTATIS)

| author1 = Herve Abdi

| author2 = Joseph P. Dunlop

| author3 = Lynne J. Williams

| journal = NeuroImage

| volume = 45

| year = 2009

| pages = 89–95

| doi = 10.1016/j.neuroimage.2008.11.008

| pmid = 19084072

| issue = 1

}}

Definitions

The definition of the RV-coefficient makes use of ideas{{Cite journal

| author1 = Escoufier, Y.

| title = Le Traitement des Variables Vectorielles

| journal = Biometrics

| year = 1973

| volume = 29

| issue = 4

| pages = 751–760

| doi = 10.2307/2529140

| publisher = International Biometric Society

| jstor = 2529140

}}

concerning the definition of scalar-valued quantities which are called the "variance" and "covariance" of vector-valued random variables. Note that standard usage is to have matrices for the variances and covariances of vector random variables.

Given these innovative definitions, the RV-coefficient is then just the correlation coefficient defined in the usual way.

Suppose that X and Y are matrices of centered random vectors (column vectors) with covariance matrix given by

:\Sigma_{XY}=\operatorname{E}( XY^\top) \,,

then the scalar-valued covariance (denoted by COVV) is defined by

:\operatorname{COVV}(X,Y)= \operatorname{Tr}(\Sigma_{XY}\Sigma_{YX}) \, .

The scalar-valued variance is defined correspondingly:

:\operatorname{VAV}(X)= \operatorname{Tr}(\Sigma_{XX}^2) \, .

With these definitions, the variance and covariance have certain additive properties in relation to the formation of new vector quantities by extending an existing vector with the elements of another.

Then the RV-coefficient is defined by

:\mathrm{RV}(X,Y) = \frac

{ \operatorname{COVV}(X,Y) }

{ \sqrt{ \operatorname{VAV}(X) \operatorname{VAV}(Y) } } \, .

Shortcoming of the coefficient and adjusted version

Even though the coefficient takes values between 0 and 1 by construction, it seldom attains values close to 1 as the denominator is often too large with respect to the maximal attainable value of the denominator.{{Cite journal

| author1 = Pucetti, G.

| title = Measuring Linear Correlation Between Random Vectors

| journal = SSRN

| year = 2019

| url = http://dx.doi.org/10.2139/ssrn.3116066s

}}

Given known diagonal blocks \Sigma_{XX} and \Sigma_{YY} of dimensions p\times p and q\times q respectively, assuming that p \le q without loss of generality, it has been proved{{Cite journal

| author1 = Mordant Gilles

| author2 = Segers Johan

| title = Measuring dependence between random vectors via optimal transport,

| journal = Journal of Multivariate Analysis

| volume = 189

| year = 2022

}}

that the maximal attainable numerator is \operatorname{Tr}(\Lambda_X \Pi \Lambda_Y),

where \Lambda_X (resp. \Lambda_Y ) denotes the diagonal matrix of the eigenvalues of \Sigma_{XX} (resp. \Sigma_{YY} ) sorted decreasingly from the upper leftmost corner to the lower rightmost corner and \Pi is the p \times q matrix (I_p \ 0_{p\times (q-p)} ).

In light of this, Mordant and Segers proposed an adjusted version of the RV coefficient in which the denominator is the maximal value attainable by the numerator. It reads

:\bar{\operatorname{RV}}(X,Y) = \frac{\operatorname{Tr}(\Sigma_{XY}\Sigma_{YX})}{\operatorname{Tr}(\Lambda_X \Pi \Lambda_Y)} = \frac{\operatorname{Tr}(\Sigma_{XY}\Sigma_{YX})}{\sum_{j=1}^{min(p,q)} (\Lambda_X)_{j,j} (\Lambda_Y)_{j,j}}.

The impact of this adjustment is clearly visible in practice.

See also

References