library("testthat") library("mapi") test_that("MAPI_EstimateHalfwidth", { data("samples") hw50 <- MAPI_EstimateHalfwidth(samples, crs=3857, beta=0.5) expect_equal(hw50, 286.835888252914) hw25 <- MAPI_EstimateHalfwidth(samples, crs=3857, beta=0.25) expect_equal(hw25, 143.417944126457) })