Forward-backward algorithm called in fit.MSAR.
Estep.MSAR(data, theta, smth = FALSE,
verbose = FALSE,
covar.emis = covar.emis, covar.trans = covar.trans)array of univariate or multivariate series with dimension T*N.samples*d. T: number of time steps of each sample, N.samples: number of realisations of the same stationary process, d: dimension.
model's parameter; object of class MSAR. See also init.theta.MSAR. .
If smth=FALSE, only the forward step is computed for forecasting probabilities. If smth=TRUE, the smoothing probabilities are computed too.
if verbose=TRUE some results are printed at each iteration.
covariables for emission probabilities.
covariables for transition probabilities.
A list including
log likelihood
smoothing probabilities: \(P(S_t=s|y_0,\cdots,y_T)\)
one step smoothing probabilities: \(P(S_t=s,S_{t+1}|y_0,\cdots,y_T)\)
Ailliot P., Monbet V., (2012), Markov switching autoregressive models for wind time series. Environmental Modelling & Software, 30, pp 92-101.
fit.MSAR, Mstep.hh.MSAR