r - R: 计算积分累计变化？

Question

有谁知道一个用一个清晰​​的例子计算 CUSUM 的包？我找到了绘制 CUSUM 的“qcc”。但是，我想要提交的数据范围的 CUSUM 值而不是图表。另一个新的好包“changepoint”，它针对不同的分布计算 CUSUM。我需要一个简单的包或函数，它可以根据 CUSUM- 的基本定义计算 CUSUM- 而无需任何参数：

As its name implies, CUSUM involves the calculation of a cumulative sum 
which is what makes it "sequential"). 
Samples from a process are assigned weights, 
and summed as follows:
For samples i = 0 ->infinity
   C0 = 0 :
   Ci+1 = max(0;Ci + Xi - mean(Xi) )

在Wiki中定义

任何计算基本 CUSUM 的包或函数？

Answer 1

为什么不简单for-loop？像这样，例如：

set.seed(20)
x <- rnorm(10, 1, 3)
mx <- mean(x)
CUMSUM <- function(x) {
    res <- numeric(length(x))
    for (i in seq_along(x) + 1) {
        res[i] <- max(0, res[i-1] + x[i-1] - mx)
    }
    res
}
CUMSUM(x)