R/simdist.R

In degreenet: Models for Skewed Count Distributions Relevant to Networks

#  File degreenet/R/simdist.R
#  Part of the statnet package, http://statnet.org
#
#  This software is distributed under the GPL-3 license.  It is free,
#  open source, and has the attribution requirements (GPL Section 7) in
#    http://statnet.org/attribution
#
# Copyright 2003 Mark S. Handcock, University of California-Los Angeles
# Copyright 2007 The statnet Development Team
######################################################################
#
# These are the basic simulation routines
#
simyule <- function(n=100, rho=4, maxdeg=10000){
  sample(x=1:maxdeg, size=n, replace=TRUE, prob=dyule(v=rho,x=1:maxdeg))
}
simnb <- function(n=100, v=c(5,0.2), maxdeg=10000){
  sample(x=1:maxdeg, size=n, replace=TRUE,
         prob=dnbinom(x=1:maxdeg, size=v[2]*v[1], prob=v[2]))
}
simwar <- function(n=100, v=c(3.5,0.1), maxdeg=10000){
  sample(x=1:maxdeg, size=n, replace=TRUE,
         prob=dwar(v=v,x=1:maxdeg))
}
simdqe <- function(n=100, v=c(3.5,1), maxdeg=10000){
  sample(x=1:maxdeg, size=n, replace=TRUE,
         prob=ddqe(v=v,x=1:maxdeg))
}
simdp <- function(n=100, v=3.5, maxdeg=10000){
  sample(x=1:maxdeg, size=n, replace=TRUE,
         prob=ddp(v=v,x=1:maxdeg))
}
simcmp <- function(n=100, v=c(7,2.6), maxdeg=10000){
  sample(x=1:maxdeg, size=n, replace=TRUE,
         prob=dcmp_mu(v=v,x=1:maxdeg))
}