methylKit icon indicating copy to clipboard operation
methylKit copied to clipboard

dbdir= argument from tilemethylCounts(), does not work with methylRawListDB data

Open AnjaCaRa opened this issue 3 years ago • 1 comments

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)
) 

AnjaCaRa avatar Oct 10 '22 08:10 AnjaCaRa

Hi @AnjaCaRa,

thanks for reporting this issue and providing the example code! I will have a closer look at the problem.

Best, Alex

alexg9010 avatar Oct 10 '22 09:10 alexg9010