Displaying 1 result from an estimated 1 matches for "corsp".
Did you mean:
corp
2013 Sep 23
0
Correlate rows of 2 matrices
Hi,
You may try:
set.seed(49)
m1 = matrix(rnorm(30), nrow = 3)
m2 = matrix(rnorm(30), nrow = 3)
?corsP<-vector()
? for(i in 1:3) corsP[i] =? cor(m1[i,], m2[i,])
?corsP
#[1]? 0.58411274 -0.02382329? 0.03760757
diag(cor(t(m1),t(m2)))
#[1]? 0.58411274 -0.02382329? 0.03760757
#or
mNew<- rbind(m1,m2)
?indx<-rep(seq(nrow(mNew)/2),2)
?sapply(split(seq_len(nrow(mNew)),indx),function(x) cor(t(mNew...