[R] Correlation matrix from a vector of pairwise correlations
Ben Bolker
bolker at zoo.ufl.edu
Sun Dec 4 03:12:52 CET 2005
Serguei Kaniovski <kaniovsk <at> wifo.ac.at> writes:
>
> I've a vector of pairwise correlations in the order low-index element
> precedes the high-index element, say:
>
> corr(1,2)=0.1, corr(1,3)=0.2, corr(2,3)=0.3, corr(3,4)=0.4
>
> How can I construct the corresponding correlation matrix?
Not absolutely sure what you want to do, but guessing a little
bit (should your list above include corr(1,4) and corr(2,4) as well?)
If you have
corrs <- c(c12=0.1,c13=0.2,c14=0.3,c23=0.4,c24=0.5,c34=0.6)
## names unnecessary but included for clarity
m = diag(4) ## 4x4 diagonal matrix
m[lower.tri(m)] = corrs ## fill in lower triangle
## two ways to make it symmetric:
m = t(m) ## flip around matrix
m[lower.tri(m)] = corrs ## fill in lower triangle
## OR
m[row(m)>col(m)] = t(m)[row(m)<col(m)]
don't know which is more efficient, but either should work ---
if that's what you were trying for
Ben Bolker
More information about the R-help
mailing list