summary method for objects of class "gnm"## S3 method for class 'gnm':
summary(object, dispersion = NULL, correlation = FALSE,
symbolic.cor = FALSE, with.eliminate = FALSE, ...)## S3 method for class 'summary.gnm':
print(x, digits = max(3, getOption("digits") - 3),
signif.stars = getOption("show.signif.stars"),
symbolic.cor = x$symbolic.cor, ...)
"gnm"."summary.gnm".object.TRUE, the correlation matrix of
the estimated parameters is returned.TRUE, the correlations are
printed in a symbolic form rather than numbers (see
symnum).TRUE, "significance stars" are
printed for each coefficient.TRUE, any eliminated
coefficients are included in the summary.summary.gnm returns an object of class "summary.gnm",
which is a list with components"call" component from object."ofInterest" component from object."family" component from object."deviance" component from object."aic" component from object."df.residual" component from object."iter" component from object.residuals.glm.with.eliminate = TRUE a matrix of
eliminated coefficients, standard errors, z-values and p-values.NULL.dispersion (see vcov.gnm for more details).correlation is TRUE) the
estimated correlations of the estimated coefficients.correlation is TRUE) the value
of the argument symbolic.cor.print.summary.gnm prints the original call to gnm; a
summary of the deviance residuals from the model fit; the coefficients
of the model; the residual deviance; the Akaike's Information
Criterion value, and the number of main iterations performed. Standard errors, z-values and p-values are printed alongside the
coefficients, with "significance stars" if signif.stars is
TRUE.
When the "summary.gnm" object has a "correlation"
component, the lower triangle of this matrix is also printed, to two
decimal places (or symbolically); to see the full matrix of
correlations print summary(object, correlation =
TRUE)$correlation directly.
The standard errors returned by summary.gnm are scaled by
sqrt(dispersion). If the dispersion is not specified, it is
taken as 1 for the binomial and Poisson families,
and otherwise estimated by the residual Chi-squared statistic divided
by the residual degrees of freedom. For coefficients that have been
constrained or are not estimable, the standard error is returned as
NA.
gnm, summary## Following on from example(gnm)
data(cautres)
set.seed(1)
## Fit model as before
doubleUnidiff <- gnm(Freq ~ election:vote + election:class:religion +
Mult(Exp(election), religion:vote) +
Mult(Exp(election), class:vote), family = poisson,
data = cautres)
## Summarize results
summary(doubleUnidiff)Run the code above in your browser using DataLab