Gaussian EM-step with random initialisation.
em_rinit(y, order, partempl)etk2tau(etk)
for em_rinit
, an object from class "MixARGaussian"
for etk2tau
, a matrix representing tau
(i-th row
contains probabilities corresponding to the i-th observation)
time series.
MixAR order, vector of length the number of components.
parameter template, a list containing one element for each mixture
component, see randomArCoefficients
.
MixAR component residuals, a matrix.
Georgi N. Boshnakov
em_rinit
generates random MAR residuals, performs a non-distributional
E-step, and a Gaussian M-step.
etk2tau
estimates tau
from component residuals
only. Note that this is unlike em_tau
, which also needs
the noise pdf's, as well as estimates of the mixture probabilities.
em_rinit
uses etk2tau
to start the EM algorithm.