Centre of the Parameter Space for GLMs

Computes the COPS for one or more parameters in a fitted generalized linear model. The object has "glm". The functions comprise an S3 generic and method.

Usage

cops3(object, ...)
# Default S3 method
cops3(object, ...)
# S3 method for class 'glm'
cops3(object, level = 0.999,
  beta.range = confint.default(object, level = level),
  muxrange = 4, iter.max = 8, tol = 1e-5,
  subset = NULL, do1 = TRUE, ...)

Arguments

object: a fitted GLM, created by glm.
level,beta.range: Same as copsvglm.
muxrange,iter.max: Same as copsvglm.
tol, subset, do1: Same as copsvglm.
...: additional argument(s) for methods.

Value

See copsvglm.

Details

While VGAM is written in S4, this cops3 is an S3 generic function so that the COPS can be computed for glm models. Since vglm and glm have different convergence criteria, computing the COPS for effectively the same model may differ.

Note that not all GLMs have a finite COPS, for example, a poisson regression.

Examples

if (FALSE) { # \dontrun{
fit1 <- glm(cbind(dead, n-dead) ~ logdose, binomial, fbeetle, tr = TRUE)
coef(fit1)
cops3(fit1)
fit2 <- vglm(cbind(dead, n-dead) ~ logdose, binomialff, fbeetle, tr = TRUE)
coef(fit2)
cops(fit2)
} # }