0

我在使用 GGfortify 包的 2 个时间序列的单个图上遇到问题。

我在最后一行遇到了问题。是否可以进行更正?

library(Ecdat)
inflation <- as.numeric(Mishkin[, 1])

data_ts <- ts(inflation,start=c(1950,2), frequency = 12)

AR_data <- arima(data_ts, order = c(1,0,0))

AR_data_fitted <- data_ts - residuals(AR_data)

combined <- cbind(data_ts,AR_data_fitted)

gg <- autoplot(combined, columns = c(1,2),cmain = "Inflation Data")
4

1 回答 1

1

使用以下代码

library(ggfortify)

combined <- cbind(data_ts,AR_data_fitted)

colnames(combined) <- c("Original Data", "AR fitted data")

autoplot(combined,main = "Inflation Data",facets = FALSE)

在此处输入图像描述

更新

autoplot(combined,main = "Inflation Data",facets = FALSE)+
  scale_color_manual(values=c("blue", "black"))

在此处输入图像描述

于 2020-04-15T07:12:12.730 回答