The Relationship of Causal Factors Affecting the Future Equilibrium Change of Total Final Energy Consumption in Thailand’s Construction Sector under a Sustainable Development Goal: Enriching the SE-VAR X Model

: This study aims to analyze the inﬂuence of the relationship between causal factors that affect the future equilibrium of the total ﬁnal energy consumption in the construction sector of Thailand under the sustainable development policy for the period of 10 years (2019–2028). This analysis was achieved with the application of the Structure Equilibrium-Vector Autoregressive with Exogenous Variables model (SE-VAR X model). This model was developed to ﬁll research gaps and differs from those of previous studies. In the selection of variables, the study focused on Sustainable Development (SD)-based variables available through the lens of Thailand. The exogenous variables included real GDP, population growth, urbanization rate, industrial structure, oil price, foreign direct investment, international tourist arrivals, and total exports and imports. Every variable had a co-integration at level (1) and was used to structure the SE-VAR X model. This particular model can effectively analyze the inﬂuence of the direct relationship and meet the criteria of goodness of ﬁt without spuriousness. This SE-VAR X model allowed us to discover that every variable in the model had an inﬂuence on the equilibrium change, where the real GDP is the fastest variable to adjust to the equilibrium while the total ﬁnal energy consumption has the slowest adjustment ability. The SE-VAR X model can be used to project the total ﬁnal energy consumption, as veriﬁed by the performance test result. The test was measured based on the Mean Absolute Percentage Error ( MAPE ) and Root Mean Square Error ( RMSE ), and their results were 1.09% and 1.01%, respectively. This performance result had the highest value compared to other models in the past. Thus, the SE-VAR X model is suitable for forecasting over the next 10 years (2019–2038). The results of this study reveal that the total ﬁnal energy consumption in the construction sector of Thailand will exhibit a continuously increasing growth rate from 2019 to 2028, amounting to about 144.29% or equivalent to 364.01 ktoe. In addition, the study also found that future government plans may be difﬁcult to achieve as planned. Therefore, the introduced model should be integrated into national development planning and strategies to achieve sustainable development in the future and to enable its application to other sectors.

