## ----消息= false -------------------------------------------------------------------------------------------------------------库(HireadsProcessor)##为避免在Windows机器上避免使用软件包构建故障## if(.platform $ os.type ==“ Windows”){注册(SerialParam())} ## -------------------------------------------------------------------------------------------------------------------------------------- runData <- system.file("extdata/FLX_sample_run/",package =“ hireadsprocessor”)list.files(rundata,recursive = true)## ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------读。seqfolder(rundata,seqfilepattern =“。------------------------------------------------------------------------------------------------------------------------------------------------------------》 seqprops $ soceors $ $”1“ $样本## ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ seqProps <- findBarcodes(seqProps, sector = "all", showStats = TRUE)seqprops seqprops $部门## ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------加载(file.path(system.file(“数据”,package =)“ hireadsprocessor”),“ flx_seqprops.rdata”))## --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- # seqProps <- findPrimers(seqProps, showStats=TRUE) ## ---- eval=FALSE-------------------------------------------------------------- # seqProps <- findLTRs(seqProps, showStats=TRUE) ## ---- eval=FALSE-------------------------------------------------------------- # seqProps <- findVector(seqProps, showStats = TRUE) ## ---- eval=FALSE-------------------------------------------------------------- # seqProps <- findLinkers(seqProps, showStats=TRUE, doRC=TRUE) ## ---- eval=FALSE-------------------------------------------------------------- # seqProps <- findIntegrations(seqProps, # genomeIndices = c("hg18" = "/usr/local/genomeIndexes/hg18.noRandom.2bit"), # numServers = 2) ## ----------------------------------------------------------------------------- sampleSummary(seqProps) ## ----------------------------------------------------------------------------- sessionInfo()