outlier.test

0th

Percentile

Bonferroni Outlier Test

Reports the Bonferroni p-value for the most extreme observation. At present, there are methods for studentized residuals in linear and generalized linear models.

Keywords
regression, htest
Usage
outlier.test(model, ...)

## S3 method for class 'lm':
outlier.test(model, labels=names(rstud), ...)

## S3 method for class 'glm':
outlier.test(model, labels=names(rstud), ...)

## S3 method for class 'outlier.test':
print(x, digits=options("digits")[[1]], ...)
Arguments
model
a suitable model object.
labels
an optional vector of observation names.
...
arguments passed down to methods functions.
x
outlier.test object.
digits
number of digits for printed output.
Details

For a linear model, the p-value reported is for the largest absolute studentized residual, using the $t$ distribution with degrees of freedom one less than the residual df for the model. For a generalized linear model, the largest absolute studentized residual is also used, but with the standard-normal distribution. The Bonferroni adjustment multiplies the usual two-sided p-value by the number of observations.

Value

  • an object of class outlier.test, which is normally just printed.

References

Cook, R. D. and Weisberg, S. (1984) Residuals and Influence in Regression. Wiley. Fox, J. (1997) Applied Regression, Linear Models, and Related Methods. Sage. Williams, D. A. (1987) Generalized linear model diagnostics using the deviance and single case deletions. Applied Statistics 36, 181--191.

Aliases
  • outlier.test
  • outlier.test.lm
  • outlier.test.glm
  • print.outlier.test
Examples
outlier.test(lm(prestige~income+education, data=Duncan))
##  max|rstudent| df unadjusted p Bonferroni p
##       3.134519 41  0.003177202    0.1429741
## 
## Observation: minister
Documentation reproduced from package car, version 1.2-10, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.