R/AIC_HMM.R

AIC_HMM <- 
function(logL, m, k) 
{ 
  foo <- try(-2 * logL + 2 * (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.