Airline Sustainability Modeling: A New Framework with Application of Bayesian Structural Equation Modeling

: There are many factors which could inﬂuence the sustainability of airlines. The main purpose of this study is to introduce a framework for a ﬁnancial sustainability index and model it based on structural equation modeling (SEM) with maximum likelihood and Bayesian predictors. The introduced framework includes economic performance, operational performance, cost performance, and ﬁnancial performance. Based on both Bayesian SEM (Bayesian-SEM) and Classical SEM (Classical-SEM), it was found that economic performance with both operational performance and cost performance are signiﬁcantly related to the ﬁnancial performance index. The four mathematical indices employed are root mean square error, coefﬁcient of determination, mean absolute error, and mean absolute percentage error to compare the efﬁciency of Bayesian-SEM and Classical-SEM in predicting the airline ﬁnancial performance. The outputs conﬁrmed that the framework with Bayesian prediction delivered a good ﬁt with the data, although the framework predicted with a Classical-SEM approach did not prepare a well-ﬁtting model. The reasons for this discrepancy between Classical and Bayesian predictions, as well as the potential advantages and caveats with the application of Bayesian approach in airline sustainability studies, are debated.


Introduction
Measuring, predicting, and estimating the sustainability indices of airline industries has always been of great value to airline directors and researchers. In this regard, some researchers in their sustainability modelling focused on financial indicators [1][2][3], some only dealt with operational indicators [4][5][6][7], while few of them concentrated estimating modelling on both financial and operational performance [8] indices. Moreover, some sustainability modelling has been focused on the cost indicators. In these types of the studies, generally, researchers consider a cost indicator as a function of operational indicators with [9,10] or without [7,11] another cost indicator. Most of these types of studies do not consider country economic indicators in their research framework, and they focus on internal indicators of the company.
Accounting and financial indices have been the focus of much research across many industries. The concerned indices stand for one of the most essential communicational means applicable to senior senior management [12]. Therefore, assessment of the performance is needed, particularly in the financial arena, and, as expected, considerable capital is vital for the sustainability of these airline companies. Financial performance indices have a particularly critical role in the survival of an airline. Consequently, the airline needs to evaluate and assess the financial performance indices to determine its financial situation between the competing companies and firms. Therefore, the first purpose of the current study is to introduce a new framework which is able to fill previous gaps of airline sustainability modelling by considering economic performance, operational performance, and cost performance for estimating financial performance.
Linear and nonlinear regression modeling analyses have become the basic techniques for airline sustainability modeling, however, individual regression analysis for each dependent variable is hardly challenged as a realistic approach in the situations where the outcomes are logically and naturally related. Furthermore, some research frameworks are difficult to analyze by a regression model when an outcome is determined not only by direct impacts of the predictor variables but also by their unobserved common cause. Classical structural equation model (Classical-SEM) is a suitable technique that can address the above limitations, providing a robust method for studying interdependencies among a set of correlated variables.
In recent years, Classical-SEM has attracted the attention of many researchers as a commonly adopted method used for tasks such as data analysis in airline disciplines including sustainability [32], low cost [33], job satisfaction [34], and service quality [35]. This application presents an advanced version of linear regression with the main goal of examining the hypothesis that observes a covariance matrix for a set of measured indicators that is equal to the covariance matrix described by the hypothesized model. In linear analysis, normal distribution of residuals is a vital assumption. Otherwise, it is possible to determine the sample covariance matrix with a standard approach. Therefore, to overcome this setback, Yao, Xu [36] suggested using Bayesian-SEM for superior estimation. For the parameters of interest, Bayesian-SEM allows researchers to apply prior information to update current information. This involves utilizing the Gibbs sampler [37] to obtain samples of arbitrary size to summarize the posterior distribution for describing the parameters of interest. With these samples, the point estimates, standard deviations and interval estimates can be To analyze airline sustainability, different statistical methods are employed, such as analysis of variance (ANOVA) [18,19], panel data modeling [20][21][22], time series [4,23,24], data envelopment analysis (DEA) [25][26][27][28], neural networks [29], neuro-fuzzy systems [30], and Classical-SEM [8,31].
Linear and nonlinear regression modeling analyses have become the basic techniques for airline sustainability modeling, however, individual regression analysis for each dependent variable is hardly challenged as a realistic approach in the situations where the outcomes are logically and naturally related. Furthermore, some research frameworks are difficult to analyze by a regression model when an outcome is determined not only by direct impacts of the predictor variables but also by their unobserved common cause. Classical structural equation model (Classical-SEM) is a suitable technique that can address the above limitations, providing a robust method for studying interdependencies among a set of correlated variables.
In recent years, Classical-SEM has attracted the attention of many researchers as a commonly adopted method used for tasks such as data analysis in airline disciplines including sustainability [32], low cost [33], job satisfaction [34], and service quality [35]. This application presents an advanced version of linear regression with the main goal of examining the hypothesis that observes a covariance matrix for a set of measured indicators that is equal to the covariance matrix described by the hypothesized model. In linear analysis, normal distribution of residuals is a vital assumption. Otherwise, it is possible to determine the sample covariance matrix with a standard approach. Therefore, to overcome this setback, Yao, Xu [36] suggested using Bayesian-SEM for superior estimation. For the parameters of interest, Bayesian-SEM allows researchers to apply prior information to update current information. This involves utilizing the Gibbs sampler [37] to obtain samples of arbitrary size to summarize the posterior distribution for describing the parameters of interest. With these samples, the point estimates, standard deviations and interval estimates can be computed for the purpose of making an inference. The Bayesian approach is attractive, as it allows for the use of prior information to update current information regarding the parameters of interest.
Lee [38], provides some advantages of Bayesian-SEM prediction: • Mainly first moment properties of the raw individual observations are used for statistical methods, which make improvements of analyses much simpler compared to the second moment properties of the sample covariance matrix. Hence, it is easier to use in more complex states. • Direct impact of the latent variables (construct) is possible which makes obtaining factor score estimates simpler compared to that of the classical regression techniques.

