Entering edit mode
Hi,
I am working with EPIC v2 methylation data and trying to read the data using champ.import function. But i am getting the below error message.
data = champ.import(directory = "~/Documents/EPIC_may/IDATs/", arraytype = "EPICv2")
[===========================]
[<<<< ChAMP.IMPORT START >>>>>]
-----------------------------
[ Section 1: Read PD Files Start ]
CSV Directory: /Users/medha/Documents/EPIC_may/IDATs//pd_MS.csv
Find CSV Success
Reading CSV File
Replace Sentrix_Position into Array
Replace Sentrix_ID into Slide
[ Section 1: Read PD file Done ]
[ Section 2: Read IDAT files Start ]
Loading:/Users/medha/Documents/EPIC_may/IDATs//207218120055_R01C01_Grn.idat ---- (1/96)
Loading:/Users/medha/Documents/EPIC_may/IDATs//207218120055_R02C01_Grn.idat ---- (2/96)
.
.
.
Loading:/Users/medha/Documents/EPIC_may/IDATs//207215830010_R08C01_Red.idat ---- (96/96)
Extract Mean value for Green and Red Channel Success
Your Red Green Channel contains 1105209 probes.
[ Section 2: Read IDAT Files Done ]
[ Section 3: Use Annotation Start ]
Reading EPIC Annotation >>
Fetching NEGATIVE ControlProbe.
Totally, there are 411 control probes in Annotation.
Your data set contains 411 control probes.
Generating Meth and UnMeth Matrix
Extracting Meth Matrix...
Totally there are 865918 Meth probes in EPIC Annotation.
Your data set contains 865918 Meth probes.
Extracting UnMeth Matrix...
Totally there are 865918 UnMeth probes in EPIC Annotation.
Your data set contains 865918 UnMeth probes.
Error in champ.import(directory = "~/Documents/EPIC_may/IDATs/", arraytype = "EPICv2") :
Meth Matrix and UnMeth Matrix seems not paried correctly.
sessionInfo( )
R version 4.2.3 (2023-03-15)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS Monterey 12.0.1
Any Idea/suggestions on how to solve this.
Thanks, Medha.
Have you checked the answers here ? : Problem at champ.load
Hi, I have also the same problem. Any hints on how to resolve it? thanks
Did you read the post that Basti linked to?
I have exactly the same problem and I never had this problem with Epic V1. The answer posted in the link is not the same problem as the one we are having. That problem is that there were different numbers of Unmeth and meth probes. This problem is that the numbers are the same but their annotation is different. This is likely down to EPIC v2 probe annotation and it may be that we can no longer use ChAMP pipeline on it without a fix
The
ChAMP
pipeline doesn't recognize the difference between the V1 and V2 arrays. This is also true forminfi
, but it's simple to fix that situation. You will have to install the annotation packages that Zuguang Gu generated usingAnd then change the annotation of your
RGset
after you read in the data.Thank you very much for this comment. The annotation packages that Zuguang Gu generated worked well for me in the minfi pipeline, but not in ChAMP. Have someone managed to run it in ChAMP?
I'm having same issues too. Did anyone manage to resolve it?