Hi, thanks for your excellent work on graph generation.

I have some problems during the implementation of the baseline method ERGM in your paper. I followed your parameter settings in R, but the simulation results are always bad (such as the explosion of the number of edges).

This is my R code for graph simulation by ERGM:
result.network <- simulate( ergm( input.network ~ edges + degcor + degree1.5 + gwesp + density, control=control.ergm( MCMLE.maxit=20, MCMLE.density.guard=100000, parallel=2, parallel.type="PSOCK" )) )
where 'input.network' is the graph to be simulated, and the function 'simulate' is from package 'ergm'.
I wonder how to have the right hyper-parameters to have the best simulation result like your paper's.