gof {ergm} | R Documentation |
gof
calculates p-values for geodesic
distance, degree, and reachability summaries to
diagnose the goodness-of-fit of exponential family random graph
models. See ergm
for more information on these models.
## Default S3 method: gof(object,...) ## S3 method for class 'formula' gof(formula, ..., theta0=NULL, nsim=100, burnin=10000, interval=1000, GOF=NULL, constraints=~., control=control.gof.formula(), seed=NULL, verbose=FALSE) ## S3 method for class 'ergm' gof(object, ..., nsim=100, GOF=NULL, burnin=10000, interval=1000, constraints=NULL, control=control.gof.ergm(), seed=NULL, theta0=NULL, verbose=FALSE)
object |
an R object. Either a formula or an |
formula |
formula; An R formula object, of the form
|
theta0 |
When given either a formula or an object of class ergm, theta0 are the parameters from which the sample is drawn. By default set to a vector of 0. |
nsim |
The number of simulations to use for the MCMC p-values. This is the size of the sample of graphs to be randomly drawn from the distribution specified by the object on the set of all graphs. |
GOF |
formula; an R formula object, of the form
|
burnin |
Number of proposed edge toggles before any MCMC sampling is done. If the model is correct this can be 0. Currently, there is no support for any check of the Markov chain mixing, so burnin should be set to a fairly large number. |
interval |
Number of proposed edge toggles between sampled statistics. The program prints a warning if too few proposed toggles are being accepted (if the number of proposed toggles between sampled observations ever equals an integral multiple of 100*(1+the number of toggles accepted)). |
constraints |
A one-sided formula specifying one or more constraints
on the support of the distribution of the networks being
modeled. See the help for similarly-named argument in
|
control |
A list to control parameters, constructed using
|
seed |
integer; random number integer seed. Defaults to |
verbose |
Provide verbose information on the progress of the simulation. |
... |
Additional arguments, to be passed to lower-level functions in the future. |
A sample of graphs is randomly drawn from the specified model.
The first argument is typically
the output of a call to ergm
and the model
used for that call is the one fit.
A plot of the summary measures is plotted.
More information can be found by looking at the documentation of
ergm
.
For gof.ergm
and gof.formula
, default behavior depends on
the directedness of the network involved; if undirected then degree,
espartners, and distance are used as default properties to examine. If
the network in question is directed, “degree” in the above is replaced
by idegree and odegree.
gof
, gof.ergm
, and gof.formula
return an object of class gofobject
.
This is a list of the tables of statistics and p-values.
This is typically plotted using plot.gofobject
.
ergm, network, simulate.ergm, summary.ergm, plot.gofobject
data(florentine) gest <- ergm(flomarriage ~ edges + kstar(2)) gest summary(gest) # test the gof.ergm function gofflo <- gof(gest) gofflo # Plot all three on the same page # with nice margins par(mfrow=c(1,3)) par(oma=c(0.5,2,1,0.5)) plot(gofflo) # And now the log-odds plot(gofflo, plotlogodds=TRUE) # Use the formula version of gof gofflo2 <-gof(flomarriage ~ edges + kstar(2), theta0=c(-1.6339, 0.0049)) plot(gofflo2)