I have a left censored dependent variable where many of the observations have a value of zero. The data is clustered (multiple measurements over time for each person). I initially decided to use a Tobit model with random effects to model my outcome variable. But I am currently thinking of using a multilevel linear model with random effects instead (for various reasons that I don't want to discuss in this post).

Normally, a multilevel linear model is not suitable for modeling censored outcome variables. However, I thought of a trick that might fix the issue, but I wanted to seek the community's feedback to check whether my rationale is correct. Before I discuss the proposed solution, I'll give some background about the study's setting. The main variable of interest represents a specific component of service time in a large service organization. This component of service time is optional (i.e., workers are supposed to perform it only when needed). This is why a lot of the observations had a zero value for this component. Let's call the component $x_1$. There are two more components of service time that are performed by workers on regular basis: $x_2$ and $x_3$. These components are not censored.

My idea is to construct a new variable $z_1$ that represents total service time (i.e., $z_1 = x_1 + x_2 + x_3$). This new variable is not censored, so I can use multilevel linear models to regress $z_1$ on the independent variables of interest. But I am interested in the effects of those IVs on $x_1$ rather than $z_1$. So I'll control for $x_2$ and $x_3$ in the model to examine the effects of the IVs on $x_1$. In other words, my model would be: $z_1 = \beta_0 + \beta_1 x_2 + \beta_2 x_3 + W \Gamma + \epsilon$ where $W$ is a vector of IVs of interest. My rationale is that if I control for $x_2$ and $x_3$, then the only remaining service component in $z_1$ is $x_1$, which means that $\Gamma$ represents the effects of the IVs on $x_1$ rather than $z_1$. At the same time, I can use the multilevel linear model instead of a Tobit model since $z_1$ is not censored.

Is this approach correct? Did I miss anything?

Edit: Made a correction to the subscripts of the model parameters.

  • $\begingroup$ I know you said that you did not want to get into the matter, but unless there is a very strong reason not to use a generalized linear model, I do not view the proposed approach as valid. If this were to be put in a manuscript for review or in a dissertation/thesis proposal, you would be called on it. In fact, given the data you describe, it sounds like a Poission model for the x_1 outcome might be more appropriate rather than tobit. I say this because the outcome sounds like a count of the number of times something happened. With lots of 0s, then a negative binomial model might be better. $\endgroup$ – Erik Ruzek Jan 11 at 19:23
  • $\begingroup$ In short, if I use a Tobit model I would need to calculate the marginal effects for the betas in order to be able to interpret the data properly. This can be easily achieved using major statistical software packages. However, from what I understood, I cannot use those marginal effect procedures with interaction or quadratic terms. And my main hypotheses depend on those non-linear terms. Thus, while the Tobit model results are helpful in showing the significance and direction of the effects, I cannot interpret the results properly. [1 of 2] $\endgroup$ – GXR200 Jan 11 at 19:42
  • $\begingroup$ This is why I thought of using a linear model which doesn't suffer from this issue, given of course that the method I described is not flawed, which is why I asked the question (otherwise I'll omit from using the method). As for the Poisson model suggestion, my outcome is actually a time measurement (time it takes to complete service), so it is a continuous variable. In the case of zero observed value, it simply means that the worker didn't engage in the optional service stage. [2 of 2] $\endgroup$ – GXR200 Jan 11 at 19:46
  • $\begingroup$ BTW, why do view the proposed approach as not valid? I am interested in hearing you thoughts. $\endgroup$ – GXR200 Jan 11 at 20:18

You could consider using a two-part mixed-effects model for semi-continuous data. This combines a mixed-effects logistic regression for the dichotomous outcome $\mbox{I}(x_1 = 0)$ with a linear mixed model for the logarithm positive responses $\log(x_1)$ for $x_1 > 0$.

This model can be fitted, for example, in the GLMMadaptive package in R. For an example, check here.


