R/Mstep.binom.R

In HiddenMarkov: Hidden Markov Models

Defines functions Mstep.binom

Documented in Mstep.binom

Mstep.binom <- function(x, cond, pm, pn){
    if (is.null(pn)) stop("Variable \"size\" must be specified in pn")
    prob <- as.numeric(matrix(x, nrow=1) %*% cond$u)/
                as.numeric(matrix(pn$size, nrow=1) %*% cond$u)
    return(list(prob=prob))
}

Try the HiddenMarkov package in your browser