To install click the Add extension button. That's it.

The source code for the WIKI 2 extension is being checked by specialists of the Mozilla Foundation, Google, and Apple. You could also do it yourself at any point in time.

4,5
Kelly Slayton
Congratulations on this excellent venture… what a great idea!
Alexander Grigorievskiy
I use WIKI 2 every day and almost forgot how the original Wikipedia looks like.
Live Statistics
English Articles
Improved in 24 Hours
Added in 24 Hours
Languages
Recent
Show all languages
What we do. Every page goes through several hundred of perfecting techniques; in live mode. Quite the same Wikipedia. Just better.
.
Leo
Newton
Brights
Milds

# Studentized residual

In statistics, a studentized residual is the quotient resulting from the division of a residual by an estimate of its standard deviation. Typically the standard deviations of residuals in a sample vary greatly from one data point to another even when the errors all have the same standard deviation, particularly in regression analysis; thus it does not make sense to compare residuals at different data points without first studentizing. It is a form of a Student's t-statistic, with the estimate of error varying between points.

This is an important technique in the detection of outliers. It is among several named in honor of William Sealey Gosset, who wrote under the pseudonym Student. Dividing a statistic by a sample standard deviation is called studentizing, in analogy with standardizing and normalizing.

• 1/5
Views:
3 802
7 563
5 423
4 510
17 780
• ✪ Lecture21 (Data2Decision) Leverage in Regression
• ✪ Outlier analysis in linear regression
• ✪ SPSS regression residuals - unstandardized, standardized, studentized
• ✪ Grubbs Outlier Test - Introduced and Demonstrated
• ✪ Statistics with R: Example of outlier and leverage analysis part 1 of 3

## Motivation

The key reason for studentizing is that, in regression analysis of a multivariate distribution, the variances of the residuals at different input variable values may differ, even if the variances of the errors at these different input variable values are equal. The issue is the difference between errors and residuals in statistics, particularly the behavior of residuals in regressions.

Consider the simple linear regression model

${\displaystyle Y=\alpha _{0}+\alpha _{1}X+\varepsilon .\,}$

Given a random sample (XiYi), i = 1, ..., n, each pair (XiYi) satisfies

${\displaystyle Y_{i}=\alpha _{0}+\alpha _{1}X_{i}+\varepsilon _{i},\,}$

where the errors ${\displaystyle \varepsilon _{i}}$, are independent and all have the same variance ${\displaystyle \sigma ^{2}}$. The residuals are not the true errors, but estimates, based on the observable data. When the method of least squares is used to estimate ${\displaystyle \alpha _{0}}$ and ${\displaystyle \alpha _{1}}$, then the residuals ${\displaystyle {\widehat {\varepsilon \,}}}$, unlike the errors ${\displaystyle \varepsilon }$, cannot be independent since they satisfy the two constraints

${\displaystyle \sum _{i=1}^{n}{\widehat {\varepsilon \,}}_{i}=0}$

and

${\displaystyle \sum _{i=1}^{n}{\widehat {\varepsilon \,}}_{i}x_{i}=0.}$

(Here εi is the ith error, and ${\displaystyle \scriptstyle {\widehat {\varepsilon \,}}_{i}}$ is the ith residual.)

The residuals, unlike the errors, do not all have the same variance: the variance decreases as the corresponding x-value gets farther from the average x-value. This is not a feature of the data itself, but of the regression better fitting values at the ends of the domain. It is also reflected in the influence functions of various data points on the regression coefficients: endpoints have more influence. This can also be seen because the residuals at endpoints depend greatly on the slope of a fitted line, while the residuals at the middle are relatively insensitive to the slope. The fact that the variances of the residuals differ, even though the variances of the true errors are all equal to each other, is the principal reason for the need for studentization.

It is not simply a matter of the population parameters (mean and standard deviation) being unknown – it is that regressions yield different residual distributions at different data points, unlike point estimators of univariate distributions, which share a common distribution for residuals.

## Background

For this simple model, the design matrix is

${\displaystyle X=\left[{\begin{matrix}1&x_{1}\\\vdots &\vdots \\1&x_{n}\end{matrix}}\right]}$

and the hat matrix H is the matrix of the orthogonal projection onto the column space of the design matrix:

${\displaystyle H=X(X^{T}X)^{-1}X^{T}.\,}$

The leverage hii is the ith diagonal entry in the hat matrix. The variance of the ith residual is

${\displaystyle \operatorname {var} ({\widehat {\varepsilon \,}}_{i})=\sigma ^{2}(1-h_{ii}).}$

In case the design matrix X has only two columns (as in the example above), this is equal to

${\displaystyle \operatorname {var} ({\widehat {\varepsilon \,}}_{i})=\sigma ^{2}\left(1-{\frac {1}{n}}-{\frac {(x_{i}-{\bar {x}})^{2}}{\sum _{j=1}^{n}(x_{j}-{\bar {x}})^{2}}}\right).}$

## Calculation

Given the definitions above, the Studentized residual is then

