Running ‘combine.networks.R’ [2s/4s]
  Running ‘degree.mean.age.R’ [0s/1s]
  Running ‘dynamic_EGMME.R’ [0s/1s]
  Running ‘dynamic_EGMME.simple.R’ [0s/1s]
  Running ‘dynamic_MLE.R’ [2s/4s]
  Running ‘dynamic_MLE_2.R’ [0s/1s]
  Running ‘dynamic_MLE_blockdiag.R’ [0s/1s]
  Running ‘dynamic_MLE_blockdiag.bipartite.R’ [0s/1s]
  Running ‘sim_gf_sum.R’ [2s/4s]
  Running ‘simulate_networkDynamic.R’ [2s/3s]
  Running ‘tergm_offset_tests.R’ [0s/1s]
  Running ‘tergm_parallel.R’ [0s/1s]
Running the tests in ‘tests/combine.networks.R’ failed.
Complete output:
  > #  File tests/combine.networks.R in package tergm, part of the Statnet suite
  > #  of packages for network analysis, https://statnet.org .
  > #
  > #  This software is distributed under the GPL-3 license.  It is free,
  > #  open source, and has the attribution requirements (GPL Section 7) at
  > #  https://statnet.org/attribution
  > #
  > #  Copyright 2008-2020 Statnet Commons
  > #######################################################################
  > # This tests the internal function combine.networks.
  > 
  > n <- 10
  > T <- 10
  > 
  > library(tergm)
  Loading required package: ergm
  Loading required package: network
  network: Classes for Relational Data
  Version 1.16.1 created on 2020-10-06.
  copyright (c) 2005, Carter T. Butts, University of California-Irvine
                      Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Martina Morris, University of Washington
                      Skye Bender-deMoll, University of Washington
   For citation information, type citation("network").
   Type help("network-package") to get started.
  
  
  ergm: version 3.11.0, created on 2020-10-14
  Copyright (c) 2020, Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Carter T. Butts, University of California -- Irvine
                      Steven M. Goodreau, University of Washington
                      Pavel N. Krivitsky, UNSW Sydney
                      Martina Morris, University of Washington
                      with contributions from
                      Li Wang
                      Kirk Li, University of Washington
                      Skye Bender-deMoll, University of Washington
                      Chad Klumb
                      Michał Bojanowski, Kozminski University
                      Ben Bolker
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("ergm").
  
  NOTE: Versions before 3.6.1 had a bug in the implementation of the bd()
  constraint which distorted the sampled distribution somewhat. In
  addition, Sampson's Monks datasets had mislabeled vertices. See the
  NEWS and the documentation for more details.
  
  NOTE: Some common term arguments pertaining to vertex attribute and
  level selection have changed in 3.10.0. See terms help for more
  details. Use 'options(ergm.term=list(version="3.9.4"))' to use old
  behavior.
  
  Loading required package: networkDynamic
  
  networkDynamic: version 0.10.1, created on 2020-01-16
  Copyright (c) 2020, Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Pavel N. Krivitsky, University of Wollongong
                      Skye Bender-deMoll, University of Washington
                      with contributions from
                      Zack Almquist, University of California -- Irvine
                      David R. Hunter, Penn State University
                      Li Wang
                      Kirk Li, University of Washington
                      Steven M. Goodreau, University of Washington
                      Jeffrey Horner
                      Martina Morris, University of Washington
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("networkDynamic").
  
  
  tergm: version 3.7.0, created on 2020-10-15
  Copyright (c) 2020, Pavel N. Krivitsky, UNSW Sydney
                      Mark S. Handcock, University of California -- Los Angeles
                      with contributions from
                      David R. Hunter, Penn State University
                      Steven M. Goodreau, University of Washington
                      Martina Morris, University of Washington
                      Nicole Bohme Carnegie, New York University
                      Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Skye Bender-deMoll
                      Li Wang
                      Kirk Li, University of Washington
                      Chad Klumb
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("tergm").
  
  > yl <- replicate(T,
  +                 {
  +                   y <- network.initialize(n,dir=FALSE)
  +                   y <- simulate(y~edges, coef=-1)
  +                   y %n% "x" <- matrix(runif(n*n),n,n)
  +                   y %v% "v" <- runif(n)
  +                   y %e% "e" <- runif(network.edgecount(y))
  +                   y
  +                 },
  +                 simplify=FALSE)
  Error in set.vertex.attribute(g, "na", rep(FALSE, n), 1:n) : 
    attempt to set index 0/0 in SET_VECTOR_ELT
  Calls: replicate ... lapply -> FUN -> network.initialize -> set.vertex.attribute
  Execution halted
