logLik.ergm | R Documentation |
logLik
method for ergm
fits.A function to return the log-likelihood associated with an
ergm
fit, evaluating it if
necessary. If the log-likelihood was not computed for
object
, produces an error unless eval.loglik=TRUE
.
## S3 method for class 'ergm' logLik( object, add = FALSE, force.reeval = FALSE, eval.loglik = add || force.reeval, control = control.logLik.ergm(), ... ) ## S3 method for class 'ergm' deviance(object, ...) ## S3 method for class 'ergm' AIC(object, ..., k = 2) ## S3 method for class 'ergm' BIC(object, ...)
object |
An |
add |
Logical: If |
force.reeval |
Logical: If |
eval.loglik |
Logical: If |
control |
A list of control parameters for algorithm tuning,
typically constructed with |
... |
Other arguments to the likelihood functions. |
k |
see help for |
The form of the output of logLik.ergm
depends on
add
: add=FALSE
(the default), a
logLik
object. If add=TRUE
(the default), an
ergm
object with the log-likelihood
set.
As of version 3.1, all likelihoods for which logLikNull
is
not implemented are computed relative to the reference
measure. (I.e., a null model, with no terms, is defined to have
likelihood of 0, and all other models are defined relative to
that.)
deviance.ergm
: A deviance()
method.
AIC.ergm
: An AIC()
method.
BIC.ergm
: A BIC()
method.
Hunter, D. R. and Handcock, M. S. (2006) Inference in curved exponential family models for networks, Journal of Computational and Graphical Statistics.
logLik
, logLikNull
, ergm.bridge.llr
,
ergm.bridge.dindstart.llk
# See help(ergm) for a description of this model. The likelihood will # not be evaluated. data(florentine) ## Not run: # The default maximum number of iterations is currently 20. We'll only # use 2 here for speed's sake. gest <- ergm(flomarriage ~ kstar(1:2) + absdiff("wealth") + triangle, eval.loglik=FALSE) gest <- ergm(flomarriage ~ kstar(1:2) + absdiff("wealth") + triangle, eval.loglik=FALSE, control=control.ergm(MCMLE.maxit=2)) # Log-likelihood is not evaluated, so no deviance, AIC, or BIC: summary(gest) # Evaluate the log-likelihood and attach it to the object. # The default number of bridges is currently 20. We'll only use 3 here # for speed's sake. gest.logLik <- logLik(gest, add=TRUE) gest.logLik <- logLik(gest, add=TRUE, control=control.logLik.ergm(bridge.nsteps=3)) # Deviances, AIC, and BIC are now shown: summary(gest.logLik) # Null model likelihood can also be evaluated, but not for all constraints: logLikNull(gest) # == network.dyadcount(flomarriage)*log(1/2) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.