ergmAPL: Adjustment of ERGM pseudolikelihood

Description

Function to estimate the transformation parameters for adjusting the pseudolikelihood function.

Usage

ergmAPL(
  formula,
  aux.iters = NULL,
  n.aux.draws = NULL,
  aux.thin = NULL,
  ladder = NULL,
  estimate = c("MLE", "CD"),
  seed = 1,
  ...
)

Arguments

formula: formula; an ergm formula object, of the form <network> ~ <model terms> where <network> is a network object and <model terms> are ergm-terms.
aux.iters: count; number of auxiliary iterations used for drawing the first network from the ERGM likelihood. See control.simulate.formula.
n.aux.draws: count; Number of auxiliary networks drawn from the ERGM likelihood. See control.simulate.formula.
aux.thin: count; Number of auxiliary iterations between network draws after the first network is drawn. See control.simulate.formula.
ladder: count; Length of temperature ladder (>=3).
estimate: If "MLE" (the default), then an approximate maximum likelihood estimator is returned. If "CD" , the Monte-Carlo contrastive divergence estimate is returned. See ergm.
seed: integer; seed for the random number generator. See set.seed.
...: Additional arguments, to be passed to the ergm function. See ergm.

References

Bouranis, L., Friel, N., & Maire, F. (2018). Bayesian model selection for exponential random graph models via adjusted pseudolikelihoods. Journal of Computational and Graphical Statistics, 27(3), 516-528. https://arxiv.org/abs/1706.06344