N <- 500
zeta <- 1
# Generate the spherical coordinates
z <- rvMF(N, c(0, 0, 0))
# Genetate the gregariousness
nu <- rnorm(N, -1.35, 0.37)
# Generate degrees
d <- runif(N, 0, 45)
dist <- sim.dnetwork(nu, d, zeta, z)
Run the code above in your browser using DataLab