${\displaystyle t_{i}={{\widehat {\varepsilon \,}}_{i} \over {\widehat {\sigma }}{\sqrt {1-h_{ii}\ }}}}$

where ${\displaystyle {\widehat {\sigma }}}$ is an appropriate estimate of σ (see below).

## Internal and external studentization

The usual estimate of σ2 is the internally studentized residual

${\displaystyle {\widehat {\sigma }}^{2}={1 \over n-m}\sum _{j=1}^{n}{\widehat {\varepsilon \,}}_{j}^{\,2}.}$

where m is the number of parameters in the model (2 in our example).

But if the i th case is suspected of being improbably large, then it would also not be normally distributed. Hence it is prudent to exclude the i th observation from the process of estimating the variance when one is considering whether the i th case may be an outlier, and instead use the externally studentized residual, which is

${\displaystyle {\widehat {\sigma }}_{(i)}^{2}={1 \over n-m-1}\sum _{\begin{smallmatrix}j=1\\j\neq i\end{smallmatrix}}^{n}{\widehat {\varepsilon \,}}_{j}^{\,2},}$

based on all the residuals except the suspect i th residual. Here is to emphasize that ${\displaystyle {\widehat {\varepsilon \,}}_{j}^{\,2}(j\neq i)}$ for suspect i are computed with i th case excluded.

If the estimate σ2 includes the i th case, then it is called the internally studentized residual, ${\displaystyle t_{i}}$ (also known as the standardized residual [1]). If the estimate ${\displaystyle {\widehat {\sigma }}_{(i)}^{2}}$ is used instead, excluding the i th case, then it is called the externally studentized, ${\displaystyle t_{i(i)}}$.

## Distribution

If the errors are independent and normally distributed with expected value 0 and variance σ2, then the probability distribution of the ith externally studentized residual ${\displaystyle t_{i(i)}}$ is a Student's t-distribution with n − m − 1 degrees of freedom, and can range from ${\displaystyle \scriptstyle -\infty }$ to ${\displaystyle \scriptstyle +\infty }$.

On the other hand, the internally studentized residuals are in the range ${\displaystyle \scriptstyle 0\,\pm \,{\sqrt {\nu }}}$, where ν = n − m is the number of residual degrees of freedom. If ti represents the internally studentized residual, and again assuming that the errors are independent identically distributed Gaussian variables, then:[2]

${\displaystyle t_{i}\sim {\sqrt {\nu }}{t \over {\sqrt {t^{2}+\nu -1}}}}$

where t is a random variable distributed as Student's t-distribution with ν − 1 degrees of freedom. In fact, this implies that ti /ν follows the beta distribution B(1/2,(ν − 1)/2). The distribution above is sometimes referred to as the tau distribution;[2] it was first derived by Thompson in 1935.[3]

When ν = 3, the internally studentized residuals are uniformly distributed between ${\displaystyle \scriptstyle -{\sqrt {3}}}$ and ${\displaystyle \scriptstyle +{\sqrt {3}}}$. If there is only one residual degree of freedom, the above formula for the distribution of internally studentized residuals doesn't apply. In this case, the ti are all either +1 or −1, with 50% chance for each.

The standard deviation of the distribution of internally studentized residuals is always 1, but this does not imply that the standard deviation of all the ti of a particular experiment is 1. For instance, the internally studentized residuals when fitting a straight line going through (0, 0) to the points (1, 4), (2, −1), (2, −1) are ${\displaystyle {\sqrt {2}},\ -{\sqrt {5}}/5,\ -{\sqrt {5}}/5}$, and the standard deviation of these is not 1.

Note that any pair of studentized residual ti and tj (where ${\displaystyle i\neq j}$), are NOT i.i.d. They have the same distribution, but are not independent due to constraints on the residuals having to sum to 0 and to have them be orthogonal to the design matrix.

## Software implementations

Many programs and statistics packages, such as R, Python, etc., include implementations of Studentized residual.

Language/Program Function Notes
R rstandard(model, ...) internally studentized. See [2]
R rstudent(model, ...) externally studentized. See [3]

## References

1. ^
2. ^ a b Allen J. Pope (1976), "The statistics of residuals and the detection of outliers", U.S. Dept. of Commerce, National Oceanic and Atmospheric Administration, National Ocean Survey, Geodetic Research and Development Laboratory, 136 pages, [1], eq.(6)
3. ^ Thompson, William R. (1935). "On a Criterion for the Rejection of Observations and the Distribution of the Ratio of Deviation to Sample Standard Deviation". The Annals of Mathematical Statistics. 6 (4): 214–219. doi:10.1214/aoms/1177732567.
This page was last edited on 15 July 2019, at 14:22
Basis of this page is in Wikipedia. Text is available under the CC BY-SA 3.0 Unported License. Non-text media are available under their specified licenses. Wikipedia® is a registered trademark of the Wikimedia Foundation, Inc. WIKI 2 is an independent company and has no affiliation with Wikimedia Foundation.