Hello,
If I have two correlation matrices (e.g. one for each of two treatments) and
then perform cor() on those two correlation matrices is this third
correlation matrix interpreted as the correlation between the two
treatments?
In my sample below I would interpret that the treatments are 0.28
correlated. Is this correct?
> var1<- c(.000000000008, .09, .1234, .5670008, .00110011002200, 8, 9,
12.34, 56.7, 1015245698745859)
> var2<- c(11, 222, 333, 4444, 55555, 1.1, 2.2, 3.3, 4.4, .55)
> var3 <-c(22.2, 66.7, 99.9, 1000008, 123123, .1, .2, .3, .4, .5)
> var4<- c(.000001,.00001,.0001, .001, .1, .12345, .56789, .67890, .78901,
.89012)
> dat <- cbind(var1,var2,var3,var4)
> dat.d <- data.frame(dat)
> treatment1 <- dat.d[1:5,]
> treatment2 <-dat.d[6:10,]
> t1.d.cor <- cor(treatment1)
> t2.d.cor <- cor(treatment2)
> I <-lower.tri(t1.d.cor)
> t1.t2 <- cor(cbind(T1 = t1.d.cor[I], T2 = t2.d.cor[I]))
> t1.t2
T1 T2
T1 1.0000000 0.2802750
T2 0.2802750 1.0000000
My code may be unpolished,
Thanks,
Michael
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
