if (FALSE) {
library(BSgenome.Hsapiens.UCSC.hg19)
RegionStats(
object = rownames(atac_small),
genome = BSgenome.Hsapiens.UCSC.hg19
)
}
if (FALSE) {
library(BSgenome.Hsapiens.UCSC.hg19)
RegionStats(
object = atac_small[['peaks']],
genome = BSgenome.Hsapiens.UCSC.hg19
)
}
if (FALSE) {
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