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