RV coefficient

From testwiki
Jump to navigation Jump to search

In statistics, the RV coefficient[1] is a multivariate generalization of the squared Pearson correlation coefficient (because the RV coefficient takes values between 0 and 1).[2] 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:[1]

One application of the RV coefficient is in functional neuroimaging where it can measure the similarity between two subjects' series of brain scans[3] or between different scans of a same subject.[4]

Definitions

The definition of the RV-coefficient makes use of ideas[5] 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

ΣXY=E⁡(XY⊤),

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

COVV⁡(X,Y)=Tr⁡(ΣXYΣYX).

The scalar-valued variance is defined correspondingly:

VAV⁡(X)=Tr⁡(ΣXX2).

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.[5]

Then the RV-coefficient is defined by[5]

RV(X,Y)=COVV⁡(X,Y)VAV⁡(X)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.[6]

Given known diagonal blocks ΣXX and ΣYY of dimensions p×p and q×q respectively, assuming that p≤q without loss of generality, it has been proved[7] that the maximal attainable numerator is Tr⁡(ΛXΠΛY), where ΛX (resp. ΛY) denotes the diagonal matrix of the eigenvalues of ΣXX(resp. ΣYY) sorted decreasingly from the upper leftmost corner to the lower rightmost corner and Π is the p×q matrix (Ip 0p×(q−p)).

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

RV¯(X,Y)=Tr⁡(ΣXYΣYX)Tr⁡(ΛXΠΛY)=Tr⁡(ΣXYΣYX)∑j=1min(p,q)(ΛX)j,j(ΛY)j,j.

The impact of this adjustment is clearly visible in practice.[7]

See also

References

Template:Reflist