inst/1-B12/ire/Rsource/EMiterate_1-B12.R

In jlivsey/EMsignal:

#' single iteration of EM algorithm
#'
#' @param Sig list with each element of list being an NxN sigma matrix
#' @param lMS list with lMS[[1]] = M an [N,N,T,T] array and lMS[[2]] = S[T,N]
#' matrix
#'
#' @return updated values of all inputs
#' @export
#'

EMiterate_1_B12 = function(Sig, lMS, data, mdl){
  J = length(Sig)
  for(j in 1:J){
    Sig[[j]] = EMcritical_1_B12(j, Sig, lMS)
  }
  param = sig2param(Sig, mdl, data)
  lMS = sigexRun_1_B12(param, data, mdl)
  return(list(Sig, lMS))
}