R/likelihood.R

Defines functions 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 scio package in your browser

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

scio documentation built on March 3, 2020, 1:07 a.m.