•
As it directly models manifest variables with their latent variables through the familiar regression functions, it provides a more direct interpretation and enables the use the common techniques in regression modeling such as residual and outlier analyses in conducting statistical analysis.
With Bayesian predictors, as pointed out by Scheines, Hoijtink [39], Lee and Song [40], and Dunson [41], this technique allows the researchers to use prior experts' beliefs in addition to the sample information to produce better outputs and deliver valuable statistics and indices including the mean and percentiles of the posterior distribution of the unknown parameters. In conclusion, more reliable results for small samples can be achieved. In contrast, the Bayesian approach has much more flexibility in handling complex situations. Even though many studies have been done on determining the financial performance index, not much is done on modeling of this index using SEM, particularly when information on economic performance, operational performance, and cost performance are considered. Therefore, the second purpose of the study is to illustrate the value of Classical-SEM and Bayesian-SEM for developing a model that describes the sustainability index of an airline established in the Asia-Pacific region. The interrelationships among the latent variables, such as economic performance, cost performance, and operational performance, as well as between the latent variables and their respective manifest variables, are determined using panel data obtained from an Asia-Pacific airlines.
The structure of the current paper consists of seven sections. The first section outlines the main gaps in the existing airline sustainability frameworks. Then, the different types of statistical methods in airline sustainability modeling and the application of Classical-SEM in current studies is explained. Moreover, we mention some limitations of Classical-SEM and introduce Bayesian-SEM for better estimating of the research parameters. In this section, the main advantages of Bayesian-SEM in comparison to previous types of linear and nonlinear classical modeling is presented. Section 2 explains the literature review and research hypotheses. In this section the main trend of airline sustainability modeling is explained. This section also presents the research framework that is based on six main hypotheses. Section 3 of the paper is about the main theories of Classical-SEM and Bayesian-SEM. This section shows the procedures of dealing with prior and posterior distribution functions based on our research data structures. Section 4 presents materials and methods and explains sampling procedure and data collection. Results of the study are presented in Section 5. The outputs of Classical-SEM, Bayesian-SEM, and a comparison study between them based on familiar statistical indices are discussed in Section 5. Finally, Sections 6 and 7 are the discussion and conclusion of the study.

