Hello,

I am trying to conduct a DMR analyses using the DMRcate package, and all seems to work well until the point where I attempt to extract the results. The error says:

snapshotDate(): 2020-04-27 see ?DMRcatedata and browseVignettes('DMRcatedata') for documentation loading from cache Error in h(simpleError(msg, call)) : error in evaluating the argument 'x' in selecting a method for function 'as.data.frame': object 'NCList_find_overlaps_in_groups' not found

I have tried to look on github what it means and my DMR results look fine (1300 regions) I just can't figure it out why this error keeps on showing.

The code for my analyses is bellow:

library(DMRcate) library(IRanges) fit1 <- lmFit(object = Combined_First_Second_M, design = design) #run with lm() fit2 <- eBayes(fit1) #do not worry about this -- in stats we trust

fdr=1 #you want the full table results_zscorefull=topTable(fit2, coef = "zscore", number = nrow(Combined_First_Second_M), adjust.method = "BH", p.value = fdr)

fit1beta <- lmFit(object = Combined_First_Second_beta, design = design) #run with lm() fit2beta <- eBayes(fit1) #do not worry about this -- in stats we trust

fdr=1 #you want the full table results_zscorefullbeta=topTable(fit2beta, coef = "zscore", number = nrow(Combined_First_Second_beta), adjust.method = "BH", p.value = fdr)

CpGs <- rownames(results_zscorefull)

rownames(results_zscorefull)=as.character(CpGs) annotation_overlap_only=annotation[is.element(annotation$probeID,intersect(annotation$probeID,CpGs)),] annotation_overlap_only_2=arrange(annotation_overlap_only,probeID)

annotated <- GRanges(as.character(annotation_overlap_only_2[CpGs,"CpG_chrm" ]), #chromosome IRanges(start=c(annotation_overlap_only_2[CpGs,"CpG_beg"]),end=c(annotation_overlap_only_2[CpGs,"CpG_end"])), #add location #put the same location twice so it recognises as a region stat = results_zscorefull[CpGs,"t"], #t-statistic diff = results_zscorefullbeta[CpGs,"logFC"], #effect size ind.fdr = results_zscorefull$adj.P.Val, #adjusted p-value is.sig = results_zscorefull$adj.P.Val < 0.005) #p-value threshold #logical operator TRUE of FALSE so it recognises later on

names(annotated) <- annotation_overlap_only_2$probeID #name each line with the CpG annotated <- sort(annotated) #sorting the CpGs by position on the DNA annotated <- new("CpGannotated", ranges = annotated) #create a "CpGannotated" object

DMR <- dmrcate(annotated, lambda=1000, C=2, min.cpgs = 2) #C is an statistical factor, optimal at 2 #minimal number of CpGs that wou would consider to have a DMR default=2 DMR

obtain results
resultsRanges <- extractRanges(DMR, genome = "hg38")

sessionInfo( )

snapshotDate(): 2020-04-27 see ?DMRcatedata and browseVignettes('DMRcatedata') for documentation loading from cache Error in h(simpleError(msg, call)) : error in evaluating the argument 'x' in selecting a method for function 'as.data.frame': object 'NCList_find_overlaps_in_groups' not found

Can you include the results of your sessionInfo() ? Please also try doing a BiocManager::valid() and BiocManager::install() to update any out of date packages. In some searches it appears this may be related to an out-dated GenomicRanges version.