no doubt about the ongoing impact on fossil fuel consumption due to such growth. This has led to crisis in its quantity and price. Most importantly, it continues to impact the environment negatively. If we look at a global scale, China and parts of Europe, particularly those in the Organization for Economic Co-Operation and Development (OECD) [1], their economies are growing quickly, while their growth rate is soaring upward at the same time. In addition, for the OECD group members, the economies of these countries are growing rapidly with an annual growth rate of 2% [1,2]. One noticeable example is that of India where the population continues to increase [2]. Therefore, the growth in the economy and population will affect the global energy consumption at a growth rate of 0.9% (2017/2016), in which the oil consumption rate will become dominant and high [1,2]. The ongoing increase in global energy demand has had a large impact on the ecosystem, contributing to the unwelcoming phenomenon of climate change. As the CO 2 emissions keep increasing, so does the greenhouse effect at the same time [3]. However, although many countries have implemented various policies and collaborated between one another to reduce greenhouse emissions, the issue remains unresolved and adds no value in terms of sustainability. As for Thailand, it is one of many other countries willing to cooperate and take initiative to address the problem. One such action is the country's participation and involvement in the United Nations Framework Convention on Climate Change, Kyoto Protocol, and Paris Agreement [1,3] in order to provide a better framework for achieving sustainable development.
Today, Thailand has put policies in place to accelerate economic development, aiming towards developing a country where the economic, social, and environmental development is emphasized based on national sustainable development policies and strategies. This serves for the short term of 1 to 5 years, the medium term of 6 to 10 years, and the long term of 11 to 20 years. The above development has been set forth through the promotion and acceleration of overall development in all aspects [4], including the promotion of foreign direct investment, urbanization, and industrial structure in all major areas. There is also foreign encouragement in the tourism industry and a focus on export activities. Additionally, the government promotes the increment of wages and intervenes in securing the oil and diesel prices, while a special fund has been established to support various industries. Comparing 2017 to 2016, with the abovementioned polices, Thailand has continued to grow with an economic growth rate of 4.5%. Moreover, the social growth rate was seen to rise to 2.59%. This success has been reported from 2000 to 2017 [4,5].
However, the economic and social development of Thailand does have an effect on the environment as well; the amount of CO 2 emissions from the national energy consumption increased by 1.5% in 2017 compared to 2016 [1,5]. In particular, the construction sector emits a high volume of CO 2 emissions. This reflects that the ecological system is affected, and greenhouse gas emissions continue to increase, mainly contributed by the construction sector. In addition, methane (CH 4 ) is ranked as having the second highest emission levels [3,5]. This impact deteriorates the environmental system, though the government has tried to implement measures and strategies to minimize the impact on the environment. In the context of Thailand, the government has undertaken many initiatives in terms of supportive policies to reduce energy consumption, be it through promotional or penalty policies. However, energy consumption and CO 2 emissions continue to increase regardless, especially energy consumption in the construction sector on what appears to be a daily basis. This reflects the historical absence in Thailand of an important and effective tool and measure to formulate accurate and effective policies. Based on the relevant studies available in Thailand, spurious models continue to be applied, and such applications result in the failing of the policies implemented in all application terms. Concerning the sustainable development goal of 1 to 5 years, to support the medium term of 6 to 10 years and the long term of 11 to 20 year in Thailand, the efforts to meet this goal seem unsuccessful due to limited implementation tools for effective forecasting that can be universally applied to different sectors. Thus, this study aims to fill this gap, so Thailand can perform better in national strategic planning. Therefore, we attempted to develop a Structure Equilibrium-Vector Autoregressive with Exogenous Variables (SE-VAR X ) model for such planning. After reviewing various relevant studies, we found that the SE-VAR X model is a model that differs from other models in the past, and it is an ideal model to be deployed in different contexts and sectors.
In this section, a number of well-known studies are revisited in order to investigate the causal relationships between various factors through different research methods and techniques. Also, in the establishment of this paper's model, many literature sources have been discussed to identify and support the model with relevant contributions in terms of the Structure Equation Model. Ashley and Tsang [6] proposed a testing procedure using cross-sample validation (CSV) to analyze the credible Granger-causality. The findings reveal that the price in regards to power, based on the in-sample Granger-causality F test, is controllable. As investigated by Yu and Yang [7], the causality between both linear and nonlinear process variables was studied by proposing a transfer entropy with an optimization of binary alarm series. The results of this study's analysis provide evidence of method suitability for industrial situations that are interfered with by noise. Ghafele and Gibert [8] analyzed the impact of fair use policy relative to the growth in private copying technology and the copyright markets of Singapore by referring to a method of differences-in-differences (DiD). This study affirmed a correlation between fair use policy with faster growth rates in the studied industries, but it did not show such a correlation with higher growth in the copyright industries.
Liu et al. [9] attempted to structure a method of Economy-Society-Ecology (ESE) subsystems by using a fuzzy analytical hierarchy process (FAHP), and entropy method (EM) to support the urbanization process of Chengdu, China. The study found that such a subsystem has a smaller effect on the coupling coordinated system while it is being enhanced with its ideal outputs for healthy urban development. Wang et al. [10] tried to justify the impact of changing urban forms on CO 2 emissions from commuting workers in Beijing by proposing an integrated analytical framework. The study's results illustrated the great impact of such a change in terms of increased transport CO 2 emissions. In addition, the change relative to a more decentralized urban area has a greater impact in the above matter, resulting in a large rise in CO 2 emissions.
Also, [11] investigated the connection between water use change, urbanization and economic growth in China from 1997 to 2013 by applying a co-integration test and a VECM (vector error correction model). This study provided evidence of the existence of a long-run equilibrium relationship among those indicators, while the short-run effects and Granger causal relationship were found to be insignificant. This implies that a water crisis is likely to arise as a result of urbanization and economic growth. In the same case study of China, [12] explored the link between economic, construction land use and ecology subsystems in Zhejiang from 2000 to 2012 by optimizing a Granger test and a coupling coordination model. The results of the study demonstrate that the changes in value of the economic subsystem are the cause of changes in the ecology and construction land use (CLU) subsystems, while the changes of value of the ecology and CLU subsystems are the root cause of each other. In addition, Liu et al. [13] explored the trends between Rural-to-Urban Migrants (RUMs) and Urban Land Expansion (ULE) in Hubei, China from 2009 to 2014 by applying a decoupling model. The study concluded that an irregularly distributed large percentage of the rural population was moving to urban areas, and that urban land use is both rising rapidly and is inefficient.
Ameyaw and Yao [14] analyzed the relationship between Gross Domestic Product (GDP) and CO 2 emissions from 2007 to 2014 in five West African countries. They also proposed a non-assumption driven forecasting technique by using a bidirectional long short-term memory (BiLSTM) sequential algorithm formulation to forecast the long-term total CO 2 emissions from 2015 to 2020. Their results indicated the presence of a unidirectional relationship between GDP and CO 2 emissions. Also, they were able to suggest sustainable policy targets based on CO 2 emissions projections. Zhixin and Xin [15] investigated the causal link between energy consumption and economic growth of Shandong, China, from 1980 to 2008, by implementing a Generalized Least Square (GLS) method, which concluded that a long-term relationship between these two factors exists, implying that the economic growth of the studied area is dependent on energy consumption. Du et al. [16] applied the Bernoulli-Nash Social Welfare Function (BNSWF) to examine the effects of the emission cap on emission-dependent manufacturers' decision-making. The study shows that the system-wide and manufacturer's profits are enhanced with such an emission cap, while the permit supplier is found to decline. However, it is suggested that the manufacturer and permit supplier can increase the profit through coordination in the supply chain.
In the USA, Yildirim et al. [17] explored the casual relationship between energy consumption and economic growth by using the Toda-Yamamoto procedure and bootstrap-corrected causality test. The study illustrated that there is a causal relationship between biomass-waste-derived energy consumption and real GDP, suggesting that the USA should focus on energy production from waste as an additional energy resource. Madlener and Sunak [18] applied different processes and mechanisms of urbanization to investigate the influences of urbanization on urban structures and energy demand. They found that urban energy demand has risen significantly, and the fuel mix tends to change because of different mechanisms of urbanization within various economic sectors. Al-mulali [19] analyzed the influence of nuclear energy consumption on GDP growth and CO 2 emissions by conducting a panel data analysis in 30 major nuclear energy-based countries from 1990 to 2010. The study concluded that there is a positive long-run effect of nuclear energy consumption on GDP growth, but not on CO 2 emissions. Furthermore, Cutlip and Fath [20] investigated causality between carbon emissions and economic development through a study of case studies in six countries. This investigation showed a connection between those variables.
Bernstein and Madlener [21] forecasted both short-and-long term electricity demand elasticities for the German manufacturing industry at different sub-sectoral levels through a co-integrated VAR method. As a result, they found a long-term relationship for five studied subsectors. Additionally, they concluded that a sub-sectoral approach is an ideal method for the estimation of electricity demand elasticities. In Lebanon, [22] assessed the connection between energy consumption and economic growth from 1980 until 2009 by deploying various causality tests, including the Hsiao, Toda-Yamamoto, and vector error correction-based Granger causality tests. Their study revealed a bidirectional short-and-long term relationship. This implies that energy is a vital source for the economic growth of Lebanon. In assessing the impact of energy consumption and CO 2 emissions on the development of economy and finance, Al-mulali and Che Sab [23] selected 19 countries with an application of a panel model covering the period of 1980 to 2008. The assessment provides evidence of energy consumption contributing to the studied development. However, the study found that the high levels of said development have increased CO 2 emissions.
In an attempt to comprehend the connection between urbanization and economic development in China, Chen et al. [24] took a sample of 31 provinces through a method of quadrant mapping. This study showed that the relationship between those two variables is similar to the Matthew effect. Also, at a provincial level, such a relationship exists with a significant difference. Jiao and Yang [25] analyzed the relationship between the economic development and carbon emissions in western China by implementing an econometric FGLS (Feasible generalized least squares) analysis. The study found the presence of a positive correlation between carbon emissions and industrial transformation and a rise of consumption as parts of economic factors. In particular, the significance of such a correlation is found after the implementation of the Western development strategy. Also, Liu and Lathrop [26] tried to detect newly urbanized areas by applying an artificial neural network. This case study produced a result of 20-30% accuracy with its method while reducing the computational cost. Williams, Hahs and Vesk [27] tried to discover the relationship between plant traits and the urban environment by doing a synthetization of study results and an examination of plant traits or niche indicators of urban flora. Their study found an increase of niche indicators for nutrients, temperature and alkalinity across many studies.
In the next line of studies, different studies applied different techniques, and structural equation modeling is an ideal method among them, enabling the success of this paper. Ye et al. [28] studied transfer commuting attitudes in Nanjing, China, by applying a method of market segmentation consisting of K-means clustering, factor analysis and structural equation modelling (SEM). The study found that sensitivity to time is found to be the vital factor affecting transfer willingness. Liu and Yang [29] tried to identify factors that affect behavioral intention (BI) toward the sharing economy by expanding on the technology acceptance model. Throughout their study, they found that perceived usefulness (PU) and perceived ease of use (PEOU) are observed to be the key factors affecting BI. Shih [30] did a study to identify the determinants of different factors: enterprises' strategic orientations, competitive advantages, radical innovation and business performance with a sample of 86 firms via regression analysis using the PLS SEM model. This study concluded that a firm's radical innovation capability and adoption can be enhanced through either market orientation factor or entrepreneurial orientation factor, while brand advantage has a significant impact on firms' performance. Also, [31] attempted to determine the determinants of effective Lean Manufacturing by applying the method of Covariance-based Structural Equation Modeling (CB-SEM) analysis. Their findings showed that those determinants can be used as guidelines to attain LM effectiveness, including a provision of necessary steps for successful lean implementation and achievement of lean cost reduction and time-saving. Grewal et al. [32] did a study to detect multicollinearity and measurement error in structural equation models in the theory testing by conducting a Monte Carlo simulation. Their analysis explained multicollinearity may become an error when the multicollinearity is extreme, between 0.6 and 0.8 and between 0.4 and 0.5. In Malaysia, Habidin et al. [33] conducted a confirmatory factor analysis for sustainable manufacturing practices in the automotive industry of Malaysia by applying the SEM technique. At the end of the research, they concluded that those automotive manufacturers are found to be intrigued by and involved in sustainability. They were observed to have a clear integration between sustainability and their business strategy.
Jiao et al. [34] deployed a novel approach for the estimation of sustainable urbanization performance in 31 provinces of China using SEM. Their findings suggest that effective policies and strategies can be improved through the assessments of the study in order to improve urban sustainability. Bandalos [35] examined the preciseness and utility of overall error and error estimators in the structural equation models by using the Monte Carlo method. This study showed that the rescaled non-centrality parameter and a sample estimator of overall error (EFo) produced a highly precise estimate of the approximation error and overall error amount, respectively. Shen et al. [36] produced a ranked list of the sustainable performance of urbanization in the prospect of global view between 111 countries, with data retrieved from the World Bank database and the United Nation database. Their study provides a referential outcome in which a difference in urbanization will contribute to a varied level of sustainable performance. Woo [37] revisited the international economic policy of China after three decades of its transformation and opening regionally and globally. His analysis shows evidence of negative spillovers due to large trade surpluses and greenhouse gases. Also, Liang et al. [38] developed a principal component analysis (PCA) and used Grey TOPSIS methodology to measure the sustainable urban development capacity of 13 cities of Jiangsu, China. In the study, an index system is established, and its measurement result has shown that all kinds of index weights in the measurement are able to increase in the comprehensive and persuasive assessment results.
Meulemans et al. [39] studied the occurrence of exercise dependence and eating disorders in the different genders of young adults by a Path analysis optimizing a Structural Equation Modeling (SEM). The findings have shown a small direct impact of gender on eating disorders. Han et al. [40] explored the factors affecting travel mode choice behavior in Jinan city of China by deploying a SEM-Logit Integration model. Their study found that gender, monthly income, trip purpose, travel distance, safety and convenience service are found to have a great impact on the selected traveling mode. Cheng et al. [41] proposed a structural equation model (SEM) to evaluate factors that influence passenger satisfaction towards the bus transfer process in China, and found bus fare and transfer distance are found to be the most influential factors for passenger satisfaction. Furthermore, Wang et al. [42] aimed to discover potential factor mechanisms of the long-distance travel mode choice by adapting a Structural Equation Model with a variance (ANOVA) approach. The study's findings indicate that all factors except gender, service demand for safety, and departure time have a significant impact on the travel mode choice. Furthermore, they also illustrated that travel distance is found to be the top most influential variable on the choice of travel mode, while the service demand, performance evaluation, and personal attributes are the next most influential factors. Geng et al. [43] examined the effect of multiple motivations, government procedures, and demographic features on residents' travel mode choice in Jiangsu Province, China, by structuring a multinomial logistic regression model with 1244 urban residents in total. Their analysis revealed that a pro-environmental motivation plays a great role in the promotion of green travel mode choices compared to car use, yet it is driven by self-interested motivations rather than pro-environmental motivation. Other variables, like gender, age, vehicle ownership, travel distance, income and government instruments have a great impact on the travelling mode choices.
Dobrovolskienė et al. [44] attempted to identify what enables project managers to achieve sustainability in their decision-making in the construction and automotive industry of Lithuania. The study shows that these managers are not that concerned with sustainability, and few sustainability criterions are put in place. However, their concern is greater than their team members. Prause and Atari [45] investigated the connections between networking, organizational development, structural frame conditions and sustainability under the concept of Industry 4.0 by applying a method of semi-structured expert interviews among a set of samples of high-tech manufacturing corporations in Estonia. Based on their findings, a poor organizational environment causes a negative effect on the performance of a networked construction, apart from technical matters. Additionally, Vegera et al. [46] conducted a study to determine the stages of the technological cycle and their accounting objects of the industrial waste. In their study, "Appearance" is found to be the first stage, while "Collection", "Accumulation", "Waste Preparation for Use", "Storage", "Use" and "Burial" are the stages that follow. In addition, the waste of ferrous and non-ferrous metals, construction waste, waste due to reservoirs cleaned and costs are the accounting objects found in the cycle. Yevhen et al. [47] studied the role of industrial ecology and logistics to achieve regional sustainable development goals. Upon their analysis, the integration of the environmental vision with economic content becomes the foundation of policy formulation and implementation to attain the above goals. Fomina et al [48] sought to understand the essence of industrial cooperation, and examined its influence on the growth of a sustainable economy by using a log-linear model. The investigation's results showed that an additional growth in the economy of about 2.3-3.0% is provided by the industrial cooperation. Meanwhile, such a system of cooperation helps boost the factor impact of economic growth determinants. Lastly, Oganisjana et al [49] analyzed three-stage research based on a social innovation project urging people to get involved in tackling social problems. As for the study's results, the intrapersonal, interpersonal and external factors were found to be motivational drivers for people to get involved in solving the problems.
The review of relevant research revealed that most modeling structures are no different from one another, and some models are buggy and produce potential errors if applied to other sectors or contexts. Thus, we introduced the SE-VAR X model to fill the above gap as well as to provide maximum efficiency and effectiveness in application. The SE-VAR X model was developed for medium-and long-term forecasting, yet was found to be the best model. We used the time series data for the period 1990-2017, and applied linear structural relations (LISREL, version 9.2, Chicago, IL, USA) software in later estimations. The research process is shown below.