Literature Review and Research Hypotheses
There is a vast amount of literature concerning airline sustainability modelling using a variety of approaches. Early studies of Caves, Christensen [42] and Sickles [43,44] tend to employ energy, material, capital, and labor. A couple of years later, computerized reservation systems and related indicators including number of computers for ticket selling and number of agencies were considered by some researchers like Borenstein [45], Banker and Johnston [46], and Duliba, Kauffman [47]. Since the early 2000s, country economic indicators have been considered as vital indicators for estimating performance for many airline sustainability modeling studies. Previous studies confirmed that Gross domestic product [48], human development index [8,14], and foreign direct investment [49] are the main country economic key indicators which affect airline performance. Therefore, in this study, the combination of those indicators were defined as the economic performance latent variable. In the current paper, we define financial performance as a grouping of familiar financial indicators. Total assets [6,50], operating profit [14], and total revenue [51] are the most commonly used performance indicators in airline sustainability modelling. In this study we define the combination of total assets, operating profit, and total revenue as the financial performance latent variable. Studies by Moon, Lee [52] and Ismail and Jenatabadi [8] confirmed the impact of country economic indicators on airline financial performances. Therefore, we consider our first hypothesis of study with the following statement: • H 1 : Economic performance has a significant impact on financial performance.
Operational performance measures have been broadly applied by a good number of corporations since the early 1990s to measure current performance, identify requirements needed to enhance performance, and make the achievement of far-fetched strategic goals possible [53]. Recently, operational performance measures have been able to gain a global prevalence as myriad organizations and companies around the world have shifted their attention and reliance from the traditional method based mainly on financial performance measures to a range of non-traditional value indices [54]. Revenue passenger kilometer, revenue tone kilometer [55,56], and number of departures [6] are the main operational performance indicators in the airline industry. In our study, we define operational performance with a combination of these three indicators. Logically and empirically operational performance has an impact on financial performance and the relationship between economic performance and operational performance is confirmed by previous studies [57]. Therefore, we considered these relations in the research model and tested the following hypotheses: • H 2 : Operational performance has a significant impact on financial performance. • H 3 : There is a significant relationship between economic performance and operational performance.
Cost function is another type of airline sustainability study. In this type of study, researchers consider cost indicators of two types. The first type cost indicator is a function of operational performance: Cost indicator = f (operational indicators) Zuidberg [7] and Hansen, Gillen [11] have done their modeling based on the above function. The second type is a financial indicator and it is a function of cost and operational indicators; Financial indicator = f (cost indicators, operational indicators) Johnston and Ozment [9] and Oum and Zhang [10] have done their modeling based on the above function. However, the combination of two types of modeling, especially with the leveraging of country economic performance, is rare. Therefore, this study considers cost performance as the fourth latent variable with a combination of operating cost, labor cost, and fuel cost indicators based on the Zuidberg [7] study. Considering this latent variable lead to the development in our study of the following hypotheses: • H 4 : Cost performance has significant impact on financial performance. • H 5 : There is a significant relationship between cost performance and operational performance. • H 6 : There is a significant relationship between economic performance and cost performance. Figure 2 shows the hypothesized research model with the latent variables, while their indicators serve to show the impact of economic performance with both cost performance and operational performance on financial performance. The figure illustrates that the first three constructs are interrelated. As a result, the present research model includes four constructs and twelve observed variables.

Classical-SEM & Bayesian-SEM Theories
To perform statistical analysis, classical as well as Bayesian paradigms are initially used. According to traditional principles, supposing the parameter of interest is constant (nonstochastic), inferential subjects about are handled based on likelihood/log-likelihood. If the likelihood is denoted by ( ), it is assumed that the information about can be obtained only through sample ∈ X, and the likelihood is a function of conditional on the observed value of . However, by adopting the Bayesian paradigm, it can be assumed that Θ is stochastic and it can be incorporated in the model as a random variable. In this regard, has a probability measure ( ), or a prior distribution that gives more information about the parameter than likelihood. The information may be from different sources such as physical reasoning or expert views. Then the information in ( ) about can be updated using the likelihood information, yielding the posterior distribution denoted by ( | = ). From a Bayesian point of view, an inference about can only be made through posterior distribution.
Note that the likelihood can be considered as the distribution of the data given the parameter value. Based on Figure 3, the major portion of the prior distribution has a lower parameter value than that at the peak of the likelihood. The posterior is obtained as a compromise between the prior and the likelihood. There are three main types of prior probability distributions (informative, uninformative, and weakly informative) that vary in their degree of (un)certainty about the parameter of interest in the research model [59]. Applying informative priors means employing data from theories, literature, expert opinion, or previous experiments. Informative priors can have a significant influence on the final parameter estimates. Uninformative priors are considered by the researchers when there is no

