From 99ec89fd85c7e1bafb859c2fcdd2e45dca1310f7 Mon Sep 17 00:00:00 2001 From: Waldir Leoncio Date: Fri, 27 Sep 2024 13:18:10 +0200 Subject: [PATCH] Fixed calculation of distances --- R/writeMixtureInfo.R | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/R/writeMixtureInfo.R b/R/writeMixtureInfo.R index 0b80583..6514454 100644 --- a/R/writeMixtureInfo.R +++ b/R/writeMixtureInfo.R @@ -196,7 +196,10 @@ writeMixtureInfo <- function( prior[1, nollia] <- 1 for (pop1 in 1:npops) { squeezed_COUNTS_prior <- squeeze(globals$COUNTS[, , pop1]) + prior - d[, , pop1] <- squeezed_COUNTS_prior / sum(squeezed_COUNTS_prior) + repmat_squeezed_COUNTS_prior <- matlab2r::repmat( + colSums(squeezed_COUNTS_prior), c(maxnoalle, 1L) + ) + d[, , pop1] <- squeezed_COUNTS_prior / repmat_squeezed_COUNTS_prior } ekarivi <- as.character(npops) if (verbose) cat(ekarivi)