context("aggregateSoilDepth")
data(sp1)
depths(sp1) <- id ~ top + bottom
site(sp1) <- ~ group
# set horizon designation in SPC
hzdesgnname(sp1) <- 'name'
## tests
test_that("works as expected", {
a <- aggregateSoilDepth(sp1, groups = 'group', crit.prob = 0.9)
# structure
expect_true(inherits(a, 'data.frame'))
expect_true(ncol(a) == 3)
expect_true(nrow(a) == length(unique(sp1$group)))
# values
expect_true(all(a$group == c(1, 2)))
## TODO: cross-reference with getSoilDepthClass
## these don't currently give the same results
# site(sp1) <- getSoilDepthClass(sp1)
# tapply(sp1$depth, sp1$group, quantile, probs = c(0.9))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.