Skip to main content
1 de 3
David
  • 195
  • 10

Dudas con Forecast y auto.arima predicción

soy nuevo en la comunidad y con R. Estoy intentando hacer una predicción sobre el número de casos de covid en mi municipio en base a la fecha, pero me está dando problemas y no acabo de entender por qué. Espero puedan ayudarme.

Adjunto el código que hasta ahora hice:

# Carga datos:
datos_covid<-read.delim(file.choose(),header = T, sep = ";", fill=TRUE)
View(datos_covid)
datos2<-datos_covid


# Hacemos plot(), para ver el grafico:
plot(datos2)

# Cargamos paquete de datos ggplot:
install.packages("ggplot2")
library(ggplot2)

# Anadimos mas detalles al grafico:
ggplot(data=datos2,
       aes(x=FECHA,
           y=NUMCASOS, group=1))+
  geom_line(colour="red")+
  labs(x="FECHA", y="NUMERO CASOS",
       title="CASOS COVID MUNICIPIO POR FECHA")

class(datos2)

# Cargamos libreria Forecast:
install.packages("forecast")
library(forecast)

# Ajustamos los datos al eje de las X. Es decir, especificamos la columna del NUMERO DE CASOS,   que es la que queremos predecir:
datos3<-datos2$NUMCASOS
View(datos3)

# La data sera el eje de las Y. Es decir, convertir unos datos en clase temporal:
install.packages("zoo") 
require(zoo)

# Indicamos el inicio y el final de los datos en la grafica, para predecir:
start=as.Date("2020-02-01")
start=as.yearmon(start)
end = as.yearmon(as.Date("2020-12-31"));end

datos4<-ts(datos3,start=start,end=end,frequency=12)

# Ver que aspecto tiene el objeto
datos4
ts.plot(datos4)

# Predicciones:
mod1<-auto.arima(datos4)

# Se coge el 75% de los primeros datos (DATOS DE ENTRENAMIENTO) 
train_datos<-datos4[1:573]
View(train_datos)

prediccion <- forecast(train_datos, h=10, model=mod1)
plot(prediccion)

# Cogemos los registros restantes.
lines(x=seq(573, 839, 1), y=datos4[573:839], col="red") 

Al realizar este último paso, hay datos que me salen en blanco (NA). Hice la prueba con otros datos y una menor cantidad de observaciones (unas 30, más o menos) y no tuve este problema.

Así es como se muestran los datos csv, hasta 839 obs.

FECHA SEXODESCRIPCION SEXOCODIGO NUMCASOS 09/03/2020 Hombre 0 1 10/03/2020 Mujer 1 3 10/03/2020 Hombre 0 4 11/03/2020 Hombre 0 2 11/03/2020 Mujer 1 2 12/03/2020 Mujer 1 4

Espero que puedan ayudarme. Gracias por adelantado. Saludos.

David
  • 195
  • 10