Entering edit mode
5.5 years ago
Researcher
▴
130
Hi All, I have ranked my chip-seq peaks and calculated the probable super enhancers using ROSE for my samples from different subtypes of cancer. Now I want to perform unsupervised clustering based on enriched SE regions to see whether the samples from same cancer subtype cluster together or not. Can anybody please guide me how can I do this? Have someone done this before? Any help or recommendations will be highly appreciated. Thanks
Nothing special to be done. Follow the usual steps: Get a count matrix for the entire ChIP-seq dataset (all peaks), normalize e.g. with TMM from
edgeR
(please not RPKM), filter for the regions of interest, transform to Z-scale and then plug intohclust
. I would start withward.D2
as the agglomeration method and see if things look meaningful. I recommendComplexHeatmap
which has many convenience functions for clustering/visualization/plotting.Thanks ATpoint for your reply. But I am wondering about the sites that have different start and end loci but are overlapping. How should I deal with that?