The E-step of the EM algorithm.
This is for ML (not REML), but it seems to be good enough to use for REML starting values.
Parameters: | fe_params : 1d ndarray
cov_re : 2d ndarray
scale : positive scalar
|
---|---|
Returns: | m1x : 1d ndarray
m1y : scalar
m2 : 2d ndarray
m2xx : 2d ndarray
|