tests/testthat/test_total_biomass.R

In harryganz/rvc: A statistical package for the Reef Visual Census

context("STRAT total_biomass function")

load("./../test_data.Rdata")

strat = subset(fk2012, STRAT == "FDLR" & PROT == 0 & SPECIES_CD == "MYC BONA")
sbiom = strat_total_biomass(psu_biomass(ssu_biomass(strat, list(WLEN_A = 2e-4, WLEN_B =3))), ntot2012)

test_that("returns correct total biomass",
          expect_equal(signif(sbiom$total_biomass, 4), 2.274e6)
          )
test_that("returns correct variance",
          expect_equal(signif(sbiom$var, 4), 4.380e11)
          )

context("DOMAIN total_biomass function")

domain = subset(fk2012, SPECIES_CD == "CEP FULV")
dbiom = domain_total_biomass(strat_biomass(psu_biomass(ssu_biomass(domain, list(WLEN_A = 2e-4, WLEN_B =2.8))),
                                           ntot2012), ntot2012)

test_that("returns correct total biomass",
          expect_equal(signif(dbiom$total_biomass, 4), 1.124e5)
          )
test_that("returns correct variance",
          expect_equal(signif(dbiom$var, 4), 1.094e9)
          )