Match DNA sequence characteristics

Return a vector if genomic regions that match the distribution of a set of query regions for any given set of characteristics, specified in the input meta.feature dataframe.

MatchRegionStats(
  meta.feature,
  regions,
  features.match = c("GC.percent"),
  n = 10000,
  verbose = TRUE,
  ...
)

Arguments

meta.feature	A dataframe containing DNA sequence information
regions	Set of query regions. Must be present in rownames.
features.match	Which features of the query to match when selecting a set of regions. A vector of column names present in the feature metadata can be supplied to match multiple characteristics at once. Default is GC content.
n	Number of regions to select, with characteristics matching the query
verbose	Display messages
...	Arguments passed to other functions

Value

Returns a character vector

Examples

metafeatures <- Seurat::GetAssayData(
  object = atac_small[['peaks']], slot = 'meta.features'
)
MatchRegionStats(
  meta.feature = metafeatures,
  regions = head(rownames(metafeatures), 10),
  features.match = "percentile",
  n = 10
)
#> Matching percentile distribution
#>  [1] "chr1:1166366-1168282" "chr1:873704-873830"   "chr1:870713-871075"  
#>  [4] "chr1:954372-958413"   "chr1:1549446-1552535" "chr1:994056-995654"  
#>  [7] "chr1:1071616-1073471" "chr1:841866-842572"   "chr1:948133-951142"  
#> [10] "chr1:1309645-1311492"