如何在R中更改ACF图的滞后间隔？实现整数滞后显示

问题根源 #

你在创建时间序列对象时设置的frequency = 184是错误的。这个参数的作用是定义单位时间内的观测数量：当设为184时，R会将时间单位识别为“年”，1天的滞后就等于1/184≈0.005年，因此ACF图的横轴会显示以年为单位的小数，而非你需要的整数天数。

两种解决方案 #

方案一：修正ts对象的frequency参数

对于每日降水数据，应将frequency设为一年的天数（通常用365），这样R会把最小时间间隔视为1天，滞后值就会显示为整数天数：

precip.dat <- read.csv("Precipitation.csv") # 导入CSV文件
precip <- precip.dat$Precipitation # 提取降水数据
# 修正frequency为365，匹配每日观测的时间单位
precip.ts <- ts(precip, start = c(2010,5), end = c(2022,10), frequency = 365)
acf(precip.ts, lag.max = 7) # 绘制ACF图

方案二：直接使用原始向量计算ACF

如果不需要时间序列的时间标签，直接对原始降水向量计算ACF即可，此时滞后代表观测点的间隔（即天数），自然显示为整数：

precip.dat <- read.csv("Precipitation.csv")
precip <- precip.dat$Precipitation
# 直接对向量计算ACF，无需转换为ts对象
acf(precip, lag.max = 7)

注意：如果数据存在缺失值，方案二的处理更直接，因为ts对象会自动填充缺失值为NA，但这不会影响ACF的计算结果。

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