# Bayesian Approach to Zero-Inflated Bivariate Ordered Probit Regression Model, with an Application to Tobacco Use

Authors: Shiferaw Gurmu; Getachew A. Dagne
Journal: Journal of Probability and Statistics
 (2012)

---

## Abstract

This paper presents a Bayesian analysis of bivariate ordered probit regression model with excess of zeros. Specifically, in the context of joint modeling of two ordered outcomes, we develop zero-inflated bivariate ordered probit model and carry out estimation using Markov Chain Monte Carlo techniques. Using household tobacco survey data with substantial proportion of zeros, we analyze the socioeconomic determinants of individual problem of smoking and chewing tobacco. In our illustration, we find strong evidence that accounting for excess zeros provides good fit to the data. The example shows that the use of a model that ignores zero-inflation masks differential effects
of covariates on nonusers and users.


---

## Body

## 1. Introduction

This paper is concerned with joint modeling of two ordered data outcomes allowing for excess zeros. Economic, biological, and social science studies often yield data on two ordered categorical variables that are jointly dependent. Examples include the relationship between desired and excess fertility [1, 2], helmet use and motorcycle injuries [3], ownership of dogs and televisions [4], severity of diabetic retinopathy of the left and right eyes [5], and self-assessed health status and wealth [6]. The underlying response variables could be measured on an ordinal scale. It is also common in the literature to generate a categorical or grouped variable from an underlying quantitative variable and then use ordinal response regression model (e.g., [4, 5, 7]). The ensuing model is usually analyzed using the bivariate ordered probit model.Many ordered discrete data sets are characterized by excess of zeros, both in terms of the proportion of nonusers and relative to the basic ordered probit or logit model. The zeros may be attributed to either corner solution to consumer optimization problem or errors in recording. In the case of individual smoking behavior, for example, the zeros may be recorded for individuals who never smoke cigarettes or for those who either used tobacco in the past or are potential smokers. In the context of individual patents applied for by scientists during a period of five years, zero patents may be recorded for scientists who either never made patent applications or for those who do but not during the reporting period [8]. Ignoring the two types of zeros for nonusers or nonparticipants leads to model misspecification.The univariate as well as bivariate zero-inflated count data models are well established in the literature for example, Lambert [9], Gurmu and Trivedi [10], Mullahy [11], and Gurmu and Elder [12]. The recent literature presents a Bayesian treatment of zero-inflated Poisson models in both cross-sectional and panel data settings (see [13, 14], and references there in). By contrast, little attention has been given to the problem of excess zeros in the ordered discrete choice models. Recently, an important paper by Harris and Zhao [15] developed a zero-inflated univariate ordered probit model. However, the problem of excess zeros in ordered probit models has not been analyzed in the Bayesian framework. Despite recent applications and advances in estimation of bivariate ordered probit models [1–6], we know of no studies that model excess zeros in bivariate ordered probit models.This paper presents a Bayesian analysis of bivariate ordered probit model with excess of zeros. Specifically, we develop a zero-inflated ordered probit model and carry out the analysis using the Bayesian approach. The Bayesian analysis is carried out using Markov Chain Monte Carlo (MCMC) techniques to approximate the posterior distribution of the parameters. Bayesian analysis of the univariate  zero-inflated ordered probit will be treated as a special case of the zero-inflated bivariate order probit model. The proposed models are illustrated by analyzing the socioeconomic determinants of individual choice problem of bivariate ordered outcomes on smoking  and chewing tobacco. We use household tobacco prevalence survey data from Bangladesh. The observed proportion of zeros (those identifying themselves as nonusers of tobacco) is about 76% for smoking and 87% for chewing tobacco.The proposed approach is useful  for the analysis of ordinal data with natural zeros. The empirical analysis clearly shows the importance of accounting for excess zeros in ordinal qualitative response models. Accounting for excess zeros provides good fit to the data. In terms of both the signs and magnitudes of marginal effects, various covariates have differential impacts on the probabilities associated with the two types of zeros, nonparticipants and zero-consumption. The usual analysis that ignores excess of zeros masks these differential effects, by just focusing on observed zeros. The empirical results also show the importance of taking into account the uncertainty in the parameter estimates. Another advantage of the Bayesian approach to modeling excess zeros is the flexibility, particularly computational, of generalizing to multivariate ordered response models.The rest of the paper is organized as follows. Section2 describes the proposed zero-inflated bivariate probit model. Section 3 presents the MCMC algorithm and model selection procedure for the model. An illustrative application using household tobacco consumption data is given in Section 4. Section 5 concludes the paper.


## 2. Zero-Inflated Bivariate Ordered Probit Model

### 2.1. The Basic Model

