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 _______________________________________________ Bioc-sig-sequencing mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
