#acl ChrisSeidel:read,write,delete,revert All:read

Describe R/Notes here.

data smoothing:

loess(), runmed()

track <- import("foo.wig")

scores <- score(track)

chr1 <- track["chr1"] # limit to single chromosome

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

x <- as.integer(unlist(ranges(chr1)))

lo <- loess(y ~ x)