Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • ERANGE to call peaks from ChIP-Seq data

    Hello!
    I try to call peaks from a bigger data set (~6 million reads) with erange.
    Running it with the default parameters works but only results into ~7000 peaks. With my factor I expect around 20000 binding sites and I try to get 40000 to analyze in a less restrictive way.
    However I fail generating that many peaks with erange. Either I get not enough peaks or the server gets frozen! The latter happened with the following command:


    python ERANGE3.2/commoncode/findall.py er chip.rds ErangePeaks_min1_minPeak0.000001_spacing200_fold 1_shift -control mock.rds ‐minimum 1 -ratio 2 ‐minPeak 0.000001 ‐revbackground -shift 80 -spacing 200

    I set -minPeak very low to get a lot of peaks...

    Can anyone help???

  • #2
    If letting it learn the shift value instead of specifying it to be 80 doesn't help, then 7000 is probably what you have in that library.

    Where does the number 20000 come from? Another peak caller or simply expectations?

    Comment


    • #3
      Thanks for replying!
      I tried what you said and I still do not get more peaks.
      Also it doesnt matter if I put the parameter ‐minPeak to 0.00001 or to 0.01

      The 20,000 is what I get with the default options of other peak callers. And in total, this data set contains more than 100,000 clusters (this is what I can reach if I am very little restrictive).

      Does anyone have experience with this? Is there anything I could try?
      Thank you!

      Comment


      • #4
        In general, the settings you are using are very liberal, I personally wouldn't feel comfortable working with such peak calls - in my view it's better to have a few false negatives than a lot of false positives when looking at ChIP-Seq data and "the more peaks the better" approach isn't exactly the best approach to understanding the biology of the protein.

        That said, if you don't have any high confidence prior expectation on the number of calls you should get (say you're working with something that has never been ChIP-ed before and not with NRSF or CTCF) a few things would be helpful to know:

        - How many of the peaks contain an explanatory motif (although this wouldn't help if you're not dealing with a sequence-specific transcription factor, of course). How many of the peaks you are getting with a not so conservative peak caller contain it?

        - what is the size of the genome you're working with?

        - what is your IP efficiency as calculated by ERANGE? (second to last line in the output file)

        - How many peaks do you get with more conservative settings? (-minimum 2, -ratio 3, -minimum 4 -ratio 4, etc.)

        Comment


        • #5
          The reason why I need so many peaks is because I want to compare the top 500, 1000, ... , 40000 peaks of different methods in terms of %motif (fraction of peaks with motif).

          Basically I am interested in the scoring system of different methods. If there is no possibility to increase the number of peaks I have to use another method for my comparisons..

          I am working on hg18.

          The IP efficiency:
          #stats: 77936.5 RPM in 7313 regions

          I didnt try to use more conservative settings because I want to be as little restrictive as possible. As I said I want to explore the scoring system.

          BTW, may I ask who you are ? Did you help developing the program or are you a user?
          thanks!

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Strategies for Sequencing Challenging Samples
            by seqadmin


            Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
            03-22-2024, 06:39 AM
          • seqadmin
            Techniques and Challenges in Conservation Genomics
            by seqadmin



            The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

            Avian Conservation
            Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
            03-08-2024, 10:41 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Yesterday, 06:37 PM
          0 responses
          10 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, Yesterday, 06:07 PM
          0 responses
          9 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-22-2024, 10:03 AM
          0 responses
          49 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-21-2024, 07:32 AM
          0 responses
          67 views
          0 likes
          Last Post seqadmin  
          Working...
          X