Open alexg9010 opened 1 year ago
library(methylKit)
data(methylKit)
mdb <- makeMethylDB(methylRawList.obj)
# works for "base" -resolution
unite(mdb,destrand=TRUE)
mdb # is fine
# breaks for "region" -resolution
mdbt <- tileMethylCounts(mdb)
# works
unite(mdbt)
# deletes tabix files of mdbt
unite(mdbt,destrand=TRUE)
mdbt # breaks !!
https://github.com/al2na/methylKit/blob/4735b6fca820c8bda1f84b03e12d19bada2df1fc/R/methylDBFunctions.R#L747