如果问题是要解释代码,那么它就是取向量 x 的长度 n 的移动平均值。例如,如果没有 NA 且 n=2,则输出的前几个元素为 (x[1] + x[2])/2、(x[2] + x[3])/2 等.
n <- 2
x <- c(1, 3, 4, 7, 9)
cx <- c(0, cumsum(ifelse(is.na(x), 0, x))) # 0 1 4 8 15 24
cn <- c(0, cumsum(ifelse(is.na(x), 0, 1))) # 0 1 2 3 4 5
rx <- cx[(n+1):length(cx)] - cx[1:(length(cx) - n)] # 4 7 11 16
rn <- cn[(n+1):length(cx)] - cn[1:(length(cx) - n)] # 2 2 2 2
rsum <- rx / rn # 2.0 3.5 5.5 8.0
cx 为 0,后跟 x 的累积和,但在计算累积和时将 NA 替换为 0。
cn 为 0,后跟非 NA 的累积数。
rx 是累积总和减去 n 个位置的累积总和。
rn 是非 NA 的数量减去 n 个位置的非 NA 的数量。
rsum 是最后两个的比率。