max/ourMELONS
1
0
Fork 0
Code Issues 1 Pull requests Projects Releases Packages Wiki Activity Actions

translated laskeLoggis()

Browse source
This commit is contained in:
Waldir Leoncio 2020-07-14 15:33:00 +02:00
parent b215a663ac
commit dfc9659180
2 changed files with 10 additions and 11 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
R/greedyMix.R
View file

@ -978,14 +978,4 @@ greedyMix <- function(
# else
# popDiff = diff([0 pops npops+1]);
# emptyPop = min(find(popDiff > 1));
# end
# %------------------------------------------------------
# function loggis = laskeLoggis(counts, sumcounts, adjprior)
# npops = size(counts,3);
# logml2 = sum(sum(sum(gammaln(counts+repmat(adjprior,[1 1 npops]))))) ...
# - npops*sum(sum(gammaln(adjprior))) - ...
# sum(sum(gammaln(1+sumcounts)));
# loggis = logml2;
# end

9
R/laskeLoggis.R Normal file
View file

@ -0,0 +1,9 @@
laskeLoggis <- function(counts, sumcounts, adjprior) {
npops <- size(counts, 3)
sum1 <- sum(sum(sum(gammaln(counts + repmat(adjprior, c(1, 1, npops))))))
sum3 <- sum(sum(lgamma(adjprior))) - sum(sum(lgamma(1 + sumcounts)))
logml2 <- sum1 - npops * sum3
loggis <- logml2
return(loggis)
}