Почему ggfortify терпит неудачу с этим графиком временного ряда? - PullRequest
0 голосов
/ 15 апреля 2020

У меня проблемы с одним графиком из двух временных рядов с использованием пакета GGfortify.

Я столкнулся с проблемой в последней строке. Есть ли исправление, которое можно сделать?

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")

1 Ответ

1 голос
/ 15 апреля 2020

Используйте следующий код

library(ggfortify)

combined <- cbind(data_ts,AR_data_fitted)

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

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

enter image description here

Обновление

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

enter image description here

...