CommencerCommencer gratuitement

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.

Cet exercice fait partie du cours

ChIP-seq with Bioconductor in R

Afficher le cours

Exercice interactif pratique

Essayez cet exercice en complétant cet exemple de code.

# Create Venn diagram
___(ar_diff, mask=1:4)
Modifier et exécuter le code