Trazando series de tiempo con ggplot2 y ggfortify

Instalación

Primero instalamos ggfortify

install.packages("ggfortify")

Traza ts objetos

ggfortify y ggplot2 sabe cómo interpretar ts objetos. Después de cargar ggfortify, puede usar la función ggplot2 :: autoplot para ts objetos.

library(ggfortify)

## Loading required package: ggplot2

## Warning: package 'ggplot2' was built under R version 3.4.2

autoplot(AirPassengers)

Para cambiar el color de línea y el tipo de línea, use las opciones ts.colour y ts.linetype. Use help (autoplot.ts) (o help (auto plot. *) Para cualquier otro objeto) para verificar las opciones disponibles.

autoplot(AirPassengers, ts.colour = "red", ts.linetype = "dashed")

Las series de tiempo multivariantes se dibujarán con facetas.

autoplot(AirPassengers)

autoplot(AirPassengers, ts.geom = "bar", fill = "blue")

autoplot(AirPassengers, ts.geom = "point", shape = 3, ts.colour = "green")

ggfortify admite las siguientes estadísticas relacionadas con series de tiempo en el paquete de estadísticas:

• stl, decomposed.ts
• acf, pacf, ccf
• spec.ar, spec.pgram
• cpgram (cubierto porggcpgram)

autoplot(stl(AirPassengers, s.window = "periodic"), ts.colour="blue")

NOTA Con acf y spec. *, Especifique plot = FALSE para suprimir el trazado predeterminado.

autoplot(acf(AirPassengers,plot = FALSE))

autoplot(pacf(AirPassengers, plot = FALSE))

Puede pasar algunas opciones al trazar acf a través de autoplot. El acf incorporado calcula el intervalo de confianza en el momento del trazado y no contiene el resultado; las opciones equivalentes se pueden pasar al trazado automático. El siguiente ejemplo muestra cómo cambiar el valor del intervalo de confianza y el método (use ma asumiendo que la entrada sigue al modelo MA).

autoplot(acf(AirPassengers, plot = FALSE), conf.int.fill = "#0000FF", conf.int.value = 0.8,
         conf.int.type = "ma")

autoplot(spec.ar(AirPassengers, plot = FALSE))

ggcpgram debe generar el periodograma acumulado como lo mismo que cpgram. Debido a que cpgram no tiene valor de retorno, no podemos usar autoplot (cpgram (...)).

Periodograma acumulado

ggcpgram(arima.sim(list(ar = c(0.7, -0.5)), n = 50))

ggtsdiag debería generar el diagrama similar a tsdiag.

library(forecast)

## 
## Attaching package: 'forecast'

## The following object is masked from 'package:ggfortify':
## 
##     gglagplot

ggtsdiag(auto.arima(AirPassengers))

gglagplot es para lag.plot.

gglagplot(AirPassengers, lags = 4)

ggfreqplot es un moth.plot generalizado. Puede pasar la freq si así lo desea, de lo contrario se usará la frecuencia de la serie temporal.

ggfreqplot(AirPassengers)

ggfreqplot(AirPassengers, freq = 4)

Trazado con paquete forecast

library(forecast)
arima1<-auto.arima(AirPassengers)
forecast1<-forecast(arima1,level = c(95), h = 50)
autoplot(forecast1)

### Opciones para cambiar la configuración básica.

autoplot(forecast1, ts.colour = "firebrick1", predict.colour = "red",
         predict.linetype = "dashed", conf.int = FALSE)

Plotting with vars package

library(vars)

## Warning: package 'vars' was built under R version 3.4.2

## Warning: package 'strucchange' was built under R version 3.4.2

## Warning: package 'urca' was built under R version 3.4.2

## Warning: package 'lmtest' was built under R version 3.4.2

data("Canada")
vselect1<-VARselect(Canada,lag.max = 5, type = "const")$selection[1]
var1<-VAR(Canada, p = vselect1, type = "const")

autoplot(predict(var1, n.ahead = 50), ts.colour = "dodgerblue4",
         predict.colour = "blue", predict.linetype = "dashed")

Trazado con el paquete changepoint

library(changepoint)

## Warning: package 'changepoint' was built under R version 3.4.2

## Successfully loaded changepoint package version 2.2.2
##  NOTE: Predefined penalty values changed in version 2.2.  Previous penalty values with a postfix 1 i.e. SIC1 are now without i.e. SIC and previous penalties without a postfix i.e. SIC are now with a postfix 0 i.e. SIC0. See NEWS and help files for further details.

autoplot(cpt.meanvar(AirPassengers))

autoplot(cpt.meanvar(AirPassengers), cpt.colour = "blue", cpt.linetype = "solid")

Trazado con el paquete strucchange

library(strucchange)
autoplot(breakpoints(Nile ~ 1), ts.colour = 'blue', ts.linetype = 'dashed',
         cpt.colour = 'dodgerblue3', cpt.linetype = 'solid')