1.
Determine all variables for the SE-VAR X model so as to analyze the influence of the relationship between causal factors affecting the equilibrium adjustment of the total final energy consumption based on the sustainable development policy for 10 years (2019-2028). This process applies Augmented Dickey-Fuller theory [50,51] along with the concept of walk drift and time trend at the first difference using data from 1990 to 2017.

2.
Take stationary variables at the first difference for a co-integration test based on the Johansen-Juselius concept [52,53].

3.
Construct the SE-VAR X model to analyze the influence of the relationship in both the long term and short term. Forecast the total final energy consumption by using the SE-VAR X model for the years 2019-2028, totaling 10 years. The flowchart of the SE-VAR X model is shown in Figure 1.  The remainder of this paper is as follows: Section 2 discusses the materials and methods. Section 3 shows the results. Section 4 summarizes the conclusion and discussion.

Stationary Features
The data used to structure the SE-VARX model have to be selective and identified as stationary. This selection is illustrated below.
If t Y is a stochastic variable time series and stationary, it must consist of three properties as shown in Equations (1)-(3). The remainder of this paper is as follows: Section 2 discusses the materials and methods. Section 3 shows the results. Section 4 summarizes the conclusion and discussion.

Stationary Features
The data used to structure the SE-VAR X model have to be selective and identified as stationary. This selection is illustrated below.
If Y t is a stochastic variable time series and stationary, it must consist of three properties as shown in Equations (1)-(3).
The equation with a mean has to be as follows: The equation with a variance becomes as follows: The equation with a covariance can be displayed as below: The height or spacing of covariance values is the same. The distance between the two Y t values does not depend on time. These three properties are referred to as stationary.
The stationary stochastic process or stationary is a time series of mean or expected values, as well as the values of a variance and a covariance (constant overtime). It does not depend on time, but rather on the distance or lag.
The Dickey-Fuller concept has been used to solve non-stationary data problems in the form of a difference regression called the first-order autoregressive process, via one of the following three equations: where ∆Y t represents the variables to be studied, δ is the coefficient of the lag, and ε t is the error term with normal distribution. The mean = 0, the variance = σ 2 for the hypothesis of the unit root test, and it is denoted as H 0 :δ = 0, non-stationary, H a :δ < 1, stationary. ∆Y t will be qualified as on-stationary when accepting H 0 and the variance δ is increased exponentially or explosively when the time increases. Details of each model can be shown as follows: Equation (7) is a simple model with a mean of zero suitable for data with a mean of only zero. The model was modified by adding a drift term in the following equation: where α 1 is the drift term. Equation (8) can also improve the unit root test as a component of Trend Stationary (TS) and Difference Stationary (DS) processes, as in the following equation: where T is the time trend, α 2 is the coefficient of the trend. ε t has a stationary feature with a mean of 0, and the variance is equal to σ 2 . However, causing the time series to make the first difference or higher difference will make the serial data available as a stationary through the difference stationary method, as follows: In Equation (10), the amount of time increases with the increasing level. If rejecting H a or accepting H 0 indicates that the data are non-stationary, then the coefficient δ = 0 is significant. This shows that the tau statistics of the coefficients calculated in the absolute term are less than the Dickey-Fuller (DF) critical in the absolute term. In this case, it shows that ε t has no white noise feature and it shows that the analyzed data have an autocorrelation feature where the Augmented Dickey-Fuller (ADF) Test must be used with the assumptions and goodness of fit, as also used in similar tests. The equation was modified by increasing the number of lagged variables of the dependent variable in a higher order to eliminate the problem of autocorrelation shown in the format as follows: Hypothesis of unit root test: H o :δ = 0, non-stationary, H a :δ < 1, stationary. At level, if reject H o that the data are stationary. The White Noise or time series of variables is stationary, which shows that the ∆Y t integrated at d or written in the symbol ∆Y t ∼ I(d).
where ∆ is the lagged values of the first difference of the variable. Equations (11)- (13) show that the most appropriate form of equation is the Augmented Dickey-Fuller, as shown in Equation (14).
However, the DF and ADF equations do not differ significantly, as the aim of the ADF is to solve the problem of error term and turn it into the feature of white noise or to have an error term of zero mean.

