R/Mstep.exp.R

In HiddenMarkov: Hidden Markov Models

Defines functions Mstep.exp

Documented in Mstep.exp

Mstep.exp <- function(x, cond, pm, pn){
    rate <- apply(cond$u, MARGIN=2, FUN=sum)/
          as.numeric(matrix(x, nrow=1) %*% cond$u)
    return(list(rate=rate))
}

Try the HiddenMarkov package in your browser