simulate.network  R Documentation 
The simulate.network
and simulate.networkDynamic
wrappers
are provided for backwards compatibility. It is recommended that new
code make use of the simulate_formula.network
and
simulate_formula.networkDynamic
functions instead. See
simulate.tergm
for details on these new functions.
## S3 method for class 'network'
simulate(
object,
nsim = 1,
seed = NULL,
formation,
dissolution,
coef.form,
coef.diss,
constraints = ~.,
monitor = NULL,
time.slices = 1,
time.start = NULL,
time.burnin = 0,
time.interval = 1,
time.offset = 1,
control = control.simulate.network(),
output = c("networkDynamic", "stats", "changes", "final", "ergm_state"),
stats.form = FALSE,
stats.diss = FALSE,
verbose = FALSE,
...
)
## S3 method for class 'networkDynamic'
simulate(
object,
nsim = 1,
seed = NULL,
formation,
dissolution,
coef.form = attr(object, "coef.form"),
coef.diss = attr(object, "coef.diss"),
constraints = ~.,
monitor = NULL,
time.slices = 1,
time.start = NULL,
time.burnin = 0,
time.interval = 1,
time.offset = 1,
control = control.simulate.network(),
output = c("networkDynamic", "stats", "changes", "final", "ergm_state"),
stats.form = FALSE,
stats.diss = FALSE,
verbose = FALSE,
...
)
object 
an object of type 
nsim 
Number of replications (separate chains of networks) of the
process to run and return. The 
seed 
Seed value (integer) for the random number generator. See

formation, dissolution 
Onesided 
coef.form 
Parameters for the formation model. 
coef.diss 
Parameters for the dissolution (persistence) model. 
constraints 
A formula specifying one or more constraints
on the support of the distribution of the networks being modeled. Multiple constraints
may be given, separated by “+” and “” operators. See
The default is to have no constraints except those provided through
the Together with the model terms in the formula and the reference measure, the constraints define the distribution of networks being modeled. It is also possible to specify a proposal function directly either
by passing a string with the function's name (in which case,
arguments to the proposal should be specified through the
Note that not all possible combinations of constraints and reference measures are supported. However, for relatively simple constraints (i.e., those that simply permit or forbid specific dyads or sets of dyads from changing), arbitrary combinations should be possible. 
monitor 
A onesided formula specifying one or more terms whose
value is to be monitored. If 
time.slices 
Number of time slices (or statistics) to return from each
replication of the dynamic process. See below for return types. Defaults to
1, which, if 
time.start 
An optional argument specifying the time point at which the simulation is to start. See Details for further information. 
time.burnin 
Number of time steps to discard before starting to collect network statistics. 
time.interval 
Number of time steps between successive recordings of network statistics. 
time.offset 
Argument specifying the offset between the point when the
state of the network is sampled ( 
control 
A list of control parameters for algorithm tuning,
constructed using

output 
A character vector specifying output type: one of

stats.form, stats.diss 
Logical: Whether to return
formation/dissolution model statistics. This is not the recommended method:
use the 
verbose 
A logical or an integer to control the amount of
progress and diagnostic information to be printed. 
... 
Further arguments passed to or used by methods. 
Note that return values may be structured differently than in past versions.
Remember that in stergm
, the dissolution formula is parameterized in
terms of tie persistence: negative coefficients imply lower rates of persistence
and postive coefficients imply higher rates. The dissolution effects are simply the
negation of these coefficients.
Because the old dissolution
formula in stergm
represents
tie persistence, it maps to the new Persist()
operator
in the tergm
function, NOT the Diss()
operator
Depends on the output
argument. See simulate.tergm
for details. Note that some formation/dissolution separated
information is also attached to the return value for calls made through
simulate.network
and simulate.networkDynamic
in
an attempt to increase backwards compatibility.
logit<function(p)log(p/(1p))
coef.form.f<function(coef.diss,density) log(((1+exp(coef.diss))/(density/(1density)))1)
# Construct a network with 20 nodes and 20 edges
n<20
target.stats<edges<20
g0<network.initialize(n,dir=TRUE)
g1<san(g0~edges,target.stats=target.stats,verbose=TRUE)
S<10
# To get an average duration of 10...
duration<10
coef.diss<logit(11/duration)
# To get an average of 20 edges...
dyads<network.dyadcount(g1)
density<edges/dyads
coef.form<coef.form.f(coef.diss,density)
# ... coefficients.
print(coef.form)
print(coef.diss)
# Simulate a networkDynamic
dynsim<simulate(g1,formation=~edges,dissolution=~edges,
coef.form=coef.form,coef.diss=coef.diss,
time.slices=S,verbose=TRUE)
# "Resume" the simulation.
dynsim2<simulate(dynsim,formation=~edges,dissolution=~edges,time.slices=S,verbose=TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.