Visualising differences in protein binding

Both PCA plot and dendrogram indicate that the samples from primary tumors and treatment-resistant tumors from one cluster each. This is encouraging but tells you very little about what these differences look like. In this exercise, you will create a heatmap that compares peak intensity across samples. This can help to highlight patterns in protein binding that distinguish between groups of samples.

The set of peaks, consolidated across samples, is available as peaks. To create a heatmap of peaks from this, you'll have to know how many peaks there are in the dataset. Details of the merged peak set are available in the merged entry of peaks.

Cet exercice fait partie du cours

ChIP-seq with Bioconductor in R

Afficher le cours

Instructions

Print the peaks object.
Obtain the coordinates of the merged peaks.
Extract the number of peaks present in the data.
Create a heatmap using the dba.plotHeatmap function.

Exercice interactif pratique

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

# Print the `peaks` object
print(___)

# Obtain the coordinates of the merged peaks
merged_peaks <- peaks$___

# Extract the number of peaks present in the data
peak_count <- nrow(___)

# Create a heatmap using the `dba.plotHeatmap()` function
___(peaks, maxSites = ___, correlations = FALSE)

Modifier et exécuter le code