How to generate a pie chart of genomic annotations of ChIP peaks

Question

[cbio]

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?

[h.mon]

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.

Answer 1

[vivekbhr]

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.

Answer 2

[ATpoint]

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.