Bayesian Estimation of Archimedean Copula-Based SUR Quantile Models

We propose a high-dimensional copula to model the dependence structure of the seemingly unrelated quantile regression. As the conventional model faces with the strong assumption of the multivariate normal distribution and the linear dependence structure, thus, we apply the multivariate exchangeable copula function to relax this assumption. As there are many parameters to be estimated, we consider the Bayesian Markov chain Monte Carlo approach to estimate the parameter interests in the model. Four simulation studies are conducted to assess the performance of our proposed model and Bayesian estimation. Satisfactory results from simulation studies are obtained suggesting the good performance and reliability of the Bayesian method used in our proposed model. The real data analysis is also provided, and the empirical comparison indicates our proposed model outperforms the conventional models in all considered quantile levels.


Introduction
Typically, the seemingly unrelated regression (SUR) model, developed by Zellner [1], is a system of many structural equations, each having its own dependent variable and probably different sets of independent variables. However, these equations may be correlated through their error terms assuming the existence of the multivariate normal distribution. is model has been applied in many works such as White and Hewings [2]; Adelegan [3]; and Frankel and Poonawala [4].
Recently, the performance of this model has been questioned by many scholars, such as Jun and Pinkse [5] and Waldmann and Kneib [6], as the regression lines are fitted through the means of the independent variables leaving aside the outliers which might carry meaningful information. In general, the data may exhibit a heavy tail. Hence, the normal distribution assumption may fail to provide reliable estimated results. To deal with these issues, Koenker and Bassett [7] introduced the quantile regression. Such an extended linear model has several advantages. For example, the parameter estimates are more robust against outliers, and the extensive analysis of the relationship between the dependent and independent variables is uncovered [6]. Recognizing these advantages of the quantile regression approach, some studies made use of the seemingly unrelated quantile regression (SUQR) model (see [5,6]). ey extended the conventional SUR model by allowing the error term of each quantile equation to be correlated at all quantiles but still with the same normal distribution assumption regarding the marginals. However, the assumption of the multivariate normal distribution may not be appropriate for joining the equation errors which actually have an asymmetric distribution. To relax this assumption, the copula functions are employed to capture the nonlinear dependence structure as well as join the asymmetric error of equations in the multivariate SUQR model.
In the literature, the copula-based models have already been introduced by many studies of which findings demonstrated a higher accuracy in parameter estimation. e pioneering work of Wichitaksorn and Choy [8] suggested using the copula function to join the error terms of the linear regression equation and binary choice equation. ey showed that their proposed model is superior to the conventional model in terms of estimation accuracy. Pastpipatkul et al. [9] used the copula to join the errors of the SUR model. Louzada and Ferreira [10] also suggested using the Clayton copula to join the error terms of the bivariate seemingly unrelated regression (SUN) tobit model. ey mentioned that their model has the ability to capture the lower tail dependence of the SUN model. Ivanov et al. [11] extended a copula approach to model the dependence of unobserved multivariate factors in the dynamic factor models. ey confirmed that the copula-based approach is more general and applicable to several factor models. Recently, Zou et al. [12] suggested using the copula to establish the link between the occurrence of wildlife-vehicle collisions and the underreporting probability. ey revealed that the Gaussian copula-based empirical Bayes method is superior to the traditional EB method.
Among various studies using the copula-based models to successfully confirm the role of copula functions in improving the efficiency of the system equation, none of them used the copulas to join the errors of the SUQR model except for Tansuchat et al. [13] who used multivariate elliptical copulas consisting of Gaussian and Student-t to construct the joint equations of the SUQR model. is paper, therefore, attempts to extend the approach of Tansuchat et al. [13] by applying the multivariate Archimedean copulas to join the model errors. Using this copula class has several additional advantages over the elliptical copulas. First, this copula class can model the asymmetric structure of extreme dependence between the errors, which are particularly important in financial modeling. Second, it provides more flexibility to model the dependence structure of the errors as there are many copula families introduced in the class.
Consequently, in this study, we apply various multivariate exchangeable Archimedean copulas to join the equations in the SUQR model. We believe that our model will become more flexible and applicable to investigate the entire conditional distribution of the dependent variable and becomes more robust against outliers. In the estimation aspect, as our proposed model contains large parameter estimates and the full likelihood function is quite complicated, thus, the Bayesian estimation is employed in this study. To confirm the accuracy and reliability of our model and estimation, we conduct the simulation study and real data analysis to evaluate the performance of our model and the Bayesian estimation. To the best of our knowledge, no research introduced the multivariate exchangeable Archimedean copulas to join the error of equations in the SUQR model. To this end, we suggest a Bayesian approach to the multivariate exchangeable Archimedean copula-based SUQR.
e outline of the remaining sections is as follows. In Section 2, we explain the copula-based SUQR model and multivariate copula functions. e posterior of the Bayesian is provided in Section 3. Four simulation studies are provided in Section 4. In Section 5, the real data example is used to show the performance of our model. Finally, Section 6 is the conclusion.

