[BioC] DESeq multi-factor export results
Simon Anders
anders at embl.de
Wed Mar 21 16:20:24 CET 2012
Hi Ingrid
On 03/21/2012 04:08 PM, Ingrid Lindquist wrote:
> See below. If I'm not mistaken the padjGLM file has the adjusted pvals
> in the order that is in fit1, so I've been manually combining them. Let
> me know though if this is not the correct assumption..
Yes, this is correct.
You can also combine the results into a data frame:
> write.table(fit1, file="fit_c_ctrl")
> write.table(padjGLM, file="padjGLM_c_ctrl")
becomes
write.table(
data.frame(
geneID = row.names(counts(cds)),
pval = pvalsGLM,
padj = padjGLM,
fit1 ),
file = "result.txt" )
Simon
More information about the Bioconductor
mailing list