Research Article  Open Access
Jianxu Liu, Duangthip Sirikanchanarak, Songsak Sriboonchitta, Jiachun Xie, "Analysis of Household Consumption Behavior and Indebted SelfSelection Effects: Case Study of Thailand", Mathematical Problems in Engineering, vol. 2018, Article ID 5486185, 12 pages, 2018. https://doi.org/10.1155/2018/5486185
Analysis of Household Consumption Behavior and Indebted SelfSelection Effects: Case Study of Thailand
Abstract
Copula is deemed to be a good approach for relaxing bivariate or multivariate distributions in econometric models. This paper combines static and dynamic copula functions with endogenous switching to study selfselection effects and interdependence between error terms. This technique, copulabased models, is applied to analyze household consumption behavior and indebted selfselection effects in Thailand. The independent, Gaussian, Frank, Clayton, Gumbel, and Joe copula functions and the relatively rotated copula functions were employed in the empirical work. The best model was selected by the information criterion, AIC. We separated the households into four groups, indebted households, debtfree households, households with housing/land loans, and households without housing/land loans, which favors the examination of the treatment effects of indebted households or households with housing debts. The main results indicate that dynamic copulabased models offer better performance than others, such as classical endogenous switching models or all static copulabased models. Also, the I–I and the G–G models underestimate the treatment effects relative to the best models. Additionally, importantly, the traditional normal bivariate distribution or the static copula function could characterize the relationship as regards errors between household debt choice and household consumption and can lead to very misleading implications about the treatment effects.
1. Introduction
With rapid economic development and urbanization, household consumption plays an important role in the economic growth of a nation. Household consumption contributes toward dampening private consumption and, thus, aggregate demand and economic activity. Thailand’s economic growth has been mainly driven by household consumption. From Figure 1, we can see that household consumption expenditure accounts for 58% of Thailand’s GDP (gross domestic product) in terms of final demand and that the same tends to increase after 2008. Likewise, we can see that the growth of GDP and the household consumption expenditure seem to move with each other closely, which is shown in Figure 2. This is because there have been several government policies oriented to improve economy by boosting household spending, such as ricefarming subsidies, business development loans in rural areas, and tax breaks on new car purchases, and all of them have been used after widespread flooding which devastated the country’s economy in 2011. Therefore, this can be taken to imply that household consumption will become a key driving force in Thailand’s economy. Obviously, boosting household consumption is an effective measure to retain economic growth. In 2015, the Thai government launched multiple stimulus measures focusing on the rural population and small and mediumsized enterprises (SMEs). For rural population, the first measure involves giving out no and low interest loans worth a total of THB59 billion. For SMEs, the stimulus measure includes low interest loans worth THB100 billion, credit guarantee of THB100 billion, tax exemptions, and tax cuts, as well as a venture capital fund worth THB6 billion. The purpose of the abovementioned measures is to create shortterm consumption in the form of pushing up household debt, thereby boosting the economic growth. Moreover, household debt could not only help people smoothen their consumption over their lifetimes but also help households build up their physical assets such as land and house. As a result, it is worth considering how household debts quantitatively affect consumption.
In recent years, there have been many studies that have investigated the effects of household debt on household consumption expenditure. For example, Ogawa and Wan [1] examined the effects of household debt on consumption by estimating the consumption functions of total expenditures and of individual expenditure classified by type using Ordinary Least Square (OLS) regression. Some empirical studies suggested that the ratio of the debt associated with land and housing purchases to the market value of tangible asset had significantly negative effects on expenditure after the bubble burst. Moreover, they found that the debtasset ratio had negative effects on expenditures on semidurable and nondurable goods. Ekici and Dunn [2] investigated the relationship between credit card debt and consumption using Tobit model for the US. Meniago et al. [3] used the Vector Error Correction Model (VECM) to estimate the relationship between household debt and other macroeconomic variables. He found the existence of a longrun cointegrating relationship between household debt and other macroeconomic determinants. Increase in household debt is affected by positive changes in consumer price index, gross domestic product, and household consumption, and household prices and household saving were found to positively contribute to a rise in household debt. Ni and Seol [4] examined how nondurable consumption expenditures in households headed by government employees respond to predictable income changes. They found that moderately excess sensitivity in consumption can be attributed to about nine percent of the households. Kukk [5] used the OLS regression to investigate the estimations and showed that indebtedness measured by the debttoincome ratio hampers consumption as a whole during the recession rather than during the precrisis and the postcrisis periods in a new EU member country. Moreover, the findings suggest that household indebtedness is amplifying the recession, and the debt repayment burden indicates the mechanism. Baker [6] studied the consumption response to household income shocks. They found that highly indebted households are more sensitive to these income fluctuations and that a one standard deviation increase in debttoasset ratios increases the elasticity of consumption by approximately 25%.
In general, an indebted household is one with specific assets, much higher household income, high education levels, large household expenditure, and so forth. On the other hand, a debtfree household is one with low assets, no fixed household income, less household expenditure, and so forth. Assume that the indebted household consumption is much higher than the debtfree household consumption. Then, the difference in the consumption between the debtfree household and the indebted household is attributed to households selfselecting themselves into loans based on consumption desires. We argue in all articles that household consumption should not be viewed in isolation from household debt choice. The statistical models ignoring this type of selectivity can induce statistical biases into model estimates. Also, it will indirectly have an effect on the provision of information and policy advice.
To address this drawback, a copulabased endogenous switching model is applied for investigating the effects of household debt on household consumption expenditure and measuring the treatment effect of household debt. Copulabased endogenous switching models allow for better flexibility in joint distributions than the bivariate normal distribution. Also, many copula functions are used for providing a more general dependence structure, as opposed to treating it as a simple linear correlation. This implies that copulabased models can better fit the data than standard endogenous switching models. Smith [7] firstly proposed a copulabased endogenous switching model which was applied to child labor. Thereafter, this model has been widely employed in many empirical analyses. For example, Bhat and Eluru [8] applied the copulabased approach to accommodate residential selfselection effects in travel behavior modeling. Eberth and Smith [9] modeled the participation decision and duration of sporting activity in Scotland. Wichian and Sriboonchitta [10] used it to analyze private and public wage determination for older workers in Thailand.
The contribution of this paper can be listed as follows: first, we propose dynamic copulabased endogenous switching models to measure the dynamic dependence structure between household debt and household consumption in Thailand. Dynamic copulabased endogenous switching models can be used to measure the possibly asymmetric and tail dependence between household debt and household consumption. At the same time, dynamic copulabased models can further explain the dependence structure between household debt and household consumption, which is better than the capability of the static copulabased model that was proposed by [7]. Second, this paper is the first study that considers a nonrandom selection of household debt for household consumption behavior. It favors the distinguishing of the differences between debtfree household consumption and indebted household consumption.
The remainder of this paper is organized as follows: In the next section, we introduce the endogenous switching model and the copula functions. In Section 3, we construct the copulabased endogenous switching model. Section 4 describes the data and the variables. In Section 5, we report the empirical results. Section 6 gives us some conclusions.
2. Background
As previously mentioned, the endogenous switching model and copulas are the two building blocks of the copulabased model. Therefore, we firstly summarize the endogenous switching model in Section 2.1, and then review the copulas in Section 2.2.
2.1. Endogenous Switching Model
By following Bhat and Eluru [8], the endogenous switching model by two regimes can be expressed as where the first two equations represent binary choice model and represents a latent variable, such as or 1, implying that a worker is a public or a private sector employee. and , , denote the dependent variable to be explained in each of the two regimes. and represent the independent variables corresponding to each regime. , , and are the parameter vectors. is the error term of the binary choice model, and and are the error terms of the two regimes, respectively. In general, the standard endogenous switching model assumes that there exists a potential linear correlation between the pairs of error terms: and, and , and that both the pairs are bivariate normal distributions, and a full information maximum likelihood method is used to estimate the endogenous switching model. In addition, there are also some researchers and academicians who extended this model by using nonnormal bivariate distributions, twostep estimation methods, and so forth (see [8]).
In endogenous switching models, measurement of the treatment effect is very interesting. A “treatment effect” is the average causal effect of a binary (01) variable on an outcome variable of scientific or policy interest. The term “treatment effect” originated in medical literature concerned with the causal effects of binary, yesorno “treatments,” such as an experimental drug or a new surgical procedure. Average treatment effect (ATE) and average treatment effect on the treated (ATT) are two common methods of measure for causal effects [11]. ATE and ATT can be expressed as where and are the estimated ATE and ATT, respectively. ATE means the difference between the average outcome () for observations that receive treatment and the average outcome for observations that receive control.
2.2. Copulas
Copulas are popular in highdimensional statistical applications as they allow one to easily model and estimate the distribution of random vectors by estimating margins and dependence separately. There are many parametric copula families available, which usually have parameters that control the strength of dependence. In this study, we focus on the presentation of the bivariate copula, which will be used later. Given a joint distribution of two random variables and , the function can be defined on by the following equation:is a copula, where and are marginal distributions of and , respectively. satisfies the following properties:(1) for all and in ;(2) and for all and in ;(3)For all and ,Fréchet–Hoeffding bounds for copula functions are crucial to select an appropriate copula that covers the sample space between the lower and the upper bounds and that as a parameter of copula approaches the lower (upper) bound of its permissible range. The Fréchet–Hoeffding bounds for copula functions can be expressed as where is the Fréchet–Hoeffding lower bound and is the Fréchet–Hoeffding upper bound (see [12, 13]).
If the random vector has a joint density , it can be expressed as a function of the copula density bywhere and are the marginal densities, , and . The small letter represents the density function of the copula.
The different copula functions have different characteristics. For example, Gaussian copula cannot capture tail dependence; Clayton copula can capture lower tail dependence, while Gumbel copula is used to measure upper tail dependence of random variables. Figure 3 demonstrates several copula contour plots under standard normal distribution. The contour plots are generated based on the value of Kendall’s tau equaling 0.5. These plots show the evidence that different copula functions have different characteristics in terms of tail dependence, and it is necessary that this study employs the copula model to provide a flexible method of constructing multivariate distributions, given the marginal distributions and the dependence structures separately. However, these copulas can only capture positive dependence except the Gaussian and the Frank copulas. Fortunately, these copulas may then be “rotated” and applied again. A thorough review of rotated copulas may be found in [14, 15]. There are three rotated forms: rotated 90 degrees, rotated 180 degrees, and rotated 270 degrees. When rotating the copulas by 180 degrees, one obtains the corresponding survival copulas. Nonrotated Clayton copula can capture lower tail dependence, while survival Clayton copula can measure upper tail dependence. Also, rotation by 90 degrees and 270 degrees allows for the modeling of negative dependence. The distribution functions of copulas rotated by 90 degrees, 180 degrees, and 270 degrees are given as follows:
The copula approaches the Fréchet–Hoeffding lower (upper) bound. Therefore, we can make use of some copula families that cover a large dependence as much as possible, such as the (rotated) Clayton, the (rotated) Gumbel, and the (rotated) Joe copula families, which can represent both negative and positive dependence.
3. CopulaBased Endogenous Switching Model
The classical endogenous switching model assumes that the two pairs of error terms have bivariate normal distribution. Smith [7] proposed to relax this assumption and to model the dependence between and by using Archimedean copulas. In this study, we use many copula families to construct a copulabased endogenous switching model and select the best fit in terms of the Akaike information criterion (AIC). It is different in that we will assume the copula parameter depends on some forcing variables corresponding to some evolution equations. We introduce the full information maximum likelihood approach to estimate the copulabased endogenous switching model in this section.
We assume that the error term has standard normal distribution, while both and have normal distribution, but they have scale parameters and , respectively. Since the marginals are known, the joint distribution of and , , can be constructed by copula , where is the copula parameter for each individual. Similarly, the joint distribution of and , , corresponds to . Following [8], the likelihood function of the copulabased endogenous switching model can be written aswhere , and Different copulas have different characteristics, such as upper tail dependence, lower dependence, positive dependence, and negative dependence. Therefore, it is an excellent choice that we apply more copula families and relevant rotated copulas that can potentially capture the appropriate dependence between two random variables. In this study, we systematically consider several copula families including the Gaussian, Frank, Clayton, Gumbel, and Joe families as well as their relevant rotated versions. The model with the best fitcomplexity tradeoff is selected using the AIC criterion [16]. Other popular copula families, such as the Farlie–Gumbel–Morgenstern (FGM) copula and the Ali–Mikhail–Haq (AMH) copula, have been discarded because they cannot achieve Fréchet bounds and can only accommodate relatively weak dependence between the margins. The ranges of dependence of Kendall’s tau and Spearman’s rho for the FGM copula are and , respectively; the Kendall’s tau of the AMH copula is bounded to the interval , and the range of Spearman’s rho is (see [17, 18]).
Bhat and Eluru [8] showed that standard endogenous switching models may mislead the implications about the treatment effect by comparing them using copulabased models. However, Smith [7] and Bhat and Eluru [8] did not provide any interpretation for dependence parameters. We consider that the parameters vary according to some evolution equation. We propose the following evolution equations for the Gaussian and the other copulas:where or 1 denotes the different regimes and and represent Pearson’s correlation of the Gaussian copula and Kendall’s tau of the other copulas. and are the logistic transformation, which can guarantee the dependence with the interval (−1, 1). is a vector of forcing variables, which can be determined by practical issues. is a constant term, and is the parameters of the forcing variables.
4. Data
The household dataset used in this study is from the Socioeconomic Surveys (SES) in 2013. The sampling was conducted through random proportional sampling, and the dataset contains 21,214 observations. This dataset was gathered by the National Statistical Office workers who visit the sample households every two years.
4.1. Data Descriptive
For the purpose of characterizing the indebted and the debtfree household consumption behavior, we classified all households into four groups: the indebted household, the debtfree household, the household with housing/land loans, and the household without housing/land loans groups. Households that have debts and households that have housing/land debts are expressed by the dummy variables 0 and 1, respectively. The loans of indebted households consist of business, education, consumption, housing, land, and so forth. Table 1 summarizes the descriptive statistics of the main socioeconomic variables that include monthly consumption per household, monthly income per household, monthly expenditure per household, total asset of the household, age of the head of the household, number of household members, and number of children in the families. First, there are 10,961 households with debt and 10,253 households without debt in our total sample. At the same time, the number of households with housing/land debt is 1,633, which represents 15 percent of all households with debt. The leftover is no housing/land debt, which comprises 19,581 households. Second, the head of the household with debt is younger than the head of the household without debt at any quantile level, and the levels of family size, income, consumption, expenditure, and asset are larger in households with debt than in households without debt at any quantile level, which is consistent with the requirements of household indebtedness. In addition, the consumption, income, expenditure, and assets in households with housing/land debt are about double those of households without housing/land debt at any quantile level as well. Upon comparing the characteristics of households between households with total debt and households with land, we can conclude that households with housing/land debt have higher consumption and income as well as assets and liabilities.