Running the tests in ‘tests/dynamic_MLE.R’ failed.
Complete output:
  > #  File tests/dynamic_MLE.R in package tergm, part of the Statnet suite
  > #  of packages for network analysis, https://statnet.org .
  > #
  > #  This software is distributed under the GPL-3 license.  It is free,
  > #  open source, and has the attribution requirements (GPL Section 7) at
  > #  https://statnet.org/attribution
  > #
  > #  Copyright 2008-2020 Statnet Commons
  > #######################################################################
  > library(tergm)
  Loading required package: ergm
  Loading required package: network
  network: Classes for Relational Data
  Version 1.16.1 created on 2020-10-06.
  copyright (c) 2005, Carter T. Butts, University of California-Irvine
                      Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Martina Morris, University of Washington
                      Skye Bender-deMoll, University of Washington
   For citation information, type citation("network").
   Type help("network-package") to get started.
  
  
  ergm: version 3.11.0, created on 2020-10-14
  Copyright (c) 2020, Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Carter T. Butts, University of California -- Irvine
                      Steven M. Goodreau, University of Washington
                      Pavel N. Krivitsky, UNSW Sydney
                      Martina Morris, University of Washington
                      with contributions from
                      Li Wang
                      Kirk Li, University of Washington
                      Skye Bender-deMoll, University of Washington
                      Chad Klumb
                      Michał Bojanowski, Kozminski University
                      Ben Bolker
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("ergm").
  
  NOTE: Versions before 3.6.1 had a bug in the implementation of the bd()
  constraint which distorted the sampled distribution somewhat. In
  addition, Sampson's Monks datasets had mislabeled vertices. See the
  NEWS and the documentation for more details.
  
  NOTE: Some common term arguments pertaining to vertex attribute and
  level selection have changed in 3.10.0. See terms help for more
  details. Use 'options(ergm.term=list(version="3.9.4"))' to use old
  behavior.
  
  Loading required package: networkDynamic
  
  networkDynamic: version 0.10.1, created on 2020-01-16
  Copyright (c) 2020, Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Pavel N. Krivitsky, University of Wollongong
                      Skye Bender-deMoll, University of Washington
                      with contributions from
                      Zack Almquist, University of California -- Irvine
                      David R. Hunter, Penn State University
                      Li Wang
                      Kirk Li, University of Washington
                      Steven M. Goodreau, University of Washington
                      Jeffrey Horner
                      Martina Morris, University of Washington
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("networkDynamic").
  
  
  tergm: version 3.7.0, created on 2020-10-15
  Copyright (c) 2020, Pavel N. Krivitsky, UNSW Sydney
                      Mark S. Handcock, University of California -- Los Angeles
                      with contributions from
                      David R. Hunter, Penn State University
                      Steven M. Goodreau, University of Washington
                      Martina Morris, University of Washington
                      Nicole Bohme Carnegie, New York University
                      Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Skye Bender-deMoll
                      Li Wang
                      Kirk Li, University of Washington
                      Chad Klumb
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("tergm").
  
  > options(tergm.eval.loglik=FALSE)
  > 
  > tolerance<-0.05
  > n<-10
  > m<-6
  > theta<--1.5
  > 
  > z.error <- function(truth, fit){
  +   if(truth==coef(fit)) 0 # Infinite case
  +   else abs(truth-coef(fit))/sqrt(diag(vcov(fit)))
  + }
  > 
  > logit<-function(p) log(p/(1-p))
  > 
  > form.mle<-function(y0,y1){
  +   logit(network.edgecount(y1-y0,na.omit=TRUE)/(network.dyadcount(y1)-network.edgecount(y0-is.na(y1))))
  + }
  > 
  > diss.mle<-function(y0,y1){
  +   -logit(network.edgecount(y0-y1,na.omit=TRUE)/(network.edgecount(y0-is.na(y1))))
  + }
  > 
  > do.run <- function(dir, bip=FALSE, prop.weights="default"){
  + if(bip){ # Extreme theta creates networks with too few ties to properly test.
  +   theta <- theta/2
  + }
  +   
  + y0<-network.initialize(n,dir=dir,bipartite=bip)
  + set.seed(321)
  + y0<-simulate(y0~edges, coef=theta, control=control.simulate(MCMC.burnin=n^2*2))
  + 
  + cat("Complete data:\n")
  + 
  + set.seed(123)
  + y1<-simulate(y0~edges, coef=theta, control=control.simulate(MCMC.burnin=n^2*2))
  + 
  + # Force CMPLE
  + set.seed(543)
  + fit<-stergm(list(y0,y1), formation=~edges, dissolution=~edges, estimate="CMPLE", times=c(1,2), eval.loglik=FALSE)
  + 
  + stopifnot(fit$estimate=="CMPLE", fit$formation.fit$estimate=="MPLE", fit$dissolution.fit$estimate=="MPLE")
  + stopifnot(z.error(form.mle(y0,y1), fit$formation.fit) <= tolerance)
  + stopifnot(z.error(diss.mle(y0,y1), fit$dissolution.fit) <= tolerance)
  + 
  + # Autodetected CMPLE
  + set.seed(543)
  + fit<-stergm(list(y0,y1), formation=~edges, dissolution=~edges, estimate="CMLE", times=c(1,2), eval.loglik=FALSE)
  + 
  + stopifnot(fit$estimate=="CMLE", is.null(fit$formation.fit$sample), is.null(fit$dissolution.fit$sample))
  + stopifnot(z.error(form.mle(y0,y1), fit$formation.fit) <= tolerance)
  + stopifnot(z.error(diss.mle(y0,y1), fit$dissolution.fit) <= tolerance)
  + 
  + # Force CMLE
  + for(prop.weight in prop.weights){
  + cat("====",prop.weight,"====\n")
  + set.seed(544)
  + fit<-stergm(list(y0,y1), formation=~edges, dissolution=~edges, estimate="CMLE", control=control.stergm(CMLE.control=control.ergm(force.main=TRUE, MCMC.prop.weights=prop.weight)), times=c(1,2), eval.loglik=FALSE)
  + 
  + stopifnot(fit$estimate=="CMLE", fit$formation.fit$estimate=="MLE", fit$dissolution.fit$estimate=="MLE")
  + stopifnot(z.error(form.mle(y0,y1), fit$formation.fit) <= tolerance)
  + stopifnot(z.error(diss.mle(y0,y1), fit$dissolution.fit) <= tolerance)
  + }
  + 
  + cat("Missing data:\n")
  + 
  + y1m<-network.copy(y1)
  + set.seed(765)
  + e <- as.edgelist(y1)[1,]
  + y1m[e[1], e[2]] <- NA
  + y1m[1,m+1] <- NA
  + 
  + # Force CMPLE
  + set.seed(765)
  + fit<-stergm(list(y0,y1m), formation=~edges, dissolution=~edges, estimate="CMPLE", times=c(1,2), eval.loglik=FALSE)
  + 
  + stopifnot(fit$estimate=="CMPLE", fit$formation.fit$estimate=="MPLE", fit$dissolution.fit$estimate=="MPLE")
  + stopifnot(z.error(form.mle(y0,y1m), fit$formation.fit) <= tolerance)
  + stopifnot(z.error(diss.mle(y0,y1m), fit$dissolution.fit) <= tolerance)
  + 
  + # Autodetected CMPLE
  + set.seed(765)
  + fit<-stergm(list(y0,y1m), formation=~edges, dissolution=~edges, estimate="CMLE", times=c(1,2), eval.loglik=FALSE)
  + 
  + stopifnot(fit$estimate=="CMLE", is.null(fit$formation.fit$sample), is.null(fit$dissolution.fit$sample))
  + stopifnot(z.error(form.mle(y0,y1m), fit$formation.fit) <= tolerance)
  + stopifnot(z.error(diss.mle(y0,y1m), fit$dissolution.fit) <= tolerance)
  + 
  + # Force CMLE
  + for(prop.weight in prop.weights){
  + cat("====",prop.weight,"====\n")
  + set.seed(234)
  + fit<-stergm(list(y0,y1m), formation=~edges, dissolution=~edges, estimate="CMLE", control=control.stergm(CMLE.control=control.ergm(force.main=TRUE, MCMC.prop.weights=prop.weight)), times=c(1,2), eval.loglik=FALSE)
  + 
  + stopifnot(fit$estimate=="CMLE", fit$formation.fit$estimate=="MLE", fit$dissolution.fit$estimate=="MLE")
  + stopifnot(z.error(form.mle(y0,y1m), fit$formation.fit) <= tolerance)
  + stopifnot(z.error(diss.mle(y0,y1m), fit$dissolution.fit) <= tolerance)
  + }
  + }
  > 
  > cat("=========== Directed test ===========\n")
  =========== Directed test ===========
  > do.run(TRUE, prop.weights=c("default","random"))
  Error in set.vertex.attribute(g, "na", rep(FALSE, n), 1:n) : 
    attempt to set index 0/0 in SET_VECTOR_ELT
  Calls: do.run -> network.initialize -> set.vertex.attribute
  Execution halted
