outlier.test {car}R Documentation

Bonferroni Outlier Test

Description

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

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.

Author(s)

John Fox jfox@mcmaster.ca

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.

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 

[Package car version 1.1-0 Index]