Seemingly Unrelated Quantile Regression under Asymmetric Laplace
Distribution. Jun and Pinkse [5] and Waldmann and Kneib [6] extended the SUR model of Zellner [1] to quantile in feature to gain more robustness against outliers in the response measurements and allow explaining the entire conditional distribution of the outcome variables in the system equation. us, the asymmetric Laplace distribution (ALD) is generally specified as the likelihood of the model. Tu et al. [14] mentioned that the ALD has several mixture representations, for example, a scale mixture of normal with exponential distribution (2015) and a scale mixture of uniform with Gamma distribution according to Wichitaksorn et al. [15]. In this study, we consider the ALD as a mixture of normal as it was proved to be more efficient than as a scale mixture of uniform.
Let y i,t and x i,t be the dependent and independent variables of equation i at time t for i � 1, . . . , n and t � 1, . . . , T. e model is formulated as where β α i is the vector of parameters of equation i at given quantile α. ε i,t is the unobserved error term of equation i at time t which is assumed to follow the ALD with mean zero, variance σ α , and quantile level or quantile α. We note that the errors are allowed to be correlated across equations to gain more efficiency. e quantile level α has the range 0 < α < 1.
us, α th conditional quantile of y i given x i is simply Typically, ε 1,t , . . . , ε n,t are joined through a multivariate distribution, especially the multivariate ALD. Nevertheless, this joint distribution is based on a linear relationship and the same error distributions. To relax this restriction, the copulas are suggested to model the nonnormal and nonlinear dependence structure in the multivariate SUQR model.
us, the n-dimensional joint distribution function H(·) can be defined as where z 1 , . . . , z n are the realization of random variables (or the standardized residuals). If the marginals F i , i � 1, 2, . . . , n, are continuous, then the copula C associated to H(·) is unique, and the expression in equation (3) can be rewritten as 2 Complexity where u ∈ [0, 1] n is the uniform distribution. To construct the copula density distribution, it can be obtained by c u 1 , . . . , u n � zC u 1 , . . . , u n zu 1 , . . . , zu n . (5) In this study, we consider two copula classes, namely, elliptical copulas and Archimedean copulas. e explicit form of these copula classes is presented in the next section.

Posterior Distribution of Multivariate Copula Families.
Two classes of copulas, namely, elliptical and Archimedean copulas, are presented. Elliptical copulas consist of two families, i.e., Gaussian and Student-t copulas. ese copulas are of symmetric dependence structure. Another class is Archimedean copulas consisting of many families, but in this study, we consider only Clayton, Gumbel, Joe, and Frank copulas, which are prominent families and mostly used in many previous studies. To construct the posterior distribution of the copula, we multiply the copula density with the prior distribution.

