sartorlab / methylSig

R package for DNA methylation analysis
17 stars 5 forks source link

merging percent methylation with methylation differences output #19

Closed shawpa closed 4 years ago

shawpa commented 6 years ago

This is not a bug but I was hoping you could assist me. I am trying to get the percent methylation values appended to the methylSigDiff output. I know I can generate the the percent methylation with the command x = meth[,"numCs"]/meth[, "coverage"] colnames(x) = meth@sample.ids

However, when I attempt to append the chromosomal positions to that, I think they are somehow being sorted because they no longer match with the coordinates from the methylSigDiff output. I am sure the error is on my end and not on yours. I am not experienced dealing with S4 data objects so it's difficult for me to figure this out. Any assistance would be greatly appreciated.