#install.packages("forecast")
library(forecast)
myvector = sample(1:45, 13, replace=F)
mydata <- ts(myvector, start=c(2020, 1, 1), end=c(2022, 12, 10), frequency=12)
trend_mydata = ma(mydata, order = 4, centre = T)
plot(as.ts(mydata))
lines(trend_mydata)
plot(as.ts(mydata))
detrend_mydata = mydata - trend_mydata
plot(as.ts(detrend_mydata))
m_mydata = t(matrix(data= detrend_mydata, nrow=4))
seasonal_mydata = colMeans(m_mydata, na.rm = T)
plot(as.ts(rep(seasonal_mydata, 16)))
random_mydata = mydata - trend_mydata - seasonal_mydata
plot(as.ts(random_mydata))
# Dawid Biały, Bartłomiej Bartoszek, Daniel Sobczak