Re: ChIPpeakAnno: Invalid row.names length in peaksNearBDP
"Zhu, Lihua (Julie)" <Julie.Zhu-3WprALB+yeL2fBVCVOL8/[email protected]> Fri, 29 Aug 2014 20:12:31 +0000
| Newsgroups | gmane.science.biology.informatics.conductor |
|---|---|
| Message-ID | <D02656F9.BE48%[email protected]> |
Jens, Thanks for reporting the bug in function peaksNearBDP! Jianhong has fixed it and the new version is 2.14.1 (revision 93700). Best regards, Julie On 8/26/14 2:42 AM, "Preussner, Jens" <Jens.Preussner-6GLbxAiBP3Zo1qOY/[email protected]> wrote: Dear Julie, I noticed an unwanted breaking of the peaksNearBDP function in your package. When searching for bidirectional promoters in annotationData that contains just one entry, the function gives an “invalid row.names error”. Is this intended? Just for completeness and reproducability, I used the following code: library("ChIPpeakAnno") peaks = BED2RangedData(file.path("test.bed"),header=F) annotation = GFF2RangedData(file.path("ref.gtf"),header=F,comment.char="#",sep="\t") ovl_data = peaksNearBDP(peaks, AnnotationData = annotation, MaxDistance = 3000, PeakLocForDistance = "middle", FeatureLocForDistance = "TSS") The two files are attached. Thanks for your time and best regards, Jens [[alternative HTML version deleted]] _______________________________________________ Bioconductor mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/bioconductor Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor