[BioC] tukey test

David martin vilanew at gmail.com
Wed Dec 16 17:53:28 CET 2009


Hi,
One quick question.

I have an expression data set mydata.eset.
I build up the design and contrast matrix as follows.


contrast.matrix <- design.pairs(c(group))  # design pairs is a quick 
function to do it automatically
fit2 <- contrasts.fit(fit, contrast.matrix)
fit2 <- eBayes(fit2)
results <- decideTests(fit2)

Ok so far so good, it works

Now i want to compute a tukyHSD test between the same groups.
It is not possible to use the fit object directly , how would you do 
that from your eset object and design matrix ??

thanks,
david



More information about the Bioconductor mailing list