我有这个矩阵:
矩阵 =
我想要的：
roll_aver=
我知道我可以用这个来做滚动平均：
roll_aver <- t(apply(matrix, 1, function(x) rollmean(x, k=4, fill = 0, align ="right")))
但通过这种方式，它得到了一个包含零的滚动平均值。
想象一下，它相当于做一个 excel：averageif（行中最后 4 个的范围，
我们可以通过在行上循环来从左侧第一个非零元素应用＜code＞rollapplyrollapply使用FUN，在获得均值时排除零元素
library(zoo)
t(apply(matrix, 1, FUN = function(x) {
i1 <- cumsum(x != 0) > 0
x[i1] <- rollapply(x[i1], width = 4, fill = 0, align = "right",
FUN = function(u) mean(u[u!=0]), partial = TRUE)
x}))
-输出
value1 value2 value3 value4 value5 value6 value7 value8 value9 value10
1 0 0 45 25.5 19.0 16.25 5.75 5.666667 5 4.0
2 0 0 0 2.0 5.5 4.00 5.00 5.250000 4 5.5
matrix <- structure(c(0L, 0L, 0L, 0L, 45L, 0L, 6L, 2L, 6L, 9L, 8L, 1L,
3L, 8L, 0L, 3L, 4L, 0L, 5L, 0L), .Dim = c(2L, 10L), .Dimnames = list(
c("1", "2"), c("value1", "value2", "value3", "value4", "value5",
"value6", "value7", "value8", "value9", "value10")))
