Dear all,

I was running MEDIPS. I use the following command lines. But error occurs, which is showed below. I really appreciate it if someone can help me to solve it. Thanks.

library(MEDIPS)

bam.file.HIV.Rep1.MeDIP = "Jurkat_HIV.sort.rmdup.bam"

bam.file.HIV.Rep1.MeDIP = "Jurkat_HIV.sort.rmdup.bam"

library(BSgenome.Hsapiens.UCSC.hg38)

HIV_MeDIP = MEDIPS.createSet(file=bam.file.HIV.Rep1.MeDIP, BSgenome="BSgenome.Hsapiens.UCSC.hg38", uniq=1e-3, extend=300, shift=0, window_size=100)

HIV5aza_MeDIP = MEDIPS.createSet(file=bam.file.HIV5aza.Rep1.MeDIP, BSgenome="BSgenome.Hsapiens.UCSC.hg38", uniq=1e-3, extend=300, shift=0, window_size=100)

CS = MEDIPS.couplingVector(pattern="CG", refObj=HIV5aza_MeDIP)

mr.edgeR = MEDIPS.meth(MSet1=HIV_MeDIP, MSet2=HIV5aza_MeDIP, CSet=CS, p.adj="holm", diff.method="edgeR", MeDIP="F", CNV=F, minRowSum=10)

Calculating genomic coordinates... Creating Granges object for genome wide windows... Preprocessing MEDIPS SET 1 in MSet1... Preprocessing MEDIPS SET 1 in MSet2... Differential coverage analysis... Extracting count windows with at least 10 reads... Execute edgeR for count data of 0 windows... (Neglecting parameter 'type') Creating a DGEList object... Apply trimmed mean of M-values (TMM) for library sizes normalization... There is no replication, setting dispersion to bcv^2 where bcv=0.01. Please consider http://www.bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf section 2.9. Calculating differential coverage... *Error in rownames<-(*tmp*, value = c("tag.1", "tag.0")) : length of 'dimnames' [1] not equal to array extent

In addition: Warning message: In edgeR::DGEList(counts = values[filter, ], group = edRObj.group) : Zero library size detected.*