data(rsgcc)
x <- rnaseq[1:4,]
#compute correlation between 1th and 4th genes
#significance level of the computed correlation
#is calcuated with 200 permutation tests.
corpair <- cor.pair(c(1,4), GEMatrix = x, rowORcol = "row",
cormethod = "GCC", pernum = 200,
sigmethod = "two.sided")
#get the final correlation and p-value of GCC method
gcc.corfinal(corpair)
Run the code above in your browser using DataLab