5-Automatic_MAPPA_model.Rmd

---
title: "Automatic_MAPPA"
author: "Elias Mayer"
date: "21 8 2021"
output:
  word_document: default
  html_document: default
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

1. Simple batch MAPPA model evaluated against different other forecasting models (CV - snaive) for more details see report: https://github.com/EliDerDeli/Forecasting-M3-experiments/blob/main/Forecast_report.pdf 

```{r MAPA}
library(Mcomp)
library(forecast)
library(dplyr)
library(tsibble)
library(ggplot2)
library(ggthemes)
library(stringr)
library(MAPA)

#------------------------------------------------------------


arr_ids <- array(NA,c(70,1))

count = 1 


for (ex in 1:length(M3)){     #figure our correct indexes     DATA Batch
    

   if ((M3[[ex]]$period == "QUARTERLY") & (ex >= 701) & (ex <= 1400) & (ex %% 10 == 2)) {
     
    arr_ids[count,1] <- ex
       
    count = count + 1
    
   }
  
}


arr_all <- array(NA,c(70,13))

count =1 

colnames(arr_all) <- c("Series_ID", "n_observations", "type", "period", "MAPE_MAPA", "MASE_MAPA", "not in use", "not in us", "MAPE_NAIVE", "MAPE_MEAN", "MAPE_SNAIVE", "Trend","Seasonality")

for (nameId in arr_ids) {

  y_training <- M3[[nameId]]$x  
  y_test <- M3[[nameId]]$xx
  
  #Standard data
  
  y_id <- M3[[nameId]]$sn
  y_type <- M3[[nameId]]$type
  y_period<- M3[[nameId]]$period
  y_num <- M3[[nameId]]$n
  
  arr_all[count,1] <- y_id
  arr_all[count,2] <- y_num
  arr_all[count,3] <- y_type
  arr_all[count,4] <- y_period
  
  
#Forecast sepcific 
  
map <- mapafitX <- mapaest(y_training,paral=2)

map_f <- mapafor(y_training,mapafitX,fh=8, outplot=0)   

outFc <- map_f$outfor

acMAPE <- forecast::accuracy(outFc, y_test)
  

  #Fill errors
  
  arr_all[count,5] <- acMAPE[1,5] #MAPE
  
  arr_all[count,6] <- "not applied"
  
  arr_all[count,7] <-  "not applied " 
  
  arr_all[count,8] <- "not applied "
  
  
  
  fx1 <- forecast::naive(y_training)
  
  arr_all[count,9] <- forecast::accuracy(fx1, y_test)[2,5]
  
  fx2 <- forecast::meanf(y_training)
  
  arr_all[count,10] <- forecast::accuracy(fx2, y_test)[2,5]
  
  #SNAIVE 
  
  try( # if seasonality 
  fx3 <- forecast::snaive(y_training))
  try(arr_all[count,11] <- forecast::accuracy(fx3, y_test)[2,5])
 
  #Trend
  
  ac <- forecast::mstl(y_training) 


  value = ac[[1,2]] - ac[[length(y_training),2]] 
  
    
  arr_all[count,12] <- value
  
  
  arr_all[count,13] <- forecast::findfrequency(y_training)
  
  
  count = count + 1
}
```



```{r examination 1 - a, echo=FALSE}


arr_tsb <- arr_all %>% as_tibble() 


arr_tsb

```