I have been trying out your package on a set of bam files covering a relatively small region that should contain a fair number of aberrant splicing events. I have been following the vignette closely and have successfully derived a "counts" object from which I can extract gene, bin, and junction info. The counts, features, bam, and target objects all appear to be working correctly. However, when I run:
as <- AsDiscover( counts, targets, features, bam, readLength=150L,threshold = 5)
I get the following error message:
Error in aggregate.data.frame(mf[1L], mf[-1L], FUN = FUN, ...) :
no rows to aggregate
I am going to look at a genomic region that is a bit larger to see if that works for me but in the meantime, if you can offer any advice it would be much appreciated.
I am doing pretty much the exact same thing as Jon and get the same error message. Do you have any idea of how to fix the error? I can provide more details if necessary.
Hello!
I have been trying out your package on a set of bam files covering a relatively small region that should contain a fair number of aberrant splicing events. I have been following the vignette closely and have successfully derived a "counts" object from which I can extract gene, bin, and junction info. The counts, features, bam, and target objects all appear to be working correctly. However, when I run:
as <- AsDiscover( counts, targets, features, bam, readLength=150L,threshold = 5)
I get the following error message:
Error in aggregate.data.frame(mf[1L], mf[-1L], FUN = FUN, ...) : no rows to aggregate
I am going to look at a genomic region that is a bit larger to see if that works for me but in the meantime, if you can offer any advice it would be much appreciated.
Thanks very much in advance!
Jon