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. The term was coined by Arthur Goldberger in reference to James Tobin, who developed the model in 1958 to mitigate the problem of zero-inflated data for observations of household expenditure on durable goods. 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. 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. 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.
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 y L {\displaystyle y_{L}} when the latent variable y j ∗ ≤ y L {\displaystyle y_{j}^{*}\leq y_{L}} . In writing out the likelihood function, we first define an indicator function I {\displaystyle I} :
I ( y ) = { 0 if y ≤ y L , 1 if y > y L . {\displaystyle I(y)={\begin{cases}0&{\text{if }}y\leq y_{L},\\1&{\text{if }}y>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
L ( β , σ ) = ∏ j = 1 N ( 1 σ φ ( y j − X j β σ ) ) I ( y j ) ( 1 − Φ ( X j β − y L σ ) ) 1 − I ( y j ) {\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(y_{j})}\left(1-\Phi \left({\frac {X_{j}\beta -y_{L}}{\sigma }}\right)\right)^{1-I(y_{j})}}
and the log likelihood is given by
… excerpt ends here. Continue reading the full article.
