BaitCapture {SingleMoleculeFootprinting}R Documentation

Bait capture efficiency

Description

check bait capture efficiency. Expected to be ~70

Usage

BaitCapture(sampleSheet, genome, baits, clObj = NULL)

Arguments

sampleSheet

QuasR sample sheet

genome

BS genome

baits

Full path to bed file containing bait coordinates. If chromosome names are in e.g. "1" format, they'll be temporarily converted to "chr1"

clObj

cluster object to emply for parallel processing created using the parallel::makeCluster function. Defaults to NULL

Value

bait capture efficiency

Examples

Qinput = paste0(tempdir(), "/NRF1Pair_Qinput.txt")
library(BSgenome.Mmusculus.UCSC.mm10)

if(file.exists(Qinput)){
    # DO NOT RUN
    # clObj = parallel::makeCluster(5)
    # BaitRegions = SingleMoleculeFootprintingData::EnrichmentRegions_mm10.rds()
    # BaitCaptureEfficiency = BaitCapture(sampleSheet = Qinput, genome = BSgenome.Mmusculus.UCSC.mm10, baits = BaitRegions, clObj=clObj)
    # parallel::stopCluster(clObj)
}


[Package SingleMoleculeFootprinting version 1.2.0 Index]