ANOVA for ERGM Fits
Compute an analysis of variance table for one or more ERGM fits.
## S3 method for class 'ergm' anova(object, ..., eval.loglik = FALSE) ## S3 method for class 'ergmlist' anova(object, ..., eval.loglik = FALSE, scale = 0, test = "F")
object, ... |
|
eval.loglik |
a logical specifying whether the log-likelihood will be evaluated if missing. |
scale |
numeric. An estimate of the noise variance sigma^2. If zero this will be estimated from the largest model considered. |
test |
a character string specifying the test statistic to be used. Can
be one of |
Specifying a single object gives a sequential analysis of variance table for that fit. That is, the reductions in the residual sum of squares as each term of the formula is added in turn are given in the rows of a table, plus the residual sum of squares.
The table will contain F statistics (and P values) comparing the mean square for the row to the residual mean square.
If more than one object is specified, the table has a row for the residual degrees of freedom and sum of squares for each model. For all but the first model, the change in degrees of freedom and sum of squares is also given. (This only make statistical sense if the models are nested.) It is conventional to list the models from smallest to largest, but this is up to the user.
Optionally the table can include test statistics. Normally the F statistic
is most appropriate, which compares the mean square for a row to the
residual sum of squares for the largest model considered. If scale
is specified chi-squared tests can be used. Mallows' Cp statistic
is the residual sum of squares plus twice the estimate of
sigma^2 times the residual degrees of freedom.
If any of the objects do not have estimated log-likelihoods, produces an
error, unless eval.loglik=TRUE.
An object of class "anova" inheriting from class
"data.frame".
The comparison between two or more models will only be
valid if they are fitted to the same dataset. This may be a problem if there
are missing values and 's default of na.action = na.omit is used, and
anova.ergmlist will detect this with an error.
The model fitting function ergm, anova,
logLik.ergm for adding the log-likelihood to an existing
ergm object.
data(molecule)
molecule %v% "atomic type" <- c(1,1,1,1,1,1,2,2,2,2,2,2,2,3,3,3,3,3,3,3)
fit0 <- ergm(molecule ~ edges)
anova(fit0)
fit1 <- ergm(molecule ~ edges + nodefactor("atomic type"))
anova(fit1)
fit2 <- ergm(molecule ~ edges + nodefactor("atomic type") + gwesp(0.5,
fixed=TRUE), eval.loglik=TRUE) # Note the eval.loglik argument.
anova(fit0, fit1)
anova(fit0, fit1, fit2)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.