The support.bioconductor.org editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: limma: get list of genes after classifyTestsF
0
gravatar for Bandyopadhyay, Somnath Som
12.1 years ago by

I am using limma to analyze affy data. i made some contrasts and using classifyTestsF(object,p.value=0.01) I get the following results:

  E.OvsC.O I.OvsC.O IM.OvsC.O   E.TvsC.T I.TvsC.T
 IM.TvsC.T C.TvsC.O  E.TvsE.O IM.TvsIM.O M.TvsM.O
-1      34     229     798     369     469     430      26     146    1938     820
 0   54582   54142   53302   54105   53955   53971   54611   54417   50907   53367
 1      59     304     575     201     251     274      38     112    1830     488

Now I wanna make a list of the probests which change in every contrast (increases and decreases separately). how do I make a list of these probesets. Any help would be appreciated and thanks in advance.
Som.

limma classifytestsf • 368 views
ADD COMMENTlink modified 3.9 years ago by Gordon Smyth36k • written 12.1 years ago by Bandyopadhyay, Somnath Som60
Answer: limma: get list of genes after classifyTestsF
1
gravatar for Gordon Smyth
3.9 years ago by
Gordon Smyth36k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth36k wrote:

The output you show is from summary and decideTests:

results <- classifyTestsF(object,p.value=0.01)
summary(decideTests(results))

Suppose you want to see genes that are significantly up for every contrast. The following code will give you a topTable of just those genes:

i <- apply(results>0, 1, all)
topTable(fit[i,], n=Inf)

To see those that are significantly down for every contrast:

i <- apply(results<0, 1, all)
topTable(fit[i,], n=Inf)
ADD COMMENTlink modified 3.9 years ago • written 3.9 years ago by Gordon Smyth36k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 202 users visited in the last hour