tests/testthat/test-ps_pca.R

In benmarwick/karon: Compositional Data Analysis of Archaeological Artefacts

# test validity of ps_pca()

data(ObsidianSources)
analyticVars<-c("Rb","Sr","Y","Zr","Nb")
pca <- ps_pca(data=ObsidianSources, ID="ID", GroupVar="Code",
       Groups="All", AnalyticVars=analyticVars, pcPlot=FALSE,PlotPoints=FALSE,PlotEllipses=FALSE)

test_that("ps_pca output is a list", {
  expect_type(pca, "list")
})

#  use saved valid output from the function
data(test_pca)

# algorithm can change sign of a principal component

test_that("ps_pca Predicted pvalues are correct: first column is character", {
  expect_equal(abs(pca$Predicted[,-1]), abs(test_pca$Predicted[,-1]))
})