[BioC] limma, TopTable help

Gordon Smyth smyth at wehi.edu.au
Sat Nov 6 04:20:44 CET 2004

>Date: Thu, 4 Nov 2004 15:18:17 +1100 (EST)
>From: Anthony Bosco <unbelieve_04 at yahoo.com.au>
>Subject: [BioC] limma, TopTable help
>To: bioconductor at stat.math.ethz.ch
>Message-ID: <20041104041817.15180.qmail at web53708.mail.yahoo.com>
>Content-Type: text/plain; charset=iso-8859-1
>When I go through the limma examples with Hgu133plus2
>data I can't get the M, A, and probe set IDs in the
>toptable result.

That is because you are not using the same function as in the examples. Use 
topTable() not toptable().

When I updated the older functions lm.series(), ebayes() and toptable() to 
be more object-orientated, I kept the original function because some people 
prefer them for programming purposes. The newer more oop user-interface 
functions are lmFit(), eBayes() and topTable(). If you want all the bells 
and whistles of the examples, you must use these.


>my code below.
>fit<-lmFit(eset, design)
>fit2<-contrasts.fit(fit, cont.matrix)
>toptable(fit2, coef=2, adjust.method="holm",number=10,

More information about the Bioconductor mailing list