calcARDL <- function(ModFitado, horizonte = 1){
  ytplus1 <- predict(ModFitado, horizonte)
  return(ytplus1)
}


calcSerie <- function(serieOriginal, MaxDefasagens, horizonte = 1){
  resultSet <- serieOriginal[1:(length(serieOriginal)-MaxDefasagens)]
  for(i in 0:MaxDefasagens){
    serieX <- serieOriginal[1:(length(serieOriginal)-MaxDefasagens+i)]
    predicX <- as.numeric(calcARDL(serieX, horizonte)$pred)
    resultSet <- append(resultSet, predicX)
  }
  return(ts(resultSet,start = start(serieOriginal), end = end(serieOriginal), frequency = 4))
}


calcRes <- function(serieOri, serieCortada, defasagens){
  errosTot <- vector()
  for(i in 1:defasagens){
    erroX <- serieOri[length(serieOri)-defasagens+i] - serieCortada[length(serieOri)-defasagens+i]
    erroXX <- erroX^2
    errosTot <- append(errosTot, erroXX)
  }
  sumErros <- sum(errosTot)
  return(sumErros)
}