# Tobit model

The Tobit model is a statistical model proposed by James Tobin (1958)[1] to describe the relationship between a non-negative dependent variable ${\displaystyle y_{i}}$ and an independent variable (or vector) ${\displaystyle x_{i}}$. The term Tobit was derived from Tobin's name by truncating and adding -it by analogy with the probit model.[2] The Tobit model is distinct from the truncated regression model, which is in general different and requires a different estimator.[3]

The model supposes that there is a latent (i.e. unobservable) variable ${\displaystyle y_{i}^{*}}$. This variable linearly depends on ${\displaystyle x_{i}}$ via a parameter (vector) ${\displaystyle \beta }$ which determines the relationship between the independent variable (or vector) ${\displaystyle x_{i}}$ and the latent variable ${\displaystyle y_{i}^{*}}$ (just as in a linear model). In addition, there is a normally distributed error term ${\displaystyle u_{i}}$ to capture random influences on this relationship. The observable variable ${\displaystyle y_{i}}$ is defined as the ramp function: equal to the latent variable whenever the latent variable is above zero, and zero otherwise.

${\displaystyle y_{i}={\begin{cases}y_{i}^{*}&{\text{if }}y_{i}^{*}>0,\\0&{\text{if }}y_{i}^{*}\leq 0,\end{cases}}}$

where ${\displaystyle y_{i}^{*}}$ is a latent variable:

${\displaystyle y_{i}^{*}=\beta x_{i}+u_{i},\qquad u_{i}\sim N(0,\sigma ^{2}).\,}$

## Etymology

When asked why it was called the "Tobit" model, instead of Tobin, James Tobin explained that this term was introduced by Arthur Goldberger, either as a contraction of "Tobin probit", or as a reference to the novel The Caine Mutiny, a novel by Tobin's friend Herman Wouk, in which Tobin makes a cameo as "Mr Tobit". Tobin reports having actually asked Goldberger which it was, and the man refused to say.[4]

## Consistency

If the relationship parameter ${\displaystyle \beta }$ is estimated by regressing the observed ${\displaystyle y_{i}}$ on ${\displaystyle x_{i}}$, the resulting ordinary least squares regression estimator is inconsistent. It will yield a downwards-biased estimate of the slope coefficient and an upward-biased estimate of the intercept. Takeshi Amemiya (1973) has proven that the maximum likelihood estimator suggested by Tobin for this model is consistent.[5]

## Interpretation

The ${\displaystyle \beta }$ coefficient should not be interpreted as the effect of ${\displaystyle x_{i}}$ on ${\displaystyle y_{i}}$, as one would with a linear regression model; this is a common error. Instead, it should be interpreted as the combination of (1) the change in ${\displaystyle y_{i}}$ of those above the limit, weighted by the probability of being above the limit; and (2) the change in the probability of being above the limit, weighted by the expected value of ${\displaystyle y_{i}}$ if above. [6]

## Variations of the Tobit model

Variations of the Tobit model can be produced by changing where and when censoring occurs. Amemiya (1985, p. 384) 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.[7]

### Type I

The Tobit model is a special case of a censored regression model, because the latent variable ${\displaystyle y_{i}^{*}}$ cannot always be observed while the independent variable ${\displaystyle x_{i}}$ is observable. A common variation of the Tobit model is censoring at a value ${\displaystyle y_{L}}$ different from zero:

${\displaystyle y_{i}={\begin{cases}y_{i}^{*}&{\text{if }}y_{i}^{*}>y_{L},\\y_{L}&{\text{if }}y_{i}^{*}\leq y_{L}.\end{cases}}}$

Another example is censoring of values above ${\displaystyle y_{U}}$.

