[BioC] Most stable gene pairs in array experiment

David martin vilanew at gmail.com
Mon Oct 19 10:37:08 CEST 2009


Hi,
I have the following matrix with normalized log2 values:
CondA    CondB    CondC    CondD    CondE
geneA    -6.19    -5.74    -5.82    -5    -5.59
geneB    -6.33    -5.32    -5.6    -4.88    -5.39
geneC    -6.15    -6.07    -5.6    -4.88    -5.9
geneD    -6.57    -6.11    -6.36    -5.36    -5.96
geneD    -6.74    -6.2    -5.49    -5.35    -5.95
geneE    -6.75    -6.24    -5.73    -5.63    -6.02


Created as follows:
geneA<-c(-6.19,   -5.74,   -5.82,   -5,  -5.59)
geneB<-c(-6.33,   -5.32,   -5.6,    -4.88,   -5.39)
geneC<-c(-6.15, -6.07, -5.6, -4.88, -5.9)
geneD<-c(-6.57,   -6.11,   -6.36,   -5.36,   -5.96)
geneD<-c(-6.74,   -6.2,    -5.49,   -5.35,   -5.95)
geneE<-c(-6.75,   -6.24,   -5.73,   -5.63,   -6.02)

mygenes<-rbind(geneA, geneB, geneC, geneD, geneE)
colnames(mygenes)<-c("CondA",   "CondB",   "CondC",   "CondD",
"CondE")

I'm looking for most stable pair genes across conditions. I'm not 
looking for individual gene variance but really for most stable pairs 
ratios.
For e.g What is the variance of geneA vs geneB across all conditions. 
What is the most stable pair ?

Any help would be appreciated.

david



More information about the Bioconductor mailing list