Elliptical
Copulas. Following Smith [17] and Smith et al. [18], the uniform prior is assumed for the copula parameter; thus, the posterior of elliptical copulas can be formulated as in the following: (1) Posterior Gaussian copula: Let Φ(·) be the standard normal cumulative distribution and Φ − 1 (·) be the quantile function of the standard normal distribution; the posterior distribution of the Gaussian copula can be written as where θ α G is the n × n dependence matrix of the Gaussian copula on quantile α. e priorπ(θ α G ) is assumed to have uniform distribution [0, 1]. us, the posterior is solely dependent on the copula density. As the Bayesian estimation is used to draw the parameter, we thus draw the updated dependence parameters from the truncated normal distribution [− 1, 1] interval.
(2) Posterior Student-t copula: where π(v α ) and π(θ α T ) are the uniform prior density for θ α T and the exponential prior density for v α , respectively. f(θ α T ) is the copula density of the cumulative Student-t-distributed random vector. According to Demarta and McNeil [19], this copula density is where θ α T is the n × n dependence matrix of the Studentt copula on quantile α and v is the degree of freedom of is the vector of the inverse Student-t distribution function. Γ is the gamma distribution. In the proposal distribution, we randomly select θ α T and v α from the truncated normal distribution [− 1, 1] interval and distribution [0, ∞] interval, respectively.

Archimedean
Copulas. Different from the elliptical copula case, we apply the uninformative prior to these density functions; thus, the posterior density for each of the Archimedean copula is as follows: e density of this class' copulas, namely, Frank, Clayton, Gumbel, and Joe, varies as proposed in Hofert et al. [20]. We can write these copulas as follows: (1) Frank copula: where denotes the polylogarithm of order n − 1 at z.
To draw the updated dependence parameters of these copula functions, we randomly select the dependency parameter ( α θ ) from the proposal truncated normal distribution [0, ∞] for Clayton, [1, ∞] for Gumbel and Joe, and [− ∞, ∞] for Frank.

Bayesian Inference
We consider the Bayesian approach for estimating all unknown parameters in our proposed model. e Bayesian estimation requires the specification of the likelihood function and the prior distribution for all the estimated parameters. Hence, the posterior density of our proposed model is constructed by multiplying the full likelihood function of the model (ALD densities and copula density) with a prior density of the parameters. Let us consider the first part of the likelihood; according to Yu and Moyeed [21], the density of the asymmetric Laplace distribution is given by where ρ α (·) is the so-called check function defined by , with Ι · { } denoting the usual indicator function. en, the posterior distribution for β α and σ α can be written as where π A (β α , σ α ) is the prior distribution for β α and σ α . Note that the copula function is used to join the errors in SUQR equations; hence, the full conditional posterior distribution of this model can be formulated as follows.
For elliptical copulas, For Archimedean copulas, where Ψ α � β α , σ α , θ α is the parameter set in the model. In the estimation aspect, all parameters are drawn by an iterative Gibbs sampler with the Metropolis-Hastings algorithm over a partition of parameter blocks: (i) the unknown parameter β α ; (ii) the variance of the model σ α ; and (iii) the copula dependence parameter θ α . Concerning the prior specification in these three groups, we assume the following priors: where B is the vector of prior means, s is the vector of hyperparameters for variance of B, and a and b are the positive hyperparameters for σ α . We select these three priors since the sign of β α can be either positive or negative and asymmetric. For the case of copula dependence parameter θ α , we assume it to be uniformly distributed with minimum zero and maximum one as there are various copula families considered in this model, and the prior information of the copula parameter is generally unknown as well as difficult to specify. To simplify our 4 Complexity conditional posterior distribution, the uniform prior for the copula parameter is assumed. As there are large parameter estimates in our model, the block Gibbs sampler with the MH sampling is considered to sample the parameters in the chain. e adaptive sampler algorithm can be explained as follows [22]: (1) Starting at an initial parameter value Ψ α 0 � β α 0 , σ α 0 , θ α 0 : in this step, the initial value of β α 0 is computed from the traditional quantile regression estimation while σ α 0 � T t�1 (ε t /(T − k)) and θ α 0 � 0, θ α 0 � 2 for elliptical copulas and Archimedean copulas, respectively.
(2) Updating the candidate parameter based on the proposal function; the proposal function and p is the number of coefficients in each equation.
where Gamma(·) is the random generation for the gamma distribution.
(iv) Simulate θ α from the truncated uniform distribution, U(L, U), where L and U are the lower bound and upper bound of the copula parameter, respectively. (3) At the j-th iteration, the acceptance function is employed to strike a balance between the following two constraints: (a) the sampler should tend to approach higher probability areas under the full posterior distribution and (b) the sampler should explore the space and avoid getting stuck at one site. Each candidate ψ α j in each iteration is considered to be a proposal with acceptance probability equal to one, with the "proposal function" selected appropriately. In this computation aspect, the acceptance function can be defined as the ratio between posterior based on candidate parameters ψ α j and posterior based on previously updated parameters ψ α j− 1 . us, we can calculate the acceptance probability: Repeat steps 2 and 3 for j � 1, . . . , J in order to obtain samples Ψ α 1 , . . . , Ψ α J . In the sampling method, we specify the number iteration to be 50,000, whereas the first 20,000 iterations are discarded as burn-in. en, we can obtain the estimated parameter by averaging the remaining 30,000 simulated sets of parameters Γ α 20001 , ..., Γ α 50000 . As we consider many copula families, we employ the deviance information criterion (DIC) to compare the performance of different copula families.

