Displaying 1 result from an estimated 1 matches for "644730".
Did you mean:
644430
2012 Oct 16
2
sliding window analysis with rollapply
Dear List members
I want to do the sliding window analysis of some specific values. Here is my code:
require(zoo)
dat <- read.table("chr1.txt", header = TRUE, sep="\t")
dat2 <- cbind(dat[1,3]) #The first column is also important. It represents the position of the site on the chromosome.
TS <- zoo(c(dat2))
a <- rollapply(TS, width=1000000, by=200000, FUN=mean,