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.