methylKit
methylKit copied to clipboard
dbdir= argument from tilemethylCounts(), does not work with methylRawListDB data
Hi Alex,
during my analysis I came across the issue that the tileMethylCounts() function's argument "dbdir=" does not work in combination with a methylRawListDB object. The output will rather be produced to the working directory of the script. Below you find a reproducible example:
library(methylKit)
data(methylKit)
output_path_issue <- "/local/arathge/"
getwd()
tiles.issue <- tileMethylCounts(
object = methylRawList.obj,
win.size = 500,
step.size = 500,
save.db = TRUE,
suffix = "tiled",
dbdir = file.path(output_path_issue)
)
Hi @AnjaCaRa,
thanks for reporting this issue and providing the example code! I will have a closer look at the problem.
Best, Alex