${\displaystyle y_{i}={\begin{cases}y_{i}^{*}&{\text{if }}y_{i}^{*}

Yet another model results when ${\displaystyle y_{i}}$ is censored from above and below at the same time.

${\displaystyle y_{i}={\begin{cases}y_{i}^{*}&{\text{if }}y_{L}

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

Type II Tobit models introduce a second latent variable.

${\displaystyle y_{2i}={\begin{cases}y_{2i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$

Heckman (1987) falls into the Type II Tobit. In Type I Tobit, the latent variable absorb both the process of participation and 'outcome' of interest. Type II Tobit allows the process of participation/selection and the process of 'outcome' to be independent, conditional on x.

### Type III

Type III introduces a second observed dependent variable.

${\displaystyle y_{1i}={\begin{cases}y_{1i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$
${\displaystyle y_{2i}={\begin{cases}y_{2i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$

The Heckman model falls into this type.

### Type IV

Type IV introduces a third observed dependent variable and a third latent variable.

${\displaystyle y_{1i}={\begin{cases}y_{1i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$
${\displaystyle y_{2i}={\begin{cases}y_{2i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$
${\displaystyle y_{3i}={\begin{cases}y_{3i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$

### Type V

Similar to Type II, in Type V only the sign of ${\displaystyle y_{1i}^{*}}$ is observed.

${\displaystyle y_{2i}={\begin{cases}y_{2i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$
${\displaystyle y_{3i}={\begin{cases}y_{3i}^{*}&{\text{if }}y_{1i}^{*}>0,\\0&{\text{if }}y_{1i}^{*}\leq 0.\end{cases}}}$

## The likelihood function

Below are the likelihood and log likelihood functions for a type I Tobit. This is a Tobit that is censored from below at ${\displaystyle y_{L}}$ when the latent variable ${\displaystyle y_{j}^{*}\leq y_{L}}$. In writing out the likelihood function, we first define an indicator function ${\displaystyle I(y_{j})}$ where:

${\displaystyle I(y_{j})={\begin{cases}0&{\text{if }}y_{j}\leq y_{L},\\1&{\text{if }}y_{j}>y_{L}.\end{cases}}}$

Next, let ${\displaystyle \Phi }$ be the standard normal cumulative distribution function and ${\displaystyle \varphi }$ to be the standard normal probability density function. For a data set with N observations the likelihood function for a type I Tobit is

${\displaystyle {\mathcal {L}}(\beta ,\sigma )=\prod _{j=1}^{N}\left({\frac {1}{\sigma }}\varphi \left({\frac {y_{j}-X_{j}\beta }{\sigma }}\right)\right)^{I\left(y_{j}\right)}\left(1-\Phi \left({\frac {X_{j}\beta -y_{L}}{\sigma }}\right)\right)^{1-I\left(y_{j}\right)}}$

and the log likelihood is given by

${\displaystyle \log {\mathcal {L}}(\beta ,\sigma )=\sum _{j=1}^{n}I(y_{j})\log \left({\frac {1}{\sigma }}\varphi \left({\frac {y_{j}-X_{j}\beta }{\sigma }}\right)\right)+(1-I(y_{j}))\log \left(1-\Phi \left({\frac {X_{j}\beta -y_{L}}{\sigma }}\right)\right)}$

Note that this is different from the likelihood function of the truncated regression model.[8]

## Non-parametric version

If the underlying latent variable ${\displaystyle y_{i}^{*}}$ is not normally distributed, one must use quantiles instead of moments to analyze the observable variable ${\displaystyle y_{i}}$. Powell's CLAD estimator offers a possible way to achieve this.[9]

## Dynamic unobserved effects Tobit model

In a panel data Tobit model,[10][11] if the outcome ${\displaystyle Y_{i,t}}$ partially depends on the previous outcome history ${\displaystyle Y_{i,0},\ldots ,Y_{t-1}}$ this Tobit model is called "dynamic". For instance, taking a person who finds a job with a high salary this year, it will be easier for her to find a job with a high salary next year because the fact that she has a high-wage job this year will be a very positive signal for the potential employers. The essence of this type of dynamic effect is the state dependence of the outcome. The "unobservable effects" here refers to the factor which partially determines the outcome of individual but cannot be observed in the data. For instance, the ability of a person is very important in job-hunting, but it is not observable for researchers. A typical dynamic unobserved effects Tobit model can be represented as

${\displaystyle Y_{i,t}=Y_{i,t}^{1}[Y_{i,t}>0];}$
${\displaystyle Y_{i,t}=z_{i,t}\delta +\rho y_{i,t-1}+c_{i}+u_{i,t};}$
${\displaystyle c_{i}\mid y_{i,0},\ldots ,y_{i,t-1}\sim F(y_{i,0}x_{i});}$
${\displaystyle u_{i,t}\mid z_{i,t},y_{i,0},\ldots ,y_{i,t-1}\backsim _{N}(0,1).}$

In this specific model, ${\displaystyle \rho y_{i,t-1}}$ is the dynamic effect part and ${\displaystyle c_{i}}$ is the unobserved effect part whose distribution is determined by the initial outcome of individual i and some exogenous features of individual i.

Based on this setup, the likelihood function conditional on ${\displaystyle \{y_{i,0}\}_{i-1}^{N}}$ can be given as

${\displaystyle \prod _{i=1}^{N}\int f_{\theta }(c_{i}\mid y_{i,0},x_{i})\left[\prod _{t=1}^{T}{\Bigl (}1[y_{i,t}=0][1-\Phi (z_{i,t}\delta +\rho y_{i,t-1}>0]{\frac {\varphi (z_{i,t}\delta +\rho y_{i,t-1}+c_{i})}{\Phi (z_{i,t}\delta +\rho y_{i,t-1}+c_{i})}}{\biggr )}\right]\,dc_{i}}$

For the initial values ${\displaystyle \{y_{i,0}\}_{i-1}^{N}}$ ,there are two different ways to treat them in the construction of the likelihood function: treating them as constant, or imposing a distribution on them and calculate out the unconditional likelihood function. But whichever way is chosen to treat the initial values in the likelihood function, we cannot get rid of the integration inside the likelihood function when estimating the model by Maximum Likelihood Estimation (MLE). Expectation Maximum (EM) algorithm is usually a good solution for this computation issue.[12] Based on the consistent point estimates from MLE, Average Partial Effect (APE)[13] can be calculated correspondingly.[14]

## Applications

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)[15] analyse a sample of 115 municipalities (42 of which received a grant). Dubois and Fattore (2011)[16] use a Tobit model to investigate the role of various factors in European Union fund receipt by applying Polish sub-national governments. 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.[17]

## References

1. ^ Tobin, James (1958). "Estimation of relationships for limited dependent variables". Econometrica. 26 (1): 24–36. doi:10.2307/1907382. JSTOR 1907382.
2. ^
3. ^ Park, B.U., L. Simar, and V. Zelenyuk (2008). "Local likelihood estimation of truncated regression and its partial derivatives: Theory and application," Journal of Econometrics 146(1), pages 185-198.
4. ^ The ET Interview: Professor James Tobin
5. ^ Amemiya, Takeshi (1973). "Regression analysis when the dependent variable is truncated normal". Econometrica. 41 (6): 997–1016. doi:10.2307/1914031. JSTOR 1914031.
6. ^ McDonald, John F.; Moffit, Robert A. (1980). "The Uses of Tobit Analysis". The Review of Economics and Statistics. The MIT Press. 62 (2): 318–321. doi:10.2307/1924766. JSTOR 1924766.
7. ^ Schnedler, Wendelin (2005). "Likelihood estimation for censored random vectors". Econometric Reviews. 24 (2): 195–217. doi:10.1081/ETC-200067925.
8. ^ Park, B.U., L. Simar, and V. Zelenyuk (2008). "Local likelihood estimation of truncated regression and its partial derivatives: Theory and application," Journal of Econometrics 146(1), pages 185-198.
9. ^ Powell, James L (1 July 1984). "Least absolute deviations estimation for the censored regression model". Journal of Econometrics. 25 (3): 303–325. doi:10.1016/0304-4076(84)90004-6.
10. ^ Greene, W. H. (2003), Econometric Analysis , Prentice Hall , Upper Saddle River, NJ.
11. ^ The model framework comes from Wooldridge, J. (2002): Econometric Analysis of Cross Section and Panel Data, MIT Press, Cambridge, Mass, pp 542. But the author revises the model more general here.
12. ^ For more details, refer to: Cappé, O., Moulines, E. and Ryden, T. (2005): “Inference in Hidden Markov Models”, Springer-Verlag New York, Chapter 2
13. ^ Wooldridge, J. (2002): Econometric Analysis of Cross Section and Panel Data, MIT Press, Cambridge, Mass, pp 22.
14. ^ For more details, refer to: Takeshi Amemiya (1984),” Tobit models: A survey”, Journal of Econometrics, 24, pp 3–6
15. ^ Dahlberg, Matz; Johansson, Eva (2002-03-01). "On the Vote-Purchasing Behavior of Incumbent Governments". American Political Science Review. null (01): 27–40. doi:10.1017/S0003055402004215. ISSN 1537-5943.
16. ^ Dubois, Hans F. W.; Fattore, Giovanni (2011-07-01). "Public Fund Assignment through Project Evaluation". Regional & Federal Studies. 21 (3): 355–374. doi:10.1080/13597566.2011.578827. ISSN 1359-7566.
17. ^ Baltas, George (2001). "Utility-consistent Brand Demand Systems with Endogenous Category Consumption: Principles and Marketing Applications". Decision Sciences. 32 (3): 399–422. doi:10.1111/j.1540-5915.2001.tb00965.x. ISSN 0011-7315.