library(readxl)
beer <- read_excel("Documents/Data Analysis and econometrics/beer.xlsx")
summary(beer)
## Month Monthly beer production
## Length:476 Min. : 64.8
## Class :character 1st Qu.:112.9
## Mode :character Median :139.2
## Mean :136.4
## 3rd Qu.:158.8
## Max. :217.8
##libraries needed##
library(forecast)
## Registered S3 method overwritten by 'quantmod':
## method from
## as.zoo.data.frame zoo
library(ggplot2)
##times series plot##
beerTS=ts(beer$`Monthly beer production`)
autoplot(beerTS)
##decompisition Additive##
beerTS=ts(beer$`Monthly beer production`,frequency = 12)
decompADD=decompose(beerTS,type="additive")
plot(decompADD)
##decomposition Multiplicative##
decompMult=decompose(beerTS,type = "multiplicative")
plot(decompMult)