We consider the basic Bayesian approach to a bivariate latent variable regression model with excess of zeros. To develop notation, letỹ1i* and ỹ2i*  denote the bivariate latent variables. We consider two observed ordered response variables ỹ1i and ỹ2i taking on values 0,1,…,Jr, for r=1,2. Define two sets of cut-off parameters αr=(αr2,αr3,…,αrJr), r=1,2, where the restrictions αr0=-∞, αrJr+1=∞, and αr1=0 have been imposed. We assume that (ỹ1i*,ỹ2i*)′≡ỹi* follows a bivariate regression model(2.1)ỹri*=xri′βr+ɛri,r=1,2,
where xri is a Kr-variate of regressors for the ith individual (i=1,…,N) and  ɛri are the error terms. For subsequent analysis, let β=(β1′,β2′)′, ϵi=(ϵ1i,ϵ2i)′, and (2.2)Xi=(x1i′0′0′x2i′).
Analogous to the univariate case, the observed bivariate-dependent variables are defined as(2.3)ỹri={0ifỹri*≤0,1if0<ỹri*≤αr2,jifαrj<ỹri*≤αrj+1,j=2,3,…,Jr-1,Jrifỹi*≤αrJr,
where r=1,2. Let ỹi=(ỹ1i,ỹ2i)′.We introduce inflation at the point(ỹ1i=0,ỹ2i=0), called the zero-zero state. As in the univariate case, define the participation model:(2.4)si*=zi′γ+μi,si=I(si*>0).
In the context of the zero-inflation model, the observed response random vector yi=(y1i,y2i)′ takes the form(2.5)yi=siỹi.
We observe yi=0 when either the individual is a non-participant (si=0) or the individual is a zero-consumption participant (si=1 and ỹi=0).  Likewise,  we observe positive outcome (consumption) when the individual is a positive consumption participant for at least one good (si=1 and ỹi≠0).LetΦ(a) and ϕ(a) denote the respective cumulative distribution and probability density functions of standardized normal evaluated at a. Assuming normality and that μi is uncorrelated with (ɛ1i,ɛ2i), but corr(ɛ1i,ɛ2i)=ρ12≠0, and each component with unit variance, the zero-inflated bivariate ordered probit (ZIBOP) distribution is(2.6)fb(yi*,yi,si*,si∣Xi,zi,Ψ)={Pr⁡(si=0)+(1-Pr⁡(si=0))Pr⁡(ỹ1i=0,ỹ2i=0),for(ỹ1i,ỹ2i)=(0,0)(1-Pr⁡(si=0))Pr⁡(ỹ1i=j,ỹ2i=l),for(ỹ1i,ỹ2i)≠(0,0),
where j=0,1,…,J1,l=0,1,…,J2,Pr(si=0)=Φ(-zi′γ),Pr(si=1)=Φ(-zi′γ). Further, for (ỹ1i,ỹ2i)=(0,0) in (2.6), we have αr0=-∞,αr1=0 for r=1,2 so that(2.7)Pr(ỹ1i=0,ỹ2i=0)=Φ2(-x1i′β1,-x2i′β2,ρ12),
where Φ2(·) is the cdf for the standardized bivariate normal. Likewise, Pr(ỹ1i=j,ỹ2i=l) in (2.6) are given by(2.8)Pr(ỹ1i=j,ỹ2i=l)=Φ2(α1j+1-x1i′β1,α2l+1-x2i′β2;ρ12)-Φ2(α1j-x1i′β1,α2l-x2i′β2,ρ12)forj=1,…,J1-1;l=1,…,J2-1;Pr(ỹ1i=J1,ỹ2i=J2)=1-Φ2(α1J1-x1i′β1,α2J2-x2i′β2,ρ12).
The ensuing likelihood contribution for N-independent observations is(2.9)Lb(y*,y,s*,s∣X,z,Ψb)=∏i=1N∏(j,l)=(0,0)[Pr(si=0)+(1-Pr(si=0))Pr(ỹ1i=0,ỹ2i=0)]dijl×∏i=1N∏(j,l)≠(0,0)[(1-Pr(si=0))Pr(ỹ1i=j,ỹ2i=l)]dijl,
where dijl=1 if ỹ1i=j and ỹ2i=l, and dijl=0 otherwise. Here, the vector Ψb consists of β,γ,α1,α2, and the parameters associated with the trivariate distribution of (ϵ,μ).Regarding identification of the parameters in the model defined by (2.1) through (2.5) with normality assumption, we note that the mean parameter (joint choice probability associate with the observed response vector yi) depends nonlinearly on the probability of zero inflation (Φ(-zi′γ)) and choice probability (Pr(ỹ1i=j,ỹ2i=l)) coming from the BOP submodel. Since the likelihood function for ZIBOP depends separately on the two regression components, the parameters of ZIBOP model with covariates are identified as long as the model is estimated by full maximum likelihood method. The same or different sets of covariates can affect the two components via zi and xri. When using quasi-likelihood estimation or generalized estimating equations methods rather than full ML, the class of identifiable zero-inflated count and ordered data models is generally more restricted; see, for example, Hall and Shen [16] and references there in. Although the parameters in the ZIBOP model above are identified through a nonlinear functional form estimated by ML, for more robust identification we can use traditional exclusion restrictions by including instrumental variables in the inflation equation, but excluding them from the ordered choice submodel. We follow this strategy in the empirical section.About 2/3 of the observations in our tobacco application below have a double-zero-state,(y1=0,y2=0). Consequently, we focused on a mixture constructed from a point mass at (0,0) and a bivariate ordered probit. In addition to allowing for inflation in the double-zero-state, our approach can be extended to allow for zero-inflation in each component.


### 2.2. Marginal Effects

It is common to use marginal or partial effects to interpret covariate effects in nonlinear models; see, for example, Liu et al. [17]. Due to the nonlinearity in zero-inflated ordered response models and in addition to estimation of regression parameters, it is essential to obtain the marginal effects of changes in covariates on various probabilities of interest. These include the effects of covariates on probability of nonparticipation (zero-inflation), probability of participation, and joint and/or marginal probabilities of choice associated with different levels of consumption.From a practical point of view, we are less interested in the marginal effects of explanatory variables on the joint probabilities of choice from ZIBOP. Instead, we focus on the marginal effects associated with the marginal distributions ofyri for r=1,2. Define a generic (scalar) covariate wi that can be a binary or approximately continuous variable. We obtain the marginal effects of a generic covariate wi on various probabilities assuming that the regression  results are based on ZIBOP. If wi is a binary regressor, then the marginal effect of wi on probability, say P, is the difference in the probability evaluated at 1 and 0, conditional on observable values of covariates: P(wi=1)-P(wi=0). For continuous explanatory variables, the marginal effect is given by the partial derivative of the probability of interest with respect to wi, ∂P(·)/∂wi.Regressorwi can be a common covariate in vectors of regressors xri and zi or appears in either xri or zi. Focusing on the continuous regressor case, the marginal effects of wi in each of the three cases are presented below. First, consider the case of common covariate in participation and main parts of  the model, that is, wi in both xri and zi. The marginal effect on the probability of participation is given by(2.10)Mi(si=1)=∂Pr(si=1)∂wi=ϕ(zi′γ)γwi,
where again ϕ(·) is the probability density function (pdf) of the standard normal distribution and γwi is the coefficient in the inflation part associated with variable wi. In terms of  the zeros category, the effect on the probability of nonparticipation (zero inflation) is(2.11)Mi(si=0)=∂Pr(si=0)∂wi=-ϕ(-zi′γ)γwi,
while(2.12)Mi(s=1,ỹri=0)=∂Pr(si=1)Pr(ỹri=0)∂wi=Φ(-xri′βr)ϕ(zi′γ)γwi-Φ(zi′γ)ϕ(-xri′βr)βrwi,r=1,2,
represents the marginal effect on the probability of zero-consumption. Here the scalar βrwi is the coefficient in the main part of the model associated with wi.Continuing with the case of common covariate, the marginal effects ofwi on the probabilities of choice are given as follows. First, the total marginal effect on the probability of observing zero-consumption is obtained as a sum of the marginal effects in (2.11) and (2.12); that is,(2.13)Mi(yri=0)=[Φ(-xri′βr)-1]ϕ(zi′γ)γwi-Φ(zi′γ)ϕ(-xri′βr)βrwi.
The effects for the remaining choices for outcomes r=1,2 are as follows:(2.14)Mi(yri=1)=[Φ(αr2-xri′βr)-Φ(-xri′βr)]ϕ(zi′γ)γwi-Φ(zi′γ)[ϕ(αr2-xri′βr)-ϕ(-xri′βr)]βrwi;Mi(yri=j)=[Φ(αr,j+1-xri′βr)-Φ(αrj-xri′βr)]ϕ(zi′γ)γwi-Φ(zi′γ)[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]βrwi,forj=2,…,Jr-1;Mi(yri=Jr)=[1-Φ(αr,Jr-xri′βr)]ϕ(zi′γ)γwi+Φ(zi′γ)ϕ(αr,Jr-xri′βr)βrwi.Now consider case 2, where a generic independent variablewi is included only in xri, the main part of the model. In this case, covariate wi has obviously no direct effect on the inflation part. The marginal effects of wi on various choice probabilities can be presented as follows:(2.15)Mi(yri=j)=∂Pr(yri=j)∂wi=-Φ(zi′γ)[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]βrwi,forj=0,1,…,Jr,
with αr0=-∞, αr1=0, and αr,Jr+1=∞. The marginal effects in (2.15) can be obtained by simply setting γwi=0 in (2.13) and (2.14).For case 3, wherewi appears only in zi, its marginal effects on participation components given in (2.10) and (2.11) will not change. Since βrwi=0 in case 3, the partial effects of wi on various choice probabilities take the form:(2.16)Mi(yri=j)=[Φ(αr,j+1-xri′βr)-Φ(αrj-xri′βr)]ϕ(zi′γ)γwiforj=0,1,…,Jr.
Again, we impose the restrictions αr0=-∞, αr1=0 and αr,Jr+1=∞.As noted by a referee, it is important to understand the sources of covariate effects and the relationship between the marginal effects and the coefficient estimates. Since(2.17)Pr(yri=j)=[Pr(si=1)Pr(ỹri=j)]
for j=0,1,…,Jr, the total effect of a generic covariate wi on probability of consumption at level j comes from two (weighted) sources: the participation part (Pr(si=1)) and the main ordered probit part (Pr(ỹri=j)) such that(2.18)∂Pr(si=1)∂wi=ϕ(zi′γ)γwi;(2.19)∂Pr(ỹri=j)∂wi=-[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]βrwi
with αr0=-∞, αr1=0,s and αr,Jr+1=∞. This shows that sign(γwi) is the same as sign(∂Pr(si=1)/∂wi)—the participation effect in (2.18)—but sign(βrwi) is not necessarily the same as the sign of (∂Pr(ỹri=j)/∂wi). The latter is particularly true in the left tail of the distribution, where the coefficient (βrwi) and the main (unweighted) effect in (2.19) have opposite signs because (2.20){-[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]}≡ϖ
is negative. In this case, a positive effect coming from the main part requires βrwi to be negative. By contrast, ϖ is positive in the right tail, but can be positive or negative when the terms (αr,j-xri′βr) and (αr,j+1-xri′βr) are on the opposite sides of the mode of the distribution. This shows that a given covariate can have opposite effects in the participation and main models. Since the total effect of an explanatory variable on probability of choice is a weighted average of (2.18) and (2.19), interpretation of results should focus on marginal effects of covariates rather than the signs of estimated coefficients. This is the strategy adopted in the empirical analysis below.


### 2.3. A Special Case

Since the zero-inflated univariate ordered probit (ZIOP) model has not been analyzed previously in the Bayesian framework, we provide a brief sketch of the basic framework for ZIOP. The univariate ordered probit model with excess of zeros can be obtained as a special case of the ZIBOP model presented previously. To achieve this, letρ12=0 in the ZIBOP model and focus on the first ordered outcome with r=1. In the standard ordered response approach, the model for the latent variable ỹ1i* is given by (2.1) with r=1. The observed ordered variable ỹ1i can be presented compactly as(2.21)ỹ1i=∑j=0JjI(α1j<ỹ1i*≤α1j+1),
where I(w∈A) is the indicator function equal to 1 or 0 according to whether w∈A or not. Again α10,α11,…,α1J1 are unknown threshold parameters, where we set α10=-∞, α11=0, and α1J1+1=∞.Zero-inflation is now introduced at pointỹ1i=0. Using the latent variable model (2.4) for the zero inflation, the observed binary variable is given by si=I(si*>0), where I(si*>0)=1 if si*>0, and 0 otherwise. In regime 1, si=1  or si*>0 for participants (e.g., smokers), while, in regime 0, si=0  or si*≤0 for nonparticipants. In the context of the zero-inflation model, the observed response variable takes the form y1i=siỹ1i. We observe y1i=0 when either the individual is a non-participant (si=0) or the individual is a zero-consumption participant (si=1 and ỹ1i=0). Likewise, we observe positive outcome (consumption) when the individual is a positive consumption participant (si=1andỹ1i*>0).Assume thatϵ1 and μ are independently distributed. Harris and Zhao [15] also consider the case where ϵ1 and μ are correlated. In the context of our application, the correlated model did not provide improvements over the uncorrelated ZIOP in terms of deviance information criterion. The zero-inflated ordered multinomial distribution, say Pr(y1i), arises as a mixture of a degenerate distribution at zero and the assumed distribution of the response variable ỹ1i as follows:(2.22)f1(y1i*,y1i,si*,si∣x1i,zi,Ψ1)={Pr(si=0)+Pr(si=1)Pr(ỹ1i=0),forj=0Pr(si=1)Pr(ỹ1i=j),forj=1,2,…,J1,
where, for any parameter vector Ω10 associated with the distribution of (ϵ1,μ), Ψ1=(β1,γ,α1,Ω10) with α1=(α12,…,α1J1). For simplicity, dependence on latent variables, covariates, and parameters has been suppressed on the right-hand side of (2.22). The likelihood based on N-independent observations takes the form(2.23)L1(y1*,y1,s*,s∣x1,z,Ψ1)=∏i=1N∏j=0J1[Pr(y1i=j∣x1i,zi,Ψ1)]dij=∏i=1N∏j=0[Pr(si=0)+Pr(si=1)Pr(ỹ1i=j)]dij×∏i=1N∏j>0[Pr(si=1)Pr(ỹ1i=j)]dij,
where, for example, y1*=(y1*,…,yN*)′, and dij=1 if individual i chooses outcome j, or dij=0 otherwise.Different choices of the specification of the joint distribution of(ϵ1i,μi) give rise to various zero-inflated ordered response models. For example, if the disturbance terms in the latent variable equations are normally distributed, we get the zero-inflated ordered probit model of  Harris and Zhao [15]. The zero-inflated ordered logit model can be obtained by assuming that ϵ1i and μi are independent, each of the random variables following the logistic distribution with cumulative distribution function defined as Λ(a)=ea/(1+ea). Unlike the ordered probit framework, the ordered logit cannot lend itself easily to allow for correlation between bivariate discrete response outcomes. Henceforth, we focus on the ordered probit paradigm in both univariate and bivariate settings.Assuming thatϵ1i and μi are independently normally distributed, each with mean 0 and variance 1, the required components in (2.22) and consequently (2.23) are given by:(2.24)Pr(si=0)=Φ(-zi′γ),Pr(ỹ1i=0)=Φ(-x1i′β1),Pr(ỹ1i=j)=Φ(α1j+1-x1i′β1)-Φ(α1j-x1i′β1),forj=1,…,J1-1withα10=0,Pr(ỹ1i=J1)=1-Φ(α1J1-x1i′β1).
The marginal effects for the univariate ZIOP are given by Harris and Zhao [15]. Bayesian analysis of the univariate ZIOP will be treated as a special case of the zero-inflated bivariate order probit model in the next section.


## 2.1. The Basic Model

We consider the basic Bayesian approach to a bivariate latent variable regression model with excess of zeros. To develop notation, letỹ1i* and ỹ2i*  denote the bivariate latent variables. We consider two observed ordered response variables ỹ1i and ỹ2i taking on values 0,1,…,Jr, for r=1,2. Define two sets of cut-off parameters αr=(αr2,αr3,…,αrJr), r=1,2, where the restrictions αr0=-∞, αrJr+1=∞, and αr1=0 have been imposed. We assume that (ỹ1i*,ỹ2i*)′≡ỹi* follows a bivariate regression model(2.1)ỹri*=xri′βr+ɛri,r=1,2,
where xri is a Kr-variate of regressors for the ith individual (i=1,…,N) and  ɛri are the error terms. For subsequent analysis, let β=(β1′,β2′)′, ϵi=(ϵ1i,ϵ2i)′, and (2.2)Xi=(x1i′0′0′x2i′).
Analogous to the univariate case, the observed bivariate-dependent variables are defined as(2.3)ỹri={0ifỹri*≤0,1if0<ỹri*≤αr2,jifαrj<ỹri*≤αrj+1,j=2,3,…,Jr-1,Jrifỹi*≤αrJr,
where r=1,2. Let ỹi=(ỹ1i,ỹ2i)′.We introduce inflation at the point(ỹ1i=0,ỹ2i=0), called the zero-zero state. As in the univariate case, define the participation model:(2.4)si*=zi′γ+μi,si=I(si*>0).
In the context of the zero-inflation model, the observed response random vector yi=(y1i,y2i)′ takes the form(2.5)yi=siỹi.
We observe yi=0 when either the individual is a non-participant (si=0) or the individual is a zero-consumption participant (si=1 and ỹi=0).  Likewise,  we observe positive outcome (consumption) when the individual is a positive consumption participant for at least one good (si=1 and ỹi≠0).LetΦ(a) and ϕ(a) denote the respective cumulative distribution and probability density functions of standardized normal evaluated at a. Assuming normality and that μi is uncorrelated with (ɛ1i,ɛ2i), but corr(ɛ1i,ɛ2i)=ρ12≠0, and each component with unit variance, the zero-inflated bivariate ordered probit (ZIBOP) distribution is(2.6)fb(yi*,yi,si*,si∣Xi,zi,Ψ)={Pr⁡(si=0)+(1-Pr⁡(si=0))Pr⁡(ỹ1i=0,ỹ2i=0),for(ỹ1i,ỹ2i)=(0,0)(1-Pr⁡(si=0))Pr⁡(ỹ1i=j,ỹ2i=l),for(ỹ1i,ỹ2i)≠(0,0),
where j=0,1,…,J1,l=0,1,…,J2,Pr(si=0)=Φ(-zi′γ),Pr(si=1)=Φ(-zi′γ). Further, for (ỹ1i,ỹ2i)=(0,0) in (2.6), we have αr0=-∞,αr1=0 for r=1,2 so that(2.7)Pr(ỹ1i=0,ỹ2i=0)=Φ2(-x1i′β1,-x2i′β2,ρ12),
where Φ2(·) is the cdf for the standardized bivariate normal. Likewise, Pr(ỹ1i=j,ỹ2i=l) in (2.6) are given by(2.8)Pr(ỹ1i=j,ỹ2i=l)=Φ2(α1j+1-x1i′β1,α2l+1-x2i′β2;ρ12)-Φ2(α1j-x1i′β1,α2l-x2i′β2,ρ12)forj=1,…,J1-1;l=1,…,J2-1;Pr(ỹ1i=J1,ỹ2i=J2)=1-Φ2(α1J1-x1i′β1,α2J2-x2i′β2,ρ12).
The ensuing likelihood contribution for N-independent observations is(2.9)Lb(y*,y,s*,s∣X,z,Ψb)=∏i=1N∏(j,l)=(0,0)[Pr(si=0)+(1-Pr(si=0))Pr(ỹ1i=0,ỹ2i=0)]dijl×∏i=1N∏(j,l)≠(0,0)[(1-Pr(si=0))Pr(ỹ1i=j,ỹ2i=l)]dijl,
where dijl=1 if ỹ1i=j and ỹ2i=l, and dijl=0 otherwise. Here, the vector Ψb consists of β,γ,α1,α2, and the parameters associated with the trivariate distribution of (ϵ,μ).Regarding identification of the parameters in the model defined by (2.1) through (2.5) with normality assumption, we note that the mean parameter (joint choice probability associate with the observed response vector yi) depends nonlinearly on the probability of zero inflation (Φ(-zi′γ)) and choice probability (Pr(ỹ1i=j,ỹ2i=l)) coming from the BOP submodel. Since the likelihood function for ZIBOP depends separately on the two regression components, the parameters of ZIBOP model with covariates are identified as long as the model is estimated by full maximum likelihood method. The same or different sets of covariates can affect the two components via zi and xri. When using quasi-likelihood estimation or generalized estimating equations methods rather than full ML, the class of identifiable zero-inflated count and ordered data models is generally more restricted; see, for example, Hall and Shen [16] and references there in. Although the parameters in the ZIBOP model above are identified through a nonlinear functional form estimated by ML, for more robust identification we can use traditional exclusion restrictions by including instrumental variables in the inflation equation, but excluding them from the ordered choice submodel. We follow this strategy in the empirical section.About 2/3 of the observations in our tobacco application below have a double-zero-state,(y1=0,y2=0). Consequently, we focused on a mixture constructed from a point mass at (0,0) and a bivariate ordered probit. In addition to allowing for inflation in the double-zero-state, our approach can be extended to allow for zero-inflation in each component.


## 2.2. Marginal Effects

It is common to use marginal or partial effects to interpret covariate effects in nonlinear models; see, for example, Liu et al. [17]. Due to the nonlinearity in zero-inflated ordered response models and in addition to estimation of regression parameters, it is essential to obtain the marginal effects of changes in covariates on various probabilities of interest. These include the effects of covariates on probability of nonparticipation (zero-inflation), probability of participation, and joint and/or marginal probabilities of choice associated with different levels of consumption.From a practical point of view, we are less interested in the marginal effects of explanatory variables on the joint probabilities of choice from ZIBOP. Instead, we focus on the marginal effects associated with the marginal distributions ofyri for r=1,2. Define a generic (scalar) covariate wi that can be a binary or approximately continuous variable. We obtain the marginal effects of a generic covariate wi on various probabilities assuming that the regression  results are based on ZIBOP. If wi is a binary regressor, then the marginal effect of wi on probability, say P, is the difference in the probability evaluated at 1 and 0, conditional on observable values of covariates: P(wi=1)-P(wi=0). For continuous explanatory variables, the marginal effect is given by the partial derivative of the probability of interest with respect to wi, ∂P(·)/∂wi.Regressorwi can be a common covariate in vectors of regressors xri and zi or appears in either xri or zi. Focusing on the continuous regressor case, the marginal effects of wi in each of the three cases are presented below. First, consider the case of common covariate in participation and main parts of  the model, that is, wi in both xri and zi. The marginal effect on the probability of participation is given by(2.10)Mi(si=1)=∂Pr(si=1)∂wi=ϕ(zi′γ)γwi,
where again ϕ(·) is the probability density function (pdf) of the standard normal distribution and γwi is the coefficient in the inflation part associated with variable wi. In terms of  the zeros category, the effect on the probability of nonparticipation (zero inflation) is(2.11)Mi(si=0)=∂Pr(si=0)∂wi=-ϕ(-zi′γ)γwi,
while(2.12)Mi(s=1,ỹri=0)=∂Pr(si=1)Pr(ỹri=0)∂wi=Φ(-xri′βr)ϕ(zi′γ)γwi-Φ(zi′γ)ϕ(-xri′βr)βrwi,r=1,2,
represents the marginal effect on the probability of zero-consumption. Here the scalar βrwi is the coefficient in the main part of the model associated with wi.Continuing with the case of common covariate, the marginal effects ofwi on the probabilities of choice are given as follows. First, the total marginal effect on the probability of observing zero-consumption is obtained as a sum of the marginal effects in (2.11) and (2.12); that is,(2.13)Mi(yri=0)=[Φ(-xri′βr)-1]ϕ(zi′γ)γwi-Φ(zi′γ)ϕ(-xri′βr)βrwi.
The effects for the remaining choices for outcomes r=1,2 are as follows:(2.14)Mi(yri=1)=[Φ(αr2-xri′βr)-Φ(-xri′βr)]ϕ(zi′γ)γwi-Φ(zi′γ)[ϕ(αr2-xri′βr)-ϕ(-xri′βr)]βrwi;Mi(yri=j)=[Φ(αr,j+1-xri′βr)-Φ(αrj-xri′βr)]ϕ(zi′γ)γwi-Φ(zi′γ)[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]βrwi,forj=2,…,Jr-1;Mi(yri=Jr)=[1-Φ(αr,Jr-xri′βr)]ϕ(zi′γ)γwi+Φ(zi′γ)ϕ(αr,Jr-xri′βr)βrwi.Now consider case 2, where a generic independent variablewi is included only in xri, the main part of the model. In this case, covariate wi has obviously no direct effect on the inflation part. The marginal effects of wi on various choice probabilities can be presented as follows:(2.15)Mi(yri=j)=∂Pr(yri=j)∂wi=-Φ(zi′γ)[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]βrwi,forj=0,1,…,Jr,
with αr0=-∞, αr1=0, and αr,Jr+1=∞. The marginal effects in (2.15) can be obtained by simply setting γwi=0 in (2.13) and (2.14).For case 3, wherewi appears only in zi, its marginal effects on participation components given in (2.10) and (2.11) will not change. Since βrwi=0 in case 3, the partial effects of wi on various choice probabilities take the form:(2.16)Mi(yri=j)=[Φ(αr,j+1-xri′βr)-Φ(αrj-xri′βr)]ϕ(zi′γ)γwiforj=0,1,…,Jr.
Again, we impose the restrictions αr0=-∞, αr1=0 and αr,Jr+1=∞.As noted by a referee, it is important to understand the sources of covariate effects and the relationship between the marginal effects and the coefficient estimates. Since(2.17)Pr(yri=j)=[Pr(si=1)Pr(ỹri=j)]
for j=0,1,…,Jr, the total effect of a generic covariate wi on probability of consumption at level j comes from two (weighted) sources: the participation part (Pr(si=1)) and the main ordered probit part (Pr(ỹri=j)) such that(2.18)∂Pr(si=1)∂wi=ϕ(zi′γ)γwi;(2.19)∂Pr(ỹri=j)∂wi=-[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]βrwi
with αr0=-∞, αr1=0,s and αr,Jr+1=∞. This shows that sign(γwi) is the same as sign(∂Pr(si=1)/∂wi)—the participation effect in (2.18)—but sign(βrwi) is not necessarily the same as the sign of (∂Pr(ỹri=j)/∂wi). The latter is particularly true in the left tail of the distribution, where the coefficient (βrwi) and the main (unweighted) effect in (2.19) have opposite signs because (2.20){-[ϕ(αr,j+1-xri′βr)-ϕ(αrj-xri′βr)]}≡ϖ
is negative. In this case, a positive effect coming from the main part requires βrwi to be negative. By contrast, ϖ is positive in the right tail, but can be positive or negative when the terms (αr,j-xri′βr) and (αr,j+1-xri′βr) are on the opposite sides of the mode of the distribution. This shows that a given covariate can have opposite effects in the participation and main models. Since the total effect of an explanatory variable on probability of choice is a weighted average of (2.18) and (2.19), interpretation of results should focus on marginal effects of covariates rather than the signs of estimated coefficients. This is the strategy adopted in the empirical analysis below.


## 2.3. A Special Case

Since the zero-inflated univariate ordered probit (ZIOP) model has not been analyzed previously in the Bayesian framework, we provide a brief sketch of the basic framework for ZIOP. The univariate ordered probit model with excess of zeros can be obtained as a special case of the ZIBOP model presented previously. To achieve this, letρ12=0 in the ZIBOP model and focus on the first ordered outcome with r=1. In the standard ordered response approach, the model for the latent variable ỹ1i* is given by (2.1) with r=1. The observed ordered variable ỹ1i can be presented compactly as(2.21)ỹ1i=∑j=0JjI(α1j<ỹ1i*≤α1j+1),
where I(w∈A) is the indicator function equal to 1 or 0 according to whether w∈A or not. Again α10,α11,…,α1J1 are unknown threshold parameters, where we set α10=-∞, α11=0, and α1J1+1=∞.Zero-inflation is now introduced at pointỹ1i=0. Using the latent variable model (2.4) for the zero inflation, the observed binary variable is given by si=I(si*>0), where I(si*>0)=1 if si*>0, and 0 otherwise. In regime 1, si=1  or si*>0 for participants (e.g., smokers), while, in regime 0, si=0  or si*≤0 for nonparticipants. In the context of the zero-inflation model, the observed response variable takes the form y1i=siỹ1i. We observe y1i=0 when either the individual is a non-participant (si=0) or the individual is a zero-consumption participant (si=1 and ỹ1i=0). Likewise, we observe positive outcome (consumption) when the individual is a positive consumption participant (si=1andỹ1i*>0).Assume thatϵ1 and μ are independently distributed. Harris and Zhao [15] also consider the case where ϵ1 and μ are correlated. In the context of our application, the correlated model did not provide improvements over the uncorrelated ZIOP in terms of deviance information criterion. The zero-inflated ordered multinomial distribution, say Pr(y1i), arises as a mixture of a degenerate distribution at zero and the assumed distribution of the response variable ỹ1i as follows:(2.22)f1(y1i*,y1i,si*,si∣x1i,zi,Ψ1)={Pr(si=0)+Pr(si=1)Pr(ỹ1i=0),forj=0Pr(si=1)Pr(ỹ1i=j),forj=1,2,…,J1,
where, for any parameter vector Ω10 associated with the distribution of (ϵ1,μ), Ψ1=(β1,γ,α1,Ω10) with α1=(α12,…,α1J1). For simplicity, dependence on latent variables, covariates, and parameters has been suppressed on the right-hand side of (2.22). The likelihood based on N-independent observations takes the form(2.23)L1(y1*,y1,s*,s∣x1,z,Ψ1)=∏i=1N∏j=0J1[Pr(y1i=j∣x1i,zi,Ψ1)]dij=∏i=1N∏j=0[Pr(si=0)+Pr(si=1)Pr(ỹ1i=j)]dij×∏i=1N∏j>0[Pr(si=1)Pr(ỹ1i=j)]dij,
where, for example, y1*=(y1*,…,yN*)′, and dij=1 if individual i chooses outcome j, or dij=0 otherwise.Different choices of the specification of the joint distribution of(ϵ1i,μi) give rise to various zero-inflated ordered response models. For example, if the disturbance terms in the latent variable equations are normally distributed, we get the zero-inflated ordered probit model of  Harris and Zhao [15]. The zero-inflated ordered logit model can be obtained by assuming that ϵ1i and μi are independent, each of the random variables following the logistic distribution with cumulative distribution function defined as Λ(a)=ea/(1+ea). Unlike the ordered probit framework, the ordered logit cannot lend itself easily to allow for correlation between bivariate discrete response outcomes. Henceforth, we focus on the ordered probit paradigm in both univariate and bivariate settings.Assuming thatϵ1i and μi are independently normally distributed, each with mean 0 and variance 1, the required components in (2.22) and consequently (2.23) are given by:(2.24)Pr(si=0)=Φ(-zi′γ),Pr(ỹ1i=0)=Φ(-x1i′β1),Pr(ỹ1i=j)=Φ(α1j+1-x1i′β1)-Φ(α1j-x1i′β1),forj=1,…,J1-1withα10=0,Pr(ỹ1i=J1)=1-Φ(α1J1-x1i′β1).
The marginal effects for the univariate ZIOP are given by Harris and Zhao [15]. Bayesian analysis of the univariate ZIOP will be treated as a special case of the zero-inflated bivariate order probit model in the next section.


## 3. Bayesian Analysis

### 3.1. Prior Distributions

The Bayesian hierarchical model requires prior distributions for each parameter in the model.  For this purpose, we can use noninformative conjugate priors. There are two reasons for adopting noninformative conjugate priors. First, we prefer to let the data dictate the inference about the parameters with little or no influence from prior distributions. Secondly, the noninformative priors facilitate resampling using Markov Chain Monte Carlo algorithm (MCMC) and have nice convergence properties. We assume noninformative (vague or diffuse) normal priors for regression coefficientsβ, with mean β* and variance Ωβ which are chosen to make the distribution proper but diffuse with large variances. Similarly, γ~N(γ*,Ωγ).In choosing prior distributions for the threshold parameters,α’s, caution is needed because of the order restriction on them. One way to avoid the order restriction is to reparameterize them. Following Chib and Hamilton [18] treatment in the univariate ordered probit case, we reparameterize the ordered threshold parameters(3.1)τr2=log⁡(αr2);τrj=log⁡(αrj-αrj-1),j=3,…,Jr;r=1,2
with the inverse map(3.2)αrj=∑m=2jexp⁡(τrm),j=2,…,Jr;r=1,2.
For r=1,2, let τr=(τr2,τr3,…,τrJ)′ so that τ=(τ1,τ2). We choose normal prior τ~N(τ*,Ωτ) without order restrictions for τr’s.The only unknown parameter associate with the distribution of(ϵ,μ) in (2.1) and (2.4) is ρ12, the correlation between ϵ1 and ϵ2. The values of ρ12 by definition are restricted to be in the −1 to 1 interval. Therefore, the choice for prior distribution for ρ12 can be uniform(-1,1) or a proper distribution based on reparameterization. Let ν denote the hyperbolic arc-tangent transformation of ρ12, that is, (3.3)ν=atanh(ρ12),
and taking hyperbolic tangent transformation of ν gives us back ρ12=tanh(ν). Then parameter ν is asymptotically normal distributed with stabilized variance, 1/(N-3), where N is the sample size. We may also assume that ν~N(ν*,σν2).


### 3.2. Bayesian Analysis via MCMC

For carrying out a Bayesian inference, the joint posterior distribution of the parameters of the ZIBOP model in (2.6) conditional on the data is obtained by combining the likelihood function given in (2.9) and the above-specified prior distributions via Bayes’ theorem, as:(3.4)f(Ψb∣x,z)∝∏i=1N∏(j,l)=(0,0)[Φ(-zi′γ)+Φ(zi′γ)Φ2(-x1i′β1,-x2i′β2,ρ12)]dijl×∏i=1N∏(j,l)≠(0,0)[Φ(zi′γ)[Φ2(α1j+1-x1i′β1,α2l+1-x2i′β2;ρ12)-Φ2(α1j-x1i′β1,α2l-x2i′β2,ρ12)]]dijl×f(Ψb),
where f(Ψb)∝f(β)f(γ)f(τ)f(ν)  and the  parameter vector Ψb  now consists of β=(β1′,β2′)′, γ, τ=(τ1,τ2),s and ν=atanh(ρ12). Here f(β)∝|Ωβ|-1/2exp{-1/2(β-β*)′Ωβ-1(β-β*)};f(γ)∝|Ωγ|-1/2exp{-1/2(γ-γ*)′Ωγ-1(γ-γ*)};f(τ)∝|Ωτ|-1/2exp{-1/2(τ-τ*)′Ωτ-1(τ-τ*)};τrj are defined in (3.1),  and αrj are given via the inverse map (3.2).Full conditional posterior distributions are required to implement the MCMC algorithm [19–22], and they are given as follows:      (1)
fixed effects:
(a)
zero state:(3.5)f(γ∣x,z,Ψ-γ)∝|Ωγ|-1/2exp⁡{-12(γ-γ*)′Ωγ-1(γ-γ*)}×f(Ψb∣x,z);(b)
nonzero state:(3.6)f(β∣x,z,Ψ-β)∝|Ωβ|-1/2exp⁡{-12(β-β*)′Ωβ-1(β-β*)}×f(Ψb∣x,z).(2)
thresholds:(3.7)f(τ∣x,z,Ψ-τ)∝|Ωτ|-1/2exp⁡{-12(τ-τ*)′Ωτ-1(τ-τ*)}×∏i=1N∏(j,l)≠(0,0)[Φ(zi′γ)[Φ2(α1j+1-x1i′β1,α2l+1-x2i′β2;ρ12)-Φ2(α1j-x1i′β1,α2l-x2i′β2,ρ12)]]dijl.(3)
bivariate correlation:(3.8)f(ν∣x,z,Ψ-ν)∝σν-1exp⁡{-(ν-ν*)22σν2}×f(Ψb∣x,z).The MCMC algorithm simulates direct draws from the above full conditionals iteratively until convergence is achieved. A single long chain [23, 24] is used for the proposed model. Geyer [23] argues that using a single longer chain is better than using a number of smaller chains with different initial values. We follow this strategy in our empirical analysis.The Bayesian analysis of the univariate ZIOP follows as a special case of that of the ZIBOP presented above. In particular, the joint posterior distribution of the parameters of the ZIOP model in (2.22) conditional on the data is obtained by combining the likelihood function given in (2.23) and the above-specified prior distributions (with modified notations) via Bayes' theorem, as follows: (3.9)f(Ψ∣x,z,)∝∏i=1N∏j=0[Φ(-zi′γ)+Φ(zi′γ)Φ(-xi′β)]dij×∏i=1N∏j>0[Φ(zi′γ){Φ(αj+1-xi′β)-Φ(αj-xi′β)}]dij×f(β)f(γ)f(τ),
where, using notation of Section 2.3 for β and the other parameter vectors, f(β)∝|Ωβ|-1/2exp{-1/2(β-β*)′Ωβ-1(β-β*)};f(γ)∝|Ωγ|-1/2exp{-1/2(γ-γ*)′Ωγ-1(γ-γ*)};   f(τ)∝|Ωτ|-1/2exp{-1/2(τ-τ*)′Ωτ-1(τ-τ*)}, τ2=log(α2) and τj=log(αj-αj-1),j=3,…,J.   Apart from dropping the bivariate correlation, we basically replace the bivariate normal  cumulative distribution Φ2(·,·;ρ12) by the univariate counterpart Φ(·). Details are available upon request from the authors.Apart from Bayesian estimation of the regression parameters, the posterior distributions of other quantities of interest can be obtained. These include posteriors for marginal effects and probabilities for nonparticipation, zero-consumption, and joint outcomes of interest. These will be considered in the application section. Next, we summarize model selection procedure.The commonly used criteria for model selection like BIC and AIC are not appropriate for the multilevel models (in the presence of random effects), which complicates the counting of the true number of free parameters. To overcome such a hurdle, Spiegelhalter et al. [25] proposed a Bayesian model comparison criterion, called Deviance Information Criterion (DIC).  It is given as (3.10)DIC=goodness-of-fit+penaltyforcomplexity,
where the “goodness-of-fit” is measured by the deviance for θ=(β,γ,α)(3.11)D(θ)=-2log⁡L(data∣θ)
and complexity is measured by the “effective number of parameters”:(3.12)pD=Eθ|y[D(θ)]-D(Eθ|y[θ])=D¯-D(θ¯);
that is, posterior mean deviance minus deviance evaluated at the posterior mean of the parameters. The DIC is then defined analogously to AIC as(3.13)DIC=D(θ¯)+2pD=D¯+pD.
The idea here is that models with smaller DIC should be preferred to models with larger DIC. Models are penalized both by the value of D¯, which favors a good fit, but also (similar to AIC and BIC) by the effective number of parameters pD. The advantage of DIC over other criteria, for Bayesian model selection, is that the DIC is easily calculated from the MCMC samples. In contrast, AIC and BIC require calculating the likelihood at its maximum values, which are not easily available from the MCMC simulation.


## 3.1. Prior Distributions

The Bayesian hierarchical model requires prior distributions for each parameter in the model.  For this purpose, we can use noninformative conjugate priors. There are two reasons for adopting noninformative conjugate priors. First, we prefer to let the data dictate the inference about the parameters with little or no influence from prior distributions. Secondly, the noninformative priors facilitate resampling using Markov Chain Monte Carlo algorithm (MCMC) and have nice convergence properties. We assume noninformative (vague or diffuse) normal priors for regression coefficientsβ, with mean β* and variance Ωβ which are chosen to make the distribution proper but diffuse with large variances. Similarly, γ~N(γ*,Ωγ).In choosing prior distributions for the threshold parameters,α’s, caution is needed because of the order restriction on them. One way to avoid the order restriction is to reparameterize them. Following Chib and Hamilton [18] treatment in the univariate ordered probit case, we reparameterize the ordered threshold parameters(3.1)τr2=log⁡(αr2);τrj=log⁡(αrj-αrj-1),j=3,…,Jr;r=1,2
with the inverse map(3.2)αrj=∑m=2jexp⁡(τrm),j=2,…,Jr;r=1,2.
For r=1,2, let τr=(τr2,τr3,…,τrJ)′ so that τ=(τ1,τ2). We choose normal prior τ~N(τ*,Ωτ) without order restrictions for τr’s.The only unknown parameter associate with the distribution of(ϵ,μ) in (2.1) and (2.4) is ρ12, the correlation between ϵ1 and ϵ2. The values of ρ12 by definition are restricted to be in the −1 to 1 interval. Therefore, the choice for prior distribution for ρ12 can be uniform(-1,1) or a proper distribution based on reparameterization. Let ν denote the hyperbolic arc-tangent transformation of ρ12, that is, (3.3)ν=atanh(ρ12),
and taking hyperbolic tangent transformation of ν gives us back ρ12=tanh(ν). Then parameter ν is asymptotically normal distributed with stabilized variance, 1/(N-3), where N is the sample size. We may also assume that ν~N(ν*,σν2).


## 3.2. Bayesian Analysis via MCMC

For carrying out a Bayesian inference, the joint posterior distribution of the parameters of the ZIBOP model in (2.6) conditional on the data is obtained by combining the likelihood function given in (2.9) and the above-specified prior distributions via Bayes’ theorem, as:(3.4)f(Ψb∣x,z)∝∏i=1N∏(j,l)=(0,0)[Φ(-zi′γ)+Φ(zi′γ)Φ2(-x1i′β1,-x2i′β2,ρ12)]dijl×∏i=1N∏(j,l)≠(0,0)[Φ(zi′γ)[Φ2(α1j+1-x1i′β1,α2l+1-x2i′β2;ρ12)-Φ2(α1j-x1i′β1,α2l-x2i′β2,ρ12)]]dijl×f(Ψb),
where f(Ψb)∝f(β)f(γ)f(τ)f(ν)  and the  parameter vector Ψb  now consists of β=(β1′,β2′)′, γ, τ=(τ1,τ2),s and ν=atanh(ρ12). Here f(β)∝|Ωβ|-1/2exp{-1/2(β-β*)′Ωβ-1(β-β*)};f(γ)∝|Ωγ|-1/2exp{-1/2(γ-γ*)′Ωγ-1(γ-γ*)};f(τ)∝|Ωτ|-1/2exp{-1/2(τ-τ*)′Ωτ-1(τ-τ*)};τrj are defined in (3.1),  and αrj are given via the inverse map (3.2).Full conditional posterior distributions are required to implement the MCMC algorithm [19–22], and they are given as follows:      (1)
fixed effects:
(a)
zero state:(3.5)f(γ∣x,z,Ψ-γ)∝|Ωγ|-1/2exp⁡{-12(γ-γ*)′Ωγ-1(γ-γ*)}×f(Ψb∣x,z);(b)
nonzero state:(3.6)f(β∣x,z,Ψ-β)∝|Ωβ|-1/2exp⁡{-12(β-β*)′Ωβ-1(β-β*)}×f(Ψb∣x,z).(2)
thresholds:(3.7)f(τ∣x,z,Ψ-τ)∝|Ωτ|-1/2exp⁡{-12(τ-τ*)′Ωτ-1(τ-τ*)}×∏i=1N∏(j,l)≠(0,0)[Φ(zi′γ)[Φ2(α1j+1-x1i′β1,α2l+1-x2i′β2;ρ12)-Φ2(α1j-x1i′β1,α2l-x2i′β2,ρ12)]]dijl.(3)
bivariate correlation:(3.8)f(ν∣x,z,Ψ-ν)∝σν-1exp⁡{-(ν-ν*)22σν2}×f(Ψb∣x,z).The MCMC algorithm simulates direct draws from the above full conditionals iteratively until convergence is achieved. A single long chain [23, 24] is used for the proposed model. Geyer [23] argues that using a single longer chain is better than using a number of smaller chains with different initial values. We follow this strategy in our empirical analysis.The Bayesian analysis of the univariate ZIOP follows as a special case of that of the ZIBOP presented above. In particular, the joint posterior distribution of the parameters of the ZIOP model in (2.22) conditional on the data is obtained by combining the likelihood function given in (2.23) and the above-specified prior distributions (with modified notations) via Bayes' theorem, as follows: (3.9)f(Ψ∣x,z,)∝∏i=1N∏j=0[Φ(-zi′γ)+Φ(zi′γ)Φ(-xi′β)]dij×∏i=1N∏j>0[Φ(zi′γ){Φ(αj+1-xi′β)-Φ(αj-xi′β)}]dij×f(β)f(γ)f(τ),
where, using notation of Section 2.3 for β and the other parameter vectors, f(β)∝|Ωβ|-1/2exp{-1/2(β-β*)′Ωβ-1(β-β*)};f(γ)∝|Ωγ|-1/2exp{-1/2(γ-γ*)′Ωγ-1(γ-γ*)};   f(τ)∝|Ωτ|-1/2exp{-1/2(τ-τ*)′Ωτ-1(τ-τ*)}, τ2=log(α2) and τj=log(αj-αj-1),j=3,…,J.   Apart from dropping the bivariate correlation, we basically replace the bivariate normal  cumulative distribution Φ2(·,·;ρ12) by the univariate counterpart Φ(·). Details are available upon request from the authors.Apart from Bayesian estimation of the regression parameters, the posterior distributions of other quantities of interest can be obtained. These include posteriors for marginal effects and probabilities for nonparticipation, zero-consumption, and joint outcomes of interest. These will be considered in the application section. Next, we summarize model selection procedure.The commonly used criteria for model selection like BIC and AIC are not appropriate for the multilevel models (in the presence of random effects), which complicates the counting of the true number of free parameters. To overcome such a hurdle, Spiegelhalter et al. [25] proposed a Bayesian model comparison criterion, called Deviance Information Criterion (DIC).  It is given as (3.10)DIC=goodness-of-fit+penaltyforcomplexity,
where the “goodness-of-fit” is measured by the deviance for θ=(β,γ,α)(3.11)D(θ)=-2log⁡L(data∣θ)
and complexity is measured by the “effective number of parameters”:(3.12)pD=Eθ|y[D(θ)]-D(Eθ|y[θ])=D¯-D(θ¯);
that is, posterior mean deviance minus deviance evaluated at the posterior mean of the parameters. The DIC is then defined analogously to AIC as(3.13)DIC=D(θ¯)+2pD=D¯+pD.
The idea here is that models with smaller DIC should be preferred to models with larger DIC. Models are penalized both by the value of D¯, which favors a good fit, but also (similar to AIC and BIC) by the effective number of parameters pD. The advantage of DIC over other criteria, for Bayesian model selection, is that the DIC is easily calculated from the MCMC samples. In contrast, AIC and BIC require calculating the likelihood at its maximum values, which are not easily available from the MCMC simulation.


## 4. Application

### 4.1. Data

We consider an application to tobacco consumption behavior of individuals based on the 2001 household Tobacco Prevalence survey data from Bangladesh. The Survey was conducted in two administrative districts of paramount interest for tobacco production and consumption in the country. Data on daily consumption of smoking and chewing tobacco along with other socioeconomic and demographic characteristics and parental tobacco consumption habits were collected from respondents of 10 years of age and above. The data set has been used previously by Gurmu and Yunus [26] in the context of binary response models. Here we focus on a sample consisting of 6000 individual respondents aged between 10 and 101 years.The ordinal outcomesyr=0,1,2,3  used in this paper correspond roughly to zero, low, moderate, and high levels of tobacco consumption in the form of smoking (y1) or chewing tobacco (y2), respectively. The first dependent variable y1 for an individual's daily cigarette smoking intensities assumes the following 4 choices: y1=0 if nonsmoker, y1=1 if smoking up to 7 cigarettes per day, y=2 if smoking between 8 and 12 cigarettes daily, and y1=3 if smoking more than 12 cigarettes daily; likewise, for the intensity of chewing tobacco, y2=0 if reported not chewing tobacco, y2=1 if uses up to 7 chewing tobacco, and y2=2 if consuming 7 or more chewing tobacco. The frequency distribution of  cigarette smoking  and tobacco chewing choices in Table 1 shows that nearly 66% of the respondents identify themselves as nonusers of tobacco. Our modeling strategy recognizes that these self-identified current nonusers of tobacco may include either individuals who never smoke or chew tobacco (genuine nonusers) or those who do, but not during the reporting period (potential users of tobacco). For example, potential tobacco users may include those who wrongly claim to be nonusers, previous  tobacco users that are currently nonusers, and those most likely to use tobacco in the future due to changes in, say, prices and income. Table 1 also shows that 76% of the respondents are non-smokers and nearly 87% identify themselves as nonusers of tobacco for chewing. Given the extremely high proportion of observed zeros coupled with sparse cells on the right tail, we employ the zero-inflated bivariate ordered probit framework.Table 1
Bivariate frequency distribution for intensity of tobacco use.
Smoke groupChew groupTotal (N)01203931302324455712651262832526353759834982935562Total (N)52203784026000Table2 gives definition of the explanatory variables as well as their means and standard deviations. The respondents are more likely to be Muslim, married, in early thirties, live in rural area, and have about 7 years of formal schooling. Although the country is mostly agrarian, only around 11% of the respondents were related to agricultural occupation in either doing agricultural operations on their own farms or working as agricultural wage laborers. About 12% of the respondents belong to the service occupation. The benchmark occupational group consists of business and other occupations. More than one-half of the fathers and slightly less than two-thirds of the mothers of the respondents currently use or have used tobacco products in the past.Table 2
Definition and summary statistics for independent variables.
NameDefinitionMeanbSt. Dev.AgeaAge in years30.35(14.9)EducationaNumber of years of formal schooling6.83(4.7)IncomeMonthly family income in 1000s of Taka7.57(10.3)Male= 1 if male54.6Married= 1 if married57.2Muslim= 1 if religion is Islam78.4Father use= 1 if father uses tobacco54.0Mother use= 1 if mother uses tobacco65.1Region= 1 if Rangpur resident,= 0 if Chittagong resident49.7Urban= 1 if  urban resident38.0Agriservice= 1 if agriculture labor  or service occupation23.2Self-employed= 1 if self-employed or household chores30.7Student= 1 if student26.8Other= 1 if business or other occupations (control)19.3aIn implementation, we also include age squared and education squared.bThe means for binary variables are in percentage.Among the variables given in Table2, the two indicators of parental use of tobacco products are included in z  as part of the participation equation (2.4). The rest of the variables are included in xr and z of (2.1) and (2.4). To allow for nonlinear effects, age and education enter all three equations using a quadratic form. Due to lack of data on prices, our analysis is limited to the study of other economic and demographic determinants of participation, smoking, and chewing tobacco.


### 4.2. Results

We estimate the standard bivariate ordered probit (BOP) and zero-inflated bivariate ordered probit regression models for smoking and chewing tobacco and report estimation results for parameters, marginal effects, and choice probabilities, along measures of model selection. An earlier version of this paper reports results from the standard ordered probit model as well as the uncorrelated and correlated versions of the univariate zero-inflated ordered probit model for smoking tobacco. Convergence of the generated samples is assessed using standard tools (such as trace plots and ACF plots) within WinBUGS software. After initial 10,000 burn-in iterations, every 10th MCMC sample thereafter was retained from the next 100,000 iterations, obtaining 10,000 samples for subsequent posterior inference of the unknown parameters. The slowest convergence is observed for some parameters in the inflation submodel. By contrast, the autocorrelations functions for most of the marginal  effects die out quickly relative to those for the associated parameters.Table3 reports the goodness-of-fit statistics for the standard bivariate ordered probit model and its zero-inflated version, ZIBOP. The ZIBOP regression model clearly dominates BOP in terms of DIC and its components; compare the DIC of 11330 for the former and 11447 for the latter model. Table 4 gives posterior means, standard deviations, medians, and the 95 percent credible intervals (in terms of the 2.5 and 97.5 percentiles) of the parameters and choice probabilities from ZIBOP model. For comparison, the corresponding results from BOP are shown in Table 6 of the appendix. Both models predict significant negative correlation between the likelihood of smoking and chewing tobacco. The posterior estimates of the cut-off points are qualitatively similar across models. In what follows, we focus on discussion of results from the preferred ZIBOP model. The 95% credible interval for the correlation parameter ρ12  from the zero-inflated model is in the range −0.25 to −0.12, indicating that smoking and chewing tobacco are generally substitutes. Results of selected predicted choice probabilities (bottom of Table 4) show that the ZIBOP regression model provides very good fit to the data. The posterior mean for the probability of (zero, zero)-inflation is about 24% while the 95% credible interval is [0.15, 0.32], indicating that a substantial proportion of zeros may be attributed to nonparticipants. These results underscore the importance of modeling excess zeros in bivariate ordered probit models.Table 3
Goodness-of-fit statistics via DIC.
ModelDbarDhatpDDICBivariate ordered probit (BOP)11417.111386.930.111447.2Zero-inflated BOP11301.111270.329.811329.9Dbar:  Posterior mean of deviance, Dhat:  Deviance evaluated at the posterior mean of the parameters, pD: Dbar-Dhat, the effective number of parameters, and DIC: Deviance information criterion.Table 4
Posterior mean, standard deviation, and 95% credible intervals of parameters from zibop for smoking and chewing tobacco.
VariableMeanSt. dev.2.50%Median97.50%Main (β1,α1): smoking (y1):Age/100.6720.1190.4440.685Age square/100−0.0700.012−0.093−0.071Education−0.0710.014−0.097−0.071Education square0.0010.001−0.0020.001Income0.0000.002−0.0050.000Male2.0920.0861.9252.091Married0.2130.0700.0740.213Muslim−0.0530.052−0.157−0.053Region−0.0070.048−0.102−0.007Urban−0.0960.051−0.198−0.097Agriservice−0.2340.056−0.345−0.233Self-employed−0.2460.087−0.414−0.247student−0.4760.137−0.742−0.478α120.2840.0170.2520.283α130.9870.0300.9280.987Main (β2,α2): chewing (y2)Age/100.6490.1330.3820.658Age square/100−0.0460.013−0.071−0.046Education−0.0200.016−0.052−0.020Education square−0.0020.001−0.005−0.002Income0.0010.003−0.0040.002Male−0.4790.081−0.641−0.479Married−0.0250.075−0.171−0.025Muslim−0.0720.056−0.181−0.072Region0.4170.0510.3170.418Urban−0.0800.058−0.194−0.079Agriservice0.0520.074−0.0960.052Self-employed0.1270.092−0.0580.126Student−0.4500.221−0.887−0.448α220.4840.0230.4390.484Inflation (γ):Age/10−0.0122.044−4.7550.253Age square/1000.5090.552−0.1970.398Education−0.2180.115−0.476−0.204Education square0.0280.0110.0100.026Income0.0060.022−0.0270.003Male0.2390.827−1.5820.417Married2.3064.478−0.4160.500Muslim−0.5280.356−1.331−0.494Mother−0.1700.267−0.716−0.164Father−0.1190.330−0.664−0.160Region0.6300.2910.0610.625Urban0.0400.357−0.7370.071Agrservice5.3125.4161.0172.674Self-employed3.7835.0250.1241.275Sstudent−0.3440.411−1.154−0.339ρ12−0.1850.033−0.249−0.186Select probabilities:P (y1=0)0.7600.0040.7520.760P (y2=0)0.8710.0040.8640.871P (y1=0, y2=0)0.6620.0050.6520.662P (zero-inflation)0.2420.0480.1510.243Results for the constant terms in the main and inflation parts have  been suppressed for brevity.To facilitate interpretation of results, we report in Tables5 and 7 the same set of posterior estimates for the marginal effects from ZIBOP and BOP models, respectively. Since age and education enter the three equations non-linearly, we report the total marginal effects coming from the linear and quadratic parts. We examine closely the marginal effects on the unconditional marginal probabilities at all levels of smoking and chewing tobacco (y1=0,1,2,3; y2=0,1,2). The marginal effects reported in Table 5 show that the results for covariates are generally plausible. Age has a negative impact on probabilities of moderate and heavy use of tobacco. For heavy smokers, education has a significant negative impact on the probability of smoking cigarettes. An additional year of schooling on average decreases probability of smoking by about 6.9% for heavy smokers. Among participants, being male or married has positive impact on probability of smoking, while the effects for being Muslim, urban resident, and student are largely negative. Male respondents are more likely to smoke cigarettes while women respondents are more likely to use chewing tobacco with heavy intensity, a result which is in line with custom of the country [26].Table 5
Posterior mean, standard deviation, and 95% credible intervals of marginal effects of covariates on probability of smoking and chewing tobacco (ZIBOP model).
VariableProbabilityMeanSt. dev.2.50%Median97.50%AgeNonparticipation−0.02590.0129−0.0556−0.0236−0.0078Zero-consumption,y10.04630.01020.02940.04530.0687All zeros,y1=00.02040.00590.00780.02130.0304y1=10.00580.00350.00090.00530.0138y1=2−0.00140.0029−0.0057−0.00190.0055y1=3−0.06900.0235−0.1223−0.0658−0.0344Zero-consumption,y20.04030.01160.01950.03860.0675All zeros,y2=00.01450.00640.00180.01490.0264y2=1−0.00340.0021−0.0071−0.00350.0008y2=2−0.00190.0014−0.0043−0.00200.0011EducationNonparticipation−0.28230.0768−0.4260−0.2837−0.1252Zero-consumption,y10.24470.07490.09170.24590.3851All zeros,y1=0−0.03770.0241−0.0853−0.03740.0094y1=10.04980.01410.02310.04940.0789y1=20.02410.01020.00450.02390.0444y1=3−0.55570.1536−0.8415−0.5588−0.2417Zero-consumption,y20.31360.07720.15610.31590.4546All zeros,y2=00.03130.0161−0.00090.03150.0618y2=1−0.01340.0080−0.0288−0.01350.0027y2=2−0.02220.0119−0.0455−0.02210.0009IncomeNonparticipation−0.00040.0015−0.0038−0.00020.0022Zero-consumption,y10.00030.0014−0.00220.00020.0035All zeros,y1=0−0.00010.0004−0.0009−0.00010.0008y1=10.00010.0003−0.00050.00000.0008y1=20.00000.0002−0.00030.00000.0004y1=3−0.00070.0030−0.0075−0.00040.0044Zero-consumption,y20.00010.0016−0.00250.00000.0036All zeros,y2=0−0.00020.0005−0.0011−0.00020.0007y2=10.00010.0002−0.00030.00010.0004y2=20.00010.0001−0.00020.00010.0003MaleNonparticipation−0.02540.0599−0.1268−0.03050.1012Zero-consumption,y1−0.35950.0611−0.4900−0.3540−0.2565All zeros,y1=0−0.38490.0116−0.4078−0.3849−0.3618y1=10.06300.00400.05550.06300.0711y1=20.15600.00650.14350.15590.1689y1=30.16590.00830.15030.16570.1829Zero-consumption,y20.10120.0623−0.03090.10640.2075All zeros,y2=00.07580.01260.05110.07590.1004y2=10.05010.00330.04380.05000.0567y2=2−0.12580.0112−0.1478−0.1258−0.1040MarriedNonparticipation−0.06800.0777−0.2274−0.04330.0346Zero-consumption,y10.02000.0705−0.0778−0.00010.1692All zeros,y1=0−0.04800.0149−0.0796−0.0472−0.0207y1=10.00560.00350.00060.00470.0132y1=20.01610.00610.00600.01540.0296y1=30.02630.00730.01160.02640.0406Zero-consumption,y20.07090.0791−0.03710.04740.2349All zeros,y2=00.00280.0119−0.02000.00260.0269y2=10.06280.00320.05660.06270.0693y2=2−0.06560.0115−0.0888−0.0654−0.0434MuslimNonparticipation0.03930.0243−0.00500.03840.0900Zero-consumption,y1−0.02390.0247−0.0752−0.02310.0216All zeros,y1=00.01540.0090−0.00160.01530.0334y1=1−0.00230.0011−0.0044−0.0022−0.0002y1=2−0.00530.0027−0.0106−0.0053−0.0001y1=3−0.00780.0060−0.0200−0.00770.0036Zero-consumption,y2−0.02600.0258−0.0797−0.02530.0222All zeros,y2=00.01330.0092−0.00460.01330.0315y2=10.06130.00300.05540.06130.0674y2=2−0.07460.0091−0.0926−0.0746−0.0569Father useNonparticipation0.01220.0187−0.02510.01240.0487Zero-consumption,y1−0.01020.0158−0.0411−0.01040.0214All zeros,y1=00.00200.0030−0.00400.00190.0082y1=1−0.00050.0008−0.0022−0.00050.0011y1=2−0.00090.0014−0.0037−0.00090.0018y1=3−0.00050.0008−0.0023−0.00050.0011Zero-consumption,y2−0.01160.0179−0.0464−0.01180.0240All zeros,y2=00.00060.0011−0.00120.00030.0033y2=1−0.00030.0006−0.0019−0.00020.0007y2=2−0.00020.0005−0.0014−0.00010.0005Mother useNonparticipation0.01290.0257−0.03430.01230.0634Zero-consumption,y1−0.01060.0215−0.0527−0.01030.0298All zeros,y1=00.00240.0043−0.00470.00200.0115y1=1−0.00060.0012−0.0031−0.00060.0014y1=2−0.00110.0019−0.0051−0.00090.0022y1=3−0.00070.0012−0.0033−0.00050.0012Zero-consumption,y2−0.01190.0242−0.0587−0.01180.0338All zeros,y2=00.00100.0016−0.00070.00040.0053y2=1−0.00060.0009−0.0030−0.00020.0005y2=2−0.00040.0007−0.0023−0.00010.0002RegionNonparticipation−0.04800.0240−0.0963−0.0470−0.0040Zero-consumption,y10.04120.0237−0.00390.04060.0889All zeros,y1=0−0.00680.0079−0.0222−0.00680.0086y1=10.00210.00110.00010.00210.0046y1=20.00330.0025−0.00160.00330.0083y1=30.00130.0052−0.00870.00140.0114Zero-consumption,y2−0.02060.0252−0.0672−0.02170.0301All zeros,y2=0−0.06860.0078−0.0840−0.0686−0.0533y2=10.07560.00380.06820.07550.0832y2=2−0.00700.0070−0.0207−0.00700.0072UrbanNonparticipation−0.00620.0261−0.0595−0.00540.0428Zero-consumption,y10.02170.0258−0.02710.02110.0733All zeros,y1=00.01550.0088−0.00180.01550.0324y1=1−0.00070.0012−0.0029−0.00080.0017y1=2−0.00420.0028−0.0096−0.00420.0014y1=3−0.01060.0056−0.0215−0.01060.0006Zero-consumption,y20.01810.0275−0.03370.01780.0739All zeros,y2=00.01190.0090−0.00620.01200.0295y2=10.05970.00360.05280.05970.0668y2=2−0.07160.0075−0.0864−0.0717−0.0566AgriserviceNonparticipation−0.19890.0521−0.3092−0.1960−0.1062Zero-consumption,y10.21020.05060.12020.20750.3161All zeros,y1=00.01130.0098−0.00840.01150.0297y1=10.00580.00180.00260.00570.0097y1=20.00230.0033−0.00390.00210.0092y1=3−0.01940.0060−0.0311−0.0194−0.0077Zero-consumption,y20.18380.05300.08710.18110.2940All zeros,y2=0−0.01510.0126−0.0400−0.01500.0091y2=10.06800.00490.05880.06780.0782y2=2−0.05290.0096−0.0716−0.0530−0.0338Self-employedNonparticipation−0.12870.0693−0.2542−0.1191−0.0122Zero-consumption,y10.15900.06860.04310.15080.2845All zeros,y1=00.03030.0166−0.00340.03050.0627y1=10.00050.0025−0.00420.00050.0058y1=2−0.00750.0060−0.0192−0.00750.0043y1=3−0.02330.0089−0.0398−0.0237−0.0046Zero-consumption,y20.10340.0704−0.01790.09410.2327All zeros, y2=0−0.02540.0147−0.0546−0.02510.0035y2=10.06840.00470.05940.06810.0781y2=2−0.04300.0118−0.0660−0.0431−0.0195StudentNonparticipation0.03050.0357−0.03120.02700.1076Zero-consumption,y10.05480.0434−0.03530.05640.1354All zeros,y1=00.08520.02060.04370.08550.1247y1=1−0.00900.0027−0.0149−0.0089−0.0041y1=2−0.02950.0079−0.0455−0.0294−0.0143y1=3−0.04680.0106−0.0657−0.0475−0.0244Zero-consumption,y20.02840.0448−0.06860.03130.1073All zeros,y2=00.05880.02390.00650.06100.0995y2=10.03900.01020.02070.03830.0604y2=2−0.09790.0142−0.1211−0.0994−0.0659Table 6
Posterior mean, standard deviation and 95% credible intervals of parameters from BOP for smoking and chewing tobacco.
VariableMeanSt. Dev.2.50%Median97.50%Smoking (y1) equation, (β1,α1)Age/101.0290.0950.8281.030Age square/100−0.1040.010−0.123−0.105Education−0.0780.014−0.105−0.078Education square0.0020.0010.0000.002Income0.0000.002−0.0040.000Male2.0660.0911.8882.067Married0.2210.0640.0930.220Muslim−0.0830.049−0.177−0.083Region0.0410.043−0.0440.041Urban−0.0910.048−0.186−0.091Agriservice−0.1210.050−0.219−0.122Self-employed−0.1490.087−0.318−0.150Sstudent−0.7200.093−0.905−0.719α120.2700.0150.2410.270α130.9560.0280.9010.956Chewing (y2) equation, (β2,α2)Age/100.7970.0910.6090.801Age square/100−0.0590.010−0.079−0.059Education−0.0230.016−0.055−0.023Education square−0.0020.001−0.005−0.002Income0.0020.003−0.0040.002Male−0.4410.074−0.586−0.441Married−0.0100.073−0.153−0.011Muslim−0.0770.056−0.187−0.077Region0.4300.0490.3340.430Urban−0.0820.056−0.193−0.081Agriservice0.0780.073−0.0670.078Self employed0.1770.0870.0100.176Student−0.7150.177−1.070−0.710α220.4800.0230.4360.480ρ12−0.1780.034−0.244−0.179Each equation includes father use and mother use variables as well as a constant term.Table 7
Posterior mean, standard deviation, and 95% credible intervals of marginal effects of covariates on probability of smoking and chewing tobacco (BOP model).
VariableProbabilityMeanSt. dev.2.50%Median97.50%AgeAll zeros,y1=00.03680.00380.02880.03690.0438y1=1−0.00040.0003−0.0009−0.00040.0002y1=2−0.00730.0010−0.0093−0.0073−0.0054y1=3−0.02920.0029−0.0345−0.0292−0.0230All zeros,y2=00.02130.00430.01250.02150.0298y2=1−0.00560.0013−0.0082−0.0057−0.0031y2=2−0.00300.0008−0.0047−0.0030−0.0014EducationAll zeros,y1=0−0.03420.0236−0.0803−0.03400.0126y1=10.00380.0025−0.00110.00390.0086y1=20.01300.0084−0.00390.01290.0293y1=30.01740.0128−0.00760.01720.0428All zeros,y2=00.03220.01560.00020.03260.0616y2=1−0.01500.0077−0.0296−0.01510.0009y2=2−0.02010.0113−0.0418−0.02030.0024IncomeAll zeros,y1=0−0.00010.0004−0.0009−0.00010.0007y1=10.00000.0000−0.00010.00000.0001y1=20.00000.0001−0.00020.00000.0002y1=30.00000.0002−0.00040.00000.0005All zeros,y2=0−0.00030.0005−0.0012−0.00030.0007y2=10.00010.0002−0.00020.00010.0004y2=20.00010.0002−0.00020.00010.0004MaleAll zeros,y1=0−0.38240.0121−0.4064−0.3826−0.3586y1=10.06410.00400.05670.06400.0722y1=20.15400.00650.14160.15400.1667y1=30.16430.00830.14870.16410.1807All zeros,y2=00.07210.01230.04810.07210.0962y2=10.05000.00320.04380.05000.0565y2=2−0.12220.0108−0.1430−0.1220−0.1014MarriedAll zeros,y1=0−0.04160.0124−0.0666−0.0415−0.0174y1=10.00390.00130.00150.00380.0067y1=20.01310.00420.00530.01300.0218y1=30.02460.00700.01060.02470.0385All zeros,y2=00.00180.0118−0.02070.00180.0254y2=10.06220.00310.05630.06220.0685y2=2−0.06400.0114−0.0873−0.0640−0.0420MuslimAll zeros,y1=00.01540.0092−0.00290.01540.0331y1=1−0.00130.0008−0.0028−0.00130.0002y1=2−0.00440.0026−0.0093−0.00440.0008y1=3−0.00970.0058−0.0211−0.00970.0018All zeros,y2=00.01260.0093−0.00510.01250.0313y2=10.06130.00310.05550.06130.0675y2=2−0.07390.0091−0.0922−0.0739−0.0562Father useAll zeros,y1=00.76040.00420.75210.76040.7684y1=10.04770.00270.04260.04770.0531y1=20.09820.00350.09150.09820.1051y1=30.09370.00320.08740.09360.1000All zeros,y2=00.87130.00390.86350.87130.8789y2=10.06230.00300.05660.06230.0684y2=20.06640.00300.06070.06640.0724Mother useAll zeros,y1=00.76040.00420.75210.76040.7684y1=10.04770.00270.04260.04770.0531y1=20.09820.00350.09150.09820.1051y1=30.09370.00320.08740.09360.1000All zeros,y2=00.87130.00390.86350.87130.8789y2=10.06230.00300.05660.06230.0684y2=20.06640.00300.06070.06640.0724RegionAll zeros,y1=0−0.00750.0079−0.0229−0.00750.0080y1=10.00060.0007−0.00070.00060.0020y1=20.00220.0023−0.00230.00220.0067y1=30.00470.0049−0.00500.00470.0144All zeros,y2=0−0.06910.0078−0.0846−0.0691−0.0539y2=10.07560.00380.06840.07550.0832y2=2−0.00650.0070−0.0200−0.00650.0072UrbanAll zeros,y1=00.01670.0087−0.00030.01670.0339y1=1−0.00140.0008−0.0030−0.00140.0000y1=2−0.00490.0026−0.0100−0.00490.0001y1=3−0.01040.0054−0.0210−0.01040.0002All zeros,y2=00.01300.0088−0.00410.01290.0303y2=10.05920.00360.05240.05910.0664y2=2−0.07210.0074−0.0866−0.0721−0.0576AgriserviceAll zeros,y1=00.02180.00880.00430.02190.0390y1=1−0.00180.0007−0.0032−0.0018−0.0004y1=2−0.00620.0025−0.0110−0.0062−0.0013y1=3−0.01380.0057−0.0250−0.0139−0.0027All zeros,y2=0−0.01270.0119−0.0366−0.01270.0106y2=10.06560.00430.05720.06550.0742y2=2−0.05280.0094−0.0711−0.0529−0.0338Self employedAll zeros,y1=00.02770.0162−0.00390.02770.0592y1=1−0.00280.0018−0.0065−0.00270.0003y1=2−0.00870.0053−0.0194−0.00860.0012y1=3−0.01630.0093−0.0335−0.01650.0025All zeros,y2=0−0.02900.0144−0.0578−0.0286−0.0017y2=10.06860.00460.06000.06850.0779y2=2−0.03960.0116−0.0617−0.0398−0.0162StudentAll zeros,y1=00.12870.01550.09800.12860.1588y1=1−0.01730.0030−0.0235−0.0171−0.0118y1=2−0.04750.0069−0.0614−0.0475−0.0343y1=3−0.06390.0063−0.0758−0.0640−0.0510All zeros,y2=00.08550.01510.05310.08660.1117y2=10.02780.00710.01540.02740.0428y2=2−0.11330.0088−0.1288−0.1140−0.0944Using (2.13), we decompose the marginal effect on probability of observing zero-consumption into two components: the effect on nonparticipation (zero inflation) and zero-consumption.  For each explanatory variable, this decomposition is shown in Table 5 in the first three rows for smoking and in rows 1, 7, and 8 for chewing tobacco.  For most variables, the effects on probabilities of nonparticipation and zero-consumption are on average opposite in sign, but this difference seems to diminish at the upper tail of the distribution. For example, looking at the posterior mean for age under smoking, getting older by one more year decreases probability of nonparticipation by about 2.6% but increases probability of zero-consumption by 4.6%, implying a net increase of 2.0% in  predicted probability of observing zero. The effect of age in the case of chewing tobacco is qualitatively similar, negative effect on genuine nonusers and positive effect on potential tobacco users, with the latter dominating in the overall effect.Income has opposite effects on probability of nonparticipation and zero-consumption, predicting on average that tobacco is an inferior good for nonparticipants and a normal good for participants. However, the 95% credible interval contains zero, suggesting that the effect of income is weak.  Generally, the opposing effects on probabilities of nonparticipation and zeroconsumption would have repercussions on both the magnitude and the statistical significance of the full effect of observing zero-consumption. Similar considerations apply to positive levels of consumption since the marginal effect on probability of observing consumption levelj(j=1,2,…) can be decomposed into the marginal effects on (i) participation P(si=1) and (ii) levels of consumption conditional on participation, P(yri=j∣si=1). These results show that policy recommendations that ignore excess zeros may lead to misleading conclusions.


## 4.1. Data

We consider an application to tobacco consumption behavior of individuals based on the 2001 household Tobacco Prevalence survey data from Bangladesh. The Survey was conducted in two administrative districts of paramount interest for tobacco production and consumption in the country. Data on daily consumption of smoking and chewing tobacco along with other socioeconomic and demographic characteristics and parental tobacco consumption habits were collected from respondents of 10 years of age and above. The data set has been used previously by Gurmu and Yunus [26] in the context of binary response models. Here we focus on a sample consisting of 6000 individual respondents aged between 10 and 101 years.The ordinal outcomesyr=0,1,2,3  used in this paper correspond roughly to zero, low, moderate, and high levels of tobacco consumption in the form of smoking (y1) or chewing tobacco (y2), respectively. The first dependent variable y1 for an individual's daily cigarette smoking intensities assumes the following 4 choices: y1=0 if nonsmoker, y1=1 if smoking up to 7 cigarettes per day, y=2 if smoking between 8 and 12 cigarettes daily, and y1=3 if smoking more than 12 cigarettes daily; likewise, for the intensity of chewing tobacco, y2=0 if reported not chewing tobacco, y2=1 if uses up to 7 chewing tobacco, and y2=2 if consuming 7 or more chewing tobacco. The frequency distribution of  cigarette smoking  and tobacco chewing choices in Table 1 shows that nearly 66% of the respondents identify themselves as nonusers of tobacco. Our modeling strategy recognizes that these self-identified current nonusers of tobacco may include either individuals who never smoke or chew tobacco (genuine nonusers) or those who do, but not during the reporting period (potential users of tobacco). For example, potential tobacco users may include those who wrongly claim to be nonusers, previous  tobacco users that are currently nonusers, and those most likely to use tobacco in the future due to changes in, say, prices and income. Table 1 also shows that 76% of the respondents are non-smokers and nearly 87% identify themselves as nonusers of tobacco for chewing. Given the extremely high proportion of observed zeros coupled with sparse cells on the right tail, we employ the zero-inflated bivariate ordered probit framework.Table 1
Bivariate frequency distribution for intensity of tobacco use.
Smoke groupChew groupTotal (N)01203931302324455712651262832526353759834982935562Total (N)52203784026000Table2 gives definition of the explanatory variables as well as their means and standard deviations. The respondents are more likely to be Muslim, married, in early thirties, live in rural area, and have about 7 years of formal schooling. Although the country is mostly agrarian, only around 11% of the respondents were related to agricultural occupation in either doing agricultural operations on their own farms or working as agricultural wage laborers. About 12% of the respondents belong to the service occupation. The benchmark occupational group consists of business and other occupations. More than one-half of the fathers and slightly less than two-thirds of the mothers of the respondents currently use or have used tobacco products in the past.Table 2
Definition and summary statistics for independent variables.
NameDefinitionMeanbSt. Dev.AgeaAge in years30.35(14.9)EducationaNumber of years of formal schooling6.83(4.7)IncomeMonthly family income in 1000s of Taka7.57(10.3)Male= 1 if male54.6Married= 1 if married57.2Muslim= 1 if religion is Islam78.4Father use= 1 if father uses tobacco54.0Mother use= 1 if mother uses tobacco65.1Region= 1 if Rangpur resident,= 0 if Chittagong resident49.7Urban= 1 if  urban resident38.0Agriservice= 1 if agriculture labor  or service occupation23.2Self-employed= 1 if self-employed or household chores30.7Student= 1 if student26.8Other= 1 if business or other occupations (control)19.3aIn implementation, we also include age squared and education squared.bThe means for binary variables are in percentage.Among the variables given in Table2, the two indicators of parental use of tobacco products are included in z  as part of the participation equation (2.4). The rest of the variables are included in xr and z of (2.1) and (2.4). To allow for nonlinear effects, age and education enter all three equations using a quadratic form. Due to lack of data on prices, our analysis is limited to the study of other economic and demographic determinants of participation, smoking, and chewing tobacco.


## 4.2. Results

We estimate the standard bivariate ordered probit (BOP) and zero-inflated bivariate ordered probit regression models for smoking and chewing tobacco and report estimation results for parameters, marginal effects, and choice probabilities, along measures of model selection. An earlier version of this paper reports results from the standard ordered probit model as well as the uncorrelated and correlated versions of the univariate zero-inflated ordered probit model for smoking tobacco. Convergence of the generated samples is assessed using standard tools (such as trace plots and ACF plots) within WinBUGS software. After initial 10,000 burn-in iterations, every 10th MCMC sample thereafter was retained from the next 100,000 iterations, obtaining 10,000 samples for subsequent posterior inference of the unknown parameters. The slowest convergence is observed for some parameters in the inflation submodel. By contrast, the autocorrelations functions for most of the marginal  effects die out quickly relative to those for the associated parameters.Table3 reports the goodness-of-fit statistics for the standard bivariate ordered probit model and its zero-inflated version, ZIBOP. The ZIBOP regression model clearly dominates BOP in terms of DIC and its components; compare the DIC of 11330 for the former and 11447 for the latter model. Table 4 gives posterior means, standard deviations, medians, and the 95 percent credible intervals (in terms of the 2.5 and 97.5 percentiles) of the parameters and choice probabilities from ZIBOP model. For comparison, the corresponding results from BOP are shown in Table 6 of the appendix. Both models predict significant negative correlation between the likelihood of smoking and chewing tobacco. The posterior estimates of the cut-off points are qualitatively similar across models. In what follows, we focus on discussion of results from the preferred ZIBOP model. The 95% credible interval for the correlation parameter ρ12  from the zero-inflated model is in the range −0.25 to −0.12, indicating that smoking and chewing tobacco are generally substitutes. Results of selected predicted choice probabilities (bottom of Table 4) show that the ZIBOP regression model provides very good fit to the data. The posterior mean for the probability of (zero, zero)-inflation is about 24% while the 95% credible interval is [0.15, 0.32], indicating that a substantial proportion of zeros may be attributed to nonparticipants. These results underscore the importance of modeling excess zeros in bivariate ordered probit models.Table 3
Goodness-of-fit statistics via DIC.
ModelDbarDhatpDDICBivariate ordered probit (BOP)11417.111386.930.111447.2Zero-inflated BOP11301.111270.329.811329.9Dbar:  Posterior mean of deviance, Dhat:  Deviance evaluated at the posterior mean of the parameters, pD: Dbar-Dhat, the effective number of parameters, and DIC: Deviance information criterion.Table 4
Posterior mean, standard deviation, and 95% credible intervals of parameters from zibop for smoking and chewing tobacco.
VariableMeanSt. dev.2.50%Median97.50%Main (β1,α1): smoking (y1):Age/100.6720.1190.4440.685Age square/100−0.0700.012−0.093−0.071Education−0.0710.014−0.097−0.071Education square0.0010.001−0.0020.001Income0.0000.002−0.0050.000Male2.0920.0861.9252.091Married0.2130.0700.0740.213Muslim−0.0530.052−0.157−0.053Region−0.0070.048−0.102−0.007Urban−0.0960.051−0.198−0.097Agriservice−0.2340.056−0.345−0.233Self-employed−0.2460.087−0.414−0.247student−0.4760.137−0.742−0.478α120.2840.0170.2520.283α130.9870.0300.9280.987Main (β2,α2): chewing (y2)Age/100.6490.1330.3820.658Age square/100−0.0460.013−0.071−0.046Education−0.0200.016−0.052−0.020Education square−0.0020.001−0.005−0.002Income0.0010.003−0.0040.002Male−0.4790.081−0.641−0.479Married−0.0250.075−0.171−0.025Muslim−0.0720.056−0.181−0.072Region0.4170.0510.3170.418Urban−0.0800.058−0.194−0.079Agriservice0.0520.074−0.0960.052Self-employed0.1270.092−0.0580.126Student−0.4500.221−0.887−0.448α220.4840.0230.4390.484Inflation (γ):Age/10−0.0122.044−4.7550.253Age square/1000.5090.552−0.1970.398Education−0.2180.115−0.476−0.204Education square0.0280.0110.0100.026Income0.0060.022−0.0270.003Male0.2390.827−1.5820.417Married2.3064.478−0.4160.500Muslim−0.5280.356−1.331−0.494Mother−0.1700.267−0.716−0.164Father−0.1190.330−0.664−0.160Region0.6300.2910.0610.625Urban0.0400.357−0.7370.071Agrservice5.3125.4161.0172.674Self-employed3.7835.0250.1241.275Sstudent−0.3440.411−1.154−0.339ρ12−0.1850.033−0.249−0.186Select probabilities:P (y1=0)0.7600.0040.7520.760P (y2=0)0.8710.0040.8640.871P (y1=0, y2=0)0.6620.0050.6520.662P (zero-inflation)0.2420.0480.1510.243Results for the constant terms in the main and inflation parts have  been suppressed for brevity.To facilitate interpretation of results, we report in Tables5 and 7 the same set of posterior estimates for the marginal effects from ZIBOP and BOP models, respectively. Since age and education enter the three equations non-linearly, we report the total marginal effects coming from the linear and quadratic parts. We examine closely the marginal effects on the unconditional marginal probabilities at all levels of smoking and chewing tobacco (y1=0,1,2,3; y2=0,1,2). The marginal effects reported in Table 5 show that the results for covariates are generally plausible. Age has a negative impact on probabilities of moderate and heavy use of tobacco. For heavy smokers, education has a significant negative impact on the probability of smoking cigarettes. An additional year of schooling on average decreases probability of smoking by about 6.9% for heavy smokers. Among participants, being male or married has positive impact on probability of smoking, while the effects for being Muslim, urban resident, and student are largely negative. Male respondents are more likely to smoke cigarettes while women respondents are more likely to use chewing tobacco with heavy intensity, a result which is in line with custom of the country [26].Table 5
Posterior mean, standard deviation, and 95% credible intervals of marginal effects of covariates on probability of smoking and chewing tobacco (ZIBOP model).
VariableProbabilityMeanSt. dev.2.50%Median97.50%AgeNonparticipation−0.02590.0129−0.0556−0.0236−0.0078Zero-consumption,y10.04630.01020.02940.04530.0687All zeros,y1=00.02040.00590.00780.02130.0304y1=10.00580.00350.00090.00530.0138y1=2−0.00140.0029−0.0057−0.00190.0055y1=3−0.06900.0235−0.1223−0.0658−0.0344Zero-consumption,y20.04030.01160.01950.03860.0675All zeros,y2=00.01450.00640.00180.01490.0264y2=1−0.00340.0021−0.0071−0.00350.0008y2=2−0.00190.0014−0.0043−0.00200.0011EducationNonparticipation−0.28230.0768−0.4260−0.2837−0.1252Zero-consumption,y10.24470.07490.09170.24590.3851All zeros,y1=0−0.03770.0241−0.0853−0.03740.0094y1=10.04980.01410.02310.04940.0789y1=20.02410.01020.00450.02390.0444y1=3−0.55570.1536−0.8415−0.5588−0.2417Zero-consumption,y20.31360.07720.15610.31590.4546All zeros,y2=00.03130.0161−0.00090.03150.0618y2=1−0.01340.0080−0.0288−0.01350.0027y2=2−0.02220.0119−0.0455−0.02210.0009IncomeNonparticipation−0.00040.0015−0.0038−0.00020.0022Zero-consumption,y10.00030.0014−0.00220.00020.0035All zeros,y1=0−0.00010.0004−0.0009−0.00010.0008y1=10.00010.0003−0.00050.00000.0008y1=20.00000.0002−0.00030.00000.0004y1=3−0.00070.0030−0.0075−0.00040.0044Zero-consumption,y20.00010.0016−0.00250.00000.0036All zeros,y2=0−0.00020.0005−0.0011−0.00020.0007y2=10.00010.0002−0.00030.00010.0004y2=20.00010.0001−0.00020.00010.0003MaleNonparticipation−0.02540.0599−0.1268−0.03050.1012Zero-consumption,y1−0.35950.0611−0.4900−0.3540−0.2565All zeros,y1=0−0.38490.0116−0.4078−0.3849−0.3618y1=10.06300.00400.05550.06300.0711y1=20.15600.00650.14350.15590.1689y1=30.16590.00830.15030.16570.1829Zero-consumption,y20.10120.0623−0.03090.10640.2075All zeros,y2=00.07580.01260.05110.07590.1004y2=10.05010.00330.04380.05000.0567y2=2−0.12580.0112−0.1478−0.1258−0.1040MarriedNonparticipation−0.06800.0777−0.2274−0.04330.0346Zero-consumption,y10.02000.0705−0.0778−0.00010.1692All zeros,y1=0−0.04800.0149−0.0796−0.0472−0.0207y1=10.00560.00350.00060.00470.0132y1=20.01610.00610.00600.01540.0296y1=30.02630.00730.01160.02640.0406Zero-consumption,y20.07090.0791−0.03710.04740.2349All zeros,y2=00.00280.0119−0.02000.00260.0269y2=10.06280.00320.05660.06270.0693y2=2−0.06560.0115−0.0888−0.0654−0.0434MuslimNonparticipation0.03930.0243−0.00500.03840.0900Zero-consumption,y1−0.02390.0247−0.0752−0.02310.0216All zeros,y1=00.01540.0090−0.00160.01530.0334y1=1−0.00230.0011−0.0044−0.0022−0.0002y1=2−0.00530.0027−0.0106−0.0053−0.0001y1=3−0.00780.0060−0.0200−0.00770.0036Zero-consumption,y2−0.02600.0258−0.0797−0.02530.0222All zeros,y2=00.01330.0092−0.00460.01330.0315y2=10.06130.00300.05540.06130.0674y2=2−0.07460.0091−0.0926−0.0746−0.0569Father useNonparticipation0.01220.0187−0.02510.01240.0487Zero-consumption,y1−0.01020.0158−0.0411−0.01040.0214All zeros,y1=00.00200.0030−0.00400.00190.0082y1=1−0.00050.0008−0.0022−0.00050.0011y1=2−0.00090.0014−0.0037−0.00090.0018y1=3−0.00050.0008−0.0023−0.00050.0011Zero-consumption,y2−0.01160.0179−0.0464−0.01180.0240All zeros,y2=00.00060.0011−0.00120.00030.0033y2=1−0.00030.0006−0.0019−0.00020.0007y2=2−0.00020.0005−0.0014−0.00010.0005Mother useNonparticipation0.01290.0257−0.03430.01230.0634Zero-consumption,y1−0.01060.0215−0.0527−0.01030.0298All zeros,y1=00.00240.0043−0.00470.00200.0115y1=1−0.00060.0012−0.0031−0.00060.0014y1=2−0.00110.0019−0.0051−0.00090.0022y1=3−0.00070.0012−0.0033−0.00050.0012Zero-consumption,y2−0.01190.0242−0.0587−0.01180.0338All zeros,y2=00.00100.0016−0.00070.00040.0053y2=1−0.00060.0009−0.0030−0.00020.0005y2=2−0.00040.0007−0.0023−0.00010.0002RegionNonparticipation−0.04800.0240−0.0963−0.0470−0.0040Zero-consumption,y10.04120.0237−0.00390.04060.0889All zeros,y1=0−0.00680.0079−0.0222−0.00680.0086y1=10.00210.00110.00010.00210.0046y1=20.00330.0025−0.00160.00330.0083y1=30.00130.0052−0.00870.00140.0114Zero-consumption,y2−0.02060.0252−0.0672−0.02170.0301All zeros,y2=0−0.06860.0078−0.0840−0.0686−0.0533y2=10.07560.00380.06820.07550.0832y2=2−0.00700.0070−0.0207−0.00700.0072UrbanNonparticipation−0.00620.0261−0.0595−0.00540.0428Zero-consumption,y10.02170.0258−0.02710.02110.0733All zeros,y1=00.01550.0088−0.00180.01550.0324y1=1−0.00070.0012−0.0029−0.00080.0017y1=2−0.00420.0028−0.0096−0.00420.0014y1=3−0.01060.0056−0.0215−0.01060.0006Zero-consumption,y20.01810.0275−0.03370.01780.0739All zeros,y2=00.01190.0090−0.00620.01200.0295y2=10.05970.00360.05280.05970.0668y2=2−0.07160.0075−0.0864−0.0717−0.0566AgriserviceNonparticipation−0.19890.0521−0.3092−0.1960−0.1062Zero-consumption,y10.21020.05060.12020.20750.3161All zeros,y1=00.01130.0098−0.00840.01150.0297y1=10.00580.00180.00260.00570.0097y1=20.00230.0033−0.00390.00210.0092y1=3−0.01940.0060−0.0311−0.0194−0.0077Zero-consumption,y20.18380.05300.08710.18110.2940All zeros,y2=0−0.01510.0126−0.0400−0.01500.0091y2=10.06800.00490.05880.06780.0782y2=2−0.05290.0096−0.0716−0.0530−0.0338Self-employedNonparticipation−0.12870.0693−0.2542−0.1191−0.0122Zero-consumption,y10.15900.06860.04310.15080.2845All zeros,y1=00.03030.0166−0.00340.03050.0627y1=10.00050.0025−0.00420.00050.0058y1=2−0.00750.0060−0.0192−0.00750.0043y1=3−0.02330.0089−0.0398−0.0237−0.0046Zero-consumption,y20.10340.0704−0.01790.09410.2327All zeros, y2=0−0.02540.0147−0.0546−0.02510.0035y2=10.06840.00470.05940.06810.0781y2=2−0.04300.0118−0.0660−0.0431−0.0195StudentNonparticipation0.03050.0357−0.03120.02700.1076Zero-consumption,y10.05480.0434−0.03530.05640.1354All zeros,y1=00.08520.02060.04370.08550.1247y1=1−0.00900.0027−0.0149−0.0089−0.0041y1=2−0.02950.0079−0.0455−0.0294−0.0143y1=3−0.04680.0106−0.0657−0.0475−0.0244Zero-consumption,y20.02840.0448−0.06860.03130.1073All zeros,y2=00.05880.02390.00650.06100.0995y2=10.03900.01020.02070.03830.0604y2=2−0.09790.0142−0.1211−0.0994−0.0659Table 6
Posterior mean, standard deviation and 95% credible intervals of parameters from BOP for smoking and chewing tobacco.
VariableMeanSt. Dev.2.50%Median97.50%Smoking (y1) equation, (β1,α1)Age/101.0290.0950.8281.030Age square/100−0.1040.010−0.123−0.105Education−0.0780.014−0.105−0.078Education square0.0020.0010.0000.002Income0.0000.002−0.0040.000Male2.0660.0911.8882.067Married0.2210.0640.0930.220Muslim−0.0830.049−0.177−0.083Region0.0410.043−0.0440.041Urban−0.0910.048−0.186−0.091Agriservice−0.1210.050−0.219−0.122Self-employed−0.1490.087−0.318−0.150Sstudent−0.7200.093−0.905−0.719α120.2700.0150.2410.270α130.9560.0280.9010.956Chewing (y2) equation, (β2,α2)Age/100.7970.0910.6090.801Age square/100−0.0590.010−0.079−0.059Education−0.0230.016−0.055−0.023Education square−0.0020.001−0.005−0.002Income0.0020.003−0.0040.002Male−0.4410.074−0.586−0.441Married−0.0100.073−0.153−0.011Muslim−0.0770.056−0.187−0.077Region0.4300.0490.3340.430Urban−0.0820.056−0.193−0.081Agriservice0.0780.073−0.0670.078Self employed0.1770.0870.0100.176Student−0.7150.177−1.070−0.710α220.4800.0230.4360.480ρ12−0.1780.034−0.244−0.179Each equation includes father use and mother use variables as well as a constant term.Table 7
Posterior mean, standard deviation, and 95% credible intervals of marginal effects of covariates on probability of smoking and chewing tobacco (BOP model).
VariableProbabilityMeanSt. dev.2.50%Median97.50%AgeAll zeros,y1=00.03680.00380.02880.03690.0438y1=1−0.00040.0003−0.0009−0.00040.0002y1=2−0.00730.0010−0.0093−0.0073−0.0054y1=3−0.02920.0029−0.0345−0.0292−0.0230All zeros,y2=00.02130.00430.01250.02150.0298y2=1−0.00560.0013−0.0082−0.0057−0.0031y2=2−0.00300.0008−0.0047−0.0030−0.0014EducationAll zeros,y1=0−0.03420.0236−0.0803−0.03400.0126y1=10.00380.0025−0.00110.00390.0086y1=20.01300.0084−0.00390.01290.0293y1=30.01740.0128−0.00760.01720.0428All zeros,y2=00.03220.01560.00020.03260.0616y2=1−0.01500.0077−0.0296−0.01510.0009y2=2−0.02010.0113−0.0418−0.02030.0024IncomeAll zeros,y1=0−0.00010.0004−0.0009−0.00010.0007y1=10.00000.0000−0.00010.00000.0001y1=20.00000.0001−0.00020.00000.0002y1=30.00000.0002−0.00040.00000.0005All zeros,y2=0−0.00030.0005−0.0012−0.00030.0007y2=10.00010.0002−0.00020.00010.0004y2=20.00010.0002−0.00020.00010.0004MaleAll zeros,y1=0−0.38240.0121−0.4064−0.3826−0.3586y1=10.06410.00400.05670.06400.0722y1=20.15400.00650.14160.15400.1667y1=30.16430.00830.14870.16410.1807All zeros,y2=00.07210.01230.04810.07210.0962y2=10.05000.00320.04380.05000.0565y2=2−0.12220.0108−0.1430−0.1220−0.1014MarriedAll zeros,y1=0−0.04160.0124−0.0666−0.0415−0.0174y1=10.00390.00130.00150.00380.0067y1=20.01310.00420.00530.01300.0218y1=30.02460.00700.01060.02470.0385All zeros,y2=00.00180.0118−0.02070.00180.0254y2=10.06220.00310.05630.06220.0685y2=2−0.06400.0114−0.0873−0.0640−0.0420MuslimAll zeros,y1=00.01540.0092−0.00290.01540.0331y1=1−0.00130.0008−0.0028−0.00130.0002y1=2−0.00440.0026−0.0093−0.00440.0008y1=3−0.00970.0058−0.0211−0.00970.0018All zeros,y2=00.01260.0093−0.00510.01250.0313y2=10.06130.00310.05550.06130.0675y2=2−0.07390.0091−0.0922−0.0739−0.0562Father useAll zeros,y1=00.76040.00420.75210.76040.7684y1=10.04770.00270.04260.04770.0531y1=20.09820.00350.09150.09820.1051y1=30.09370.00320.08740.09360.1000All zeros,y2=00.87130.00390.86350.87130.8789y2=10.06230.00300.05660.06230.0684y2=20.06640.00300.06070.06640.0724Mother useAll zeros,y1=00.76040.00420.75210.76040.7684y1=10.04770.00270.04260.04770.0531y1=20.09820.00350.09150.09820.1051y1=30.09370.00320.08740.09360.1000All zeros,y2=00.87130.00390.86350.87130.8789y2=10.06230.00300.05660.06230.0684y2=20.06640.00300.06070.06640.0724RegionAll zeros,y1=0−0.00750.0079−0.0229−0.00750.0080y1=10.00060.0007−0.00070.00060.0020y1=20.00220.0023−0.00230.00220.0067y1=30.00470.0049−0.00500.00470.0144All zeros,y2=0−0.06910.0078−0.0846−0.0691−0.0539y2=10.07560.00380.06840.07550.0832y2=2−0.00650.0070−0.0200−0.00650.0072UrbanAll zeros,y1=00.01670.0087−0.00030.01670.0339y1=1−0.00140.0008−0.0030−0.00140.0000y1=2−0.00490.0026−0.0100−0.00490.0001y1=3−0.01040.0054−0.0210−0.01040.0002All zeros,y2=00.01300.0088−0.00410.01290.0303y2=10.05920.00360.05240.05910.0664y2=2−0.07210.0074−0.0866−0.0721−0.0576AgriserviceAll zeros,y1=00.02180.00880.00430.02190.0390y1=1−0.00180.0007−0.0032−0.0018−0.0004y1=2−0.00620.0025−0.0110−0.0062−0.0013y1=3−0.01380.0057−0.0250−0.0139−0.0027All zeros,y2=0−0.01270.0119−0.0366−0.01270.0106y2=10.06560.00430.05720.06550.0742y2=2−0.05280.0094−0.0711−0.0529−0.0338Self employedAll zeros,y1=00.02770.0162−0.00390.02770.0592y1=1−0.00280.0018−0.0065−0.00270.0003y1=2−0.00870.0053−0.0194−0.00860.0012y1=3−0.01630.0093−0.0335−0.01650.0025All zeros,y2=0−0.02900.0144−0.0578−0.0286−0.0017y2=10.06860.00460.06000.06850.0779y2=2−0.03960.0116−0.0617−0.0398−0.0162StudentAll zeros,y1=00.12870.01550.09800.12860.1588y1=1−0.01730.0030−0.0235−0.0171−0.0118y1=2−0.04750.0069−0.0614−0.0475−0.0343y1=3−0.06390.0063−0.0758−0.0640−0.0510All zeros,y2=00.08550.01510.05310.08660.1117y2=10.02780.00710.01540.02740.0428y2=2−0.11330.0088−0.1288−0.1140−0.0944Using (2.13), we decompose the marginal effect on probability of observing zero-consumption into two components: the effect on nonparticipation (zero inflation) and zero-consumption.  For each explanatory variable, this decomposition is shown in Table 5 in the first three rows for smoking and in rows 1, 7, and 8 for chewing tobacco.  For most variables, the effects on probabilities of nonparticipation and zero-consumption are on average opposite in sign, but this difference seems to diminish at the upper tail of the distribution. For example, looking at the posterior mean for age under smoking, getting older by one more year decreases probability of nonparticipation by about 2.6% but increases probability of zero-consumption by 4.6%, implying a net increase of 2.0% in  predicted probability of observing zero. The effect of age in the case of chewing tobacco is qualitatively similar, negative effect on genuine nonusers and positive effect on potential tobacco users, with the latter dominating in the overall effect.Income has opposite effects on probability of nonparticipation and zero-consumption, predicting on average that tobacco is an inferior good for nonparticipants and a normal good for participants. However, the 95% credible interval contains zero, suggesting that the effect of income is weak.  Generally, the opposing effects on probabilities of nonparticipation and zeroconsumption would have repercussions on both the magnitude and the statistical significance of the full effect of observing zero-consumption. Similar considerations apply to positive levels of consumption since the marginal effect on probability of observing consumption levelj(j=1,2,…) can be decomposed into the marginal effects on (i) participation P(si=1) and (ii) levels of consumption conditional on participation, P(yri=j∣si=1). These results show that policy recommendations that ignore excess zeros may lead to misleading conclusions.


## 5. Conclusion

In this paper we analyze the zero-inflated bivariate ordered probit model in a Bayesian framework. The underlying model arises as a mixture of a point mass distribution at(0,0) for nonparticipants and the bivariate ordered probit distribution for participants. The Bayesian analysis is carried out using MCMC techniques to approximate the posterior distribution of the parameters. Using household tobacco survey data with substantial proportion of zeros, we analyze the socioeconomic determinants of individual problem of smoking and chewing tobacco. In our illustration, we find evidence that accounting for excess zeros provides very good fit to the data. The use of a model that ignores zero-inflation masks differential effects of covariates on nonusers and users at various levels of consumption, including zeros. The Bayesian approach to modeling excess zeros provides computational flexibility of generalizing to multivariate ordered response models as well as ordinal panel data models.The proposed zero-inflated bivariate model is particularly useful when most of the bivariate ordered outcomes are zero(y1=0,y2=0). In addition to allowing for inflation in the double-zero state, our approach can be extended to allow for zero inflation in each component. If needed, other states in an ordered regression model may be inflated as well. These extensions need to be justified empirically on a case-by-case basis and are beyond the scope of this paper.



---

*Source: 617678-2011-12-12.xml*