[BioC] Bayesian selection of genes and its correct template

Ian Jeffery ian.jeffery at ucd.ie
Fri Jul 16 13:45:45 CEST 2004


hi,
  I was looking at the Bayesian method of selecting genes. I used the
commands

library(limma)
fit <- lmFit(x,design=y)#x is the datafile
                        #y is the template
eb2<- eBayes(fit)
topTable(eb2,coef=,number=15,adjust="fdr")

I found that that the two groups could be defined by the template in two
ways
either
      (0,0,0,0,1,1,1,1)
or
      (-1,-1,-1,-1,1,1,1,1)
(I tryed other templates and all give one of the two gene lists)
Both methods return a gene list, so my questions are;
1. Are both methods correct?
2. If so why do they return different gene lists and what is the
significance of this?

thanks, Ian



More information about the Bioconductor mailing list