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
comp = decompose(ts.Bcase)
library(ggplot2)
library(ggfortify)
autoplot(comp)
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)
library(forecast); library(fpp2)
autoplot(ts.Bcase) + ylab("Vs") + xlab("Weeks per Year")
fc = ses(ts.Bcase, h=4)
autoplot(fc)