tests/testthat/test-linRcpp.R

In propr: Calculating Proportionality Between Vectors of Compositional Data

library(propr)
library(cccrm)

data(mail)
rho <- perb(mail, 3)
swept <- sweep(rho@logratio, 2, apply(rho@logratio, 2, mean), "-")
N <- nrow(swept)

dat <- data.frame("values" = c(swept[, 1], swept[, 2]),
                  "method" = c(rep(1, N), rep(2, N)))
ccc <- cccUst(dat, "values", "method")
names(ccc) <- NULL

lin <- propr:::linRcpp(rho@matrix, as.matrix(rho@logratio))

test_that("linRcpp performs Z-transformation correctly (alr)", {

  expect_equal(
    ccc[5],
    lin[2, 1]
  )
})

test_that("lltRcpp returns lower left triangle", {

  expect_equal(
    propr:::lltRcpp(lin),
    propr:::proprTri(lin)
  )
})

test_that("urtRcpp returns upper right triangle", {

  expect_equal(
    propr:::urtRcpp(lin),
    propr:::proprTri(t(lin))
  )
})