# ves() - Vector Exponential Smoothing

#### 2018-11-30

This vignette covers ves() function, which is a part of smooth package. In this vignette we will use data from Mcomp package, so it is advised to install it.

require(smooth)
require(Mcomp)

We will use to time series from the M3 united in a vector:

Y <- cbind(M3$N2570$x,M3$N2571$x);

ves() function allows constructing Vector Exponential Smoothing in either pure additive or pure multiplicative form. The function has several elements that can either be individual or grouped. The former means that all the time series use the same value. For example, persistence="g" means that the smoothing parameters for all the series are the same. A simple call for ves() results in estimation of VES(A,N,N) with grouped smoothing parameters, transition matrix and individual initials:

ves(Y, h=18, holdout=TRUE, silent=FALSE)
## Time elapsed: 0.02 seconds
## Model estimated: VES(ANN)
## 6 parameters were estimated for 2 time series in the process
## Cost function type: likelihood; Cost function value: 13.717
##
## Information criteria:
##      AIC     AICc      BIC     BICc
## 1906.528 1919.702 1914.283 1933.197

The output tells us how much time the estimation took, what model we estimated, how many parameters were estimated, the cost function type used and its value and finally the information criteria. Currently we do not provide error measures for the holdout, this functionality will be available with newer releases of smooth.

In some cases we may decide that the series should be connected with each other. In this case we can ask function to use “dependent” persistence. This means that along with the individual smoothing parameters, we will estimate cross-series ones. Here’s the example:

ourModel <- ves(Y, "AAN", persistence="d", h=18, holdout=TRUE, silent=FALSE)

The resulting persistence matrix contains more values than the individual one:

ourModel$persistence ## M3_N2570_x M3_N2571_x ## M3_N2570_x_level 0.65068718 0.19102137 ## M3_N2570_x_trend 0.17115192 -0.09290405 ## M3_N2571_x_level -0.07796857 0.89884351 ## M3_N2571_x_trend 0.29444941 -0.12528793 Note that some of the values of smoothing parameters are negative and the others are greater than one. This is a normal behaviour for VES model in this implementation. Currently we only have bounds derived from the stability region (bounds="admissible") and we do not do traditional restrictions yet (and not sure if we ever will). Currently we have pure additive and pure multiplicative models only, and I don’t intend introducing mixed models for VES at all, because I think that they are evil. The multiplicative model implemented in VES is in fact just an additive model applied to the data in logarithms. Let’s see how the damped trend multiplicative seasonal model with individual damping and smoothing parameters looks like: ourModel <- ves(Y, "MMdM", phi="i", persistence="i", h=18, holdout=TRUE) Number of estimated parameters in the model can be extracted via nParam() method. However, when it comes to the calculation of the number of degrees of freedom in the model, this value is divided by the number of series (Lütkepohl 2005). So both ourModel$Sigma and all the information criteria rely on the $$df = T - k_m$$, where $$T$$ is the number of observations and $$k_m = \frac{k}{m}$$ is the number of parameters $$k$$ per series ($$m$$ is the number of series).

AICc and BICc for the vector models are calculated as proposed in (Bedrick and Tsai 1994) and (Tremblay and Wallach 2004).

Currently we don’t do model selection, don’t have exogenous variables and don’t produce conditional prediction intervals. But at least it works and allows you to play around with it :).

### References

Bedrick, Edward J, and Chih-Ling Tsai. 1994. “Model Selection for Multivariate Regression in Small Samples.” Biometrics 50 (1): 226. doi:10.2307/2533213.

Lütkepohl, Helmut. 2005. New Introduction to Multiple Time Series Analysis. Springer Berlin Heidelberg. doi:10.1007/978-3-540-27752-1.

Tremblay, Marie, and Daniel Wallach. 2004. “Comparison of parameter estimation methods for crop models.” Agronomie 24 (6-7): 351–65. doi:10.1051/agro:2004033.