Nothing
context("Test MBD")
if(require("fda", quietly = TRUE)) {
test_that("Equal to fda", {
data("katowice.airpollution")
expect_equal(fda:::fMBD(t(katowice.airpollution)),
fncDepthMBD(katowice.airpollution))
})
}
context("Test MBD")
test_that("Equal interfaces", {
data("katowice.airpollution")
expect_equal(
fncDepthMBD(katowice.airpollution),
fncDepthMBD(katowice.airpollution, katowice.airpollution)
)
expect_equal(
fncDepthMBD(katowice.airpollution)[10:20],
fncDepthMBD(katowice.airpollution[10:20, ], katowice.airpollution)
)
})
test_that("Ref rank", {
x <- c(1, 1, 2, 3, 4, 4, 4, 4, 4, 5)
expect_equal(rank(x), refRank(x, x))
x <- c(5, 5, 3, 2, 11, 1, 11, 11)
expect_equal(rank(x), refRank(x, x))
x <- c(5, 5, 3, 2, 11, 1, 11, 11)
y <- c(5, 5, 3, 2, 130, 1, 12, 12)
expect_equal(refRank(y, x), rank(x))
x <- c(1, 5, 5, 10)
expect_equal(refRank(c(1, 4, 5, 10), x), c(1.0, 1.0, 2.5, 4.0))
x <- c(4, 5, 6)
expect_equal(refRank(c(4.5, 5), x), c(1, 2))
x <- c(1, 4, 4, 4, 4, 5, 6)
expect_equal(refRank(x[4:5], x), rank(x)[4:5])
x <- sort(katowice.airpollution[, 1])
rank.x <- rank(x)
expect_equal(refRank(x, x), rank.x)
expect_equal(refRank(x[10:30], x), rank.x[10:30])
expect_equal(refRank(x[70:90], x), rank.x[70:90])
})
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.