powered by
orthbasis(D)
isomLR
cenLR
data(expenditures) V <- orthbasis(ncol(expenditures)) xcen <- cenLR(expenditures)$x.clr xi <- as.matrix(xcen) %*% V xi2 <- isomLR(expenditures) all.equal(xi, xi2)
Run the code above in your browser using DataLab