summary
method for objects of class "gnm"
## S3 method for class 'gnm':
summary(object, dispersion = NULL, correlation = FALSE,
symbolic.cor = 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.summary.gnm
returns an object of class "summary.gnm"
,
which is a list with components"call"
component from object."eliminate"
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
.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