Question

Error using champ.import

0

Entering edit mode

medha.suman • 0

@5f773e94

Last seen 10 months ago

Sweden

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.

ChAMP EPICv2 • 1.7k views

ADD COMMENT • link updated 12 days ago by shilpa • 0 • written 10 months ago by medha.suman • 0

0

Entering edit mode

Have you checked the answers here ? : Problem at champ.load

ADD REPLY • link 10 months ago Basti ▴ 780

0

Entering edit mode

Hi, I have also the same problem. Any hints on how to resolve it? thanks

ADD REPLY • link 8 months ago Fulvio • 0

0

Entering edit mode

Did you read the post that Basti linked to?

ADD REPLY • link 8 months ago James W. MacDonald 65k

0

Entering edit mode

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

ADD REPLY • link 5 months ago wellitonboots • 0

1

Entering edit mode

The ChAMP pipeline doesn't recognize the difference between the V1 and V2 arrays. This is also true for minfi, but it's simple to fix that situation. You will have to install the annotation packages that Zuguang Gu generated using

library(BiocManger)
install("jokergoo/IlluminaHumanMethylationEPICv2manifest")
install("jokergoo/IlluminaHumanMethylationEPICv2anno.20a1.hg38")

And then change the annotation of your RGset after you read in the data.

ADD REPLY • link 5 months ago James W. MacDonald 65k

0

Entering edit mode

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?