library(TTR)
## Warning: package 'TTR' was built under R version 4.5.1
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.5.1
library(forecast)
## Warning: package 'forecast' was built under R version 4.5.1
## Registered S3 method overwritten by 'quantmod':
## method from
## as.zoo.data.frame zoo
library(readxl)
## Warning: package 'readxl' was built under R version 4.5.1
datos2=read_excel("datos2_examen1.xlsx")
datos3=read_excel("datos3_examen1.xlsx")
d2s=ts(datos2,start = c(2011,1),frequency = 12)
d3s=ts(datos3,start = c(2018,1),frequency = 12)
(knitr::kable(datos2)) para mostrar conjunto de datos
plot(d2s,xlab="Años",ylab="Numero de Nacimientos",col="darkblue")
interpretacion de datos 2 tiene tendencia desde el año 2014,es aditivo,tiene factor estacional,su punto max es de 30:00 y su minimo de 20:00
plot(d3s,xlab="Años",ylab="Ventas",col="darkred")
interpretacion de datos 3 tiene tendencia, es multiplicativo, tiene la parte estacional,las ventas maximas mostradas son de 104661 y su ventas mas bajas son de 1665.
d2sd=decompose(d2s)
plot(d2sd)
plot(d2sd$trend)
d3sd=decompose(d3s)
plot(d3sd)
plot(d3sd$trend)
?HoltWinters
## starting httpd help server ... done
d2se=HoltWinters(d2s,seasonal = "additive")
d2sp=predict(d2se,8)
(d2sp)
## Jan Feb Mar Apr May Jun Jul Aug
## 2025 27.30020 25.92822 29.04779 27.56001 28.88793 28.50749 30.56943 30.55133
plot(d2se,d2sp)
summary(d2s)
## datos2
## Min. :20.00
## 1st Qu.:23.28
## Median :24.96
## Mean :25.06
## 3rd Qu.:26.88
## Max. :30.00
d3se=HoltWinters(d3s,seasonal = "multiplicative")
d3sp=predict(d3se,8)
knitr::kable(d3sp)
| fit |
|---|
| 14994.74 |
| 18948.09 |
| 35767.05 |
| 28159.43 |
| 25669.26 |
| 31472.70 |
| 40496.45 |
| 41217.96 |
plot(d3se,d3sp)
summary(d3s)
## datos3
## Min. : 1665
## 1st Qu.: 5884
## Median : 8772
## Mean : 14316
## 3rd Qu.: 16889
## Max. :104661