Describe R/Notes here.

data smoothing:

loess(), runmed()

track <- import("foo.wig") scores <- score(track)

chr1 <- trackchr1 # limit to single chromosome

y <- rep(score(chr1), width(chr1))

x <- as.integer(unlist(ranges(chr1))) lo <- loess(y ~ x)