[BioC] subsetting in limma

echang4@life.uiuc.edu echang4 at life.uiuc.edu
Tue Jan 24 05:44:59 CET 2006


Hi Bioconductor users,
This is a very naive question on exporting subsets in limma.  But I am not
very familiar with pogramming.

I am trying to work through the limma user guide Sec 8.8 "Time Course
Experiments". I am able to complete the limma analysis but I am stuck
after subsetting the genes to only include genes with adj. p values < 0.05

> fit2 = contrasts.fit(fit, cont.matrix)
> fit2 = eBayes(fit2)
> results2 <- decideTests(fit2)
> selected.wt<- p.adjust(fit2$F.p.value, method="fdr") < 0.05

My question is how do I export the genes names in "selected.wt" out as a
text file?
I tried
> write.csv(selected.wt, "selected.wt.csv")
but I get a single column with TRUE and FALSE.

Thank you very much,
Edmund



More information about the Bioconductor mailing list