# NOT RUN {
metafeatures <- Seurat::GetAssayData(object = atac_small[['peaks']], slot = 'meta.features')
MatchRegionStats(
meta.feature = metafeatures,
regions = head(rownames(metafeatures), 10),
features.match = "percentile",
n = 10
)
# }
Run the code above in your browser using DataCamp Workspace