intCor {MergeMaid} | R Documentation |
Given a mergeExprSet, this function calculates the study specific correlation matrices, and, for each gene, the correlation of correlations.
intCor(x,method= c("pearson", "kendall", "spearman"),...)
x |
Object of class mergeExprSet. |
method |
Method used to calculate correlation coefficient. |
... |
Not implemented at this time |
Integrative correlation coefficients are calcualted as follows. The first step is to identify the n genes common to all studies. Within each study, we calculate the correlation coefficient between gene g, and every other common gene. This gives a vector of length n-1. For a pair of studies, S1 and S2, we calculate the correlation of correlations for gene g. When there are more than 2 studies under consideration, all pairwise correlation of correlations are calculated and averaged.
The output is an object of class mergeCor.
data(mergeData) merged <-mergeExprs(sample1,sample2,sample3) corcor <-intCor(merged,method="spearman") plot(corcor) hist(corcor) vv<-c(1,3) corcor1 <-intCor(merged[vv],method="kendall") plot(corcor1,xlab="study A",ylab="study B",main="CORRELATION OF CORRELATION",col=3,pch=4) hist(corcor1,xlab="CORRELATION OF CORRELATION")