context("Network models with formation offsets")
test_that("netsim works with standard offset models", {
nw <- network_initialize(n = 50)
nw <- set_vertex_attribute(nw, "race", rbinom(50, 1, 0.5))
est <- netest(nw, formation = ~edges + offset(nodematch("race")),
target.stats = 25, coef.form = -Inf,
coef.diss = dissolution_coefs(~offset(edges), 10, 0),
verbose = FALSE)
param <- param.net(inf.prob = 0.3, act.rate = 0.5)
init <- init.net(i.num = 10)
control <- control.net(type = "SI", nsims = 2, nsteps = 5, verbose = FALSE)
mod <- netsim(est, param, init, control)
expect_is(mod, "netsim")
plot(mod)
plot(mod, type = "formation")
plot(mod, type = "network")
test_net(mod)
})
test_that("netsim works with faux offset models", {
nw <- network_initialize(n = 50)
nw <- set_vertex_attribute(nw, "race", rbinom(50, 1, 0.5))
est <- netest(nw, formation = ~edges + nodematch("race"),
target.stats = c(25, 0),
coef.diss = dissolution_coefs(~offset(edges), 10, 0),
verbose = FALSE)
param <- param.net(inf.prob = 0.3, act.rate = 0.5)
init <- init.net(i.num = 10)
control <- control.net(type = "SI", nsims = 2, nsteps = 5, verbose = FALSE)
mod <- netsim(est, param, init, control)
expect_is(mod, "netsim")
plot(mod)
plot(mod, type = "formation")
plot(mod, type = "network")
test_net(mod)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.