j.pickl at dkfz-heidelberg.de
Mon Sep 8 19:10:38 CEST 2014
thank you very much for your help.
Which tool for DE did you use?
I used edger, however I´ve read that edgeR and DESeq2 might be overstringent for RIP-Seq (f.e. RIP-Seeker package Paper, Supplement).
Von: mailinglist.honeypot at gmail.com [mailto:mailinglist.honeypot at gmail.com] Im Auftrag von Steve Lianoglou
Gesendet: Freitag, 5. September 2014 18:35
An: Julia [guest]
Cc: bioconductor at r-project.org list; Pickl, Julia
Betreff: Re: [BioC] HTSeq-Count
On Fri, Sep 5, 2014 at 7:41 AM, Julia [guest] <guest at bioconductor.org> wrote:
> Hi all,
> I am new to the field of seq and performed a RIP-Seq experiment using HTSeq count as counter.
> I get now the following (using union, but doesnÂ´t look better for interesection_strict):
> __no_feature 1503377
> __ambiguous 490772
> __too_low_aQual 0
> __not_aligned 0
> __alignment_not_unique 5277314
> When I sum up counts for all genes, I get 3227845.
> The number for __no_feature, __ambiguous, __alignment_not_unique look very high.
> Does somebody have an idea for that?
While I haven't worked with RIP-seq data myself, I do have some experience with HITS-CLIP and PAR-CLIP, which I believe are quite similar (at least in principle) -- these experiments must incredibly difficult to pull off, however, because I'd say most of these types of datasets that came my way were notoriously/incredibly noisy.
This is just to say the problem you are seeing may not be an informatics problem, and could (quite possibly) be an experimental one.
More information about the Bioconductor