[BioC] limma: topTable
Gordon Smyth
smyth at wehi.edu.au
Fri Mar 19 01:32:32 MET 2004
At 02:32 AM 19/03/2004, Julia Engelmann wrote:
>Hi Bioconductor folks,
>
>I have really been enjoying using the limma package, but I have just come
>across a problem.
>When I use the topTable-command, the slot "result$Probe.Set.ID" does not seem
>to match the other entries I am interested in, namely M, P, t.
>I am using R 1.8.0, limma 1.5.5 on Affymetrix ATH1-121501 chips.
>
>An example would be:
>Output of the topTable-result:
>
> Probe.Set.ID M t
> P.Value B
>5598 259302_at 4.593339 38.9793 1.126260e-05 11.88238
>
>If I check on gene number 5598 it says
>geneNames(myExprSet)[5598]
>[1] "250498_at"
>
>Am I misinterpreting 5598 as the index of my ExpressionSet?
I think you've done something non-standard, anyway you can't have simply
used lmFit() and eBayes() and topTable() on our exprSet object. Please show
us enough of your code so that we can see how topTable is getting the probe
set ID's.
Gordon
>Thanks a lot for any suggestions!
>Julia
More information about the Bioconductor
mailing list