ourMELONS/R/laskeLoggis.R

laskeLoggis <- function(counts, sumcounts, adjprior) {
  npops <- size(counts, 3)
  replicated_adjprior <- array(adjprior, c(nrow(adjprior), ncol(adjprior), npops))
  sum1 <- sum(sum(sum(lgamma(counts + replicated_adjprior))))
  sum2 <- npops * sum(sum(lgamma(adjprior))) + sum(sum(lgamma(1 + sumcounts)))
  sum1 - sum2
}
translated laskeLoggis() 2020-07-14 15:33:00 +02:00			`laskeLoggis <- function(counts, sumcounts, adjprior) {`
Restyled files Ran through styler::style_dir() in the R and tests directories in preparation for #23. 2021-11-10 14:02:35 +01:00			`npops <- size(counts, 3)`
Fixed basic parsing of FASTA files (#25) 2023-09-11 12:15:30 +02:00			`replicated_adjprior <- array(adjprior, c(nrow(adjprior), ncol(adjprior), npops))`
			`sum1 <- sum(sum(sum(lgamma(counts + replicated_adjprior))))`
Fixed calculation of logml 2024-09-27 12:59:55 +02:00			`sum2 <- npops * sum(sum(lgamma(adjprior))) + sum(sum(lgamma(1 + sumcounts)))`
			`sum1 - sum2`
Restyled files Ran through styler::style_dir() in the R and tests directories in preparation for #23. 2021-11-10 14:02:35 +01:00			`}`