Simulation Study
To learn about the performance of the Bayesian estimation for fitting our proposed model, we conduct four simulation studies: (1) First, we examine and evaluate the accuracy of the Bayesian estimation on our model under Clayton, Gumbel, and Frank copulas (2) In the second part, we examine the performance of our model when misspecified copula is assumed (3) In simulation study 3, we evaluate the finite sample performance of the Bayesian estimation under various sample sizes (4) Finally, we investigate the performance of the Bayesian method in the high-dimension setting

Simulation Study 1: Accuracy in Parameter Estimation.
In the simulation study, Archimedean copula families, namely, Clayton, Gumbel, and Frank, are considered to model the dependence structure of the SUQR. In this study, the simulation is the realization of SUQR with three equations. us, our simulated equation can be written as Note that the three-dimension error terms, ε 1,t , ε 2,t , ε 3,t , are assumed to have asymmetric Laplace distribution with skewness or quantile α. In this simulation, we first simulate the uniform margins from the three-dimension copula model. We set the true value of the correlation coefficient of Gumbel, Clayton, and Frank copula parameters to be 2.5, 0.5, and 2, respectively. en, the obtained uniform data are Complexity 5 transformed to be errors u i1 , u i2 , and u i3 using the quantile function of the asymmetric Laplace distribution with variance σ 2 and skewness or quantile parameter α. e independent variables x 1,t , x 2,t , and x 3,t are randomly simulated from a standard normal distribution. e true parameters in this simulation model are provided in Tables 1-3. We generate 1,000 data sets each with T � 500. To assess the accuracy of the Bayesian estimation, we consider the average of the estimated parameters and the average of the estimated standard deviations. e simulations are conducted for three quantile levels, say α � 0.25, 0.5, and 0.75, and the results are provided, respectively, in Tables 1-3, showing the average of the estimated parameters and their standard errors. We can see that our Bayesian estimation produces the unbiased parameter estimates for the SUQR model. It is observed that the estimated parameters are close to the true values in all cases, and the average standard deviation from the parameter is reliable. is simulation result indicates the adequacy and reliability of our Bayesian estimation of all unknown parameters in the copula-based SUQR model. In this simulation study, the hyperparameters a, b, c, and d are specified as 0.1 to reflect weak prior information.