Running the tests in ‘tests/sim_gf_sum.R’ failed.
Complete output:
  > #  File tests/sim_gf_sum.R in package tergm, part of the Statnet suite
  > #  of packages for network analysis, https://statnet.org .
  > #
  > #  This software is distributed under the GPL-3 license.  It is free,
  > #  open source, and has the attribution requirements (GPL Section 7) at
  > #  https://statnet.org/attribution
  > #
  > #  Copyright 2008-2020 Statnet Commons
  > #######################################################################
  > library(tergm)
  Loading required package: ergm
  Loading required package: network
  network: Classes for Relational Data
  Version 1.16.1 created on 2020-10-06.
  copyright (c) 2005, Carter T. Butts, University of California-Irvine
                      Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Martina Morris, University of Washington
                      Skye Bender-deMoll, University of Washington
   For citation information, type citation("network").
   Type help("network-package") to get started.
  
  
  ergm: version 3.11.0, created on 2020-10-14
  Copyright (c) 2020, Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Carter T. Butts, University of California -- Irvine
                      Steven M. Goodreau, University of Washington
                      Pavel N. Krivitsky, UNSW Sydney
                      Martina Morris, University of Washington
                      with contributions from
                      Li Wang
                      Kirk Li, University of Washington
                      Skye Bender-deMoll, University of Washington
                      Chad Klumb
                      Michał Bojanowski, Kozminski University
                      Ben Bolker
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("ergm").
  
  NOTE: Versions before 3.6.1 had a bug in the implementation of the bd()
  constraint which distorted the sampled distribution somewhat. In
  addition, Sampson's Monks datasets had mislabeled vertices. See the
  NEWS and the documentation for more details.
  
  NOTE: Some common term arguments pertaining to vertex attribute and
  level selection have changed in 3.10.0. See terms help for more
  details. Use 'options(ergm.term=list(version="3.9.4"))' to use old
  behavior.
  
  Loading required package: networkDynamic
  
  networkDynamic: version 0.10.1, created on 2020-01-16
  Copyright (c) 2020, Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Pavel N. Krivitsky, University of Wollongong
                      Skye Bender-deMoll, University of Washington
                      with contributions from
                      Zack Almquist, University of California -- Irvine
                      David R. Hunter, Penn State University
                      Li Wang
                      Kirk Li, University of Washington
                      Steven M. Goodreau, University of Washington
                      Jeffrey Horner
                      Martina Morris, University of Washington
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("networkDynamic").
  
  
  tergm: version 3.7.0, created on 2020-10-15
  Copyright (c) 2020, Pavel N. Krivitsky, UNSW Sydney
                      Mark S. Handcock, University of California -- Los Angeles
                      with contributions from
                      David R. Hunter, Penn State University
                      Steven M. Goodreau, University of Washington
                      Martina Morris, University of Washington
                      Nicole Bohme Carnegie, New York University
                      Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Skye Bender-deMoll
                      Li Wang
                      Kirk Li, University of Washington
                      Chad Klumb
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("tergm").
  
  > 
  > logit<-function(p)log(p/(1-p))
  > 
  > coef.form.f<-function(coef.diss,density) -log(((1+exp(coef.diss))/(density/(1-density)))-1)
  > 
  > simtest <- function(S, edges, dur, n, dir=FALSE, bip=0){
  +   m <- if(bip) bip else n
  +   
  +   target.stats<-edges
  +   coef.diss<-logit(1-1/dur)
  +   
  +   dc <- matrix(rnorm(n*m),n,m)
  +   if(!dir && !bip) dc <- dc+t(dc)
  +   # Note that if nw1!=nw2, then
  +   # summary(nw1~edgecov(dc)+edgecov.ages(dc))!=summary(nw2~edgecov(dc)+edgecov.ages(dc))
  +   # with probability 1. Thus, they offer a sort of a checksum of a network.
  +   
  +   g0<-network.initialize(n, dir=dir, bipartite=bip)
  +   g0 %n% "dc" <- dc
  +   g0 %n% "time" <- 0
  +   g0 %n% "lasttoggle" <- -1-rgeom(network.dyadcount(g0),1/4)
  + 
  + 
  +   
  +   dyads <- network.dyadcount(g0)
  +   density<-edges/dyads
  +   
  +   coef.form<-coef.form.f(coef.diss,density)
  +   
  +   # Get a reasonably close starting network.
  +   g1<-san(g0~edges,target.stats=target.stats,verbose=TRUE)
  +   
  +   print(coef.form)
  +   print(coef.diss)
  +   
  +   # Simulate. Starting from an ordinary network:
  +   dynsim<-simulate(g1,formation=~edges,dissolution=~edges,coef.form=coef.form,coef.diss=coef.diss,time.slices=S,verbose=TRUE,monitor=~edgecov("dc")+edgecov.ages("dc"))
  +   
  +   # Resuming from a networkDynamic:
  +   dynsim2<-simulate(dynsim,formation=~edges,dissolution=~edges,coef.form=coef.form,coef.diss=coef.diss,time.slices=S,verbose=TRUE,monitor=~edgecov("dc")+edgecov.ages("dc"))
  +   
  +   # Resuming from a resumed networkDynamic:
  +   dynsim3<-simulate(dynsim2,formation=~edges,dissolution=~edges,coef.form=coef.form,coef.diss=coef.diss,time.slices=S,verbose=TRUE,monitor=~edgecov("dc")+edgecov.ages("dc"))
  + 
  +   sim.stats <- attr(dynsim3,"stats")
  +   
  +   # Replay the simulation using a networkDynamic:
  +   gf1.stats <- as.matrix(tergm.godfather(dynsim3~edgecov("dc")+edgecov.ages("dc"), start=0, end=S*3))
  +   
  +   # Replay the sim using networkDynamic with no explicit time params (should be using net.obs.period)
  +   gf1b.stats <- as.matrix(tergm.godfather(dynsim3~edgecov("dc")+edgecov.ages("dc")))
  +   
  +   # Replay the simulation using a network + list of changes:
  +   gf2.stats <- as.matrix(tergm.godfather(g1~edgecov("dc")+edgecov.ages("dc"), start=0, end=S*3, changes=attr(dynsim3,"changes")))
  +   
  +   # Evaluate summary statistics directly:
  +   summ.stats <- summary(dynsim3~edgecov("dc")+edgecov.ages("dc"), at=1:(S*3))
  +   
  +   # If they aren't all identical, we are in trouble.
  +   stopifnot(all.equal(sim.stats,gf1.stats),
  +             all.equal(sim.stats,gf1b.stats),
  +             all.equal(sim.stats,gf2.stats),
  +             all.equal(sim.stats,summ.stats))
  + }
  > 
  > set.seed(1)
  > cat("Undirected, unipartite:\n")
  Undirected, unipartite:
  > simtest(20, 10, 10, 10, FALSE, FALSE)
  Error in set.vertex.attribute(g, "na", rep(FALSE, n), 1:n) : 
    attempt to set index 0/0 in SET_VECTOR_ELT
  Calls: simtest -> network.initialize -> set.vertex.attribute
  Execution halted
