drop1Wald: Drop Single Terms of a Model and Calculate Respective Wald...

In regr0: Building regression models

Description

'drop1Wald' calculates tests for single term deletions based on the covariance matrix of estimated coefficients instead of re-fitting a reduced model. This helps in cases where re-fitting is not feasible, inappropriate or costly.

Usage

drop1Wald(object, scope=NULL, scale = 0, test = c("none", "Chisq", "F"),
  k = 2, ...)

Arguments

object	a fitted model of class 'regr'.
scope	a formula giving the terms to be considered for dropping. If 'NULL', 'drop.scope(object)' is obtained
scale	an estimate of the residual mean square to be used in computing Cp. Ignored if '0' or 'NULL'.
test	see drop1
k	the penalty constant in AIC / Cp.
...	further arguments, ignored

Details

The test statistics and Cp and AIC values are calculated on the basis of the estimated coefficients and their (unscaled) covariance matrix as provided by the fitting object. The function may be used for all model fitting objects that contain these two components as $coefficients and $cov.unscaled.

Value

An object of class 'anova' summarizing the differences in fit between the models.

Note

drop1Wald is used in regr for models of class 'lm' or 'lmrob' for preparing the 'testcoef' table.

Author(s)

Werner A. Stahel, Seminar for Statistics, ETH Zurich

See Also

drop1

Examples

  data(d.blast)
  r.blast <-
       regr(log10(tremor)~location+log10(distance)+log10(charge),
            data=d.blast)
  drop1(r.blast)
  rr.blast <-
       regr(log10(tremor)~location+log10(distance)+log10(charge),
            data=d.blast, robust=TRUE)
  drop1(rr.blast)

regr0 documentation built on May 2, 2019, 4:52 p.m.