使用rollapply计算滚动平均时遇非NA值错误求助

看起来你遇到的问题是因为直接对POSIXct格式的时间序列使用rollapply+fill='extend'计算均值导致的——底层的approx函数无法直接处理时间类型的数据，而且时间类型的均值计算需要先转换成数值格式再处理。

错误原因拆解 #

你的ks[,1]是POSIXct类型的时间数据：

[,1]
2014-02-14 14:27:00 "2014-02-14 14:27:00"
2014-02-14 14:32:00 "2014-02-14 14:32:00"
...

当你调用rollapply(ks[,1], 2, mean, align='right', fill='extend')时：

1. mean()虽然能自动把POSIXct转成数值（从1970-01-01开始的秒数）计算均值，但fill='extend'参数会触发approx函数做插值，而approx不支持直接处理POSIXct对象，因此抛出错误。
2. 直接对时间类型做滚动插值的逻辑本身也不符合数值插值的预期，必须先转成数值再处理。

解决方案 #

我们可以先把时间转成数值格式（秒数）计算滚动平均，再转回POSIXct格式，具体有两种方式：

方式一：使用rollapply分步处理

library(zoo)

# 1. 将POSIXct时间转为数值（从1970-01-01起的秒数）
ks_numeric <- as.numeric(ks[,1])

# 2. 计算滚动平均，保留fill='extend'参数
avg_numeric <- rollapply(ks_numeric, 2, mean, align='right', fill='extend')

# 3. 将数值结果转回POSIXct格式
avg <- as.POSIXct(avg_numeric, origin = "1970-01-01")

方式二：使用更高效的rollmean（zoo包专门函数）

如果只是简单的固定窗口滚动平均，rollmean比rollapply更高效，写法也更简洁：

library(zoo)

avg <- as.POSIXct(
  rollmean(as.numeric(ks[,1]), 2, align='right', fill='extend'),
  origin = "1970-01-01"
)

额外说明 #

• 如果你不需要fill='extend'来填充开头的缺失值，可以直接去掉该参数，此时结果会比原数据少1个元素（因为窗口大小为2），也不会触发approx的错误：

  avg <- as.POSIXct(rollmean(as.numeric(ks[,1]), 2, align='right'), origin = "1970-01-01")
  

• 确保你已经正确加载了zoo包，没有的话先运行install.packages("zoo")安装。

内容的提问来源于stack exchange，提问作者Liondancer