4.2. Variables Considered
In this part, we would like to introduce the variables that we used in the estimation. We used the debt decision that takes the value 1 if the household has debt or housing/land debt; otherwise, the value equals 0. The selection effects of indebtedness are conditional on other observed determinants which are the logarithms of the asset, expenditure, number of children, income, and age of the household head. When we consider the consumption equation, the dependent variable is assumed to be the ratio of consumption expenditure to disposable income so as to avoid the heteroscedasticity problem. As socioeconomic variables determine consumption, we select some variables as independent variables, which include the reciprocal of the household disposal, saving rate, number of members, number of workers, age of the household head, square of the age of the household head divided by 100, and number of workers with a stable earning job (earned from business, industry, and/or nonfarming business). Moreover, the debtasset ratio is one of the independent variables for indebted households. Thus, the endogenous switching model in this study can be given by where represents the latent variable of household debt or housing/land debt in each dataset; and are the consumption expenditures of debtfree households and indebted households, respectively; DI is the disposal income; FS is the family size; NW is the number of workers; NSEW is the number of workers with a stable income; and DAR is the debtasset ratio.
5. Empirical Results
In this study, we use six copula functions to model endogenous switching and selfselection effects. This is because of the fact that different copula functions have different characteristics. For example, Gaussian copula cannot capture tail dependence; Clayton copula can capture lower tail dependence, while Gumbel copula is used to measure upper tail dependence of random variables. Therefore, the independent, Gaussian, Clayton, Frank, Gumbel, Joe, rotated Clayton, rotated Gumbel, and rotated Joe copulas are applied to endogenous switching models. The preferred copulabased endogenous switching model is selected based on the AIC. We use the R program to estimate the model, and the starting values are obtained from “glm” and “lm” functions, and the maxLik package with the BFGS algorithm is applied to maximize the loglikelihood function of the R software.
As mentioned above, we firstly estimate all the alternative models and then select the preferred models based on the AIC. Hereafter, the copula functions among the preferable models are transformed into dynamic copula functions, and then the maximum likelihood method is used to estimate the dynamic copulabased endogenous switching models. Moreover, there are 72 models for each dataset. In the first case, the dataset is divided into two groups: indebted households and debtfree households. In the second case, we separate the dataset into households with and without housing/land debt. The results of the AIC are shown in Table 2. According to the values of AIC, we find that the Joe–Frank copulasbased endogenous switching model offers better performance than the others for the first case. The Joe and the Frank copulas are used to capture the dependence structures between the error pairs and , respectively. For the second case, the endogenous switching model based on the rotatedby90degree Joe copula and the Joe copula is the best one among all the alternative models. Therefore, the dynamic Joe–Frank copulasbased endogenous switching model was applied to the first case, while the endogenous switching model based on the dynamic rotatedby90degree Joe copula and Joe copula was used for the second case. It is well known that household income is related to household debt and household consumption. Then, household income may be explained by the relation between the error pairs and .
 
