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
                          Exploring the Dynamics of the Tumor Microenvironment
                          by seqadmin




                          The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                          07-08-2024, 03:19 PM
                        • seqadmin
                          Exploring Human Diversity Through Large-Scale Omics
                          by seqadmin


                          In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                          06-25-2024, 06:43 AM

                        ad_right_rmr

                        Collapse

                        News

                        Collapse

                        Topics Statistics Last Post
                        Started by seqadmin, 07-10-2024, 07:30 AM
                        0 responses
                        30 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 07-03-2024, 09:45 AM
                        0 responses
                        202 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 07-03-2024, 08:54 AM
                        0 responses
                        212 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 07-02-2024, 03:00 PM
                        0 responses
                        194 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X