Classical-SEM & Bayesian-SEM Theories
To perform statistical analysis, classical as well as Bayesian paradigms are initially used. According to traditional principles, supposing the parameter of interest Θ is constant (non-stochastic), inferential subjects about Θ are handled based on likelihood/log-likelihood. If the likelihood is denoted by L (Θ), it is assumed that the information about Θ can be obtained only through sample x ∈ X, and the likelihood is a function of Θ conditional on the observed value of X. However, by adopting the Bayesian paradigm, it can be assumed that Θ is stochastic and it can be incorporated in the model as a random variable. In this regard, Θ has a probability measure π (Θ), or a prior distribution that gives more information about the parameter than likelihood. The information may be from different sources such as physical reasoning or expert views. Then the information in π (Θ) about Θ can be updated using the likelihood information, yielding the posterior distribution denoted by π(Θ|X = x) . From a Bayesian point of view, an inference about Θ can only be made through posterior distribution.
Note that the likelihood can be considered as the distribution of the data given the parameter value. Based on Figure 3, the major portion of the prior distribution has a lower parameter value than that at the peak of the likelihood. The posterior is obtained as a compromise between the prior and the likelihood.

Classical-SEM & Bayesian-SEM Theories
To perform statistical analysis, classical as well as Bayesian paradigms are initially used. According to traditional principles, supposing the parameter of interest is constant (nonstochastic), inferential subjects about are handled based on likelihood/log-likelihood. If the likelihood is denoted by ( ), it is assumed that the information about can be obtained only through sample ∈ X, and the likelihood is a function of conditional on the observed value of . However, by adopting the Bayesian paradigm, it can be assumed that Θ is stochastic and it can be incorporated in the model as a random variable. In this regard, has a probability measure ( ), or a prior distribution that gives more information about the parameter than likelihood. The information may be from different sources such as physical reasoning or expert views. Then the information in ( ) about can be updated using the likelihood information, yielding the posterior distribution denoted by ( | = ). From a Bayesian point of view, an inference about can only be made through posterior distribution.
Note that the likelihood can be considered as the distribution of the data given the parameter value. Based on Figure 3, the major portion of the prior distribution has a lower parameter value than that at the peak of the likelihood. The posterior is obtained as a compromise between the prior and the likelihood. There are three main types of prior probability distributions (informative, uninformative, and weakly informative) that vary in their degree of (un)certainty about the parameter of interest in the research model [59]. Applying informative priors means employing data from theories, literature, expert opinion, or previous experiments. Informative priors can have a significant influence on the final parameter estimates. Uninformative priors are considered by the researchers when there is no There are three main types of prior probability distributions (informative, uninformative, and weakly informative) that vary in their degree of (un)certainty about the parameter of interest in the research model [59]. Applying informative priors means employing data from theories, literature, expert opinion, or previous experiments. Informative priors can have a significant influence on the final parameter estimates. Uninformative priors are considered by the researchers when there is no prior knowledge about the parameter of interest [60]. A compromise between the informative and uninformative priors is called a weakly informative prior [61][62][63]. Weakly informative priors can be recommended when the researchers want to apply a weakerer prior than what your actual knowledge would allow [62]. Weakly informative priors include some information about the parameter estimate but do not typically impact the final parameter estimate to a large extent. In our research model, we do not have any prior knowledge about the parameter of interested, therefore, the uninformative priors are specified.
In what follows, each analysis is specifically addressed. From a conventional viewpoint, the classical SEM is initially specified, and the measurement and structural relations are defined. Suppose that the measurement equation is where x i is a p × 1 vector of indicators describing the q × 1 random vector of latent variable ω i ; Λ is a p × q matrix of the loading coefficients obtained from the regressions of x i on ω i ; and ε i is p × 1 represents random vectors of the measurement errors that are summed to be the distribution according to N (0, ψ ε ) in the current setup. It is further assumed that vectors ω i , i = 1, . . . , n are independent, uncorrelated with ε i , and specifically distributed according to N (0, Φ). To accommodate the relation between endogenous and exogenous variables, ω i is partitioned as (η i , ξ i ), where η i and ξ i are the m × 1 and n × 1 vectors of the latent variables, respectively.
At this stage, the following structural equation is considered: where B is the m × m matrix of the structural parameters governing the relationship among endogenous latent variables, which is assumed to have zeros in the diagonal; Γ is the m × n regression parameter matrix for relating the endogenous with exogenous latent variables; and δ i is the m × 1 vector of disturbances, which is assumed to be distributed according to N (0, ψ δ ) where ψ δ is a diagonal covariance matrix. It is further assumed that δ i is uncorrelated with ξ i . Since only one endogenous latent variable is involved in this study, in other words, Bη i = 0, The above formula can be rewritten as η i = Γξ i + δ i for simplicity. In this study, for estimating research model parameters based on the SEM technique, the robust weighted least-squares (RWLS) estimation method is incorporated. RWLS provides parameter estimates and standard errors and computes χ 2 and the fit indices that are found using the diagonal components of the weight matrix and that are derived based on the threshold asymptotic variances and latent correlation estimates [64]. After the estimation process, model evaluation is required. In this respect, the model's goodness of fit can be checked through the related Chi-square statistic [ [65]. The judgment based on these measures is discussed in detail in the empirical study section.
For perfect SEM analysis and to improve the fit, the model can be modified using χ 2 difference, Lagrange multiplier, and Wald tests. Many programs provide modification indices that specify the fit improvement as a result of adding an extra path to the model [65].
From a Bayesian viewpoint, the prior distribution must first be specified. Beforehand, similar to Yanuar, Ibrahim [66], a threshold specification has to be identified in order to treat the ordered categorical data as manifestations of hidden continuous normal distribution. As a brief explanation about threshold specification, if adopting the parameterization by Lee [38], suppose X = (x 1 , x 2 , . . . , x n ) and Y = (y 1 , y 2 , . . . , y n ) are both latent continuous variables. The relationship between X and Y is explained using the threshold specification. The procedure for x 1 is described as an instant. More precisely, let where c is the number of categories for x 1 , and τ c − 1 and τ c denote the threshold levels associated with y 1 . For example, in this study c = 3 is considered, where τ 0 = −∞ and τ 3 = ∞. Meanwhile, the values of τ 1 and τ 2 are determined based on the proportion of cases in each category of x 1 using where Φ −1 (·) is the inverse of the standardized normal distribution, N r is the number of cases in the rth category and N is the total number of cases. It is specifically assumed that Y is distributed according to a multivariate normal. Under the Bayesian-SEM, X = (x 1 , x 2 , . . . , x n ) and Y = (y 1 , y 2 , . . . , y n ) are continuous data matrices and latent continuous variables, respectively, and Ω = (ω 1 , ω 2 , . . . , ω n ) is the matrix of latent variables. The observed data X are augmented with the latent data (Y, Ω) in the posterior analysis. The parameter space is denoted by Θ = (τ, θ, Ω), where θ = (Φ, Λ, Λ ω , Ψ δ , Ψ ε ) is the structural parameter. In line with Lee [38], the prior model is given by where due to the ordinal nature of thresholds, a diffuse prior can be adopted. Specifically, for some constant c, π (τ) = c Further, to accommodate a subjective viewpoint, a natural conjugate prior can be adopted for θ with the conditional representation π (θ) = π (Λ|Ψ ε ) π (Ψ ε ). More specifically, let where ψ εk is the kth diagonal element of ψ ε , Λ k is the kth row of Λ, and Γ denotes the gamma distribution. Finally, an inverse-Wishart distribution is adopted for Φ as follows: It is further supposed that all hyperparameters are known. Posterior distribution can be found by normalizing the product L (Θ|X = x) π (Θ).
Owing to computational difficulties in identifying the posterior distribution Θ|X = x , the Markov Chain Monte Carlo (MCMC) technique is applied to generate a sequence of random observations from Θ|X = x . Then Bayesian analysis can be performed using WinBUGS, a freely available software.
The next procedure in Bayesian-SEM is convergence testing of the research model parameters. According to Yanuar, Ibrahim [67], model diagnostics are performed by graphically designing time series diagrams to evaluate the accuracy of the research parameters with different starting values and to illustrate the diagnosis based on tracing of the diagrams [39,68].
To assess the plausibility of the proposed model, which includes measurement and structural equations, the residual estimates are plotted versus the latent variable estimates to provide information on the model fit. The residual estimates for measurement equation (ε i ) can be obtained from whereΛ andξ i are Bayesian estimates obtained via the MCMC methods. The estimates of residuals in structural equation (δ i ) can be obtained from the following estimated model: whereB,η i ,Γ andξ i are Bayesian estimates obtained from the corresponding simulated observations through MCMC.
According to Figure 2, the model hypothesized in this study consists of 12 indicator variables with three exogenous latent variables and one endogenous latent variable. The following measurement model is then formulated: where ω i = (η i , ξ i1 , ξ i2 , ξ i3 ) T . The structural part of the current SEM model has the form In the data analysis, AMOS 18 is used to estimate the parameters for Classical-SEM, while the Bayesian model is fitted to the data using winBUGS version 1.4.

