CommencezCommencez gratuitement

Quels pics sont différents?

Voyons à quel point les échantillons diffèrent sur le plan des appels de pics. En regardant les pics fusionnés avec lesquels vous avez déjà travaillé, combien de pics sont partagés entre différents échantillons?

Un type de graphique qui peut répondre à cette question est un diagramme de Venn. Vous pouvez en créer un directement à partir de l'objet DBA avec la fonction dba.plotVenn(). Une autre option est d'utiliser un graphique UpSet, que vous pouvez créer avec le paquet UpSetR. Celui-ci remplace les ellipses d'un diagramme de Venn par une série d'histogrammes qui indiquent la taille de divers sous-ensembles. La fonction upset() exige un cadre de données avec une colonne par ensemble, c.-à-d. par échantillon, qui indique l'appartenance à l'ensemble pour chaque pic. Les données nécessaires se trouvent dans le champ called de l'objet ar_diff.

Cette activité fait partie du cours

ChIP-seq avec Bioconductor en R

Voir le cours

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

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