Expectation and Maximization steps.
estep(theta, y) mstep(y, W, lambda, family = "Gaussian")
theta
y
W
lambda
family
J. Einbeck, D. Bonetti