Running the tests in ‘tests/simulate_networkDynamic.R’ failed.
Complete output:
  > #  File tests/simulate_networkDynamic.R in package tergm, part of the Statnet suite
  > #  of packages for network analysis, https://statnet.org .
  > #
  > #  This software is distributed under the GPL-3 license.  It is free,
  > #  open source, and has the attribution requirements (GPL Section 7) at
  > #  https://statnet.org/attribution
  > #
  > #  Copyright 2008-2020 Statnet Commons
  > #######################################################################
  > library(tergm)
  Loading required package: ergm
  Loading required package: network
  network: Classes for Relational Data
  Version 1.16.1 created on 2020-10-06.
  copyright (c) 2005, Carter T. Butts, University of California-Irvine
                      Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Martina Morris, University of Washington
                      Skye Bender-deMoll, University of Washington
   For citation information, type citation("network").
   Type help("network-package") to get started.
  
  
  ergm: version 3.11.0, created on 2020-10-14
  Copyright (c) 2020, Mark S. Handcock, University of California -- Los Angeles
                      David R. Hunter, Penn State University
                      Carter T. Butts, University of California -- Irvine
                      Steven M. Goodreau, University of Washington
                      Pavel N. Krivitsky, UNSW Sydney
                      Martina Morris, University of Washington
                      with contributions from
                      Li Wang
                      Kirk Li, University of Washington
                      Skye Bender-deMoll, University of Washington
                      Chad Klumb
                      Michał Bojanowski, Kozminski University
                      Ben Bolker
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("ergm").
  
  NOTE: Versions before 3.6.1 had a bug in the implementation of the bd()
  constraint which distorted the sampled distribution somewhat. In
  addition, Sampson's Monks datasets had mislabeled vertices. See the
  NEWS and the documentation for more details.
  
  NOTE: Some common term arguments pertaining to vertex attribute and
  level selection have changed in 3.10.0. See terms help for more
  details. Use 'options(ergm.term=list(version="3.9.4"))' to use old
  behavior.
  
  Loading required package: networkDynamic
  
  networkDynamic: version 0.10.1, created on 2020-01-16
  Copyright (c) 2020, Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Pavel N. Krivitsky, University of Wollongong
                      Skye Bender-deMoll, University of Washington
                      with contributions from
                      Zack Almquist, University of California -- Irvine
                      David R. Hunter, Penn State University
                      Li Wang
                      Kirk Li, University of Washington
                      Steven M. Goodreau, University of Washington
                      Jeffrey Horner
                      Martina Morris, University of Washington
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("networkDynamic").
  
  
  tergm: version 3.7.0, created on 2020-10-15
  Copyright (c) 2020, Pavel N. Krivitsky, UNSW Sydney
                      Mark S. Handcock, University of California -- Los Angeles
                      with contributions from
                      David R. Hunter, Penn State University
                      Steven M. Goodreau, University of Washington
                      Martina Morris, University of Washington
                      Nicole Bohme Carnegie, New York University
                      Carter T. Butts, University of California -- Irvine
                      Ayn Leslie-Cook, University of Washington
                      Skye Bender-deMoll
                      Li Wang
                      Kirk Li, University of Washington
                      Chad Klumb
  Based on "statnet" project software (statnet.org).
  For license and citation information see statnet.org/attribution
  or type citation("tergm").
  
  > 
  > set.seed(0)
  > 
  > logit<-function(p)log(p/(1-p))
  > coef.form.f<-function(coef.diss,density) -log(((1+exp(coef.diss))/(density/(1-density)))-1)
  > 
  > # Construct a network with 20 nodes and 20 edges
  > n<-20
  > target.stats<-edges<-20
  > g0<-network.initialize(n,dir=TRUE)
  Error in set.vertex.attribute(g, "na", rep(FALSE, n), 1:n) : 
    attempt to set index 0/0 in SET_VECTOR_ELT
  Calls: network.initialize -> set.vertex.attribute
  Execution halted