minfi::getQC - badsamplecutoff 10.5
0
0
Entering edit mode
23 days ago
June • 0

Hello,

I just started working on analyzing methylation EPIC array data using minfi. Data is generated on EPICv2 array and I'm using updated EPICv2 annotation (rgset@annotation <- c(array = "IlluminaHumanMethylationEPICv2", annotation = "20a1.hg38") Using the getQC and plotQC function. I noticed that all of samples are below 10.5 cut-off (lowest value is 9.1 and highest value is 10.48). Beta values density plot looks normal and detection p-value also looks normal. I'm confused on how should I interpret this data?

Thanks

getqc minfi methylation EPICv2 • 220 views
ADD COMMENT
0
Entering edit mode

I've recently also analyzed EPICv2 and have noticed that the log2 median intensities (and particularly the "uMed" channel), for an otherwise good QC cohort, were lower than what I had seen in the past for EPICv1 in general. But my N for this observation was N = 1 EPICv2 cohort...

In my opinion if in your plotQC scatter plot your samples cluster together, and if the minfi::qcReport control probes look OK, this is more of a "general batch" thing and would proceed (only filtering out outliers further away from the cluster of points) but it would be nice to get feedback from more people analyzing EPICv2 to see if it's typical to get less absolute signal intensities.

ADD REPLY

Login before adding your answer.

Traffic: 1817 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6