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.