context("dyad_census")
test_that("dyad_census works", {
library(igraph)
ce <- simplify(read_graph(gzfile("celegansneural.gml.gz"), format="gml"))
dc <- dyad_census(ce)
expect_that(dc, equals(list(mut=197, asym=1951, null=41808)))
expect_that(sum(which_mutual(ce)), equals(dc$mut * 2))
expect_that(ecount(as.undirected(ce, mode="collapse")) - dc$mut,
equals(dc$asym))
expect_that(sum(unlist(dc)), equals(vcount(ce) * (vcount(ce)-1) / 2))
})
Generated by dwww version 1.15 on Sun Jun 23 21:31:10 CEST 2024.