Co-Integration Test
This distance is stationary or symbolized as Y t − β 2 X t ∼ I(0), Y t − β 2 X t = u t and for time series Y t and X t has a stationary spacing model, which shows that the values of both variables are moving in one direction, not in different directions. This is because the two time series share a common stochastic trend, which we will call vector 1 −β 2 also written as [1 − β 2 ] showing a co-integrating vector.
Moreover, the equation Y t − β 2 X t = u t is referred to as a co-integrating regression. Y t and X t are in the following equations: where µ t is a random walk but is a random walk but ε Xt and ε Yt are white noise. We use Equations (15) and (16) to find the values of Y t and X t as follows: From Equation (18), when ε Yt is the White Noise, it is still White Noise. The time ε Xt series Y t − X t is a time series that is stationary where the vector showing the long-run equilibrium is 1 −1 or written as [1 − 1] . If Y t and X t have long-term equilibrium relationships or co-integration, they must have a coefficient β t which is not equal to zero, causing the following results. If the linear result of Y t and X t (β 1 Y t + β 2 X t ) is stationary, we call [β 1 − β 2 ] the co-integration vector, and we can find the error correction model as follows: where ϕ 1 (Y t−1 − X t−1 ) and Equilibrium Correction Mechanism (or Error Correction ϕ 2 (Y t−1 − X t−1 ) Mechanism) and Equations (19) and (20) are the short-term adaptive model for long-run equilibrium (Error Correction Model: ECM). Using the least squares method with the ECM model is not a spurious regression because all the dependent, independent and random variables are all stationary.

