testing/old/sbp2.R

In mcomas/coda.base: A Basic Set of Functions for Compositional Data Analysis

D = 9
n = 100
Y = matrix(exp(rnorm(D*n)), ncol = D)
x = rnorm(n)

sbp_basis(find_predictive_balance(Y,x,method=1), silent = TRUE)
sbp_basis(find_predictive_balance(Y,x,method=2), silent = TRUE)

microbenchmark(
  pb_basis(Y, method = 'exact'),
  find_principal_balance(Y))

find_principal_balance(Y)
pb_basis(Y, method = 'exact')

microbenchmark(
  sbp_basis(find_principal_balance(Y), silent = TRUE),
  sbp_basis(find_predictive_balance(Y,x,method=1), silent = TRUE),
  sbp_basis(find_predictive_balance(Y,x,method=2), silent = TRUE))

pb_basis(X, method = 'exact')[,1, drop=FALSE]

x = matrix(rnorm(n), ncol=1)
sbp_basis(sbp2_test_2(X, x), silent = TRUE)

var(coordinates(X, pb_basis(X, method = 'exact')[,1, drop=FALSE]))

library(microbenchmark)

microbenchmark(
  pb_basis(X, method = 'exact'),
  sbp2_test_1(X))