Browse code

more efficient in lengths

Ge Tan authored on 20/09/2016 11:19:20
Showing 1 changed files

... ...
@@ -156,7 +156,7 @@ read.rmskFasta <- function(fn){
156 156
   foo3 <- lapply(lapply(strsplit(as.character(seq),"") ,"%in%",
157 157
                        c("a","c","g","t")), Rle)
158 158
   foo4 <- lapply(foo3, as, "IRanges")
159
-  foo5 <- GRanges(seqnames=Rle(names(foo4), sapply(foo4, length)),
159
+  foo5 <- GRanges(seqnames=Rle(names(foo4), lengths(foo4)),
160 160
                   ranges=IRanges(start=unlist(sapply(foo4, start)),
161 161
                                  end=unlist(sapply(foo4, end))),
162 162
                   strand="+")