Nothing
context("argillic bounds and clay increase")
data(sp1, package = 'aqp')
depths(sp1) <- id ~ top + bottom
site(sp1) <- ~ group
p <- sp1[1]
hzdesgnname(p) <- "name"
hztexclname(p) <- "texture"
clay.attr <- 'prop' # clay contents %
texcl.attr <- 'texture' # class containing textural class (for finding sandy textures)
#implicitly testing the basic logic of this function :)
threshold.fun <- crit.clay.argillic
#standard for argillic
vertical.distance <- 30
test_that("get.increase.matrix() (used for getArgillicBounds())", {
m <- get.increase.matrix(p, clay.attr, threshold.fun, vertical.distance)
# correct data type?
expect_true(inherits(m, 'matrix'))
# is square matrix?
expect_true(nrow(m) == ncol(m))
# no TRUE in lower triangle?
expect_true(all(!(m * lower.tri(m))))
# correct number of eluvial horizons identified in first illuvial column?
# this is also a test of the underlying crit.clay.argillic
expect_equal(colSums(m)[4], 2)
})
test_that("argillic.clay.increase() (used for getArgillicBounds())", {
#make sure the main wrapper method for the increase functions works
d <- argillic.clay.increase.depth(p, clay.attr = 'prop')
# this is also a test of the underlying crit.clay.argillic
expect_equal(d, 49)
})
test_that("getArgillicBounds()", {
# name and texture class are determined from metadata
d <- getArgillicBounds(p, clay.attr = 'prop')
# this makes sure estimateSoilDepth() isn't broken...
expect_equivalent(d, c(49, 89))
# error when hzdesgn and texcl.attr are unknown; no guessing
expect_error(getArgillicBounds(p, hzdesgn = 'foo', clay.attr = 'prop', texcl.attr = 'bar'))
})
test_that("getArgillicBounds - error conditions", {
expect_error(getArgillicBounds(p, hzdesgn = "foo"))
expect_error(getArgillicBounds(p, texcl.attr = "foo"))
expect_error(getArgillicBounds(p, clay.attr = "foo"))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.