[BioC] NucleR: processReads/solveUserSEW0 - Error
Oscar Flores
oflores at mmb.pcb.ub.es
Thu Nov 3 13:16:39 CET 2011
Hi Stefanie,
I'm the developer of nucleR, so let's see if I can help you ;)
After the processing, processReads converts the input data to a RangedData
object for a easier manipulation later, so this error is occurs at the last
step of the call, but data can be messed in previous steps. It's
hard to tell what is happening without having a look to the input data,
which I guess is huge...
I would like to have a look to the raw data, but I know it is difficult
to send it if it's not in a public repository. Maybe you can contact me
directly about that (oflores at mmb.pcb.ub.es)...
Meanwhile, if you want to try a workaround, you can directly convert the
imported reads to RangedData format (which is the other format supported
by processReads):
(being "ar" your imported AlignedReads object)
res = RangedData(IRanges(start=position(ar),width=width(ar)),
strand=strand(ar),space=ar at chromosome)
reads_pair = processReads(res, type="paired", fragmentLen=fragment_len)
This should work, but will be nice to have a look to your data
to fix a possible problem in the AlignedReads method.
Regards,
Oscar
More information about the Bioconductor
mailing list