## Test script
library(EpiModelHIVhet)
st <- make_nw_het(part.dur = 2013)
est <- netest(st$nw,
formation = st$formation,
target.stats = st$stats,
coef.form = -Inf,
coef.diss = st$coef.diss,
constraints = ~bd(maxout = 3),
set.control.ergm = control.ergm(MCMLE.maxit = 500, MPLE.type = "penalized"))
dx <- netdx(est, nsims = 5, nsteps = 250,
set.control.ergm = control.simulate.ergm(MCMC.burnin = 1e6))
print(dx)
plot(dx)
param <- param_het()
init <- init_het(i.prev.male = 0.25, i.prev.feml = 0.25)
control <- control_het(nsteps = 2600)
sim <- netsim(est, param, init, control)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.