2020-07-14 15:33:00 +02:00
|
|
|
laskeLoggis <- function(counts, sumcounts, adjprior) {
|
|
|
|
|
npops <- size(counts, 3)
|
|
|
|
|
|
2020-11-19 07:53:36 +01:00
|
|
|
sum1 <- sum(sum(sum(lgamma(counts + repmat(adjprior, c(1, 1, npops))))))
|
2020-07-14 15:33:00 +02:00
|
|
|
sum3 <- sum(sum(lgamma(adjprior))) - sum(sum(lgamma(1 + sumcounts)))
|
|
|
|
|
logml2 <- sum1 - npops * sum3
|
|
|
|
|
loggis <- logml2
|
|
|
|
|
return(loggis)
|
|
|
|
|
}
|
