1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
context("Dot-product random graphs")
test_that("Dot product rng works", {
library(igraph)
set.seed(42)
vecs <- cbind(c(0,1,1,1,0)/3, c(0,1,1,0,1)/3, c(1,1,1,1,0)/4,
c(0,1,1,1,0))
g <- sample_dot_product(vecs)
g0 <- graph.formula(1:2:3-4)
expect_that(g[], is_equivalent_to(g0[]))
g2 <- sample_dot_product(vecs, directed=TRUE)
g20 <- graph.formula(1:2:3:4, 1-+3, 1-+4, 3-+4, 4+-1, 4+-3)
expect_that(g[], is_equivalent_to(g20[]))
vecs <- replicate(5, rep(1/2, 4))
g <- sample_dot_product(vecs)
expect_that(g[], is_equivalent_to(graph.full(5)[]))
g2 <- sample_dot_product(vecs, directed=TRUE)
expect_that(g2[], is_equivalent_to(graph.full(5, directed=TRUE)[]))
vecs <- replicate(100, rep(sqrt(1/8), 4))
g <- sample_dot_product(vecs)
expect_that(ecount(g), equals(2454))
g2 <- sample_dot_product(vecs, directed=TRUE)
expect_that(ecount(g2), equals(4938))
})
test_that("Dot product rng gives warnings", {
library(igraph)
vecs <- cbind(c(1,1,1)/3, -c(1,1,1)/3)
expect_that(g <- sample_dot_product(vecs),
gives_warning("Negative connection probability in dot-product graph"))
vecs <- cbind(c(1,1,1), c(1,1,1))
expect_that(g <- sample_dot_product(vecs),
gives_warning(paste(sep="", "Greater than 1 connection probability ",
"in dot-product graph")))
})
|
