[Bioc-devel] R: Re: Deseq2 and differentia expression

2014-07-14 Thread jarod...@libero.it
Dear Professo Thanks for the clarification. I'd like to use that method as a starting list for do some enrichment analisys using SPIA. res <-results(dds) ddsNoPrior <- DESeq(ddHTSeq, betaPrior=FALSE) #only for lessABs resGA <- results(ddsNoPrior, lfcThreshold=.5, altHypothesis="lessAbs") #great

Re: [Bioc-devel] R: Re: Deseq2 and differentia expression

2014-07-12 Thread Simon Anders
Hi Jarod Mike overlooked one point in your question On 11/07/14 12:05, jarod...@libero.it wrote: If I interested in the genes that have a foldchange more than 0.5 and 2 I need to use this comand is it right? DESeq2 reports al fold changes on a log2 scale. So your limits of 0.5 and 2 for unlo

[Bioc-devel] R: Re: Deseq2 and differentia expression

2014-07-11 Thread jarod...@libero.it
Dear Dr, Thanks so much for clarification!!! So I try the test of log fold change but I'm bit confusion on the results: If I interested in the genes that have a foldchange more than 0.5 and 2 I need to use this comand is it right? ddsNoPrior <- DESeq(ddHTSeq, betaPrior=FALSE) #only for lessABs re