which one is better EdgeR or DEGSeq?
5
5
Entering edit mode
7.4 years ago
Whoknows ▴ 880

Hi

My project is on Rice, Cufflinks could not help me, I'd like to use edgeR or DEGSeq, But edgeR has more citation.

Could you please give/introduce me a simple pipeline for analyzing RNA-SEQ from FASTQ file to end !??

Thanks a lot.

RNA-Seq RPKM • 8.5k views
10
Entering edit mode
7.4 years ago

Do you actually mean DEGseq, or do you instead mean DESeq? DEGSeq should not be used by anyone for any reason what so ever. Its statistical assumptions are wrong. DESeq (actually, do use DESeq2, it has a number of improvements), on the other hand, is quite good and what I normally use.

So then the question becomes which of edgeR and DESeq2 is better. There's really no single answer to that. DESeq2 has integrated independent filtering and per-gene outlier detection (using Cook's distance), which generally makes me favor it. edgeR, however, is also nicely written and there's no reason that the exact same features couldn't be used with it, though that'd take more work on your part. Having said that, edgeR has nicer integration with things like camera() and roast(), which is unsurprising given the overlap in authors. In general, give both a try and do some independent validation to see which is better modeling your data. That's really the most objective way to make the determination.

7
Entering edit mode
7.4 years ago
3
Entering edit mode
7.4 years ago
Josh Herr 5.7k

RNA-Seq pipelines are many -- and they are often debated between experts and users. I came up with my own for my specific project, but largely I use BWA/bowtie2 and merge the resultant mapping into Cufflinks and/or edgeR/DESeq. I would recommend cufflinks -- what did not work for you? Many of us are working on organisms with much less information than Rice -- what type of errors have you received?

Your pipeline depends what you want to do - you'll need to do mapping first before you use DESeq or edgeR.

Here's a pretty exhaustive list of tools for RNA-Seq.

I would also recommend reading some papers in your area (I study plants and I know there are hundreds of papers on plant transcriptomes) -- what are other researchers using and why?

0
Entering edit mode

That's what I was wondering, there's even an iGenomes for rice from Illumina that should play nicely with Cufflinks

http://support.illumina.com/sequencing/sequencing_software/igenome.ilmn

3
Entering edit mode
7.4 years ago
Ann ★ 2.3k

I use edgeR because I know how to use it and the documentation is extensive. The biggest challenge for people seems to be understanding how to define the design matrix and contrasts for using the linear modeling features. The documentation and many examples helps with that. But probably any of the methods mentioned above are fine. I prefer using the R libraries because then I can do everything in one language and environment, including making plots, doing exploratory data analysis, and so on. I don't like cuff* programs because they make too many files that are not well documented, and sometimes the p values don't seem to make any sense. But maybe that's because I'm not using it the right way.

3
Entering edit mode
7.4 years ago

Sound like you've already gotten a lot of feedback, but here is my two cents:

sRAP won't handle the first couple steps, but it sounds like you already have RPKM values from cufflinks and that is the starting point where sRAP takes over. I'm not saying the differential expression step in sRAP is necessary much better than DESeq, limma, etc., but you did specifically ask about a pipeline (beyond just the differential expression step) and I happened have to developed such a tool ;)