BaitCapture {SingleMoleculeFootprinting} | R Documentation |
check bait capture efficiency. Expected to be ~70
BaitCapture(sampleSheet, genome, baits, clObj = NULL)
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 |
bait capture efficiency
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) }