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

Fixed log2() on getPopDistancesByKL()

Browse source
This commit is contained in:
Waldir Leoncio 2023-02-01 13:46:58 +01:00
parent 4ba4d4b934
commit 45a41a9faa
1 changed files with 2 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
R/getPopDistancesByKL.R
View file

@ -25,9 +25,9 @@ getPopDistancesByKL <- function(adjprior) {
for (pop2 in (pop1 + 1):npops) { for (pop2 in (pop1 + 1):npops) {
dist1 <- d[, , pop1] dist1 <- d[, , pop1]
dist2 <- d[, , pop2] dist2 <- d[, , pop2]
div12 <- sum(sum(dist1 * log2((dist1 + 10^-10) / (dist2 + 10^-10)))) / div12 <- sum(sum(dist1 * base::log2((dist1 + 10^-10) / (dist2 + 10^-10)))) /
nloci nloci
div21 <- sum(sum(dist2 * log2((dist2 + 10^-10) / (dist1 + 10^-10)))) / div21 <- sum(sum(dist2 * base::log2((dist2 + 10^-10) / (dist1 + 10^-10)))) /
nloci nloci
div <- (div12 + div21) / 2 div <- (div12 + div21) / 2
distances[pointer] <- div distances[pointer] <- div