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

Description Usage Arguments Details Value Note Author(s) See Also Examples

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

1
2
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

1
2
3
4
5
6
7
8
9
  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)


Search within the regr0 package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.