Subset McPeaks by certain clusters

Usage

subset_peak_clusters(
  atac_mc,
  cluster_membership,
  clusters_to_keep,
  reverse = TRUE
)

Arguments

atac_mc

an McPeaks object

cluster_membership

which cluster each peak is a member of

clusters_to_keep

a vector of clusters to keep (or exclude, if reverse == TRUE)

reverse

(optional) - a logical/flag whether to keep (default - TRUE) or remove the clusters in clusters_to_keep

Value

the atac_mc object only with the clusters (peaks) of interest (not saved in the "ignore_..." slots)

Examples

if (FALSE) {
peak_cl_km <- gen_atac_peak_clust(atac_mc, k = 15)
atac_mc_subset <- subset_peak_clusters(atac_mc, cluster_membership = peak_cl_km, clusters_to_keep = c(4, 5, 8))
}