Я решил проблему, изменив пакет с пророка на прогноз и построив график с помощью Plotly. Я думаю, что проблема была в dyplot.prophet, но с пакетами Plotly и Forecast проблема была решена. Я использовал сюжетный код из здесь
ts_cases = ts(mex$y,ts(mex$y, start=min(mex$ds), end = max(mex$ds), frequency=1))
fit <- auto.arima(ts_cases)
f_arima<- forecast(fit, 20) # forecast
# names
f <- list(
family = "Arial",
size = 18,
color = "#7f7f7f"
)
x <- list(
title = "Tiempo en dias",
titlefont = f
)
y <- list(
title = "Numero de casos confirmados",
titlefont = f
)
library(plotly)
plot_ly() %>%
add_lines(x = time(ts_cases), y = ts_cases,hoverinfo = "text",
color = I("black"), name = "observed") %>%
add_ribbons(x = time(f_arima$mean), ymin = f_arima$lower[, 2], ymax = f_arima$upper[, 2],
color = I("gray95"), name = "95% confidence") %>%
add_ribbons(x = time(f_arima$mean), ymin = f_arima$lower[, 1], ymax = f_arima$upper[, 1],
color = I("gray80"), name = "80% confidence") %>%
add_lines(x = time(f_arima$mean), y = f_arima$mean, color = I("blue"), name = "prediction")%>% layout(xaxis = x, yaxis = y,title= "ARIMA model")