1. Learn
  2. /
  3. Courses
  4. /
  5. ChIP-seq with Bioconductor in R

Exercise

Which peaks are different?

Let's take a look at how different the samples are in terms of peak calls. Looking at the merged peaks you've worked with before, how many peaks are shared between different samples?

One plot that can answer this question is a Venn diagram. They can be created directly from the DBA object using the dba.plotVenn() function. An alternative is to use an UpSet plot, which you can create with the UpSetR package. This replaces the ellipses of a Venn diagram with a series of bar plots, indicating the size of various subsets. The upset() function requires you to provide a data frame with one column for each set, i.e. sample, indicating set membership for each peak. The required data is contained in the called field of the ar_diff object.

Instructions 1/2

undefined XP
    1
    2
  • Create a Venn diagram from the ar_diff object.