Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to plot peak density heat map

    I just start to take the bioinformatics part of my ChIP-seq experiments. and it is not too difficult for me to do mapping and peak calling. but I want to make further analysis of my data.
    Like, whether the transcription factor binding is significantly enriched in the promoters of genes belonging to the same pathway, TGFb for example. while not other pathways (JAK/STAT, Ras etc).

    So, I want to plot the binding peaks at the promoter of a specific set of genes, and rank the peaks based on the distance to TSS.

    Can anyone help to figure out, what software should I use?
    Or I have to write the script by myself, if so, should I use perl?
    I only have a bit of experience with Python.

    I attached a similar example.
    Attached Files

  • #2
    Here is another example for this kind of heat map or high density peak map.
    Attached Files

    Comment


    • #3
      Please help about this topic.

      Comment


      • #4
        looks like a custom script is needed

        Comment


        • #5
          plot

          Seqminer will be closet to do this

          Comment


          • #6
            Thank you guys.
            seqMiner seems do the job.

            Comment


            • #7
              any one has script or software to plot such heat maps

              I want revive this thread and would like to know if soemone can share the script here on forum or privatly or has any idea about a tool which can generate such heat maps of TSS distance vs tag density. Seqminer will not give thsi kind of Figure and is good mainly for Bam/ sam bed files of unprorocessed data.
              i am sure there are number of users who will be intrested in such heat maps.
              Thanks
              Last edited by mathew; 07-12-2012, 08:45 AM.

              Comment


              • #8
                Hi Mathew,

                Probably you can use the program from HOMER software.



                Originally posted by mathew View Post
                I want revive this thread and would like to know if soemone can share the script here on forum or privatly or has any idea about a tool which can generate such heat maps of TSS distance vs tag density. Seqminer will not give thsi kind of Figure and is good mainly for Bam/ sam bed files of unprorocessed data.
                i am sure there are number of users who will be intrested in such heat maps.
                Thanks

                Comment


                • #9
                  Any python scripts for this kind of heatmap?
                  Thanks!

                  Comment


                  • #10
                    You can make these kinds of plot in SeqMonk (called the aligned probes plot). You can put them over any arbitrary region of either fixed or variable sizes.

                    Comment


                    • #11
                      Thanks so much! SeqMonk looks very powerful.

                      Comment


                      • #12
                        you can plot heat map easily if you are familiar with matplotlib.

                        custom code would like this:
                        ###########
                        from matplotlib import pyplot as PLT
                        from matplotlib import cm as CM
                        from matplotlib import axes

                        ##create a figure and add a plot
                        fig = PLT.figure(facecolor='w')
                        ax1 = fig.add_subplot(2,1,1,position=[0.1,0.15,0.9,0.8])

                        ##create a cmap object. it controls the color and some style in the heatmap
                        ## you can change 'spectral' to other words. see the doc
                        cmap = CM.get_cmap('spectral', 1000)

                        ##create heatmap use the DATA and the cmap. vmin and vmax are the minimum and maximum value in your data
                        map = ax1.imshow(DATA, interpolation="nearest", cmap=cmap,aspect='auto', vmin=0,vmax=15)

                        ##then create the scale bar
                        cb = PLT.colorbar(mappable=map, cax=None, ax=None,shrink=0.5)
                        cb.set_label('Enrichment level')

                        ##plot it
                        PLT.show()
                        ###########
                        the DATA is like:
                        DATA=[
                        [1.4,2.2,4.7,6.1],
                        [4.1,13.5,15.0,1.9],
                        .....
                        ]

                        this is the website: http://matplotlib.org/
                        it is very powerful. You can plot almost everything you want using this module
                        Last edited by Luyi Tian; 02-26-2013, 05:50 AM.

                        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
                        12 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, Yesterday, 06:07 PM
                        0 responses
                        10 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-22-2024, 10:03 AM
                        0 responses
                        51 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-21-2024, 07:32 AM
                        0 responses
                        68 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X