tests/testthat/test-full.R

In mpadge/spatialcluster: R port of redcap

context ("full")

test_that ("structure", {
    set.seed (1)
    n <- 100
    xy <- matrix (runif (2 * n), ncol = 2)
    dmat <- matrix (runif (n^2), ncol = n)
    ncl <- 4
    scl <- scl_full (xy, dmat, ncl = ncl)
    expect_is (scl, "scl")
    expect_true (scl$pars$ncl == ncl)
    expect_true (all (names (scl) %in%
        c (
            "tree", "merges", "ord", "nodes",
            "pars", "statistics"
        )))
    cl <- scl$nodes$cluster [!is.na (scl$nodes$cluster)]
    expect_true (length (unique (cl)) == ncl)
})

test_that ("methods", {
    set.seed (1)
    n <- 100
    xy <- matrix (runif (2 * n), ncol = 2)
    dmat <- matrix (runif (n^2), ncol = n)
    ncl <- 8
    scl1 <- scl_full (xy, dmat, ncl = ncl, linkage = "single")
    scl2 <- scl_full (xy, dmat, ncl = ncl, linkage = "average")
    expect_true (!identical (scl1, scl2))
    cl1 <- scl1$nodes$cluster [!is.na (scl1$nodes$cluster)]
    expect_equal (length (unique (cl1)), ncl)
    cl2 <- scl2$nodes$cluster [!is.na (scl2$nodes$cluster)]
    expect_equal (length (unique (cl2)), ncl)
})

test_that ("recluster", {
    set.seed (1)
    n <- 100
    xy <- matrix (runif (2 * n), ncol = 2)
    dmat <- matrix (runif (n^2), ncol = n)
    scl <- scl_full (xy, dmat, ncl = 4)
    scl1 <- scl_full (scl, ncl = 3)
    scl2 <- scl_recluster (scl, ncl = 3)
    expect_identical (scl1, scl2)
    expect_error (
        scl3 <- scl_redcap (scl, ncl = 3),
        "scl_redcap can pass to scl_recluster only"
    )
})