Compute the GC content, region lengths, and dinucleotide base frequencies for regions in the assay and add to the feature metadata.
RegionStats(object, ...)# S3 method for default
RegionStats(object, genome, verbose = TRUE, ...)
# S3 method for ChromatinAssay
RegionStats(object, genome, verbose = TRUE, ...)
# S3 method for Seurat
RegionStats(object, genome, assay = NULL, verbose = TRUE, ...)
A Seurat object, Assay object, or set of genomic ranges
Arguments passed to other methods
A BSgenome object
Display messages
Name of assay to use
Returns a dataframe
# NOT RUN {
library(BSgenome.Hsapiens.UCSC.hg19)
RegionStats(
object = rownames(atac_small),
genome = BSgenome.Hsapiens.UCSC.hg19
)
# }
# NOT RUN {
library(BSgenome.Hsapiens.UCSC.hg19)
RegionStats(
object = atac_small[['peaks']],
genome = BSgenome.Hsapiens.UCSC.hg19
)
# }
# NOT RUN {
library(BSgenome.Hsapiens.UCSC.hg19)
RegionStats(
object = atac_small,
assay = 'bins',
genome = BSgenome.Hsapiens.UCSC.hg19
)
# }
Run the code above in your browser using DataLab