[BioC] Top gene list

daphne mouzaki (RI) daphne.mouzaki at bbsrc.ac.uk
Thu Sep 14 14:13:18 CEST 2006


Hi and thanks for your response,


actually topTable and Write.fit works if the object that I have is a
MArrayLM object but.     
The problem is that I have an object "Bgene" which is a list and
includes all the columns of interest.I am sorry that I wasn't clear.I
have already generated the results but I also tried decideTests in case
I could generated them in an other way to get all the FDRs and fold
changes for each gene.


> class(Bgene)
[1] "list"
> names(Bgene)
[1] "ttAll"      "DownSameUp"
> names(Bgene$ttAll)
[1] "T11"     "T3"      "T19"     "T19_T11"
> names(Bgene$ttAll[1])
[1] "T11"
> names(Bgene$ttAll[[1]])
 [1] "Row"        "Column"     "Block"      "Name"       "ID"        
 [6] "Reps"       "ARKcloneID" "Status"     "M"          "A"         
[11] "t"          "P.Value"    "adj.P.Val"  "B"          "FoldChange"
[16] "FDRmax"     "FDR"        "FDRup"      "FDRlow"

For example, I want to get all the FDr for time point T11 with the
respective gene name FDR11<-(Bgene$ttAll[[1]])$FDRmax 

the problem is that I don't get the respective gene names automatically
when I do that and when I assign
rownames(FDR)<-sort(unique(MA$genes$ARKcloneID)) I get an error message
that I assign it to something with no dimension which is true because it
is not a matrix but a list!! 

Cheers,

Daphne



-----Original Message-----
From: Giovanni Coppola [mailto:gcoppola at ucla.edu]
Sent: 13 September 2006 17:40
To: daphne mouzaki (RI)
Cc: Bioconductor list
Subject: Re: [BioC] topTable gene list

Hi Daphne,
perhaps ?write.fit does what you want.
Best
Giovanni

On Sep 13, 2006, at 7:13 AM, daphne mouzaki ((RI)) wrote:

> hi,
>
> I have generated a list of significant genes for 3 time points 
> separately using the function "decideTests"in Limma package.Using this

> function I generated foldchanges and FDR for every gene using method 
> Benjamini and Hochberg ("BH") and p.value=0.05.
> I was wondering if it is possible to generate all the FDR and 
> foldchanges with the respective gene name for every time point.
> I tried changing the cut-off to 1 in order not to exclude anything and

> the method set to the function "decideTests" changed to "none" but it 
> doesn't work.
> Could you recommend any other possible ways?
>
> I would appreciate any help,
>
> Daphne
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/ 
> gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list