A Forecasting Model Using the SE-VAR X Model
Concerning the SE-VARx model, we adapted the concept of the basic models, including the Structure Equilibrium-Vector Autoregressive model with exogenous variables. This SE-VAR X model employs stationary variables at the first difference with a co-integration characteristic based on Johansen's concept in order to analyze both the short-and long-term relationships of each variable in the model. Here, each causal factor has a direct effect, and is later used to construct the SE-VAR X model as explained below.
where α 0 is constant, Y t−i represents the exogenous variables, which are stationary at level, and ECT t−i represents the error correlation mechanism test. Therefore, this study made the SE-VAR X model applicable to forecasting the total final energy consumption in the construction sector as demonstrated below.

Measurement of the Forecasting Performance
In the measurement of the model performance, we employed the concepts of Mean Absolute Percentage Error (MAPE) and Root Mean Square Error (RMSE). This was done by comparing the forecasting result and accuracy of each model [54,55]. The calculation equations are displayed below.

Screening of Influencing Factors for Model Input
The variables used in the model had to be screened and specifically selected. For the SE-VAR X model, we selected variables relevant to the sustainable development policy in the construction sector by optimizing the time series data of the period of 1990 to 2017. In this study, nine variables were taken for analysis: the total final energy consumption (FEC), real GDP (GDP), population growth (P g ), urbanization rate (U r ), industrial structure (I s ), oil price (O p ), foreign direct investment (F i ), international tourist arrivals (I t ), and total exports and imports (X m ). Each of those variables was stationary at the first difference, as verified by the Augmented Dickey-Fuller unit root test shown in Table 1. Note: ln(FEC) is the natural logarithm of total final energy consumption, ln(GDP) is the natural logarithm of real GDP, ln(P g ) is the natural logarithm of population growth, ln(U r ) is the natural logarithm of urbanization rate, ln(I s ) is the natural logarithm of industrial structure, ln(O p ) is the natural logarithm of oil price, ln(F i ) is the natural logarithm of foreign direct investment, ln(I t ) is the natural logarithm of international tourist arrivals, ln(X m ) is the natural logarithm of total exports and imports, *** denotes significance at α = 0.01, and ∆ is the first difference. Sources (own research). Table 1 indicates the result of the unit root test, and this affirms that the tau test value of all variables is greater than the MacKinnon critical value. This implies a stationary character at level (1) with a significance of 1% of those variables. Consequently, we subjected those variables to a co-integration test, as suggested by Johansen and Juselius, to analyze the long-term relationship of each variable at the same level, as illustrated in Table 2. Table 2 indicates the co-integration of every variable since the values of the trace test at the first difference are 225.11 and 70.09, which are higher than the MacKinnon critical values at significance levels of 1% and 5%, respectively. The maximum eigenvalue test results are 132.70 and 49.02, which are higher than the MacKinnon critical values at the same significance levels. Therefore, we subjected all those variables to a co-integration test, so we could proceed to the construction of SE-VAR X model as illustrated below. Note: *** denotes significance at α = 0.01. Sources: own research.

