我正在尝试自定义 stl 分解图的 x 轴。我的代码顺序是:
- 我有一个矩阵,其中包含索引、每个月的月底和观察到的幅度
矩阵样本
日期 | 震级 | |
---|---|---|
1 | 2018-02-28 | 8088 |
2 | 2018-03-31 | 10830 |
3 | 2018-04-30 | 10772 |
等等等等。
- 创建时间序列对象及其分解,
m_ts <- ts(MATRIX$magnitude, start = c(year(MATRIX$date[1]), month(year(MATRIX$date[1])), frequency = 12)
decom <- stl(m_ts, s.window = "periodic", 12))
这是我发现难题的地方:
- 如果我使用plot():
plot(decom)
我得到了情节,但日期是十进制格式。
- 而,如果我使用autoplot():
decom_plot <- autoplot(decom) + scale_x_date(date_breaks = "4 month", date_labels = "%b-%Y") + theme(axis.text.x = element_text(angle = 90))
decom_plot
它允许我自定义日期标签,但绘图的其余部分具有与 plot() 不同的线型。但是,因为是余数,所以我认为最好将线型保持为广告垂直线。
那么,如果使用绘图,我如何自定义 x 轴标签,或者如果使用自动绘图,我该如何更改绘图类型。
在众多方法中,我尝试过该方法axis(1, xaxp = date)
,但“参数 xaxp 有错误的一面”。有什么建议吗?
预先感谢。