 # Correlation with constraints on pairs

An interesting question was posted on http://math.stackexchange.com/726205/…: if one knows the covariances  and , is it possible to infer ? I asked myself a question close to this one a few weeks ago (that I might also relate to a question I asked a long time ago, about possible correlations between three exchange rates, on financial markets). More precisely, if one knows the correlations  and , is it possible to say something about ?

I could not find much details (but maybe I did not look enough in the existing literature). My strategy was to consider the correlation matrix, and to use the fact that a correlation matrix is symmetric, positive semidefinite matrix (also called Gramian matrix, which is a matrix with no negative eigenvalues). given the two correlations, we should consider the function of the third correlation, which indicates whether the smallest eigenvalue is non-negative, or not. Then, I look at the range of the third correlation, to get the minimum and the maximum possible value (I guess we can prove that possible values belongs to some interval). The code to get that is simply

corrminmax=function(r1,r2){
h=function(r3){
R=matrix(c(1,r1,r2,r1,1,r3,r2,r3,1),3,3)
return(min(eigen(R)\$values)>0)}
vc=seq(-1,+1,length=1e4+1)
vr=Vectorize(h)(vc)
indx=which(vr==TRUE)
return(vc[range(indx)])
}

Using this code, it is possible to look at the smallest correlation for the third pair, as well as the maximum correlation,

x1=seq(-1,1,by=.1)
x2=seq(-1,1,by=.1)
W=M=matrix(NA,length(x1),length(x2))
for(i in 1:length(x1)){
for(j in 1:length(x2)){
C=corrminmax(x1[i],x2[j])
W[i,j]=C
M[i,j]=C
}}

If we plot those matrices, we get

par(mfrow=c(1,2))
persp(x1,x2,W,zlim=c(-1,1),col="green",
persp(x1,x2,M,zlim=c(-1,1),col="green",
shade=TRUE,theta=-30) and if we plot the difference, to get the range of the interval we clearly see that the largest range is obtained when the two correlations are null (in that case, any correlation is valid) ## 5 thoughts on “Correlation with constraints on pairs”

1. Pierre says:

Bonjour Arthur,
It might be interesting for you to look up “onion method for correlation matrices” for a precise description of the convex set of correlation matrices and the uniform measure supported by it.
Cordialement

1. Arthur Charpentier says:

🙂
ça m’amuse, oui ! ça pourrait être intéressant !

1. Pierre says:

Attention : la méthode de l’oignon, ça pique les yeux 🙂
Voilà un bon point de départ :
http://ms.mcmaster.ca/canty/seminars/Joe_vinecorr_print.pdf
J’ai hâte de lire le post de blog (forcément intéressant et instructif) que tu pondras à partir de là…

1. Arthur Charpentier says:

thanks Andy,
I have to check, but from what I understood, it is only the exchangeable case… but I look at it carefully

This site uses Akismet to reduce spam. Learn how your comment data is processed.