R/BIC_HMM.R

BIC_HMM <-
function(size, m, k, logL)  
{ 
  foo <- try(-2 * logL + log(size) * (m^2 + k * m - 1), silent = FALSE)
  if (inherits(foo, "try-error")) 
  {
    foo <- NA
  }
  return(foo)
}

Try the HMMpa package in your browser

Any scripts or data that you put into this service are public.

HMMpa documentation built on May 2, 2019, 7:58 a.m.