27 lines
670 B
R
27 lines
670 B
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)
|
|
}
|