8. 导入 rtrackalyer MACS2 输出目录中还包含一个“.narrowPeak”文件。这是一种间隔/床文件,正如我们之前所做的那样,我们可以使用 rtracklayer 包导入这些文件。 “格式”参数必须设置为“窄峰”。 library(rtracklayer)macsPeaks_GR_np<-import("data/Mel1_peaks.narrowPeak",format="narrowPeak")macsPeaks_...
### step6.7 peak calling by BayesPeak(R bioconductor package) # Bayesian Analysis of ChIP-seq data ## BayesPeak fits a Markov model to the data (the aligned reads) via Markov Chain Monte Carlo (MCMC) techniques. # 有博客里面提到I've used BayesPeak running in R. It is much easier to i...
wigpe generate wiggle files for paired reads ranger peak calling for sharp peaks ccat peak calling for broad peaks bcp peak calling for complex broad peaks ## 上面是该软件的几个用法,它直接各种格式的比对文件,我这里给的bed格式的,就是把sam转为bam再转为bed,,大家没必要那么复杂, 直接用bam格式即...