inst/tinytest/test-get_metadata.R

In chem16S: Chemical Metrics for Microbial Communities

# 20220505

info <- "Sample metrics are in same order as metadata"
# Get metrics for the Bison Pool dataset
RDPfile <- system.file("extdata/RDP/SMS+12.tab.xz", package = "chem16S")
RDP <- read_RDP(RDPfile, quiet = TRUE)
map <- map_taxa(RDP, refdb = "RefSeq_206", quiet = TRUE)
metrics <- get_metrics(RDP, map, refdb = "RefSeq_206")
# Get metadata for the Bison Pool dataset
mdatfile <- system.file("extdata/metadata/SMS+12.csv", package = "chem16S")
mdat <- get_metadata(mdatfile, metrics)
expect_true(all(mdat$metadata$Run == mdat$metrics$Run), info = info)
# This isn't a trivial test, because MG-RAST IDs are not in the sample order
expect_false(all(mdat$metadata$Run == metrics$Run), info = info)