diff --git a/R/adam.R b/R/adam.R index 78a65b77..4c673e25 100644 --- a/R/adam.R +++ b/R/adam.R @@ -9510,12 +9510,10 @@ reapply.adam <- function(object, nsim=1000, bootstrap=FALSE, heuristics=NULL, .. switch(Etype, "A"= arimaPolynomials$ariPolynomial[nonZeroARI[,1]] %*% t(profilesRecentArray[j+componentsNumberARIMA, - 1:initialArimaNumber, i]) / - tail(arimaPolynomials$ariPolynomial,1), + 1:initialArimaNumber, i]), "M"=exp(arimaPolynomials$ariPolynomial[nonZeroARI[,1]] %*% t(log(profilesRecentArray[j+componentsNumberARIMA, - 1:initialArimaNumber, i])) / - tail(arimaPolynomials$ariPolynomial,1))); + 1:initialArimaNumber, i])))); } } else{ @@ -9533,12 +9531,10 @@ reapply.adam <- function(object, nsim=1000, bootstrap=FALSE, heuristics=NULL, .. switch(Etype, "A"=arimaPolynomials$maPolynomial[nonZeroMA[,1]] %*% t(profilesRecentArray[componentsNumberETS+componentsNumberARIMA, - 1:initialArimaNumber, i]) / - tail(arimaPolynomials$maPolynomial,1), + 1:initialArimaNumber, i]), "M"=exp(arimaPolynomials$maPolynomial[nonZeroMA[,1]] %*% t(log(profilesRecentArray[componentsNumberETS+componentsNumberARIMA, - 1:initialArimaNumber, i])) / - tail(arimaPolynomials$maPolynomial,1))); + 1:initialArimaNumber, i])))); } } }