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
                          Current Approaches to Protein Sequencing
                          by seqadmin


                          Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                          04-04-2024, 04:25 PM
                        • 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

                        ad_right_rmr

                        Collapse

                        News

                        Collapse

                        Topics Statistics Last Post
                        Started by seqadmin, 04-11-2024, 12:08 PM
                        0 responses
                        13 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 04-10-2024, 10:19 PM
                        0 responses
                        17 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 04-10-2024, 09:21 AM
                        0 responses
                        14 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 04-04-2024, 09:00 AM
                        0 responses
                        43 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X