In statistics, a tobit model is any of a class of regression models in which the observed range of the dependent variable is censored in some way.[1] The term was coined by Arthur Goldberger in reference to James Tobin,[2] who developed the model in 1958 to mitigate the problem of zero-inflated data for observations of household expenditure on durable goods.[3] Because Tobin's method can be easily extended to handle truncated and other non-randomly selected samples, some authors adopt a broader definition of the tobit model that includes these cases.[4]
Tobin's idea was to modify the likelihood function so that it reflects the unequal sampling probability for each observation depending on whether the latent dependent variable fell above or below the determined threshold.[5] For a sample that, as in Tobin's original case, was censored from below at zero, the sampling probability for each non-limit observation is simply the height of the appropriate density function. For any limit observation, it is the cumulative distribution, i.e. the integral below zero of the appropriate density function. The tobit likelihood function is thus a mixture of densities and cumulative distribution functions.[6]
Below are the likelihood and log likelihood functions for a type I tobit. This is a tobit that is censored from below at
yL
* | |
y | |
j |
\leqyL
I
I(y)=\begin{cases} 0&ify\leqyL,\\ 1&ify>yL. \end{cases}
Next, let
\Phi
\varphi
l{L}(\beta,\sigma)=\prod
N | ||
\left( | ||
j=1 |
1 | |
\sigma |
\varphi\left(
yj-Xj\beta | |
\sigma |
I(yj) | |
\right)\right) |
\left(1-\Phi\left(
Xj\beta-yL | |
\sigma |
1-I(yj) | |
\right)\right) |
and the log likelihood is given by
\begin{align} logl{L}(\beta,\sigma)&=
n | |
\sum | |
j=1 |
I(yj)log\left(
1 | |
\sigma |
\varphi\left(
yj-Xj\beta | |
\sigma |
\right)\right)+(1-I(yj))log\left(1-\Phi\left(
Xj\beta-yL | |
\sigma |
\right)\right)\\ &=
\sum | |
yj>yL |
log\left(
1 | |
\sigma |
\varphi\left(
yj-Xj\beta | |
\sigma |
\right)\right)+
\sum | |
yj=yL |
log\left(\Phi\left(
yL-Xj\beta | |
\sigma |
\right)\right) \end{align}
The log-likelihood as stated above is not globally concave, which complicates the maximum likelihood estimation. Olsen suggested the simple reparametrization
\beta=\delta/\gamma
\sigma2=\gamma-2
logl{L}(\delta,\gamma)=
\sum | |
yj>yL |
\left\{log\gamma+log\left[\varphi\left(\gammayj-Xj\delta\right)\right]\right\}+
\sum | |
yj=yL |
log\left[\Phi\left(\gammayL-Xj\delta\right)\right]
which is globally concave in terms of the transformed parameters.[7]
For the truncated (tobit II) model, Orme showed that while the log-likelihood is not globally concave, it is concave at any stationary point under the above transformation.[8] [9]
If the relationship parameter
\beta
yi
xi
The
\beta
xi
yi
yi
yi
Variations of the tobit model can be produced by changing where and when censoring occurs. classifies these variations into five categories (tobit type I – tobit type V), where tobit type I stands for the first model described above. Schnedler (2005) provides a general formula to obtain consistent likelihood estimators for these and other variations of the tobit model.[12]
The tobit model is a special case of a censored regression model, because the latent variable
* | |
y | |
i |
xi
yL
yi=\begin{cases}
* | |
y | |
i |
&if
* | |
y | |
i |
>yL,\\ yL&if
* | |
y | |
i |
\leqyL. \end{cases}
Another example is censoring of values above
yU
yi=\begin{cases}
* | |
y | |
i |
&if
* | |
y | |
i |
<yU,\\ yU&if
* | |
y | |
i |
\geqyU. \end{cases}
Yet another model results when
yi
yi=\begin{cases}
* | |
y | |
i |
&ifyL<y
* | |
i |
<yU,\\ yL&if
* | |
y | |
i |
\leqyL,\\ yU&if
* | |
y | |
i |
\geqyU. \end{cases}
The rest of the models will be presented as being bounded from below at 0, though this can be generalized as done for Type I.
Type II tobit models introduce a second latent variable.[13]
y2i=\begin{cases}
* | |
y | |
2i |
&if
* | |
y | |
1i |
>0,\\ 0&if
* | |
y | |
1i |
\leq0. \end{cases}
In Type I tobit, the latent variable absorbs both the process of participation and the outcome of interest. Type II tobit allows the process of participation (selection) and the outcome of interest to be independent, conditional on observable data.
The Heckman selection model falls into the Type II tobit,[14] which is sometimes called Heckit after James Heckman.[15]
Type III introduces a second observed dependent variable.
y1i=\begin{cases}
* | |
y | |
1i |
&if
* | |
y | |
1i |
>0,\\ 0&if
* | |
y | |
1i |
\leq0. \end{cases}
y2i=\begin{cases}
* | |
y | |
2i |
&if
* | |
y | |
1i |
>0,\\ 0&if
* | |
y | |
1i |
\leq0. \end{cases}
Type IV introduces a third observed dependent variable and a third latent variable.
y1i=\begin{cases}
* | |
y | |
1i |
&if
* | |
y | |
1i |
>0,\\ 0&if
* | |
y | |
1i |
\leq0. \end{cases}
y2i=\begin{cases}
* | |
y | |
2i |
&if
* | |
y | |
1i |
>0,\\ 0&if
* | |
y | |
1i |
\leq0. \end{cases}
y3i=\begin{cases}
* | |
y | |
3i |
&if
* | |
y | |
1i |
\leq0,\\ 0&if
* | |
y | |
1i |
<0. \end{cases}
Similar to Type II, in Type V only the sign of
* | |
y | |
1i |
y2i=\begin{cases}
* | |
y | |
2i |
&if
* | |
y | |
1i |
>0,\\ 0&if
* | |
y | |
1i |
\leq0. \end{cases}
y3i=\begin{cases}
* | |
y | |
3i |
&if
* | |
y | |
1i |
\leq0,\\ 0&if
* | |
y | |
1i |
>0. \end{cases}
If the underlying latent variable
* | |
y | |
i |
yi
Tobit models have, for example, been applied to estimate factors that impact grant receipt, including financial transfers distributed to sub-national governments who may apply for these grants. In these cases, grant recipients cannot receive negative amounts, and the data is thus left-censored. For instance, Dahlberg and Johansson (2002) analyse a sample of 115 municipalities (42 of which received a grant).[17] Dubois and Fattore (2011) use a tobit model to investigate the role of various factors in European Union fund receipt by applying Polish sub-national governments.[18] The data may however be left-censored at a point higher than zero, with the risk of mis-specification. Both studies apply Probit and other models to check for robustness. Tobit models have also been applied in demand analysis to accommodate observations with zero expenditures on some goods. In a related application of tobit models, a system of nonlinear tobit regressions models has been used to jointly estimate a brand demand system with homoscedastic, heteroscedastic and generalized heteroscedastic variants.[19]