This is a test version of Biostars. For the public version, visit https://www.biostars.org.
"MethylKit" package for WGBS data

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 $<-.data.frame(*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, ...)

next-gen sequencing
• 1,001 views • link
written by shrinka.genetics • shrinka.genetics

0 answers

No answers yet.

Log in to answer this question.

More posts like this