I got the following error when running champ.DMR with DMRcate as the method. The package "IlluminaHumanMethylationEPICanno.ilm10b2.hg19" was installed successfully and can be loaded without problem. Can anyone help me resolve this? Thanks.
> res.DMRcate=champ.DMR(beta=BetaVal,
+ pheno=disease,
+ arraytype="450K",
+ method = "DMRcate",
+ minProbes=3,
+ adjPvalDmr=0.05,
+ cores=1,
+ rmSNPCH=F,
+ fdr=0.05,
+ mafcut=0.05,
+ lambda=1000,
+ C=2)
[===========================]
[<<<<< ChAMP.DMR START >>>>>]
-----------------------------
Loading required package: IlluminaHumanMethylation450kanno.ilmn12.hg19
1 cores will be used to do parallel DMRcate computing.
<< Find DMR with DMRcate Method >>
Loading required package: IlluminaHumanMethylationEPICanno.ilm10b2.hg19Error in get(ann) :
object 'IlluminaHumanMethylationEPICanno.ilm10b2.hg19' not found
Calls: champ.DMR ... cpg.annotate -> makeGenomicRatioSetFromMatrix -> get -> get
Execution halted
I just tested DMRcate function on a EPIC data on my computer, it seems works fine. Below is my code:
> myDMR <- champ.DMR(myFilter$beta,pheno=myFilter$pd$Sample_Group,method="DMRcate",arraytype="EPIC")
[===========================]
[<<<<< ChAMP.DMR START >>>>>]
-----------------------------
Loading required package: IlluminaHumanMethylationEPICanno.ilm10b2.hg19
3 cores will be used to do parallel DMRcate computing.
<< Find DMR with DMRcate Method >>
Your contrast returned 246813 individually significant probes. We recommend the default setting of pcutoff in dmrcate().
Fitting chr1...
Fitting chr10...
Fitting chr11...
Fitting chr13...
Fitting chr14...
Fitting chr16...
Fitting chr12...
Fitting chr17...
Fitting chr19...
Fitting chr15...
Fitting chr2...
Fitting chr21...
Fitting chr18...
Fitting chr4...
Fitting chr20...
Fitting chr7...
Fitting chr3...
Fitting chr22...
Fitting chr5...
Fitting chr6...
Fitting chr8...
Fitting chr9...
And seems your sessionInfo() is also correct. That's wired...Could you try directly load IlluminaHumanMethylationEPICanno.ilm10b2.hg19 package, and see if you can find object IlluminaHumanMethylationEPICanno.ilm10b2.hg19 in your R?
I just tried to load IlluminaHumanMethylationEPICanno.ilm10b2.hg19 package directly and it worked fine. But, when I tried to find the object using get("IlluminaHumanMethylationEPICanno.ilm10b2.hg19"), the same error came out saying that object 'IlluminaHumanMethylationEPICanno.ilm10b2.hg19' not found.
One thing I'd like to mention is that I got the error when running on HPC server (Linux). And everything worked fine on my laptop, where the same versions of R and ChAMP were installed.
Another thing I feel a little confused is that in the function I specified my array type to be "450K" and it seems that the annotation package for 450K was loaded successfully. I wasn't sure why it also loaded annotation package for EPIC array.
Hi Yuan,
Thanks for your reply.
I just tried to load IlluminaHumanMethylationEPICanno.ilm10b2.hg19 package directly and it worked fine. But, when I tried to find the object using get("IlluminaHumanMethylationEPICanno.ilm10b2.hg19"), the same error came out saying that object 'IlluminaHumanMethylationEPICanno.ilm10b2.hg19' not found.
One thing I'd like to mention is that I got the error when running on HPC server (Linux). And everything worked fine on my laptop, where the same versions of R and ChAMP were installed.
Another thing I feel a little confused is that in the function I specified my array type to be "450K" and it seems that the annotation package for 450K was loaded successfully. I wasn't sure why it also loaded annotation package for EPIC array.
Thanks for your help!