search

ZW-xjtlu / exomePeak2

Peak calling and differential methylation for MeRIP-Seq
25 stars 5 forks source link

exomePeak2运行出错 #48

Closed tslzs closed 6 months ago

tslzs commented 6 months ago

以下是我的运行代码：version: 1.14.3 setwd("/data2/XXX/ZJ_MeRIPseq") library(exomePeak2)

GTF_FILE <- "/data1/database/GRCH38.94/Homo_sapiens.GRCh38.94.chr.gtf"

res <- exomePeak2( bam_input = c("WT_bam/0h_Input.bam"), bam_ip = c("WT_bam/0h_IP.bam"), bam_input_treated = c("WT_bam/4h_Input.bam"), bam_ip_treated = c("WT_bam/4h_IP.bam"), genome = "hg38", gff = GTF_FILE, save_dir = "exomePeak2_WT_0h_4h" )

Extract bin features ... OK Count reads on bin features ... OK Identify background features ... Error in if (G[1] == 1) { : missing value where TRUE/FALSE needed Calls: exomePeak2 ... classifyBackground -> Mclust -> eval -> eval -> mclustBIC In addition: There were 47 warnings (use warnings() to see them) 请问可能出现的问题在哪呢？

tslzs commented 6 months ago

已经找到解决方法，是因为我们使用的bam文件里面的染色体没有chr前缀，添加之后就可以成功使用了。