context("dominance")
library(igraph)
library(magrittr)
library(Matrix)
test_that("neighborhood-inclusion is correct", {
data("dbces11")
A <- matrix(0, 11, 11)
rownames(A) <- colnames(A) <- LETTERS[1:11]
A[c(23, 35, 45, 47, 56, 67, 79, 111, 113)] <- 1
expect_equal(neighborhood_inclusion(dbces11, sparse = FALSE), A)
####
tg <- threshold_graph(10, 0)
V(tg)$name <- LETTERS[1:10]
A <- matrix(1, 10, 10)
rownames(A) <- colnames(A) <- LETTERS[1:10]
diag(A) <- 0
A[10, ] <- 0
expect_equal(neighborhood_inclusion(tg, sparse = FALSE), A)
})
test_that("positional dominance is correct", {
g <- graph.empty(n = 11, directed = FALSE)
g <- add_edges(g, c(
1, 11, 2, 4, 3, 5, 3, 11, 4, 8, 5, 9, 5, 11, 6, 7, 6, 8,
6, 10, 6, 11, 7, 9, 7, 10, 7, 11, 8, 9, 8, 10, 9, 10
))
D <- g %>%
indirect_relations() %>%
positional_dominance(map = FALSE, benefit = FALSE)
Dmap <- g %>%
indirect_relations() %>%
positional_dominance(map = TRUE, benefit = FALSE)
A <- matrix(0, 11, 11)
A[c(23, 35, 45, 47, 56, 67, 79, 111, 113)] <- 1
Amap <- matrix(0, 11, 11)
Amap[c(
2, 23, 24, 34, 35, 45, 46, 47, 48, 56, 57, 58, 59, 60,
62, 63, 64, 65, 67, 68, 69, 70, 71, 78, 79, 80, 81, 82,
84, 86, 87, 89, 90, 91, 92, 93, 95, 96, 98, 100, 101, 102,
103, 111, 112, 113, 114, 115, 117
)] <- 1
expect_equal(D, A)
expect_equal(Dmap, Amap)
Data_mat <- matrix(c(10:1, 1:10), 10, 2)
D <- positional_dominance(Data_mat, type = "two-mode")
expect_equal(sum(abs(D)), 0)
expect_error(positional_dominance("a"))
})
test_that("dominance_graph is correct", {
library(igraph)
library(magrittr)
g <- graph.empty(n = 11, directed = FALSE)
g <- add_edges(g, c(
1, 11, 2, 4, 3, 5, 3, 11, 4, 8, 5, 9, 5, 11, 6, 7, 6, 8,
6, 10, 6, 11, 7, 9, 7, 10, 7, 11, 8, 9, 8, 10, 9, 10
))
P <- neighborhood_inclusion(g, sparse = FALSE)
g_dom <- dominance_graph(P)
expect_equal(sum(P), ecount(g_dom))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.