Thank you for the additional explanation in the comments about your outcome. They help in thinking about a way to move forward. I will first focus on that issue before sharing why I find your solution to move to a linear model problematic.

  1. With regards to $x_1$ and using a tobit model, it is critical to note that tobit regression works under the assumption that the underlying variable has a continuous scale, theoretically ranging from negative to positive. However, you do not observe some piece of that range (called truncation). With your time variable, a negative value does not make sense, so tobit regression doesn't apply. See this thread. You need a different generalized linear modeling approach.

  2. For modeling $x_i$, you have a couple of possibilities given the presence of all those 0s. In both cases, you would need to move to the generalized linear modeling framework to accomodate your strictly positive dependent variable. One possibility would be a gamma GLM model. The other possibility to consider is the class of zero-inflated models combined with a gamma GLM. These models have a separate model for the 0 vs. > 0 case and then a model for the values > 0.

  3. Finally, your proposed solution for adding $x_1$, $x_2$, and $x_3$ together as a new dependent variable and then running a linear model w/ covariates $x_2$ and $x_3$ included is problematic from my perspective. You are proposing to include as independent variables in your regression the very same variables that go into forming your dependent variable. An analogous situation would be a case where the outcome was the sum of 3 conceptually related items from a Likert-scaled survey and then in a regression model predicting that sum score I include 2 of those three items as predictors. We would expect the two items to be highly correlated with the sum score, so much so that we would have trouble finding other covariates powerful enough to explain what is leftover. And what exactly is it that is leftover?

By combining the items you create a new construct, and you need to convince readers that the construct is meaningful. You are trying to use a mathematical trick to justify using a linear model and I'm not convinced that trick works as you think it does. And, at the same time, it poses all sorts of interpretation problems. I would be interested to know if others have different thoughts about this issue, but from my perspective it is problematic.

Edit: In general, my suggestion is to build a model for the variable you are actually interested in understanding rather than trying to manipulate that variable in such a way that you sort of retain its original properties. This also means that if the variable were amenable to a known transformation that maintained its properties, you could do that, allowing you to utilize a linear model for the transformed version of the outcome. For example, the log transformation for the >0 values mentioned in @Dimitris Rizopoulos' response to your post.

  • $\begingroup$ Thank you for your reply and suggestions, this is helpful. With regards to point 1 above, you suggested that a Tobit model might not be suitable for my analysis given that the Tobit model assumes that censoring might involve negative values of the outcome variable and that a negative value for my outcome variable doesn't make sense. You make a good argument, but if you look at examples of Tobit model usage discussed in three major econometrics books (Hill, Griffiths, and Lim 2017; Wooldridge 2016; Greene 2018), you'll notice that they all discussed models of married women's labor force ... $\endgroup$ – GXR200 Jan 12 at 7:15
  • $\begingroup$ participation where the outcome was the number of hours worked by the married women. The way I see it is that those models are similar to my case where the outcome variable is censored at zero (a large fraction of women in the study sample had zero hours of work), it is continuous, and it cannot assume a negative value (it doesn't make sense for someone to work negative hours). I've seen other examples of research that uses Tobit models with similar non-negative outcomes. $\endgroup$ – GXR200 Jan 12 at 7:15
  • $\begingroup$ With regards to your second point, this is very helpful. I'll check those alternatives, so thanks for sharing the resources. Finally, with regards to the third point, I agree that my approach is not conventional. To be honest, I don't recall we ever discussed such a method in any of my stat/econometrics classes. While I am worried about the validity of the method, I would love to understand why it doesn't work (if that was the case indeed). I understand your argument that there are high correlations between $x_2$/$x_3$ and the outcome variable $z_1$... $\endgroup$ – GXR200 Jan 12 at 7:30
  • $\begingroup$ Still, when I ran the aforementioned model, the parameter estimates for the other independent variables were all statistically significant and similar (but not identical) to the results of another model that uses $x_1$ as the outcome variable. So my intuition tells me that the model $𝑧_1 = \beta_0+\beta_1 𝑥_2 + \beta_2 𝑥_3+𝑊\Gamma +\epsilon$ is connected to a model of $x_1 = \mu_0 + 𝑊B + u$. The main rationale is that if you subtract $x_2$ and $x_3$ from both sides of the $z_1$, then you are left with $x_1$ on the left hand side. But still I am not sure if this is statistically correct. $\endgroup$ – GXR200 Jan 12 at 7:31
  • $\begingroup$ With regards to tobit usage, your field may accept these types of situations, in which case I would use that rather than the adding together three variables and using it as your dependent variable. Why do you think you cannot use marginal effect estimation after running tobit models. It looks possible in both R's ggeffects - strengejacke.github.io/ggeffects and in Stata's margins - timberlake.co.uk/news/…. $\endgroup$ – Erik Ruzek Jan 12 at 15:01

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.