VGAM (version 1.1-14)

cops3: Centre of the Parameter Space for GLMs

Description

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, ...)
# S3 method for default
cops3(object, ...)
# S3 method for 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, ...)

Value

See copsvglm.

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.

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.

See Also

copsvglm, confint.default, wsdm, fbeetle.

Examples

if (FALSE) {
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)
}