6.21 窗口平滑

6.21.1 问题

你想要计算序列给定窗口长度的平均数。这是平滑数据的一种方式。

6.21.2 方案

假设你有一个数值向量并且想要找出第一个连续四个数、第二个连续四个数等等如此的平均数（就像用一个窗口在序列上移动，移动一次，计算一次）。

  1. # 生成有22个（0-99）以内的随机数向量
  2. set.seed(123)
  3. x <- floor(runif(22) * 100)
  4. x
  5. #>  [1] 28 78 40 88 94  4 52 89 55 45 95 45 67 57 10 89 24
  6. #> [18]  4 32 95 88 69
  8. # 将向量长度取整到最接近的4的倍数
  9. newlength <- ceiling(length(x)/4) * 4
  10. newlength
  11. #> [1] 24
  13. # 给新增的数值取NA值
  14. x[newlength] <- NA
  15. x
  16. #>  [1] 28 78 40 88 94  4 52 89 55 45 95 45 67 57 10 89 24
  17. #> [18]  4 32 95 88 69 NA NA
  19. # 将它转换为4行的矩阵
  20. x <- matrix(x, nrow = 4)
  21. x
  22. #>      [,1] [,2] [,3] [,4] [,5] [,6]
  23. #> [1,]   28   94   55   67   24   88
  24. #> [2,]   78    4   45   57    4   69
  25. #> [3,]   40   52   95   10   32   NA
  26. #> [4,]   88   89   45   89   95   NA
  28. # 计算每一列的均值，忽略NA值
  29. colMeans(x, na.rm = TRUE)
  30. #> [1] 58.50 59.75 60.00 55.75 38.75 78.50