MeTPeak的使用

安装:

install.packages("devtools")

library("devtools")

source("https://bioconductor.org/biocLite.R")

install_github("compgenomics/MeTPeak")

使用:

library(MeTPeak)

gtf <- system.file('extdata','example.gtf',package='MeTPeak')

ip1 <- system.file('extdata','IP1.bam',package='MeTPeak')

ip2 <- system.file('extdata','IP2.bam',package='MeTPeak')

ip3 <- system.file('extdata','IP3.bam',package='MeTPeak')

input1 <- system.file('extdata','Input1.bam',package='MeTPeak')

input2 <- system.file('extdata','Input2.bam',package='MeTPeak')

input3 <- system.file('extdata','Input3.bam',package='MeTPeak')

IP_BAM <- c(ip1,ip2,ip3)

INPUT_BAM <- c(input1,input2,input3)

metpeak(GENE_ANNO_GTF=gtf,IP_BAM = IP_BAM,INPUT_BAM = INPUT_BAM,

        EXPERIMENT_NAME="example")

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容