Case kiểm tra hiệu quả của Time series; This case is cheking performance of time series

Chuẩn bị data, prepare data

Bcase = scan("C:/Users/tien.vo/Desktop/Paulus Cicero/BCase/TSBcase4years.txt")
ts.Bcase = ts(data = Bcase, start = c(2017,1), frequency = 52, deltat = 1/52)
ts.Bcase
## Time Series:
## Start = c(2017, 1) 
## End = c(2020, 5) 
## Frequency = 52 
##   [1]  4614648  4277881  6350436  8549914  6901816  6851594  6883294  7725983
##   [9]  6474413  6611395  5801667  6737096  8766446  7981651  6119259  5492393
##  [17]  6440131  6631434  5302194  5610582  7872824 10212149 10067921  9752701
##  [25] 10971478  9556546  8926904  9056203 10397286 11168614  9366451  9895920
##  [33]  9709980 11885922 12239436 11254488  7815029  7972991  8428031  7915513
##  [41]  8400323  7845750  8408058  7925964  7898994  6848700  7258329  6571322
##  [49]  6658055  6641390  6096181  6532440  7169074  6396842  7933800  6306913
##  [57]  7233875  6069921  7189243  7149594  6949889  6811420  6496747  7355842
##  [65]  8017409 10287638 12763221 14382594 11165126  9620406 10294388 13400454
##  [73] 13689181 13045490 16795020 19643486 18335780 16174950 13929293 12996372
##  [81] 14273586 12390835 17471701 18321583 16648429 13984175 15677681 13921935
##  [89]  8933402 11445061 11151125  9910043 11686773  9738036 10454577 10963451
##  [97]  9455539  8756873 11903391 12162217 11797065 11457483 11117191 11203548
## [105] 11706251 11432017  9614850  9184222  9184222 10209429 10734420 10279031
## [113] 12862144 12239622 11412318 12044220 11323803 10345438  9496123  9479539
## [121]  8960258  9791516 11783355 14455359 16462009 13707988 12163750 13550706
## [129] 10922921 15502983 12427797 11626119 12805941 14970189 14422308 14450597
## [137] 12808338  9886357  8612896 10781416  5643496  5666656  6530793  6560657
## [145]  5549700  5766594  5703512  6931988  5874779  6389493  6693185  3908297
## [153]  5720620  6060855  7458445  5911833  5231868  5903372  7462769  6298243
## [161]  8613661

Decompose time series data:

comp = decompose(ts.Bcase)
library(ggplot2)
library(ggfortify)
autoplot(comp)

Chúng ta có thể thấy chỉ số quan tâm cao dần từ tháng 1 đạt đỉnh vào cuối tháng 4, đỉnh phụ tháng 6 và xuống dần

Trend tốt dần đến cuối năm 2018 thì có xu hương giảm

Sai số ngẫu nhiên độc lập với thời gian (remainder)

Smooth data; làm mềm dữ liệu cho dễ nhìn

library(TTR)
sma4 = SMA(ts.Bcase, n=4)
sma5 = SMA(ts.Bcase, n=5)
par(mfrow=c(1,2))
plot.ts(sma4); plot.ts(sma5)

Predict the not very clearly seasonal data; Tiên lượng dữ liệu với đặc tính mùa chưa rõ ràng

library(forecast); library(fpp2) 
autoplot(ts.Bcase) + ylab("Vs") + xlab("Weeks per Year")

fc = ses(ts.Bcase, h=4)
autoplot(fc)

Hình 1 mô tả số liệu thật

Hình 2 thể hiện khoảng biến thiên của các giá trị tiên lượng

Rút ra hành động trước mắt:

+ Những tháng đầu năm tập trung phát triển nội tại

+ Những tháng giữa năm tập trung phân phối