Simulation Study 2: Robustness and Kullback-Leibler
Divergence. In this section, another simulation study is proposed to measure the performance of the model using the relative entropy, also known as Kullback-Leibler divergence (KLD) [23]. is relative entropy is the measurement of the difference between two probability distributions. Consider the continuous probability distribution, and let f ⌢ and f denote the density of probabilities F ⌢ and F; thus, the KLD is given as In this study, we define f ⌢ as an alternative or approximated posterior and define f as a true posterior function when all parameters are known. We simulate the data sets in a similar way as in the previous simulation study (Section 4.1). If the copula family is known, then the proposed estimator and model could provide an accurate result. However, this estimation will be valid only in the simulation study. In practice, we need to select the true copula function to join the error terms in the SUQR model. Noh et al. [24] discussed the copula misspecification and suggested that the selection of the wrong copula function will bring about bias in the estimation of the model. e purpose of this simulation study is to investigate the distance between the true SUQR function and its approximation when the copula is correctly specified and when the copula is misspecified. In this simulation study, we set the true posterior function to be the Student-t SUQR function for quantile levels 0.25, 0.50, and 0.75. We compare the true SUQR function and its approximation (in terms of posterior function) among the SUQR function family (i.e., Gaussian, Student-t, Joe, Clayton, Gumbel, and Frank) as well as the conventional SUQR function (noncopula-based). Figure 1 illustrates the three-panel result for quantile levels 0.25, 0.50, and 0.75. As expected, the Student-t copula-based SUQR model achieves its minimum and is close to the true SUQR function line (dashed line) at every quartile level. We also made a performance comparison between our proposed model and the conventional SUQR function (M0) and found that our proposed model performs better than the conventional model since D(f True , f Student− t ) < D(f True , f M0 ). In the case of misspecified copula function, we observed that the misspecified copula functions bring a larger deviation of the approximated SUQR function from the true function. According to these results, we can say that our proposed model is the robust model, and the incorrectly specified copula function will lead to the low accuracy of the model.

Simulation Study 3: e Finite-Sample Properties of the Bayesian Estimation.
In the third simulation study, the finitesample properties of the Bayesian estimation in our proposed model are investigated upon the calculated absolute Bias and mean squared error (MSE) of the estimator. Again, we simulate the data sets the same way as in the first simulation study. e absolute Bias and MSE can be calculated by    where R � 1,000 is the number of Monte Carlo replications and ψ r and ψ r are the estimated values and the true values, respectively. e sample sizes are fixed at 100, 500, and 1,000 d for each replication.
Tables 4-6 contain the results of the Bayesian estimation over the 1,000 simulated data sets with three different dimensions for the sample size. e most important finding from Tables 4-6 is that the Bayesian estimation provides reliable parameter estimates as the absolute Biases and MSEs are close to zero. In addition, the Biases and MSEs seem to be lower as the sample size increases. is result indicates that the Bayesian estimator is asymptotically unbiased and consistent for the estimated parameter in our model. is same pattern of convergence to zero is repeatedly obtained, considering different levels of quantile α.

Simulation Study 4: Evaluation of the Performance of the Bayesian Estimation in the High-Dimension Copula-Based SUQR.
Our model is supposed to propose a quite general model with the possibly arbitrary n-dimension copula-based model. So, the last simulation study is conducted to examine the performance of our model when the number of equations is large. We consider the sample size T � 1,000, the dimension n∈{3, 5, 7}, the dependency parameter as 3 for all dimensions, and the one-parameter (Archimedean) copula families of Gumbel, Clayton, and Frank. For each of these combinations, we generate a random sample of the corresponding size and compute θ α Gu , θ α C , θ α F . We repeat this procedure R � 1,000 times and compute the absolute Bias and MSE for each of the three quantile levels, α � (0.25, 0.50, 0.75). To simplify this simulation study, we set the marginal parameters for each equation as 1, say β α 1i � β α 2i � 1. Since the Bayesian estimation provides the similar marginal parameter estimates (parameters in each equation), we focus here on the copula parameter estimates θ α Gu , θ α C , θ α F . Once all the data are simulated, we fit the copula-based SUQR model with α � (0.25, 0.50, 0.75), and the absolute Bias and MSE are recorded. e results are shown in Figure 2. We can see a pattern of diverging from zero of the absolute Bias and MSE when n increases. As a general statement, we can say that the Bias and MSE tend to deviate from zero when the number of copula dimensions increases, indicating that the estimates based on the Bayesian estimation may not provide a good result. Our result is in line with the study of Embrechts and Hofert [25], which reveals that the estimation of the parameter in the exchangeable copula models becomes increasingly severe in higher dimensions. Another reason is that the estimation cannot be improved with the higher computation cost.
is is important in that we may not gain accurate parameter estimates in higher dimensions. However, our Bayesian estimation is still performing promisingly acceptable in high dimensions as the absolute Bias and MSE are not quite high.

