0

我正在尝试自定义 stl 分解图的 x 轴。我的代码顺序是:

  1. 我有一个矩阵,其中包含索引、每个月的月底和观察到的幅度

矩阵样本

日期 震级
1 2018-02-28 8088
2 2018-03-31 10830
3 2018-04-30 10772

等等等等。

  1. 创建时间序列对象及其分解
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))
  1. 这是我发现难题的地方:

    • 如果我使用plot()
plot(decom)

我得到了情节,但日期是十进制格式。

情节看起来就像我想要的那样,但我一直在努力改变 x 轴标签,但没有成功。

  • 而,如果我使用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 轴标签,但其余图的线类型错误

那么,如果使用绘图,我如何自定义 x 轴标签,或者如果使用自动绘图,我该如何更改绘图类型。

在众多方法中,我尝试过该方法axis(1, xaxp = date),但“参数 xaxp 有错误的一面”。有什么建议吗?

预先感谢。

4

0 回答 0