## -------------------------------------------------------------------------------------------------------------------------------＃if（！sireseenamespace（“ biocmanager”，悄悄= true））＃install.packages（“ biocmanager”）＃biocmanager ::::install（“ ssviz”）## ----加载package，消息= false ----------------------------------------------------------------------------------------------------------------------------------------------库（SSVIZ）## ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------数据（SSVIZ）## ---加载 -bam，整理= true，警告= false -----------------------------------------------------------------<-dir（system.file（“ extdata”，package =“ ssviz”），full = true，patt =“ bam $”）ctrlbam <-readbam（bam.files [1]）[2]）## ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- ctrl.count <-getCountMatrix（ctrlbam）treat.count <-getCountMatrix（treatbam）--gettingReadCounts2,warning=FALSE----------------------------------------- ctrl.计数[1，] ## ---- plotdistro1，tidy = true，eval = false ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- # plotDistro(list(ctrl.count), type="qwidth",samplenames = c（“控制”），ncounts = counts [1]）## -----------------------------------------------------------------------------------------------＃plotdistro（列表（ctrl.count，treat.count），type =“ qwidth”，samplenames = c（“控制”，“处理”），ncounts = counts）＃＃＃＃----------------------------------------------------------------------------------------------------------------------------------------------------------＃区域< - 'chr1：3015526-3080526'＃plotregion（list（ctrl.count），region = region）## ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- # pp.ctrl<-pingpong(pctrlbam.count) # plotPP(list(pp.ctrl), samplenames=c("Control")) ## ----ntfreq,eval=FALSE-------------------------------------------------------- # pctrlbam.count<-getCountMatrix(pctrlbam) # freq.ctrl<-ntfreq(pctrlbam.count, ntlength=10) # plotFreq(freq.ctrl) ## ----------------------------------------------------------------------------- sessionInfo()