代码之家 › 专栏 › 技术社区 › Will T-E

Cumsum,但数据点数量最大

moving-average r

Will T-E · 技术社区 · 7 年前

我想创造一个 cumsum() 和 TTR::runSum() 哪里 cumSum() 运行到预先指定数量的数据点,在这些点上,它的行为更像 runSum()

例如:

library(TTR)
data <- rep(1:3,2)
cumsum <- cumsum(data)
runSum <- runSum(data, n = 3)
DesiredResult <- ifelse(is.na(runSum),cumsum,runSum)

有没有一种不需要通过NAs获得资金就可以获得DesiredResult的方法?

2 回复 | 直到 7 年前

G. Grothendieck 7 年前

这就是 partial=TRUE 参数到 rollapplyr 做这里我们用 sum 还有 sd 和 IQR 。(请注意 sd公司 其中一个值是NA,我们选择IQR,因为它是可以计算标量的扩散度量,尽管在这种情况下它始终为0。)

library(zoo)

rollapplyr(data, 3, sum, partial = TRUE)
## [1] 1 3 6 6 6 6

rollapplyr(data, 3, sd, partial = TRUE)
## [1] NA 0.7071068 1.0000000 1.0000000 1.0000000 1.0000000

rollapplyr(data, 3, IQR, partial = TRUE)
## [1] 0.0 0.5 1.0 1.0 1.0 1.0

Julius Vainora 7 年前

这里有三种选择。

n <- 3

rowSums(embed(c(rep(0, n - 1), data), n)) # base R
# [1] 1 3 6 6 6 6

library(TTR)
runSum(c(rep(0, n - 1), data), n = n)
# [1] NA NA  1  3  6  6  6  6 # na.omit fixes the beginning

library(zoo)
rollsum(c(rep(0, n - 1), data), k = 3, align = "right")
# [1] 1 3 6 6 6 6

推荐文章

Amp · 使用R ggplot2删除geom_radial中axis.line和panel.border之间的空格

4 月前

Hard_Course · 用另一列中的值替换行的最后一个非NA条目

4 月前

Mark R · 使用geom_sf()删除地球仪上不需要的网格线

4 月前

Joe · 根据对工作日和本周早些时候的日期的了解,找到一个日期

4 月前

Ben · 统计向量中的单词在字符串中出现的频率

4 月前

TheCodeNovice · R中符号格式的尾随零和其他问题[重复]

4 月前

katefull06 · 在R中使用terra修改范围时,会为单独的SpatRaster重写范围

4 月前

dez93_2000 · 在R管道子功能中引用管道对象的当前状态

4 月前

accibio · 在ggplot2中为同一变量创建两个连续的颜色渐变比例

4 月前

Mankka · 如何在Ggplot2中绘制均匀的径向图

4 月前