ourMELONS/R/checkLogml.R

checkLogml <- function(priorTerm, adjprior, cliques, separators) {
  # tarkistaa logml:n

  # global CLIQCOUNTS
  # global SEPCOUNTS
  # global PARTITION

  npops <- length(unique(globals$PARTITION))
  cliqcounts <- computeCounts(cliques, separators, npops)$cliqcounts
  sepcounts <- computeCounts(cliques, separators, npops)$sepcounts

  globals$CLIQCOUNTS <- cliqcounts
  globals$SEPCOUNTS <- sepcounts

  logml <- computeLogml(adjprior, priorTerm)$logml
  spatialPrior <- computeLogml(adjprior, priorTerm)$spatialPrior

  disp(
    c(
      'logml: ',
      logml2String(logml),
      ', spatial prior: ',
      logml2String(spatialPrior)
    )
  )
  return(logml)
}