peak set before iterative merging across clusters #1301
-
|
Hi @rcorces, I used addReproduciblePeakSet function to add peaks per cluster. I know we can extract the final iteratively merged peak set across clusters from the data. But, Is it possible to get peaks that were iteratively merged to form a cluster-specific peak set and not those that were iteratively merged across the clusters (please see the image). I would like to do some model training with the peak set derived before the iterative merging process across the clusters. There must be a way to obtain this information from all the metadata stored in ArchR project? Thanks Peaks right before this step ->
|
Beta Was this translation helpful? Give feedback.
Replies: 1 comment 1 reply
-
|
On further digging, I found that under PeakCalls folder in the project there are: Cluster1-Sample-summits.rds for all the samples and Cluster1-reproduciblePeaks.gr.rds (for all the clusters). I believe Cluster1-reproduciblePeaks.gr.rds are peaks iteratively merged in a cluster and not merged across the clusters. Please let me know if my understanding is correct. Thanks |
Beta Was this translation helpful? Give feedback.
-
|
correct |
Beta Was this translation helpful? Give feedback.
On further digging, I found that under PeakCalls folder in the project there are: Cluster1-Sample-summits.rds for all the samples and Cluster1-reproduciblePeaks.gr.rds (for all the clusters). I believe Cluster1-reproduciblePeaks.gr.rds are peaks iteratively merged in a cluster and not merged across the clusters. Please let me know if my understanding is correct. Thanks