[Bioc-sig-seq] rtracklayer and wig manipulation

[Simon Anders]

Michael Lawrence wrote:
> Oops, I meant:
> y <- rep(score(chr1), width(chr1))
> 
> 
>> x <- as.integer(unlist(ranges(chr1)))
>> lo <- loess(y ~ x)


But what if the intervals in the wiggle file contain gaps? That's
perfectly legal, after all.

I also fail to see what the effect of as.integer on IRanges is supposed
to be. For my wiggle file, I simply get

> head( as.integer( unlist( ranges(w) ) ) )
[1] 2 3 4 5 6 7


  Simon