Estimation Results.
In this section, we illustrate the applicability of our proposed model and the Bayesian estimation developed in this study, using the same data set as in Tansuchat et al. [13]. e data set consists of several variables measured in 177 months. Here, we focus on three US stock returns in the NASDAQ market, consisting of  Table 7. Note that all data are presented in log-return form.  According to the data description, the mean of SMB has the largest mean value. All stock returns show a negative skewness, while SMB and HML show a positive skewness, indicating that they are more likely to have a negative return on stock markets. e kurtosis of all series is greater than 3, except for SMB. is means that the returns are very much e Jarque-Bera test is then conducted to investigate the normal distribution of these data. e result shows that data reject the null hypothesis of normality at 1% significant level. In addition, we also conduct the augmented Dicky-Fuller (ADF) test to examine the stationary our variables and found that all data are stationary at 1% significant level. Our empirical SUQR model is constructed under the Fama-French approach [26,27]. e empirical model can be shown as where y i,t is the return of asset i at time t, r f,t is the risk-free rate which is measured by the US. Treasury bill, y M,t , is the return of the market, SMB t is small minus big (a proxy for company size), and HML t is high minus low (a proxy for book-to-market values). e parameters β α 1i , β α 2i , β α 3i are referred to as beta risk, level of exposure to size risk, and value risk, respectively. e choice of hyperparameters is very delicate in our Bayesian estimation; it is important to determine an approriate prior informations. In this empirical study, we decide to suggest three priors, which are as follows: (1) Weak informative prior: the hyperparameters a, b, c, and d are specified as 0.1.
where I is the identity matrix. Santos and Bolfarine [28] suggested that this prior has not presented any problem, and the posterior distribution seems to be insensitive to minor changes. (3) Informative prior: a � 2, b � 1, c � 1, and d � 0.1.
Prior to showing the estimated results, we compare the performance of various copula-based models as well as the conventional model of Jun and Pinkse [5] and Waldmann and Kneib [6] (M0) and the multivariate elliptical copulabased SUQR of Tansuchat et al. [13]. As a sensitivity analysis, we have also investigated various hyperparameters. By using the DIC, Table 8 reports the comparison results on 0.25, 0.50, and 0.75 quantiles. Among the trial runs of several alternative copula functions as well as hyperparameter prior, we learn that Clayton presents the lowest DIC for quantile level at α � 0.25, 0.75, while Frank presents the lowest DIC for α � 0.5.
is indicates that different copula families can model different quantile models. In addition, we compare the best fit model with the conventional model of Jun and Pinkse [5] and Waldmann and Kneib [6] and with the two elliptical copula-(Gaussianand Student-t-) based models of Tansuchat et al. [13]. e results reveal that our model is more adequate for this data set at all quantiles than other symmetric copula families. We also conclude that if the same copula is used to construct the joint between errors of the SUQR model, the unreliable result may be obtained. Moreover, we also learn that the model selection result seems to be insensitive to the hyperparameter prior as similar results are obtained. However, when we compare the DIC of these three priors, we find that the informative prior seems to be more valid than the others. Table 9 shows the estimated posterior mean and standard deviation parameters for quantiles 0.25, 0.50, and 0.75. We can observe that all estimated parameters of all stocks change when the quantile changes. is indicates that our model is robust against the outliers. In addition, we also find that the values of these copula dependences are different in different quantiles. It is noticed that the estimated copula parameters are quite high for all quantiles indicating a strong positive relationship among the errors of our model. To illustrate the correlation among the 19 marginals (cumulative of the standardized residual of each equation), threedimension scatter plots among marginals for 0.25, 0.50, and 0.75 are illustrated in Figure 3. Finally, for this application study, we can see that the model estimates are reliable, and acceptable results are obtained. Moreover, as a sensitivity analysis, we have also provided the estimation results of the best fit model based on various hyperparameter priors. We find that the estimated posterior mean for each predictor is quite similar under various priors, and they do agree on the importance of the variables.
To formally check the convergence of MCMC chains, the marginal posterior distributions can be visualized by plotting the histograms of the simulated parameter draws. Figure 4 displays the posterior distributions θ α for quantiles 0.25, 0.50, and 0.75. We only present the case of copula parameters as similar results are obtained in other parameters. e plots show that the MCMC sampler converges to a normal distribution and mixes very well. We may thus conclude that the Bayesian method estimates the copula parameter quite well.