Materials and Methods
Based on an Air Transport World (ATW) report from 2013 [69], 106 airlines were listed in the Asia-Pacific region.
Nevertheless, it is notable to mention that airline companies are classified as a service-providing sector whose main task includes service provision to their customers. These types of airlines are categorized into three groups in terms of service type: airline companies specializing in transfer of passengers, airline companies specializing in cargo transfer, and airline companies specializing in both passenger and cargo transfer. This paper, however, only focuses on the airline firms specializing in passenger transfer although they also concurrently provide services for cargo transfer. The cargo transferring aspects of the case have been excluded from the present research domain (four companies). Moreover, nineteen low-cost airline carriers were eliminated from the present research domain (nineteen companies). Therefore, 23 companies that are trunk and low-cost carriers were excluded from the present research domain. In this study, 30 (36%) airline companies were selected over a 10-year period (2006 to 2015). The data were reported on an annual basis and gathered from an overall company level rather than city pairs. Therefore, 300 records were considered, seven of them were deleted due to missing information.
Mahalanobis distance is an extremely general measure that is utilized for measurement of multivariate outliers [70]. Based on Mahalanobis distance testing, ten observations (observation number ; 9, 32, 39, 86, 103, 122, 209, 252, 265, and 271) were eliminated from the list because they were considered as outliers, which could affect the model fit, R 2 , and the size and direction of parameter estimates (see Table 1). Therefore, (300 − 7 − 10 = 283) observations were considered in the final data of the study.  Figure 4 represents the results of model fitting based on the SEM approach. The values of GFI, IFI, TLI, and NFI are within the acceptable range. Therefore, the current model is fitted for our data at the 5% significance level.  Controlling for outliers and maintaining normal distribution support in adjusting the heterogeneity of the research data. The employment of the maximum likelihood estimator in this study uses the Classical-SEM procedure. The main essential assumption for the employment of the maximum likelihood is that the data are required to follow normal distribution and the scale of observed variables must be continuous. The normality testing that should be used in Classical-SEM is based on the value of skewness and kurtosis. If the absolute kurtosis value is less than 7 and the value of skewness is between −2 and +2, the endogenous variables normality is acceptable. Based on Table 2, revenue passenger kilometer, revenue ton kilometers, and number of departures are not normally distributed. Moreover, based on the output of the multivariate normality test, the kurtosis value is 18.69, and this value is not less than 10. Therefore, the multivariate normality hypothesis is rejected [71].  Controlling for outliers and maintaining normal distribution support in adjusting the heterogeneity of the research data. The employment of the maximum likelihood estimator in this study uses the Classical-SEM procedure. The main essential assumption for the employment of the maximum likelihood is that the data are required to follow normal distribution and the scale of observed variables must be continuous. The normality testing that should be used in Classical-SEM is based on the value of skewness and kurtosis. If the absolute kurtosis value is less than 7 and the value of skewness is between −2 and +2, the endogenous variables normality is acceptable. Based on Table 2, revenue passenger kilometer, revenue ton kilometers, and number of departures are not normally distributed. Moreover, based on the output of the multivariate normality test, the kurtosis value is 18.69, and this value is not less than 10. Therefore, the multivariate normality hypothesis is rejected [71].  Table 3 and Figure 5 present the output of research hypotheses regarding the relationship among the variables. The convergence statistics tests for each parameter of interest show that the R 2 values were approximately 1. In both models, the impact of economic performance (β Classical = 0.682; β Bayesian = 0.722; p − value < 0.05) and operational performance (β Classical = 0.444; β Bayesian = 0.498; p − value < 0.05) on financial performance are significant and positive. However, the impact of cost performance (β Classical = −0.522; β Bayesian = −0.607; p − value < 0.05) on financial performance is significant and negative. Moreover, the relationships between operational performance with both economic performance (β Classical = 0.418; β Bayesian = 0.439; p − value < 0.05) and cost performance (β Classical = 0.511; β Bayesian = 0.589; p − value < 0.05) is significant and positive, and the relationship between economic performance and cost performance is significant and negative (β Classical = −0.383; β Bayesian = −0.352; p − value < 0.05).    Based on Figure 5, the estimated structural equations were obtained that addressed the relationship between the performance index with economic performance, operational performance, and cost performance for the Classical-SEM and Bayesian-SEM, respectively given by Based on Figure 5, the estimated structural equations were obtained that addressed the relationship between the performance index with economic performance, operational performance, and cost performance for the Classical-SEM and Bayesian-SEM, respectively given bŷ

