d=sqrt((log(4/(4*5*6)^(1/3))-log(1/(2*3)^(1/3)))^2+(log(5/(4*5*6)^(1/3))-log(2/(2*3)^(1/3)))^2+(log(6/(4*5*6)^(1/3))-log(3/(2*3)^(1/3)))^2) test_that("AitchisonDistance", { expect_equal(AitchisonDistance(c(1,2,3), c(4,5,6)), d,tolerance=1e-6) })