## ---- 1 ------------------------------------------------------------------------------------------------------------------ library(cleanUpdTSeq) testFile <- system.file("extdata", "test.bed",package =“ clearupdtsseq”)testset <-read.table(testfile,sep =“ \ t”,header = true)peacs <-bed2grangesseq(testset,testset,perseq = false)-------------------------------------------------------------------------------------------------------峰<-Bed2Grangesseq(testset,upstream.seq.ind = 7,downstream.seq.ind = 8,fithseq = true)## ----- 3----------------------------------------------------------------------------------------------------------------------------------------- 4 -----------------------------------------------------------------------------------------------------------------------------------------库(bsgenome.drerio.ucsc.danrer7)testset.naiveBayes <-buldFeatureVector(peaks,bsgenomename = dererio,dererio,upstream = 40,下流= 30,downstream = 30,wordsize = 6,alphabet = c(“ acgt”),sampleteptys,sampleteptype types,sampleteptype=“未知”,替换量= 30,method =“ naivebayes”,zerobasedIndex = 1,fetchSeq = true)## ---- 5 ----- 5 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 数据(data.naivebayes)if(Interactive()){prectionTestSet(data.naiveBayes $ nater,data.naiveBayes $ steration.naiveSet.naiveBayes = testset.naiveSet.naiveFile =“ test-prednaiveBayes.tsv”,signtmentCutoff = 0.5)----------------------------------------------------------------------------------------------------------------数据(分类器)testResults <-predictTestSet(testset.naivebayes = testset.naiveSet.naiveBayes,classifier = classifier,upputeFile = null,null,sigitionmentCutoff = 0.5)head(testResults)## ---- sessioninfo,结果='asis'------------------------------------------------------------------------------------------------- sessioninfo()