Formation of Analysis Modeling with the SE-VAR X Model
In the structuring of the SE-VAR X model, we used all the variables mentioned above to build the model, and tested them for goodness of fit. It was found that the model met all criteria of the test. In addition, a test on the spuriousness of the model gave normal results. The fulfillment of those tests enabled the model to analyze the magnitude of the direct effect the adjustment ability toward the equilibrium, the results of which can be seen in Table 3. Note: From the above information, ECT is the error correction test, ** denotes significance at α = 0.01, * denotes significance at α = 0.05, χ 2 /df is 1.25, root mean square of error approximation (RMSEA) is 0.03, root mean squared residual (RMR) is 0.007, goodness of fit index (GFI) is 0.94, adjusted goodness of fit index (AGFI) is 0.96, R-squared is 0.93, adjusted R-squared is 0.90, the Durbin-Waston statistic is 2.01, the F-statistic is 212.05 (probability is 0.00), the Autoregressive Condition Heteroscedastic test (ARCH test) is 29.08 (probability is 0.1), and the Lagrange multiplier test (LM test) is 1.50 (probability is 0.10). Sources: own research.
According to Table 3, the SE-VAR X model was found to have both short-term and long-term relationships, with a direct effect for independent variables. This changed for dependent variables, including ∆ ln(FEC), ∆ ln(GDP), ∆ ln(P g ), ∆ ln(U r ), ∆ ln(I s ), ∆ ln(O p ), ∆ ln(F i ), ∆ ln(I t ), and ∆ ln(X m ). However, this study disclosed that when there is a change in the real GDP (∆ ln(GDP)) by 1%, it will cause the total final energy consumption (∆ ln(FEC)) to change by 6.54% at a significance level of 1%. When the population growth (∆ ln(P g )) changes by 1%, it will then cause the total final energy consumption (∆ ln(FEC)) to change by 4.02% at a significance level of 1%. Also, when the urbanization rate (∆ ln(U r )) changes by 1%, the total final energy consumption (∆ ln(FEC)) will change by 5.33% at a significance level of 1%. When the industrial structure (∆ ln(I s )) changes by 1%, the total final energy consumption (∆ ln(FEC)) will change by 5.74% at a significance level of 1%. While the oil price (∆ ln(O p )) changes by 1%, the total final energy consumption (∆ ln(FEC)) will change by 5.13% at a significance level of 1%. When the foreign direct investment (∆ ln(F i )) changes by 1%, the total final energy consumption (∆ ln(FEC)) will change by 4.91% at a significance level of 1%. Furthermore, when the international tourist arrivals (∆ ln(I t )) rate changes by 1%, the total final energy consumption (∆ ln(FEC)) will change by 5.65% at a significance level of 1%. Lastly, when the total exports and imports (∆ ln(X m )) rate changes by 1%, the total final energy consumption (∆ ln(FEC)) will change by 3.91% at a significance level of 1%. Table 3 tells us that when each of the independent variables changes, there will be an effect of change observed in the dependent variables, including ∆ ln(GDP), ∆ ln(P g ), ∆ ln(U r ), ∆ ln(I s ), ∆ ln(O p ), ∆ ln(F i ), ∆ ln(I t ), and ∆ ln(X m ), which is shown in the relevant parameters.
As for the SE-VAR X model, it indicates an adjustment ability toward the equilibrium of the dependent variables, and can be presented as ECT t−1 . Here, the analysis verifies the variables with such ability to be the following. Firstly, the real GDP (∆ ln(GDP)) is equivalent to 54.5% at a significance level of 1%. Consequently, the next variables are foreign direct investment (∆ ln(F i )), industrial structure (∆ ln(I s )), urbanization rate (∆ ln(U r )), total exports and imports (∆ ln(X m )), international tourist arrivals (∆ ln(I t )), oil price (∆ ln(O p )), population growth (∆ ln(P g )), and total final energy consumption (∆ ln(FEC)), with the ability sizes of 52.2%, 52.1%, 41.9%, 40.7%, 33.9%, 31.1%, 30.7%, and 29.1%, respectively. Those results are produced at a significance level of 1%, except for population growth and total final energy consumption which are produced at a significance level of 5%.
In addition to this, we compared the performance of the SE-VAR X model with that of other models based on MAPE and RMSE. The other models included the ML model, ANN model, BP model, Grey model, and ARIMA model.  Table 4 shows that the SE-VAR X model produced an MAPE value of 1.09%, while its RMSE value was 1.01%. This further indicates that its performance peak is better compared to other studied models. In other lines of models, the ARIMA model had an MAPE value equivalent to 3.22% and its RMSE value was 3.04%. The Grey model presented an MAPE value of 5.37% and an RMSE value of 4.41%, while the BP model had an MAPE value of 11.59% and an RMSE value of 11.26%. In addition to this, the ANN model had an MAPE value of 12.64% and an RMSE value of 12.49%. Lastly, the ML model was found to have an MAPE value of 19.16% and an RMSE value of 17.52%. Hence, through this evidence, the researchers chose the SE-VARx model for a projection of the total final energy consumption over the next 10 years (2019-2028).

