... | ... |
@@ -313,10 +313,7 @@ |
313 | 313 |
rm(p) |
314 | 314 |
|
315 | 315 |
# Pre-compute lgamma values |
316 |
- lgbeta <- lgamma((seq(0, max(.colSums( |
|
317 |
- counts, nrow(counts), |
|
318 |
- ncol(counts) |
|
319 |
- )))) + beta) |
|
316 |
+ lgbeta <- lgamma((seq(0, max(colSums(counts)))) + beta) |
|
320 | 317 |
lggamma <- lgamma(seq(0, nrow(counts) + L) + gamma) |
321 | 318 |
lgdelta <- c(NA, lgamma(seq(nrow(counts) + L) * delta)) |
322 | 319 |
|