Results
Based on the estimated structural equations, it is concluded that economic performance (ξ1) has the most influence on financial performance (η) compared with the other latent variables. In brief, cost performance and operational performance are highly linked to financial performance, showing that better leadership strategies are more able to produce a high quality performance status. Table 4 presents the values of factor loading regulated coefficients and related standard errors for every indicator variable in the measurement equations acquired through both methods. It is clear from Table 4 that both models produced nearly the same factor loading estimates. It should also be mentioned that the indicator variables speculated as predictors are remarkably associated with their specific latent factors. It is highly significant that the parameter estimate standard errors under Bayesian-SEM are lower than under Classical-SEM. Moreover, Table 4 indicates that the 95% confidence intervals related to the parameters achieved with Bayesian-SEM are mainly shorter than that of the Classical-SEM-based parameters, which is not unordinary in light of the data provided by the prior distribution. This part of the study presents the comparison analysis between the Classical-SEM and Bayesian-SEM techniques in predicting the airline financial performance index. Four mathematical indices were applied to compare the Bayesian-SEM and Classical-SEM, which are representative of the strength and correctness of the prediction analysis. Root mean square error (RMSE), mean absolute percentage error (MAPE), coefficient of determination (R 2 ), and mean absolute error (MSE) are the most familiar indices for a comparison study among different prediction techniques. Table 5 presents the formula indices and output in traditional and Bayesian approaches. Table 5. Comparative outputs of Classical-SEM and Bayesian-SEM.

