R/MCMC_master.R

In NSUM: Network Scale Up Method

nsum.simulate <- function(n, known, unknown, N, model="degree", ...)
{
  if(model=="degree") .simulate.rd(n, known, unknown, N, ...)
  else if(model=="barrier") .simulate.bar(n, known, unknown, N, ...)
  else if(model=="transmission") .simulate.trans(n, known, unknown, N, ...)
  else if(model=="combined") .simulate.comb(n, known, unknown, N, ...)
  else stop('Error: model must be one of \"degree\", \"barrier\", \"transmission\", or \"combined\"')
}  

nsum.mcmc <- function(dat, known, N, indices.k=(length(known)+1):(dim(dat)[2]), iterations=1000, burnin=100, size=iterations, model="degree", ...)
{
  if(model=="degree") .mcmc.rd(dat, known, N, indices.k, iterations, burnin, size, ...)
  else if(model=="barrier") .mcmc.bar(dat, known, N, indices.k, iterations, burnin, size, ...)
  else if(model=="transmission") .mcmc.trans(dat, known, N, indices.k, iterations, burnin, size, ...)
  else if(model=="combined") .mcmc.comb(dat, known, N, indices.k, iterations, burnin, size, ...)
  else stop('Error: model must be one of \"degree\", \"barrier\", \"transmission\", or \"combined\"')
}