R语言hts包绘制分层时间序列时x轴起始时间异常问题

问题根因 #

该问题是hts包原生plot方法的已知实现缺陷：当同时绘制多个层级序列时，底层在聚合不同层级序列的时间属性过程中，frequency = 365.25/7 带来的浮点精度误差会导致时间序列的tsp属性（起始时间、结束时间、频率）计算错误，因此出现x轴起始时间退回2016年的异常；仅绘制level 0层级时不需要做跨层级的时间属性合并，所以显示正常。

可选用的解决方案 #

• 方案1：优先使用forecast包适配的autoplot方法绘制，自动继承正确的时间属性，稳定性更高

library(forecast)
# 直接指定需要绘制的层级即可
autoplot(x, levels = c(0,2)) +
  labs(x = "时间", y = "序列值")

• 方案2：调整时间序列频率为整数值，从根源避免浮点精度误差，周频场景下取整为52的误差可忽略

# 构造ts对象时将frequency参数改为整数52
abc <- ts(abc, start=2019, frequency=52)
x <- hts(abc, characters = c(1,2))
# 此时再调用hts原生plot方法不会出现时间偏移
plot(x, levels = c(0,2))

• 方案3：如果需要保留原浮点频率参数，可手动覆盖x轴刻度

# 绘图时关闭自动生成的x轴
plot(x, levels = c(0,2), xaxt = "n")
# 手动添加正确的时间刻度
axis(1, at = seq(2019, 2022, 1), labels = seq(2019, 2022, 1))

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