Index
Formula Bayesian-SEM Classical-SEM In the formulas which are mentioned in Table 5, y i is the ith real value of the dependent variable (y) and y i is the ith predicted value. The R 2 value for the Bayesian-SEM model is bigger than the Classical-SEM, and the RMSE, MSE and MAPE values of the Bayesian-SEM are smaller than the Classical-SEM. Therefore, the performance indices with the Bayesian-SEM are better in predicting financial performance than the Classical-SEM. The main reason Bayesian-SEM performs better is the defined traditional framework, which permits simultaneous self-adjustment of parameters and effective learning of the association between inputs and outputs in causal and complex models.

Discussion
The main purpose of the present study was to demonstrate the values of the Classical-SEM and Bayesian-SEM techniques in a new airline sustainability framework with the financial performance index for the Asia-Pacific airline industry. The new framework was developed based on previous studies in airline sustainability modeling. The designed model includes four latent variables and twelve familiar indicators. Even though much research has been conducted to determine the airlines' performance index and airline cost function, not many works have addressed modeling this index using SEM, particularly with interconnection among economic performance, operational performance and cost performance. This study was determined that economic performance should be considered as an independent latent variable in relationship to both operational performance and cost performance. This latent variable includes foreign direct investment, gross domestic products, and human development index. The combination of these three indicators into an economic performance index can directly and indirectly (through cost performance and operational performance) affect the airline financial performance index. Based on the data analysis output, it was found that economic performance has a significant effect on the financial performance index. These findings are similar to studies by Jenatabadi and Ismail [14] and Ismail and Jenatabadi [8] on airline performance modeling with Classical-SEM methodology and studies by Wang and Heinonen [72] who considered effective economic indicators as including gross domestic products and foreign direct investment in their research models. The significant interconnections of three main predictors are approved. It means the relationship between economic performance and operational performance, economic performance and cost performance, and operational performance and cost performance are significant. Moreover, the impact of every one of those three predictors on financial performance are significant. Therefore, the designed framework includes the impact of three predictors and the interconnection among them on financial performance is fitted to the current data.
The data were only collected in the Asia Pacific region, therefore, the findings cannot be generalized to all airlines in the world. However, the model proposed in this study has a potential capability and ability to be applied to airline companies. In this model, in order to enhance its efficiency, all redundant measures were eliminated or modified to be as close to the requirements of the airline industry as possible. The results of data analysis verified that the model consisting of four constructs was effective in understanding their role in predicting and estimating financial performance. The final model, which has a potential to be used in airline companies, is extremely close to the needs and the requirements of the industry as all redundant measures were eliminated and the most used and proper ones were added as measures and indicators.
This information should be helpful for managers and decision-makers to distribute capital resources logically upon implementing plans to improve the overall company performance. This information can be condensed in a single measure called a performance index, which is essential for detecting the indicators that could have an impact on it.
The other part of the main objective of this study was to illustrate Bayesian-SEM for analyzing the airline performance index. Along the lines of maximum likelihood and considering the Bayesian concept, the research parameters were defined as random with a prior distribution and prior density function [38]. After gathering data, the first phase in applying the Bayes theorem entailed combining these with prior distributions. The next phase was to calculate the posterior distribution, which reveals prior knowledge and empirical research data. By performing MCMC simulation, it was possible to summarize the joint posterior distribution with regard to lower dimensional summary statistics like posterior mean and standard deviation. Therefore, Bayesian application in SEM studies is more suitable for our research data.
A computational algorithm in Classical-SEM was determined based on a normality assumption and the sample covariance matrix of the research data. However, in many studies identified, multivariate normality was not the researchers' concern or the data did not have a normal structure. Therefore, researchers including Bashir and Schilizzi [73], Ansari, Jedidi [68], and Scheines, Hoijtink [39] considered that the Bayesian approach in SEM has the capability to overcome the non-normality concern.
Unit heteroscedasticity leads to damage of the homoscedastic error assumptions [74]. Heteroscedasticity is treated quite differently in the Bayesian context relative to maximum likelihood. Since our estimate of uncertainty is considered from the posterior distribution, we only have to be concerned about fittingly modeling the process to measure that distribution correctly. In a maximum based estimator, researchers would reweight the standard errors based on group size, however, within a Bayesian framework, the inferences on each parameter fully take into consideration the uncertainty of every other parameter interested. Therefore, as long as we have heterogeneity in the contributed research model, one usually ignores the idea of heteroscedasticity.