Goodness-of-Fit Tests.
To carry out the goodness-of-fit test for our proposed models, Cramer-von Mises (CvM) method is conducted in this section. Genest et al. [29] suggested that the CvM test is the most powerful test to check the goodness of fit of copula models. In this test, significant statistics indicate that the copula models based on the data are rejected. e result is reported in Table 10. For α � 0.25, the result shows that the Clayton copula yields the highest p value, indicating Clayton copula-based SUQR offers a better fit for the data than other models at quantile 0.25. In the case of α � 0.50 and α � 0.75, Frank and Clayton copulas, respectively, provide the better fit than other copula-based models.

Conclusion
In this paper, we introduced the multivariate exchangeable Archimedean copula to join the errors of the seemingly unrelated quantile regression (SUQR). e model becomes more accurate and robust against the outlier relationship between the dependent and independent variables. We also introduced the Bayesian Markov chain Monte Carlo approach to estimate the parameter sets of our proposed model. As the posterior distribution of the copula parameter does not appear to be in any form, therefore, we employ a Bayesian estimation together with a Gibbs sampler with the Metropolis-Hastings algorithm to infer the full posterior  distribution. To examine the accuracy of our Bayesian estimation and performance of our proposed model, we present the simulation study and real data analysis. Four simulation studies are conducted. e result of the first simulation study shows the accuracy of the Bayesian estimation. e results confirm that our proposed model is well estimated as reliable estimation results are obtained for every quartile level, and the parameter estimates on average are close to their true values. In the second simulation, we adopt the Kullback-Leibler divergence (KLD) to measure the distance between the true posterior probability distribution and the approximated posterior probability when the copula function is unknown. e result confirms the robustness of our model and presents the closest distance between the correctly specified and the true probability function. e third simulation is proposed for examining the finite-sample properties of the Bayesian estimation in our proposed model. e result shows that the Bayesian estimation provides reliable parameter estimates as the absolute Biases and MSEs converge to zero when the sample size increases. Finally, the performance of the Bayesian estimation in high-dimension copula-based SUQR is investigated, and the result reveals that we may not gain the accurate parameter estimates inhigher dimensions. However, the Bayesian estimation is still performing promisingly acceptable in high dimensions as the absolute Bias and MSE are not quite high.
In the real data application, we apply our proposed model to the data set provided by Tansuchat et al. [13]. We quantify and measure the risk of the stock price through the Fama-French model with three-factor analysis. We compare the performance of our model with the Archimedean copula-based SUQR models of Tansuchat et al. [13] as well as the conventional models of Jun and Pinkse [5] and Waldmann and Kneib [6]. e result shows that our proposed model is flexible and has the potential to capture the extreme market condition. Although our proposed model presents a good performance in both simulation and real data studies, Gibbs algorithm with MH may not be the most suitable for our proposed model, given the large number of parameters to estimate. In addition, it can be difficult to choose the proposal functions in the Gibbs sampler with the MH sampling algorithm. For future study, the adaptive Gibbs sampler with the MH algorithm, which uses the history of the process to tune the appropriated proposal distribution, can be applied to our proposed model. Moreover, future research should consider applying our proposed model subject to pseudocyclical-structural changes since many financial time series exhibit behavioral change over time.

Data Availability
In this study, we use the simulated data to show the performance of our model, and the simulation processes are already explained in the paper. For the real data analysis section, we use the same data of Tansuchat et al. [13]. ese data can be freely collected from http://www. investing.com (http://www.investing.com) or from omson Reuter Datastream. By the way, the data are available from the corresponding author upon request (woraphon.econ@ gmail.com).

Conflicts of Interest
e authors declare that they have no conflicts of interest.