ergmAPL {Bergm}

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