R/likelihood.R

Defines functions likelihood

Documented in likelihood

likelihood <- function(Sigma, Omega) {
  ot <- as.numeric(unlist(determinant(Omega)))
  if (ot[2]<=0) warning("Precision matrix estimate is not positive definite!")
  tmp <- (sum(diag(Sigma%*%Omega))  - ot[1] - dim(Omega)[1])
  if(is.finite(tmp)) {
    return(tmp)
  } else {
    return(Inf)
  }
}

Try the clime package in your browser

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

clime documentation built on June 22, 2022, 5:07 p.m.