Total Final Energy Consumption Forecasting Based on the SE-VAR X Model
The SE-VAR X model was found to be suitable for such forecasting applications, and it was not spurious. Hence, it produces effective outcome with fewer errors, as shown in Figure 2. Figure 2 shows the forecasting result of the total final energy consumption in the construction sector of Thailand over the next 10 years (2019-2028). The study revealed that the total final energy consumption is increasing gradually with a growth rate of 151.04% (2028/2019). This implies that the total final energy consumption undergoes a great change in rate, which will affect the sustainable development policy of Thailand in the long run.
this evidence, the researchers chose the SE-VARx model for a projection of the total final energy consumption over the next 10 years (2019-2028).

Total Final Energy Consumption Forecasting Based on the SE-VARX Model
The SE-VARX model was found to be suitable for such forecasting applications, and it was not spurious. Hence, it produces effective outcome with fewer errors, as shown in Figure 2.  Figure 2 shows the forecasting result of the total final energy consumption in the construction sector of Thailand over the next 10 years (2019-2028). The study revealed that the total final energy consumption is increasing gradually with a growth rate of 151.04% (2028/2019). This implies that the total final energy consumption undergoes a great change in rate, which will affect the sustainable development policy of Thailand in the long run.

Conclusions and Discussion
This research focuses on the establishment of a new forecasting model, called the SE-VARX model. This model is shown to have good capacity for effective forecasting and can be used as a significant tool in policy-making to achieve sustainable development goals. In the study, each variable of the model is stationary at the first difference, yet co-integrated. The variables consist of total final energy consumption (FEC), real GDP (GDP), population growth (Pg), urbanization rate (Ur), industrial structure (Is), oil price (Op), foreign direct investment (Fi), international tourist arrivals ( t I ), and total exports and imports ( m X ). In addition, all those variables have both short-and long-term direct effects, and they are influential in the equilibrium change. Moreover, the SE-VARX model tells us that all dependent variables are adjustable to the equilibrium at different magnitudes. The real GDP ( GDP ) is the fastest variable to be adjusted the equilibrium. Meanwhile, foreign direct investment ( ln( ) ) were found to adjust to the equilibrium at slower pace. When the SE-VARX model was applied to predict a medium-term forecast of 10 years (2019-2028), the total final energy consumption was found to continuously increase with a growing rate, thus affecting the implementation of sustainable development policies in Thailand.