Conclusions
Based on the structure objective of the study we have two main conclusions. The first conclusion is about the main advantages of the current framework with previous studies in airline sustainability modeling. These advantages are:

•
This study developed Jenatabadi and Ismail [14] and Ismail and Jenatabadi [8] in estimating airline performance with application of Classical-SEM by interfering cost performance and considering operational performance separate from financial performance in different latent variables.

•
This study developed previous studies like Johnston and Ozment [9], Oum and Zhang [10], Zuidberg [7], and Hansen, Gillen [11] in airline sustainability performance in terms of cost function by leveraging country economic performance indicators.
The second conclusion is that based on the value of R 2 in Table 5, we can interpret from the Bayesian-SEM of 80.9% that variation in financial performance is related to economic performance, cost performance, and operational performance. However, this variation based on Classical-SEM is equal 67.9%. Moreover, based on the values of MPE, RMSE, and MSE, the amount of residuals of Bayesian-SEM are less than that of Classical-SEM. It means the predicted values of Bayesian-SEM are closer to the observed value than that of Classical-SEM. As a result, with higher values of R 2 and lower values of MPE, RMSE, and MSE, Bayesian-SEM have a better goodness of fit for the observations than that of the Classical-SEM.
This study does have its limitations from a practical standpoint of airline sustainability modeling structure. The main limitation is the exclusion of a few essential indicators such as ticket price, top management team, customer satisfaction, and governmental subsidy. Ticket price has been used in price modelling [75], choice modelling [76,77], and performance estimating [47]. For airline sustainability modelling, this indicator is suitable for the city pair level and not at the industrial level. For this reason, ticket price is also excluded from our research model. Information of top management team and customer satisfaction are also not available in the annual report and, therefore, are excluded from the research model. The governmental subsidy is financial support given to some airline companies by their governments [78]. The subsidy given by the government allows companies to give discounts to the ticket price for certain air travelers and has the potential to increase the load factor [79].
For future studies the following subjects are recommended: (1) Application of the framework which is introduced in Figure 2 in other areas and do comparison study between them. (2) Comparison study between Classical-SEM (as representative of parametric modeling) and Bayesian-SEM (as a representative semi-parametric modeling) was done, but comparison between these two methods with artificial neural networks (as a representative of non-parametric modeling) could be an interesting topic for future research. (3) Doing comparison analysis with the leveraging governmental subsidy as the main moderator in the research framework of Figure 2.