Source. Computation. 
In the following presentation of the empirical results, we focus on the results of the best static copulabased models and the corresponding dynamic copulabased models. All empirical results are presented in Tables 3 and 4. First, the estimation results for the dependence are discussed in Section 5.1. Second, the results of the sample selection equation are shown in Section 5.2. Third, we will discuss the results of the endogenous switching and selfselection effects in Section 5.3. Finally, the treatment effects are reported in Section 5.4.
(a) Estimation of sample selection equation  
 
(b) Estimation of consumption function for households without debt  
 
(c) Estimation of consumption function for indebted households  
 
(d) Estimation of dependence parameters  
 
Note. “” denotes significance at 1% level. Source. Computation. 

5.1. Estimation Results for Dependence
We now turn to show the estimation results of the copulabased endogenous switching model. Table 3 shows all the parameter estimates for the two cases. We can see that the parameter estimates do not have much difference between the static copulabased models and the dynamic copulabased models. Table 3(d) displays the estimated results of the dependence parameters for the two cases. First, the values of AIC in the bottom of Table 3 completely verify that the dynamic copulabased models are much better than the static copulabased models. Second, the coefficients of and for both the cases are significant and are equal to −2.265 and 2.123 in the first case, and −1.99 and 1.467 in the second case, respectively. This implies that the dependence between and is a negative correlation, while and are positively related. The negative dependency for the two cases indicates that a household that has a higher inclination to have debt or housing debt would consume more than an observationally equivalent “random” household if both these households do not have debt or housing debt. In other words, a household that makes the choice not to take loans is likely to consume less than a random household among debtfree households. While the positive dependency between and implies that a household that has a higher probability to obtain bank loans is likely to consume less than a random household in indebted households. Equivalently, if an indebted household decides not to take loans, that household would consume much more than a random household in indebted households. Also, we find that the estimated is highly statistically significant and negative for both cases, which means that when the household income is higher, the negative dependence between and is bigger. We can believe that for households without debt, and for those that have high income and high consumption, there is a better probability of taking loans. In addition, the estimated is highly statistically significant and positive, which indicates that where there is high income and high consumption, the probability of being converted into debtfree households is bigger in indebted households. Finally, Figures 4 and 5 describe the static and dynamic Kendall’s tau for the two cases. They looked more natural, as expected, because they always had a negative or a positive correlation and did not hit the upper bound or the lower bound.
5.2. Estimation Results for Sample Selection Equations
Table 3(a) reports the estimated results of sample selection equations. For the first case, all the parameter estimates are statically significant at 1% level. The estimated parameters of the asset, expenditure, and number of children are positive, which implies that households with much asset or much expenditure or more children prefer to take loans. Perhaps, these households have concrete reasons to need the money. The parameter estimates of the age of the household head and the household income are negative, which indicates that households with elderly household heads or high incomes are not likely to take loans. For the second case, the coefficients of the asset and the expenditure are positive as well, while the coefficient of income is statistically significant and positive, which is opposite to that of the first case. This may be because households with high income can apply for housing or land debt. The parameter of the number of children is insignificant, which illustrates that the number of children does not have an effect on the choice of the housing debt. Upon comparing the second case with the first case, it can be observed that the parameter estimates of the asset and the income in the second case are greater than those in the first case. This implies that household asset and household income, as expected, have a pivotal effect on making decisions regarding housing debt.
5.3. Estimation Results for Consumption Equations
Table 3(b) reports the estimate results of consumption behavior in debtfree households and households without housing debt. As expected, the effects of the reciprocal of the household disposal income, the number of members, and the number of workers having stable jobs with steady income are statistically significant and positive. It is obvious that Table 3(c) shows results similar to those of Table 3(b). While comparing the intercepts of Table 3(b) with Table 3(c), we find that the intercepts of Table 3(b) are greater than the intercepts of Table 3(c). This implies that the household disposal income of debtfree households has a higher effect on household consumption than the household disposal income of indebted households. As for the second case, the household disposal income has effects similar to household income. The parameter estimate of the debtasset ratio in the first case is negative and statistically significant, while it is positive in the second case. It is probably determined by the types of loan. For example, an excessive loan for education or consumer goods means much lower household income, and it causes a reduction in the household consumption. At the same time, a high debtasset ratio with regard to housing debt indicates a substantially high household income which may increase the household consumption.
5.4. Treatment Effects
It is apparent that there are statistically significant indebted selfselection effects, which means that the choice of the household debt affects the household consumption. To measure the magnitude of the selfselection effects, we combine the treatment effect and the bootstrap method in this study. The ATE is estimated as follows:where is the number of bootstrap draws and is the sample size. and are from the estimates of the copulabased models. Let , in this study, be a large number, 1000. Following [8], the ATT can be estimated by where is the number of households with debt or housing debt and and are given as follows:It is noticeable that and do not have a closed form. Thus, numerical integration by Simpson’s rule is used to compute and .
In this study, we calculate the point estimates of the treatment effects, the standard deviations, and the lower bound and the upper bound at 95% confidence interval. The results are presented in Table 4 for the independenceindependence (I–I) copulabased models, Gaussian and Gaussian (G–G) copulabased models, and the two dynamic copulabased models corresponding to the dynamic Frank–Joe (D–F–J) copulabased model for case 1 and the dynamic rotated Joe–Joe (D–RJ–J) copulabased model for case 2. It is not surprising that most of the treatment effects are negative values as expected. First, let us focus on ATE. In the first case, the bestfitting D–F–J model indicates that a randomly picked household will reduce household consumption by 952.82 Baht per month if it is a household with debt relative to a debtfree household. At the same time, a randomly picked household will reduce household consumption by 17,778.58 Baht per month if it is a household with housing debt relative to a household without housing debt, according to the D–RJ–J model. It is thus clear that household housing debt has a much greater impact on household consumption than other debts. Second, the estimates of TT indicate that if a randomly picked household with debt were to be a debtfree household, the household consumption would have to increase by about 4,935.51 Baht per month. If a randomly picked household with housing debt were to be a household without housing debt, the household consumption would have to increase by about 8,134.32 Baht per month. Third, the I–I model ignores the sample selection effect, while the G–G model is confined to bivariate normal assumption and linear correlation. This causes the I–I and the G–G models to underestimate the ATE and the TT relative to the best models.
6. Conclusions
The classical sample selection or endogenous switching model is restricted to a bivariate normal distribution on the error terms. Using bivariate normal distribution means there only exist linear and symmetric correlation between the error terms. These may have an impact on the sample selection effects. This study relaxed this assumption by using copula functions, and the dynamic copula was first proposed to be combined with the endogenous switching model. The dynamic copula is favorable for understanding the proper sense of the dependence structure between error terms. Moreover, the dynamic copulabased model is estimated by full information maximum likelihood method and does not entail any kind of numerical integration or simulation machinery. Thus, it is simple to apply to economic problems and will draw the attention of researchers.
This paper applied the copulabased endogenous switching approaches to model household debt choice and household consumption behavior using the data of 21,214 households in Thailand. The main results can be listed as follows: first and foremost, the dynamic copulabased models, the D–F–J copulabased model and the D–RJ–J copulabased model, have much better performance than the static copulabased models and the G–G copulabased models. On the one hand, this implies that the static dependence between error terms cannot accurately reflect selfselection effects. On the other hand, the assumption of bivariate normal distribution also affects selfselection effects. Second, we found that there exists negative dependence between debtfree household consumption and household debt choice, while it is positive dependence that exists between indebted household consumption and household debt choice. The negative dependence indicates that a household that makes the choice of not taking loans is likely to consume less than a random household in debtfree households. At the same time, the positive dependency implies that a household that has a higher probability to obtain bank loans is likely to consume less than a random household in indebted households. Also, the dynamic dependence means that if there is high income and high consumption, the probability of being converted into debtfree households is bigger in indebted households. Third, the estimates of ATE and TT of the dynamic copulabased models are very different from the estimates of the G–G copulabased models. The results showed that the classical G–G copulabased models underestimate the treatment effects. Moreover, it was observed that household housing debt has a much greater impact on household consumption than other debts, upon comparing the two cases. Last but not least, the study indicated that use of a traditional normal bivariate distribution or static copula function to characterize the relationship as regards errors between household debt choice and household consumption can lead to very misleading implications about treatment effects.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors are very grateful to Professor Vladik Kreinovich for his comments. The authors wish to thank the Puey Ungphakorn Centre of Excellence in Econometrics and the Bank of Thailand Scholarship for their financial support.
References
 K. Ogawa and J. Wan, “Household debt and consumption: A quantitative analysis based on household micro data for Japan,” Journal of Housing Economics, vol. 16, no. 2, pp. 127–142, 2007. View at: Publisher Site  Google Scholar
 T. Ekici and L. Dunn, “Credit card debt and consumption: Evidence from householdlevel data,” Applied Economics, vol. 42, no. 4, pp. 455–462, 2010. View at: Publisher Site  Google Scholar
 C. Meniago, J. MukuddemPetersen, M. A. Petersen, and I. P. Mongale, “What causes household debt to increase in South Africa?” Economic Modelling, vol. 33, pp. 482–492, 2013. View at: Publisher Site  Google Scholar
 S. Ni and Y. Seol, “New evidence on excess sensitivity of household consumption,” Journal of Monetary Economics, vol. 63, pp. 80–94, 2014. View at: Publisher Site  Google Scholar
 M. Kukk, “How did household indebtness hamper consumption during the recession? Evidence from micro data,” Journal of Comparative Economics, vol. 44, no. 3, pp. 764–786, 2016. View at: Publisher Site  Google Scholar
 S. R. Baker, “Debt and the response to household income shocks: Validation and application of linked financial account data,” Journal of Political Economy, (forthcoming), 2017. View at: Google Scholar
 M. D. Smith, “Using copulas to model endogenous switchings with an application to child labour,” Economic Record, vol. 81, no. 1, pp. S47–S57, 2005. View at: Publisher Site  Google Scholar
 C. R. Bhat and N. Eluru, “A copulabased approach to accommodate residential selfselection effects in travel behavior modeling,” Transportation Research Part B: Methodological, vol. 43, no. 7, pp. 749–765, 2009. View at: Publisher Site  Google Scholar
 B. Eberth and M. D. Smith, “Modelling the participation decision and duration of sporting activity in Scotland,” Economic Modelling, vol. 27, no. 4, pp. 822–834, 2010. View at: Publisher Site  Google Scholar
 A. Wichian and S. Sriboonchitta, “Econometric analysis of private and public wage determination for older workers using a copula and switching regression,” Thai Journal of Mathematics, pp. 111–128, 2014. View at: Google Scholar  MathSciNet
 J. Heckman, J. L. Tobias, and E. J. Vytlacil, “Four parameters of interest in the evaluation of social programs,” Southern Economic Journal, vol. 68, no. 2, pp. 210–223, 2001. View at: Publisher Site  Google Scholar
 R. B. Nelsen, An Introduction to Copulas, Springer, New York, NY, USA, 2006. View at: MathSciNet
 P. K. Trivedi and D. M. Zimmer, “Copula modeling: An introduction for practitioners,” Foundations and Trends in Econometrics, vol. 1, no. 1, pp. 1–111, 2005. View at: Publisher Site  Google Scholar
 A. Wiboonpongse, J. Liu, S. Sriboonchitta, and T. Denoeux, “Modeling dependence between error components of the stochastic frontier model using copula: application to intercrop coffee production in northern Thailand,” International Journal of Approximate Reasoning, vol. 65, pp. 34–44, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 S. Sriboonchitta, J. Liu, A. Wiboonpongse, and T. Denoeux, “A doublecopula stochastic frontier model with dependent error components and correction for sample selection,” International Journal of Approximate Reasoning, vol. 80, pp. 174–184, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 J. Du, H. Li, and Y. He, “The method of solving structural reliability with multiparameter correlation problem,” Mathematical Problems in Engineering, vol. 2017, Article ID 6976301, 12 pages, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 P. Kumar, “Probability distributions and estimation of AliMikhailHaq copula,” Applied Mathematical Sciences, vol. 4, no. 1316, pp. 657–666, 2010. View at: Google Scholar  MathSciNet
 M. D. Smith, “Stochastic frontier models with dependent error components,” Journal of Econometrics, vol. 11, no. 1, pp. 172–192, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Jianxu Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.