在看到《R數(shù)據(jù)科學(xué)》第45頁時候，突然看到一句話
image.png
盡管這里并沒有對RcppRoll包進行描述，但是經(jīng)不住好奇去查了一下這個包
install.packages("RcppRoll")
library(RcppRoll)
(x <- 1:10)
[1] 1 2 3 4 5 6 7 8 9 10
> roll_sum(x, n = 3, by = 3)
[1] 6 15 24
> roll_sum(x, n = 3, by = 2)
[1] 6 12 18 24
> roll_sum(x, n = 3, by = 1)
[1] 6 9 12 15 18 21 24 27
# 看參數(shù)說明
n 滾動求和窗口的大小
by 表示每次窗口移動的距離
這里我們可以用來計算BS-seq數(shù)據(jù)分析中每個bin里面的值们豌，比如100bp的bin, 滑動窗口為50bp
即設(shè)置為
roll_sum(x, n = 100, by = 50)
