[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