I am interested in obtaining the methylated and unmethylated signal intensities. I am using the sesame package, and obtained the sigDF object as follows:
sdfs = openSesame(baseDirectory,prep="QCDPB",func=NULL,BPPARAM = BiocParallel::MulticoreParam(4),platform = "EPICv2")
I get, for the first sample something like this:
Probe_ID MG MR UG UR col mask \
cg00000029_TC21 NA NA 2709.775 1046.9992 2 FALSE
cg00000109_TC21 NA NA 2770.632 288.9767 2 FALSE
cg00000155_BC21 NA NA 3668.775 230.5298 2 FALSE
Can't find the definitions anywhere but I am guessing that MG, MR, UG and UR are what I am looking for (methylated green, methylated red, unmethylated green and unmethylated red). However, it seems that MG and MR do not have a lot of data:
sum(!is.na(sdf[["sample1"]]$MG))
[1] 128295
sum(!is.na(sdf[["sample1"]]$UG))
[1] 937688
Is there something that we are missing or does this look normal and correspond to the methylated and unmethylated intensities?