How to generate a pie chart of genomic annotations of ChIP peaks
3
0
Entering edit mode
6.5 years ago
cbio ▴ 450

I am interested in creating a Pie charts of ChIP-Seq reads corresponding to annotated regions of the hg19 genome.  I want to count the total number of peaks at given genomic locations were counted and standardize them to each region’s occupancy of the human genome (for example, promoters only account for ~0.8% of the human genome).

NOTE: I'm aware that pie charts are often looked down by this the bioinformatics community, so if anyone has another visualization option please feel free to let me know.

How would I go about generating this type of pie chart?

ChIP-Seq • 3.5k views
ADD COMMENT
0
Entering edit mode

Most of the time, a bar chart will do the job for you. Pie charts are considered a bad way of conveying information, see this article for reasons why. The first example is compelling enough to drop pie charts entirely.

ADD REPLY
1
Entering edit mode
6.5 years ago
vivekbhr ▴ 660

You can try HOMER to annotate the peaks then plot it the way you want, or if you work in R/bioconductor, try ChIPPeakAnno. It can give you both pie and bar charts.

ADD COMMENT
1
Entering edit mode
2.3 years ago
ATpoint 62k

annotatR can annotate genomic regions based on either built-in or custom annotations. It offers various plotting functions. For a pie chart one would need to use the output table and then make a pie manually, e.g. with ggplot.

ADD COMMENT

Login before adding your answer.

Traffic: 1600 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6