Conclusions and Discussion
This research focuses on the establishment of a new forecasting model, called the SE-VAR X model. This model is shown to have good capacity for effective forecasting and can be used as a significant tool in policy-making to achieve sustainable development goals. In the study, each variable of the model is stationary at the first difference, yet co-integrated. The variables consist of total final energy consumption (FEC), real GDP (GDP), population growth (P g ), urbanization rate (U r ), industrial structure (I s ), oil price (O p ), foreign direct investment (F i ), international tourist arrivals (I t ), and total exports and imports (X m ). In addition, all those variables have both short-and long-term direct effects, and they are influential in the equilibrium change. Moreover, the SE-VAR X model tells us that all dependent variables are adjustable to the equilibrium at different magnitudes. The real GDP (GDP) is the fastest variable to be adjusted the equilibrium. Meanwhile, foreign direct investment (∆ ln(F i )), industrial structure (∆ ln(I s )), urbanization rate (∆ ln(U r )), total exports and imports (∆ ln(X m )), international tourist arrivals (∆ ln(I t )), oil price (∆ ln(O p )), population growth (∆ ln(P g )), and total final energy consumption (∆ ln(FEC)) were found to adjust to the equilibrium at slower pace. When the SE-VAR X model was applied to predict a medium-term forecast of 10 years (2019-2028), the total final energy consumption was found to continuously increase with a growing rate, thus affecting the implementation of sustainable development policies in Thailand.
In line with the importance of the present study, we can see the pertinence of the proposed hypothesis. Besides, there is no similar model as such. The SE-VAR X model has been used in policy planning to support the sustainable development goals of Thailand in the past. Most models previously used are spurious, and this has caused errors and inaccuracies as they did not meet the criterion of the "Best Model". This also has a negative impact upon the policy implementation and strategic planning of the country [5]. Thus, the SE-VAR X model was made to close the research gap left by other previous research. Through a review of relevant studies, it was shown that the research flow of this model may be similar to others, yet it is different in terms of the direct effect analysis of each causal factor. Furthermore, it provides evidence of adjustment ability to the equilibrium, thus making it a suitable model for forecasting applications. Most importantly, the SE-VAR X model is unique in its ability to be applied to different sectors and contexts. Also, it effectively supports the national policy-making and planning efforts and can be used as a guideline for knowledge discovery. In the construction of the SE-VAR X model, we employed advanced statistics and verified that it was the best model by a comparison with other previously reported models. The SE-VAR X model is an ideal model for medium-term forecasting (6-10 years) [4,5], which differs greatly from short-term models, especially in the field of policies, as these are extremely changeable over time. This factor affects future forecasts, mostly medium-and long-term forecasting [5]. Hence, it is crucial to have the right tool, such as the proposed SE-VAR X model, for the implementation of policy and strategic planning in Thailand.
Before ending this note, we affirm the methods used to build this SE-VAR X model. As it has been mentioned, only stationary causal factors at the same level were used, and the co-integration test was employed along the analysis of direct effect, which separates this model from those investigating indirect effect. No other research along these lines has been done before. Since the SE-VAR X model is a non-spurious model, it allowed us to attain a more accurate magnitude of the influence and confirmed that this model could significantly contribute to the national policy-making endeavors. It is important to confirm that the SE-VAR X model is suitable for forecasting, thus enabling Thailand to formulate proper policies in the future.
The limitation of this research is that the factors for the SE-VAR X model were determined and selected carefully in consideration of the Nationall Strategic Framework of Sustainable Development from 1990 until the present. This is also based on the Thai context, with some factors potentially differing from previous studies. Thus, results may differ if we were to use the same factors from those studies for the SE-VAR X model. Hence, researchers should select the appropriate factors for the appropriate context and sector, so as to produce the most accurate analytical results for effective policy making for national governance in the future. In addition, due to the rapid economic and social growth in Thailand resulting from the promotional policies in place, such growth inevitably leads to a continuous negative impact on the environment. Although the government is trying to tackle the issue of environmental damages, it still cannot solve such issues completely. This is because certain factors do not work in line with the market mechanism, including the factor of government intervention in the price of oil, and this will result in the distortion of real market trends. In addition, the government has been involved in the promotion of exporting energy to foreign countries for the purposes of national revenue generation. Yet, this forces Thailand to import certain energy types from abroad. Thus, this activity may influence the total exports and imports in the long term. Therefore, it is essential to consider the context of the whole country so that sustainable development can be achieved efficiently.
Concerning future contributions, the researchers wish to further this SE-VAR X model in real-world applications. By enhancing the model, forecasting at maximal efficiency can be realized. Using a method of factor analysis along with path analysis, this approach will facilitate the accurate measurement of variable selection for a given context and enable us to effectively deploy latent variables and observed variables to forecast the changes of government policies in the future.
Author Contributions: J.S. and K.K. were involved in the data collection and preprocessing phase, model constructing, empirical research, results analysis and discussion, and manuscript preparation. All authors have approved the submitted manuscript.
Funding: This research received no external funding.