glm.control {stats}R Documentation

Auxiliary for Controlling GLM Fitting

Description

Auxiliary function as user interface for glm fitting. Typically only used when calling glm or glm.fit.

Usage

glm.control(epsilon = 1e-8, maxit = 25, trace = FALSE)

Arguments

epsilon positive convergence tolerance ε; the iterations converge when |dev - devold|/(|dev| + 0.1) < ε.
maxit integer giving the maximal number of IWLS iterations.
trace logical indicating if output should be produced for each iteration.

Details

If epsilon is small, it is also used as the tolerance for the least squares solution.

When trace is true, calls to cat produce the output for each IWLS iteration. Hence, options(digits = *) can be used to increase the precision, see the example.

Value

A list with the arguments as components.

References

Hastie, T. J. and Pregibon, D. (1992) Generalized linear models. Chapter 6 of Statistical Models in S eds J. M. Chambers and T. J. Hastie, Wadsworth & Brooks/Cole.

See Also

glm.fit, the fitting procedure used by glm.

Examples

### A variation on  example(glm) :

## Annette Dobson's example ...
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
oo <- options(digits = 12) # to see more when tracing :
glm.D93X <- glm(counts ~ outcome + treatment, family=poisson(),
                trace = TRUE, epsilon = 1e-14)
options(oo)
coef(glm.D93X) # the last two are closer to 0 than in ?glm's  glm.D93
# put less so than in R < 1.8.0 when the default was 1e-4

[Package stats version 2.2.1 Index]