I am analyzing some ChIP-seq data with Diffbind, and everything works great until when I was trying to use its plotprofile funtion. The command and error messages are attached below. Since I was using Ensembl GRCh38.p13 as the reference genome, I guess this arises from the format of chromosome numbers? However, I am unexperienced with R, and have no idea how to solve the problem. I really appreciate if anyone can help!
profiles <- dba.plotProfile(Ezh2.MAPQ20)
Generating report-based DBA object... Generating profiles... Error in value[[3L]](cond) : profileplyr error: Error: BiocParallel errors element index: 1, 2, 3, 4, 5, 6 first error: BiocParallel errors element index: 1 first error: different row counts implied by arguments In addition: Warning message: In .Seqinfo.mergexy(x, y) : The 2 combined objects have no sequence levels in common. (Use suppressWarnings() to suppress this warning.)