Summarizing double generalized linear model fits
Description
These functions are all methods
for class dglm
or
summary.glm
objects.
Usage
1 2 
Arguments
object 
an object of class 
dispersion 
the dispersion parameter for the fitting family.
By default it is obtained from 
correlation 
logical; if 
... 
further arguments to be passed to 
Details
For more details, see summary.glm
.
If more than one of etastart
, start
and mustart
is specified, the first in the list will be used.
Value
summary.dglm
returns an object of class
"summary.dglm"
, a list with components
call 
the component from 
terms 
the component from 
family 
the component from 
deviance 
the component from 
aic 

constrasts 
(where relevant) the contrasts used. NOT WORKING?? 
df.residual 
the component from 
null.deviance 
the component from 
df.null 
the residual degrees of freedom for the null model. 
iter 
the component from 
deviance.resid 
the deviance residuals: see 
coefficients 
the matrix of coefficients, standard errors, zvalues and pvalues. Aliased coefficients are omitted. 
aliased 
named logical vector showing if the original coefficients are aliased. 
dispersion 
either the supplied argument or the estimated dispersion
if the latter in 
df 
a 3vector of the rank of the model and the number of residual degrees of freedom, plus number of nonaliased coefficients. 
cov.unscaled 
the unscaled ( 
cov.scaled 
ditto, scaled by 
correlation 
(only if 
dispersion.summary 
the summary of the fitted dispersion model 
outer.iter 
the number of outer iteration of the alternating iterations 
m2loglik 
minus twice the loglikelihood of the fitted model 
Note
The anova method is questionable when applied to an dglm
object with
method="reml"
(stick to method="ml"
).
Author(s)
Gordon Smyth, ported to R\ by Peter Dunn (pdunn2@usc.edu.au)
References
Smyth, G. K. (1989). Generalized linear models with varying dispersion. J. R. Statist. Soc. B, 51, 47–60.
Smyth, G. K., and Verbyla, A. P. (1999). Adjusted likelihood methods for modelling dispersion in generalized linear models. Environmetrics, 10, 696709.
Verbyla, A. P., and Smyth, G. K. (1998). Double generalized linear models: approximate residual maximum likelihood and diagnostics. Research Report, Department of Statistics, University of Adelaide.
See Also
dglm.object
, dglm.control
,
anova.dglm
,
summary.glm
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22  # Continuing the example from glm, but this time try
# fitting a Gamma double generalized linear model also.
clotting < data.frame(
u = c(5,10,15,20,30,40,60,80,100),
lot1 = c(118,58,42,35,27,25,21,19,18),
lot2 = c(69,35,26,21,18,16,13,12,12))
# The same example as in glm: the dispersion is modelled as constant
out < dglm(lot1 ~ log(u), ~1, data=clotting, family=Gamma)
summary(out)
# Try a double glm
out2 < dglm(lot1 ~ log(u), ~u, data=clotting, family=Gamma)
summary(out2)
anova(out2)
# Summarize the mean model as for a glm
summary.glm(out2)
# Summarize the dispersion model as for a glm
summary(out2$dispersion.fit)
