I am using MethyKit package for WBGS data analysis. There it has a step like
meth <- unite(myobj, destrand=FALSE) uniting...
Where it unites, if I put destrand=TRUE, then it gives this error
meth=unite(myobj, destrand=TRUE) destranding... Error in
*tmp*, "strand", value = "+") : replacement has 1 row, data has 0
My question is what does it mean? Does methylKit has any option to take destrand=TRUE? I can see that the usage is
Usage unite(object, destrand = FALSE, min.per.group = NULL, chunk.size = 1e+06, mc.cores = 1, save.db = FALSE, ...)