
Citation: Hedi Januar, Dwiyitno, Izhamil Hidayah, Irma Hermana. Seasonal heavy metals accumulation in the soft tissue of Anadara granosa mollusc form Tanjung Balai, Indonesia[J]. AIMS Environmental Science, 2019, 6(5): 356-366. doi: 10.3934/environsci.2019.5.356
[1] | Enahoro A. Iboi, Oluwaseun Sharomi, Calistus N. Ngonghala, Abba B. Gumel . Mathematical modeling and analysis of COVID-19 pandemic in Nigeria. Mathematical Biosciences and Engineering, 2020, 17(6): 7192-7220. doi: 10.3934/mbe.2020369 |
[2] | Pannathon Kreabkhontho, Watchara Teparos, Thitiya Theparod . Potential for eliminating COVID-19 in Thailand through third-dose vaccination: A modeling approach. Mathematical Biosciences and Engineering, 2024, 21(8): 6807-6828. doi: 10.3934/mbe.2024298 |
[3] | Hamed Karami, Pejman Sanaei, Alexandra Smirnova . Balancing mitigation strategies for viral outbreaks. Mathematical Biosciences and Engineering, 2024, 21(12): 7650-7687. doi: 10.3934/mbe.2024337 |
[4] | Antonios Armaou, Bryce Katch, Lucia Russo, Constantinos Siettos . Designing social distancing policies for the COVID-19 pandemic: A probabilistic model predictive control approach. Mathematical Biosciences and Engineering, 2022, 19(9): 8804-8832. doi: 10.3934/mbe.2022409 |
[5] | Weike Zhou, Aili Wang, Fan Xia, Yanni Xiao, Sanyi Tang . Effects of media reporting on mitigating spread of COVID-19 in the early phase of the outbreak. Mathematical Biosciences and Engineering, 2020, 17(3): 2693-2707. doi: 10.3934/mbe.2020147 |
[6] | Jie Bai, Xiunan Wang, Jin Wang . An epidemic-economic model for COVID-19. Mathematical Biosciences and Engineering, 2022, 19(9): 9658-9696. doi: 10.3934/mbe.2022449 |
[7] | Durward Cator, Qimin Huang, Anirban Mondal, Martial Ndeffo-Mbah, David Gurarie . Individual-based modeling of COVID-19 transmission in college communities. Mathematical Biosciences and Engineering, 2022, 19(12): 13861-13877. doi: 10.3934/mbe.2022646 |
[8] | Adil Yousif, Awad Ali . The impact of intervention strategies and prevention measurements for controlling COVID-19 outbreak in Saudi Arabia. Mathematical Biosciences and Engineering, 2020, 17(6): 8123-8137. doi: 10.3934/mbe.2020412 |
[9] | Mario Santana-Cibrian, Manuel A. Acuña-Zegarra, Jorge X. Velasco-Hernandez . Lifting mobility restrictions and the effect of superspreading events on the short-term dynamics of COVID-19. Mathematical Biosciences and Engineering, 2020, 17(5): 6240-6258. doi: 10.3934/mbe.2020330 |
[10] | Jin Guo, Aili Wang, Weike Zhou, Yinjiao Gong, Stacey R. Smith? . Discrete epidemic modelling of COVID-19 transmission in Shaanxi Province with media reporting and imported cases. Mathematical Biosciences and Engineering, 2022, 19(2): 1388-1410. doi: 10.3934/mbe.2022064 |
Energy is an integral part of every country's economic development. Electricity is an incredibly flexible energy source that fuels the efficiency of almost any sector of the economy. Pakistan's disastrous energy failure policies have left the nation with a severe power catastrophe, suffering economic damages throughout the last twenty years. Hence, understanding the energy consumption and real GDP is very important to improve its financial growth and development process from a policy perspective.
Over three decades after [1], seminal work reaffirming the causal link between energy consumption and real production, the causal path between these factors remains uncertain [2,3]. Several earlier empirical studies discovered that electricity consumption Granger-causes real GDP, such as [4,5,6], although other studies claimed that electricity consumption does not Granger-causes real GDP [7,8]. On the other hand, some studies (e.g., [9,10,11,12] specify the number of explanations why recent research on the link between electricity and economic growth has contradictory findings. According to these researches, the existing electricity demand model focused only on observed variables, economic development, and energy price, and neglects the non-observable elements. Electricity consumption often depends on non-observable variables such as energy conservation, technical advances, and customer preference. Because electricity demand is considered a derived output, energy consumption relies on electrical appliances usage and capital stock. Appropriately, such findings indicate other exogenous factors aside from economic variables such as income and price (e.g., energy consumption factors, technological change, and customer preference) should be incorporated for modeling electricity demand. argued that capturing the impact of un-observable factors through a simple timing pattern is impractical [13]. Dimitropoulos et al. (2004), Dilaver and Hunt (2011), Hunt and Ninomiya (2003) and Hunt (2003) stated that the method Structural Time Series Modeling (STSM) was established by [15,16] is an adequate approach to capture the results of non-observable variables (trends in energy usage, technical changes, and customer fondness) effectively [9,12,13,14]. Furthermore, Harvey et al. (2004) claimed that time series is broken down into non-observed patterns and other abnormal elements in the STSM technique and that these are not identified [17]. The stochastic process, also regarded as the Underlying Energy Demand Trend (UEDT), incorporates not just exogenous technological advancement but other significant social and economic implications as for customers’ expectations, resource-saving scientific changes.
The uninterrupted provision of electricity is the critical factor to flourish in the industrial sector. It plays an indispensable role, as this sector is the engine of economic growth in Pakistan. However, Pakistan has faced a severe electricity crisis, directly and indirectly, involving all economic sectors, particularly concerning the changing energy structure. In the past, energy shortage and inefficiencies’ posed damages to the economy. The government launched various plans during 2013–2018 to fix the energy crisis for ensuring the smooth distribution of energy supply. These schemes added a total power of 12,230 MW to the national energy system [18]. Electricity consumption patterns evidenced that there was no substantial change in electricity demand. The household’s electricity consumption pattern indicated a slight decrease from 51% in 2018 to 48% in 2019. The industrial sector also revealed an increasing electricity consumption trend of 25% in 2018 and 27% in 2019; however, commercial consumption remained unchanged from 2018 to 2019. The agriculture sector consumption also declined by 1% from 10% in 2018 to 9% in 2019. The surge in industrial sector electricity consumption shows a positive sign that market was recovering from the earlier economic damages but still needs improvement in other sectors. Figure 1 shows a comparison of electricity consumption patterns from July 2018 to March 2019 fiscal year.
Several past studies evaluated the indispensable role of electricity consumption in Pakistan and applied linear and nonlinear conceptual models to report agammaegate demand for electricity consumption and quantify it [3,6,19,20,21]. The past literature findings specified that real wages, real substitute and actual electricity rates, temperature, and stochastic patterns appear to be critical determinants of Pakistan's agammaegate and sectoral electricity demand. The factors mentioned above evidenced that Pakistan has to develop an appropriate and well-designed economic model linked to electricity rate, consumption, and economic growth indicators. Figure 2 indicates the rising trends in electricity consumption and GDP while prices fluctuate from 1998 to 2018 due to government subsidy for the industrial boost in Pakistan.
The study's inspiration and contribution consist of a review of Pakistan's energy demand and supply and its effects on Pakistan's economy. Comparison to prior research, the present study provides fascinating and real contributions. First, we find limited studies that address all the variables related to the industrial sector in Pakistan. Second, the current research explores the relationship between electricity consumption, price, and real GDP of Pakistan's industrial sector using data from the 1960–2018 time series. Third, the study employs the unit root test, Johansen co-integration, Vector Error Correction Model (VECM), impulse response function, and variance decomposition tests, and offers more fruitful inferences. By comparison, within a VECM system, we research the variables' relationships. It helps us determine whether there is a long-term (or co-integrating) link between the variables and investigate their existence. Fourth, the study results would help formulate energy generation and consumption policies while also articulating an ambitious power rationing strategy to minimize economic loss due to energy shortages. Finally, similar studies are based on flawed data, but specific data from Pakistan's industrial sector is considered in this study. This novelty has created a gap between earlier research.
The remainder of the study organized as follows: the next section 2 presents the literature review, and section 3 parades econometric material and methods. Section 4 shows the results and discussion, while section 5 provides conclusions and policy implications.
Several studies have been done over the last few years, examining the relationship between an economy's electricity demand and economic progress. As of now, the relationship between energy use and financial growth is recorded significantly by [4,22,23,24]. Thus, our analysis estimates electricity consumption (EC) prices and gross domestic product (GDP) from the literature context. Followed Chen and Fang (2018) investigated the co-integrating and Granger causal relationships between economic development, industrial energy use, and human resources in China using a panel of 210 prefectural cities for the 2003–2012 period [25], also Cialani and Mortazavi (2018) examined the energy demand and its determinants after the liberalization of the power markets in 29 European countries [26]. Based on panel data from 1995–2015 for these countries and using a dynamic partial adjustment model, market elasticity is calculated for both residential and industrial demand for electricity. They discovered short-run and long-run causation of EC to GDP and inferred that a lack of energy could harm economic growth (EG). The causality evidence found from energy use to EG confirms the result of [28,29,30,31], whereas the causality findings from economic progress to electricity use are aligned with [27,28,29,30,31]. Moreover, another critical study by [32] enquired about the effect of power consumption, using yearly data for the period 1971–2012, on GDP, labor, and capital stocks. The Structural Vector Self-Regression (SVAR) method analysis shows that labor demand rises with economic growth. Conversely, Nazlioglu et al. (2014) inspected the causal association between Turkey's energy and financial development between 1967 and 2007 [33]. They use three investigative econometrics methods to achieve this objective: the bounds testing method to co-integration, the linear and the nonlinear Granger causality test. The analysis of co-integration indicates that in the long run, co-integrated relationship. Further, Al-bajjali and Yacoub (2018) inspected determinants of electricity consumption from 1986–2015, focused on six independent variables by employed VECM for analysis; the results confirmed that GDP, urbanization, and agammaegate water consumption, are essential and positively related electricity consumption [34].
Similarly, Ghafoor et al. (2016) studied that energy shortages have caused economic losses of between 2% and 3% of Pakistan's GDP [35]. Further, Sidique (2014) used the Granger causality test of Hsiao and reported that the effect of EC on EG is strong and highly meaningful [36]. However, Shahbaz and Hooi (2012) utilized the Autoregressive Distributed Lag (ARDL) model and found the bidirectional causality between EC and EG [37]. Additionally, Javid and Qayyum (2014) utilized the structural time series and found a growing EC pattern in the rural, industrial, and housing sectors [38]. Jamil and Ahmad (2010) presented Johansen cointegration and VECM Granger causality tests and noticed one-way causation flowing from EG to EC [39]; Khan and Abbas (2016) implemented a panel cointegration study, and calculated electricity demand in agricultural and industrial sectors are more resilient to reform [40]. The Autoregressive Integrated Moving Average (ARIMA) models were applied by [41] and observed that EC has a clear interaction with EG. They stated that the ARIMA model was the best in the forecast and indicated that schemes' productivity would be improved by increasing energy usage.
Lin and Yousaf (2021) measured the electricity consumption in Pakistan between 1989 and 2018 using the index decomposition. The findings indicate that the economic structure's impact was the key driving force in Pakistan's increasing overall electricity usage [42]. Recently Lin and Yousaf (2020) implemented the MARKAL model and measured the supply, demand, and diversity of primary energy supplies [43]. The findings suggested that the availability of primary energy decline slightly. These analyses investigated the correlation between EG and EC, the short-term causal relationship, long-run, and regression-related tests among the parameters. The result showed that the scarcity of electricity could damage the economy of Pakistan. The outcomes are varied, according to early findings. In particular, all research on Pakistan's electricity use and economic growth relationships indicates no consent on the direction of causality among these variables. The EC, therefore, plays a vital role in Pakistan's cultural identity of economic and financial growth. Structural transition and with up-gradation might be a little more specific in this case. Therefore, the characteristics and distinctions in energy use within the multiple segments need to be addressed. Lin and Yousaf (2020) indicated that it is possible to achieve optimum production by rising capital and energy technologies [44]. The different studies' findings are inconsistent, or even substantially contradictory, partially due to variations in methods or data length. They reviewed previous studies and their findings from the causal relationship between energy consumption and economic growth reported in Table 1. Whereas we find various studies evaluating the causal association between energy and economic development in Pakistan, no research explicitly examined the relationship between energy consumption, price, and GDP growth in the industrial sector. However, these variables are used in different aspects, methodology, and duration. This research extends the relevant literature on electricity by defining the industrial sector causality path between electricity consumption, prices, and real GDP.
Author | Country (Period) | Variables | Methodology | Direction of Causality | |||
→ | ← | ↔ | ≠ | ||||
[45] | China (1971–2001) | EC, GDP | Granger Causality |
✔ | ✘ | ✘ | ✘ |
[37] | Pakistan (1972–2009) | EC, RGDP | ARDL, VECM | ✘ | ✔ | ✘ | ✘ |
[46] | Switzerland (1950–2010) | EC, RGDP | ARDL, UECM | ✘ | ✘ | ✔ | ✘ |
[46] | United Kingdom(1975–2010) | EC, RGDP | Panel bootstrap causality | ✔ | ✘ | ✔ | ✘ |
[46] | OECD(1990–2008) | EC, GDP | Panel Granger causality, | ✔ | ✘ | ✔ | ✘ |
[46] | Lebanon (1983–2014) | EC, GDP | Toda-Yamamoto | ✘ | ✘ | ✔ | ✘ |
[46] | Algeria (1971–2010) | EC, GDP | Cointegration, VECM | ✘ | ✘ | ✔ | ✘ |
[3] | Pakistan (1972–2010) | EC, RGDP | Granger Causality | ✔ | ✘ | ✘ | ✘ |
[47] | Poland (2000–2012) | EC, GDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[48] | Turkey (1970–2011) | EC, GDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[19] | Malaysia (2005–2010) | EC, RGDP | Multiplier Approach | ✔ | ✘ | ✘ | ✘ |
[24] | Algeria(1980–2012) | RE, NRE | Granger causality | ✔ | ✘ | ✔ | ✘ |
[49] | Taiwan (1998–2014) | EC, RGDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[50] | U.S.(2005–2015) | EC, RGDP | Wavelet | ✔ | ✘ | ✘ | ✘ |
[51] | Ghana (1971–2014) | EC, Industrial Growth | ARDL, ECM | ✘ | ✔ | ✘ | ✘ |
[52] | Pakistan(1981–2017) | RE, NRE | ARDL, VECM | ✘ | ✘ | ✔ | ✘ |
[53] | Pakistan (1970–2018) | EC, EP, GDP | Granger causality VECM | ✔ | ✘ | ✘ | ✘ |
Note: →, ←, ↔, and ≠ denotes unidirectional bidirectional and no causality, respectively. |
Therefore, this study aims to capture a comprehensive view, based on the most up-to-date database, of the differentiation indices of electricity consumption (EC), price, and GDP of the industrial sector in Pakistan from 1970 to 2018. Nevertheless, no consideration was given to studying the modeling techniques of electricity consumption of Pakistan in general. We may examine the efficacy of the past policies by observing the EC, price, and GDP link. In formulating future policies, a greater understanding of how the EC is evolving in Pakistan is also crucial. This analysis allows us to use the VECM in combination with the variance decomposition method to determine the influence of selected variables that affect EC in Pakistan over the time defined.
The research investigates the relationship between electricity consumption, electricity prices, and the GDP of Pakistan's industrial sector. The study used time series secondary data based on yearly observation, covering five decades approximately from 1970–2018, exhibits in Table 2. The electricity will measure in Gigawatt hour (GWh). Because electricity has been a public enterprise in Pakistan, instead of being regulated by the market, the electricity price is cross-subsidizing in all sectors. We use the average price of electricity, while real GDP uses a proxy of industrial value-added taken by the world bank's most authentic source, National Transmission & Despatch Company (NTDC) and Government of Pakistan (GOP) official websites. Figure 3 depicts the diagram of the variables.
Variables | Data Source | Scale Unit |
Industrial electricity consumption | [54] | Gigawatt hours (GWh) |
Electricity prices | [55] | Millions |
Real gross domestic product | [56] | $Millions |
The cointegration test allows the series under examination to be non-stationary in the same integration order to evaluate the long-term relationship between variables. The model forecasts that price is a significant factor in demanding electricity. Our empirical research implements the extensively employed method to maximum likelihood by Johansen for the analysis. The first step is to check all the variables for their stationarity., and two following tests are widely prominent [57,58], to detect the integration in all series, which stated in Eq 1:
Δyt=β0+δYt−1+γ1Δyt−1+γ2yt−2+.........γpΔyt−p+μt | (1) |
where, yt denotes a series and ut expresses error terms. Appropriate lags of Δyt are incorporated the whiten the errors. According to the Schwarz Bayesian criterion (SBC), the lag length is chosen, later checking for the serial correlation of first and higher-order in residuals. For the H0 null hypothesis test in Eq 1 is ‘d,’ which = 0, against the one-tailed, and alternative, which is negative. The stationarity of yt would not be rejected If δ result results significantly negative. Modeling associations among non-stationary features essentially prerequisite their differencing to make stationarity. For several years, most of the long-run economic relationship lost due to differences; therefore, extensive data needed at level variable for the conservation of the long-run. Meanwhile prevents chosen variables from being spuriously regressed. A long-run equilibrium relationship among non-stationary time series data exists as indicated by economic theories. If variables are I(1), so the co-integration approach will be applied in the long-run. Augmented Dickey-Fuller (ADF) and Phillips–Perron (PP) unit root tests are the primary step to move into the co-integration model.
Following [59] and [60], a demand-side model was used to investigate the industrial sector's reaction to unanticipated electricity consumption shock electricity prices and GDP. Impulse response function assessment ensures a systematic long-run relationship to determine the underlying predictor and does not wander much farther away from one another. The Johansen-Juselius cointegration technique is favored over the Engle and Granger two-step cointegration process for estimating long-term relationships in multivariate condition. Since if we regress Y on X in the former model, the outcomes would vary from regressing X on Y. In comparison, the Engle-Granger technique is a two-step procedure and the second step would therefore be influenced by an error happening in the first step. The technique proposed by [61] for maximum likelihood is reflected in the following VAR model:
Xt=A1Xt−1+A2Xt−2+...+ApXt−p+εt | (2) |
where Xt is an (n × 1) economic time series vector. A1, A2, ..., the (n × n) coefficient matrices are represented by Ap and εt is an (n × 1) vector with zero mean and constant variance error terms. We notice two test statistics for the estimation of cointegrating vectors, i.e., likelihood ratio or trace statistics, under the Johansen-Juselius methodology, and maximal eigenvalue statistics, as follows:
λtrace(r)=−Tn∑i=r+1ln(1−∧λr)λmax(r,r+1)=−Tln(1−∧λr+1) | (3) |
trace statistics are used to test the null hypothesis that the number of cointegrating vectors is less or equal to r, whereas max statistics evaluate the null hypothesis that cointegrating vectors are equal to r. These measures are focused on drawing parallels, and both λtrace and λmax statistics would be weak if characteristic roots are near to 0, and there will be less support for long-term series relationships.
One of the Granger interpretation theorem's ramifications is that Granger causality occurs at least in one direction if the series has a cointegrating relationship, which indicates if a conditional variable strengthens the predictor variables prediction when shown in the model [61,62]. The popular meaning of lagging independent terms is investigated under the Vector Autoregression (VAR) concept of Granger causality. The causal relationship is calculated under a VECM that differs between short-and long-term casualties while the sequence is co-integrated. Important terminology for error correction leads to long-term Granger causality, while the general sense of lagged independent variables is an example of short-term Granger causality. In addition, the corrections to the long-run equilibrium suggested in Eqs 2–4 are summarized:
ΔEC=α1+l∑i−1β1iΔECt−i+m∑i=1γ1iΔEPt−1+n∑i=1δ1iΔGDPt−1+φ1ECTr,t−1+μ1t | (4) |
ΔEP=α2+l∑i−1β2iΔECt−i+m∑i=1γ2iΔEPt−1+n∑i=1δ2iΔGDPt−1+φ2ECTr,t−1+μ2t | (5) |
ΔGDP=α3+l∑i−1β3iΔECt−i+m∑i=1γ3iΔEPt−1+n∑i=1δ3iΔGDPt−1+φ3ECTr,t−1+μ3t | (6) |
where EC,EP,andGDP represents electricity consumption, electricity price, and real GDP. Correspondingly, ‘α’ is the intercept ‘n’ is the number of lags, ‘Δ’ is the 1st difference, the joint consequence of lags β1i,γ1i,δ1i, β2i,γ2i,δ2i and β3i,γ3i,δ3i in Eqs 2–4, respectively, while ui,t For (i=1,2,3) are residuals, and the Error Correction Terms (ECT) specified by ECTr,t−1 and 'ϕ' is the adjustment of model speed towards equilibrium. For example, scale and the statistical significance of the one-period lag ECTr,t−1 coefficient determines how quick the disequilibrium in EC,EPandGDP are corrected to return to the equilibrium. All the variables are used in the natural logarithmic form. A series of short-run changes slowly rectify the divergence from the long-run equilibrium. The size and statistical significance of the ECT is a measure of the degree to which the left side variable returns to its long-run equilibrium in reaction to the random shocks of each equation for each short-run cycle, while the error correction model via the ECT provides another method for the detection of Granger causality that can be avoided in standard Granger and Sim study. The actual regression coefficients are generally hard to read in a VAR system. These are then represented by techniques connected with VAR models, i.e., impulse response functions and decomposition methods of variance. Impulse response function in earlier researches has been utilized widely for policy evolution. Energy sector reform appraisal analyses include [63,64,65,66]. Causation from EPtoEC and GDPtoEP can be tested similarly from Eqs 4–6, respectively. As explained above, whereas the existence of causality indicated by co-integration, VECM confirms the direction of causality among the variables. In addition, in [67], the error correction model makes it possible to differentiate between long-term and short-term causality in contrast to returning the missing information to the system due to the terminology of error correction. Figure 4 demonstrates the graphical representation of the study framework.
Firstly, the stationarity check of all variables to be conducted, which is necessary to prevent the spurious regression. The H0 considered as non-stationary series (has unit root). All variables should be stationery earlier containing in the model. Across all variables, the ADF and PP tests suggest the presence of unit root at a 5% significance level, indicating that any statistical inferences in levels from the series would be invalid [68]. Though the non-stationarity hypothesis was discarded in its first difference for all the series, electricity consumption, electricity prices, and real GDP were integrated with order one or I(1) as a whole shown in Table 3. All variables are integrated with order one, a necessary condition of the process variables co-integration relationship.
Industrial Variables | Philips-Perron (PP) | Augmented Dickey-Fuller (ADF) | Order of Integration | ||
Levels | First difference | Levels | First difference | ||
EC | (0.32) | (0.00)* | (0.35) | (0.00)* | I(1) |
EP | (0.49) | (0.00)* | (0.51) | (0.00)* | I(1) |
GDP | (0.31) | (0.00)* | (0.29) | (0.00)* | I(1) |
Notes: *, ** asterisk specifies significance at 1% and 5% level, respectively. p-value denotes in () parentheses. |
After the first condition has been met, we step into the next process to determine whether there is a long-run relationship among the variables. Therefore, Johansen and Juselius (1990), Johansen (1988) maximum likelihood method applied to find the co-integration, which consists of the two following estimations: trace(λ−trace)andmaximumeigenvalue(λ−max)statistics [61,62]. The H0 is rejected against the H1 at a 5% level. Hence, it is confirmed that there is a long-run relationship exists in two variables in line with [34,53]. The illustration of unrestricted co-integration outcomes shown in Table 4.
Hypothesized | r = 0 | r ≤ 1 | r ≤ 2 |
Trace Statistics | 35.40 | 12.83 | 3.29 |
P-Value | (0.01)* | (0.003)* | (0.12) |
Max-Eigen Statistics | 22.57 | 9.54 | 3.33 |
P-Value | (0.03)* | (0.002)* | (0.07) |
Notes: The numbers in the parentheses are p-values, r: indicates the number of co-integration hypothesis relationships, *, ** asterisk implies significance at 5% and 10% level, correspondingly. |
We proceed to the VECM estimate, as defined by Eqs 4–6, drawing inferences for each variable within the sample on the direction of causality and the exogenity or endogeneity. The result of the VECM is exhibited in Table 5. All variables are projected at a one-period lag. The term error-correction depends on the previous period deviation from long-run equilibrium (the error) affects the explanatory variables short-run dynamics. Hence, the coefficient of ECT, ϕ, is the speed of adjustment; it measures the speed at which explained variables return towards equilibrium after a change in the other variables [69].
Industrial Sector | VECM (Long-run relationship effects) | ECT(t-Stats) | |||||
Coefficient(p-value) | Standard Error | Coefficient(p-value) | Standard error | Coefficient (p-value) | Standard error | ||
ΔEC | 1 | - | −0.19 | −0.39 | 0.15 | −0.21 | (−2.66)* |
(0.03)* | (0.01)* | ||||||
ΔEP | −0.13 | −0.12 | 1 | - | 0.05 | −0.12 | (−3.97)** |
(0.04)* | (0.12) | ||||||
ΔGDP | 0.09 | −0.14 | −0.16 | −0.24 | 1 | - | (−3.36)* |
(0.001)* | (0.00)** | ||||||
Notes: The number in the parentheses is the p-value. *, ** asterisks indicate significance at 5% and 1% level, respectively. |
ECT's coefficient is negative in all variables and found significant at the 5% level, where a 1% increase in electricity price decreases electricity consumption by 0.13%, as well a 1% increase in GDP increase EC by 0.09%. While a 1% upsurge in electricity consumption decreases electricity price by 0.19%, if GDP increase by 1%, then electricity price decreases by 0.16%. However, a 1% increase in electricity consumption increases GDP by 0.15%, and electricity prices have an insignificant effect on GDP. The study in line with [37,70,71,72] reported similar results in Pakistan's case between energy use and economic development and energy consumers, except exogenous energy prices. Additionally, our results are similar to those of [73] for China and [74] for Turkey. These findings indicate that the key factors for Pakistan's economic progress are energy use and price. The overall vector error-correction models confirm that the industrial sector energy management could be viable. Since our study is limited to the use of electricity instead of the total use of electricity, Our findings do not merely contradict with the results in the [39] and [75] that shows energy shortages in Pakistan could affect economic growth, as causality varies from short- and long-term energy consumption to GDP. Our findings also vary from those reported in [76], which also considers the existence of unidirectional causality at the agammaegate level from electricity consumption to GDP. There could be several factors for the discrepancy in their outcomes, such as various data periods included in the research and econometric methods Granger causality form of Hsiao. Our findings show that as real economic growth rises, the electricity demand increases and distributors need to boost their capacity to maintain an adequate supply. While the increase in price decreases electricity consumption, that can decrease the productivity of the industrial sector. However, electricity consumption reflects positively in economic growth that could be considered for the policymakers.
In the vector error correction model, the IRF tracks the predictor variables' response to an unexpected shock or novelty. In this study, we calculated IRFs for sector-specific EC, EP, and GDP. IRF outcomes for the industrial sector is seen in Figure 5.
To determine how a shock affects alternative variables and how long the effect lasts, we use the generalized variable impulse responses proposed by [77] and [78] for ten years, the impulse responses from the variables given. Illustration Figure 5 shows shock in EC, EP, and GDP separately. If one standard deviation is positively shocked at electricity price, then electricity consumption is continually rising, and the GDP decline declined about two years afterward, it is continuously increasing and shows a positive long-run relationship. However, in response to electricity price positive shock increasing electricity consumption two years, then even continuously, electricity price positive shock increases GDP for three years then stable, which shows the long-run relationship. If one standard deviation shock in electricity consumption, then GDP upwards about five years after that constantly. However, the shock in electricity price declined GDP about seven years later; it seems steady.
Until now, the research has been limited to in-samples. Furthermore, to assess the comparative importance of the causality check's findings, we are now carrying out variance decomposition analyses. Because it is well renowned, the variance decomposition effects rely on ordering the variables necessary to recognize the structural VAR from the projected modified form VAR approach. In literature, the normal method is to propose the most credible order driven by economic concept and confirm the outcomes' robustness by inverting the ordering see, for example, [39,79]. For short, we are considering three years and for the long-run eight years in all models.
Table 6 depicts the findings of variance decompositions for the model (EC, EP, and GDP). In the short-run, electricity consumption analysis causes a 96.07% variation in electricity consumption shock fluctuation. Further shock to electricity price influences 3.66% in electricity consumption fluctuation, whereas GDP affects 0.27% in response to electricity consumption change. To check the long-run shock to electricity consumption causes 87.97% variation of the fluctuation is self-explained—however, electricity prices and GDP effects by 11.78% and 0.25%, respectively.
Period | Variance Decomposition of EC | |||
S.E. | EC | EP | GDP | |
1 | 0.026154 | 100 | 0 | 0 |
2 | 0.040085 | 96.92866 | 2.657988 | 0.41335 |
3 | 0.055235 | 96.07497 | 3.657699 | 0.267331 |
4 | 0.066524 | 95.30433 | 4.487052 | 0.208615 |
5 | 0.076296 | 93.71978 | 6.059193 | 0.221029 |
6 | 0.084754 | 91.77491 | 7.982465 | 0.242629 |
7 | 0.092929 | 89.84477 | 9.908605 | 0.246629 |
8 | 0.100791 | 87.9735 | 11.77734 | 0.249163 |
9 | 0.108377 | 86.54688 | 13.1834 | 0.269722 |
10 | 0.11552 | 85.69038 | 13.99968 | 0.309941 |
Table 7 illustrates the outcomes of variance decompositions for the model (EP, EC, and GDP). The short-run shock to electricity price causes a 44.97% variation of the electricity price fluctuation is generally self-explained. Further shock to electricity consumption influences 54.86% in the fluctuation of electricity prices. At the same time, GDP affects 0.17% to electricity prices. The long-run shock to electricity price causes a 24.15% variation of the fluctuation is self-explained—however, electricity consumption and GDP effects by 71.39% and 4.46%, respectively.
Period | Variance Decomposition of EP: | |||
S.E. | EC | EP | GDP | |
1 | 0.048996 | 56.37411 | 43.62589 | 0 |
2 | 0.080503 | 47.46051 | 52.45839 | 0.081095 |
3 | 0.111427 | 54.85629 | 44.96984 | 0.173865 |
4 | 0.136671 | 58.00428 | 41.5015 | 0.49422 |
5 | 0.158198 | 61.49309 | 36.99908 | 1.507833 |
6 | 0.17542 | 65.03976 | 32.20885 | 2.751394 |
7 | 0.190837 | 68.51778 | 27.79294 | 3.689283 |
8 | 0.20503 | 71.39459 | 24.14619 | 4.459221 |
9 | 0.218298 | 73.62662 | 21.30353 | 5.069856 |
10 | 0.23014 | 75.36009 | 19.18527 | 5.454636 |
Table 8 demonstrates the effects of variance decompositions for the model (GDP, EP, and EC). The short-run shock to GDP causes 95.87% variation of the fluctuation in GDP, usually is self-explained. Further shock to electricity consumption influences 1.78% in the fluctuation of GDP. However, electricity price affects 2.35% to GDP. However, the long-run shock to GDP causes 84.54% variation of the fluctuation is self-explained. But electricity consumption and electricity prices effects by 5.25% and 10.21% respectively.To summarize, the by-and-large findings of the variance decomposition endorse the results achieved by the VECM model. The primary effect that energy usage is influenced by electricity prices, whereas output and electricity consumption also positively affect real GDP.
Period | Variance Decomposition of GDP: | |||
S.E. | EC | EP | GDP | |
1 | 0.027913 | 7.08547 | 1.074455 | 91.84007 |
2 | 0.047862 | 2.562744 | 1.192117 | 96.24514 |
3 | 0.058702 | 1.780112 | 2.351816 | 95.86807 |
4 | 0.067107 | 2.679764 | 3.801723 | 93.51851 |
5 | 0.076742 | 4.208902 | 5.118915 | 90.67218 |
6 | 0.086717 | 4.775009 | 6.861578 | 88.36341 |
7 | 0.095723 | 5.065867 | 8.767336 | 86.1668 |
8 | 0.103643 | 5.248982 | 10.21453 | 84.53649 |
9 | 0.110786 | 5.330424 | 11.07375 | 83.59582 |
10 | 0.117053 | 5.316104 | 11.42037 | 83.26352 |
Pakistan's electricity generation from 1987 to 2019 was revised monthly at an average of 5,955 GWh. In 2018, the data hit an all-time peak of 16,946 GWh and a historic low of 2,313 GWh in 1989. The findings indicate that model coefficients met all predictions, and the electricity price is a reliable mediating component for the industrial sector that revealed a 1% change in electricity price decreases electricity consumption by 0.13%. The study found that as the rise in industrial electricity prices raised firms' energy costs, it could push firms to pursue technological innovation and enhance production performance. By reducing their conflicting duties, the institutional setup of energy divisions needs to be standardized. For instance, the NEPRA electricity regulator's role to advise the government on secure and reliable ways to meet electricity demands at manageable rates should be fulfilled as per law. As well a 1% change in GDP increase EC by 0.09%. Disagammaegated power use and economic development analyses would make it more useful for authorities to devise a robust energy-saving and environmental degradation strategy. The findings show the industries in which economic development over long periods is related to energy consumption. A 1% upsurge in electricity consumption decreases electricity price by 0.19%, which is quite useful for policymakers that managing industrial electricity prices can enhance manufacturing productivity. These prediction findings were correlated with related studies in the existing literature to demonstrate the accuracy of the present analysis. The findings differ from various studies focused on energy factors using different models and periods, such as [12,80,81,82,83,84,85]. In Pakistan's case, many studies concluded differently, such as [3,88–92]; however, none of them focused on these factors. The novel findings can be beneficial for the sustainable energy policy as per Pakistan's industrial sector's potential.
The current electricity provision rates in Pakistan are considered unsustainable due, in particular, to poor infrastructure. The widespread blackouts throughout the country have incredibly adverse effects on Pakistan's economic growth and social lives. The current study contributes to the emerging energy literature by investigating the causal link in the industrial sector in Pakistan between electricity consumption, prices, and GDP, by using a neoclassical development system. The empirical analysis was conducted using data for the period 1970–2018. The indication of co-integration between variables showed the long-run relationship. It ensures a short-term fluctuation in electricity consumption and GDP, which ultimately returns to long-run equilibrium. Most of the coefficients lagged in our analysis are significant at a 5% level except for GDP electricity prices. The outcomes show a causal relationship from electricity consumption to price and GDP that impact is long-term in the industrial sector, which is the foremost determinant of Pakistan's economic change.
Pakistan is a country where a lack of electricity and the electricity sector work bare capacity margin. Based on our findings obtained in this study. More interest should be given to policy recommendations as follows: (a) 1% rise in the price of energy reduces electricity usage by 0.13%, and a 1% increase in GDP raises EC by 0.09%. The result implies that the key driver is the energy price, which indicates an inverse impact on electricity usage that is visa-versa with the costs of industrial production. Also, an increase in EC was positively affecting GDP. The government should invest more in power projects to maintain energy demands and supply, which is beneficial for economic growth and boosts industrial productivity. (b) While a 1 percent rise in energy demand reduces electricity by 0.19%, the price of electricity reduces by 0.16 percent if GDP increases by 1 percent. (c) However, an increase in energy consumption by 1% raises GDP by 0.15%. The empirical results further suggest that electricity prices are the main factor in the use of industrial electricity. The findings suggest that officials should focus on relevant legislation, subsidy reforms, and tax relief in the industrial sector. Moreover, it will attract the business community; foreign investors also provide employment opportunities and boost Pakistan's GDP growth. As a whole, electricity consumption is essential for long-term economic stability and growth in the industrial sector. We suggest that infrastructure planning and investments are vital to meeting the growing electricity demand. Whereas any policies to ease the electricity market should stimulate productivity, the invention in generation and supply, resulting in lesser prices as autonomous suppliers, can stabilize the fluctuation that may positively affect the long term. Finally, we propose that planning and investing in infrastructure growth is essential to meeting the growing electricity demand.
We have some limitations concerning our studies. Based on the findings of multiple literature pieces, this study selects energy consumption, price, and GDP as useful for the industrial sector. This research was carried out in the industrial sector's specific context, so other Pakistan sectors need to be investigated.
The authors declare no conflict of interest.
[1] |
Casas S, Bacher C (2006) Modelling trace metal (Hg and Pb) bioaccumulation in the Mediterranean mussel, Mytilus galloprovincialis, applied to environmental monitoring. J Sea Res 56: 168-181. doi: 10.1016/j.seares.2006.03.006
![]() |
[2] |
Jan A, Azam M, Siddiqui K, et al. (2015) Heavy metals and human health: mechanistic insight into toxicity and counter defense system of antioxidants. Int J Mol Sci 16: 29592-29630. doi: 10.3390/ijms161226183
![]() |
[3] |
Batvari BPD, Sivakumar S, Shanthi K, et al. (2016) Heavy metals accumulation in crab and shrimps from Pulicat lake, north Chennai coastal region, southeast coast of India. Toxicol IndHealth 32: 1-6. doi: 10.5487/TR.2016.32.1.001
![]() |
[4] |
Zhang G, Bai J, Zhao Q, et al. (2016) Heavy metals in wetland soils along a wetland-forming chronosequence in the Yellow River Delta of China: Levels, sources and toxic risks. Ecol Indic 69: 331-339. doi: 10.1016/j.ecolind.2016.04.042
![]() |
[5] |
Sun Z, Mou X, Tong C, et al. (2015) Spatial variations and bioaccumulation of heavy metals in intertidal zone of the Yellow River estuary, China. Catena 126: 43-52. doi: 10.1016/j.catena.2014.10.037
![]() |
[6] |
Elturk M, Abdullah R, Rozainah MZ, et al. (2018) Evaluation of heavy metals and environmental risk assessment in the mangrove forest of Kuala Selangor estuary, Malaysia. Mar Pollut Bull 136: 1-9. doi: 10.1016/j.marpolbul.2018.08.063
![]() |
[7] |
Tjahjono A, Suwarno D (2018) The spatial distribution of heavy metals lead and cadmium pollution and coliform abundance of waters and surface sediment in Demak. J Ecol Eng 19: 43-54. doi: 10.12911/22998993/89715
![]() |
[8] |
Prartono T, Sanusi HS, Nurjaya IW (2016) Seasonal distribution and geochemical fractionation of heavy metals from surface sediment in a tropical estuary of Jeneberang River, Indonesia. Mar Pollut Bull 111: 456-462. doi: 10.1016/j.marpolbul.2016.06.106
![]() |
[9] |
Analuddin K, Sharma S, Septiana A, et al. (2017) Heavy metal bioaccumulation in mangrove ecosystem at the coral triangle ecoregion, Southeast Sulawesi, Indonesia. Mar Pollut Bull 125: 472-480. doi: 10.1016/j.marpolbul.2017.07.065
![]() |
[10] |
Rouane-Hacene O, Boutiba Z, Belhaouari B, et al. (2015) Seasonal assessment of biological indices, bioaccumulation and bioavailability of heavy metals in mussels Mytilus galloprovincialis from Algerian west coast, applied to environmental monitoring. Oceanologia 57: 362-374. doi: 10.1016/j.oceano.2015.07.004
![]() |
[11] | Azizi G, Akodad M, Baghour M, et al. (2018) The use of Mytilus spp. mussels as bioindicators of heavy metal pollution in the coastal environment, A review. J Mater Environ Sci 9: 1170-1181. |
[12] | Dabwan AH, Taufiq M (2016) Bivalves as bio-indicators for heavy metals detection in Kuala Kemaman, Terengganu, Malaysia. Indian J Sci Tech 9(9). DOI: 10.17485/ijst/2016/v9i9/88708. |
[13] | El-Moselhy KM, Saad EM, El-Shaarway RF, et al. (2016) Assessment of heavy metals pollution using sediments and bivalve Brachidontes variabilis as bioindicator in the Gulf of Suez, Egypt. Int J Mar Sci 6: 1-12. |
[14] |
Azizi G, Layachi M, Akodad M, et al. (2018) Seasonal variations of heavy metals content in mussels (Mytilus galloprovincialis) from Cala Iris offshore (Northern Morocco). Mar Pollut Bull 137: 688-694. doi: 10.1016/j.marpolbul.2018.06.052
![]() |
[15] |
Zhang G, Bai J, Xiao R, et al. (2017) Heavy metal fractions and ecological risk assessment in sediments from urban, rural and reclamation-affected rivers of the Pearl River Estuary, China. Chemosphere 184: 278-288. doi: 10.1016/j.chemosphere.2017.05.155
![]() |
[16] |
Riani E, Cordova MR, Arifin Z (2018) Heavy metal pollution and its relation to the malformation of green mussels cultured in Muara Kamal waters, Jakarta Bay, Indonesia. Mar Pollut Bull 133: 664-670. doi: 10.1016/j.marpolbul.2018.06.029
![]() |
[17] |
Xue S, Shi L, Wu C, et al. (2017) Cadmium, lead, and arsenic contamination in paddy soils of a mining area and their exposure effects on human HEPG2 and keratinocyte cell-lines. Environ Res 156: 23-30. doi: 10.1016/j.envres.2017.03.014
![]() |
[18] |
Liu Q, Liao Y, Shou L (2018) Concentration and potential health risk of heavy metals in seafoods collected from Sanmen Bay and its adjacent areas, China. Mar Pollut Bull 131: 356-364. doi: 10.1016/j.marpolbul.2018.04.041
![]() |
[19] | Mok JS, Kwon JY, Son KT, et al. (2015) Distribution of heavy metals in internal organs and tissues of Korean molluscan shellfish and potential risk to human health. J Environ Biol 36: 1161-117. |
[20] |
Sudsandee S, Tantrakarnapa K, Tharnpoophasiam P, et al. (2017) Evaluating health risks posed by heavy metals to humans consuming blood cockles (Anadara granosa) from the Upper Gulf of Thailand. Environ Sci Pollut Res 24: 14605-14615. doi: 10.1007/s11356-017-9014-5
![]() |
[21] |
Liu J, Cao L, Dou S (2017) Bioaccumulation of heavy metals and health risk assessment in three benthic bivalves along the coast of Laizhou Bay, China. Mar Pollut Bull 117: 98-110. doi: 10.1016/j.marpolbul.2017.01.062
![]() |
[22] |
Kucuk YK, Topcu A (2017) Ecological Risk Assesment and Seasonal-Spatial Distribution of Trace Elements in the Surface Sediment of Trabzon Harbour, Turkey. Open J Ecol 7: 348-363. doi: 10.4236/oje.2017.75025
![]() |
[23] |
Campillo JA, Santos-Echeandía J, Fernández, B (2019) The hydrological regime of a large Mediterranean river influences the availability of pollutants to mussels at the adjacent marine coastal area: Implications for temporal and spatial trends. Chemosphere 237: 124492. doi: 10.1016/j.chemosphere.2019.124492
![]() |
[24] |
Siregar TH, Priyanto N, Putri AK, et al. (2016) Spatial distribution and seasonal variation of the trace hazardous element contamination in Jakarta Bay, Indonesia. Mar Pollut Bull 110: 634-646. doi: 10.1016/j.marpolbul.2016.05.008
![]() |
[25] | Hammer O, Harper DAT, Ryan PD (2001) Past: palaentological statistics software package for education and data analysis. Palaentol Elec 4: 9. |
[26] |
Caccia VG, Millero FJ, Palanques A (2003) The distribution of trace metals in Florida Bay sediments. Mar Pollut Bull 46: 1420-1433. doi: 10.1016/S0025-326X(03)00288-1
![]() |
[27] |
Wang XC, Feng H, Ma HQ (2007) Assessment of metal contamination in surface sediments of Jiaozhou Bay, Qingdao, China. CLEAN Soil Air Water 35: 62-70. doi: 10.1002/clen.200600022
![]() |
[28] |
Kang M, Tian Y, Peng S (2019) Effect of dissolved oxygen and nutrient levels on heavy metal contents and fractions in river surface sediments. Sci Total Environ 648: 861-870. doi: 10.1016/j.scitotenv.2018.08.201
![]() |
[29] | Bosch AC, O'Neill B, Sigge GO, et al. (2016. Heavy metals in marine fish meat and consumer health: a review. J Sci Food Agri 96: 32-48. |
[30] |
Rajeshkumar S, Liu Y, Zhang X (2018) Studies on seasonal pollution of heavy metals in water, sediment, fish and oyster from the Meiliang Bay of Taihu Lake in China. Chemosphere 191: 626-638. doi: 10.1016/j.chemosphere.2017.10.078
![]() |
[31] |
Shakouri A, Gheytasi H (2018) Bioaccumulation of heavy metals in oyster (Saccostrea cucullata) from Chabahar bay coast in Oman Sea: Regional, seasonal and size-dependent variations. Mar Pollut Bull 126: 323-329. doi: 10.1016/j.marpolbul.2017.11.012
![]() |
[32] |
Kibria G, Hossain MM, Mallick D, et al. (2016) Trace/heavy metal pollution monitoring in estuary and coastal area of Bay of Bengal, Bangladesh and implicated impacts. Mar Pollut Bull 105: 393-402. doi: 10.1016/j.marpolbul.2016.02.021
![]() |
[33] | Le QD, Bach LG, Arai T (2015) Monitoring Heavy Metal Contamination Using Rocky Oyster (Saccostrea glomerata) in Haiphong-Halong Coastal Area, North Vietnam. Int J Environ Res 9: 1373-1378. |
[34] | Ali MM, Ali ML, Islam MS, et al. (2016) Preliminary assessment of heavy metals in water and sediment of Karnaphuli River, Bangladesh. Environ Nano, Monitor, Manag 5: 27-35. |
[35] |
Bhardwaj R, Gupta A, Garg JK (2017) Evaluation of heavy metal contamination using environmetrics and indexing approach for River Yamuna, Delhi stretch, India. Water Sci 31: 52-66. doi: 10.1016/j.wsj.2017.02.002
![]() |
[36] |
Gawad SSA (2018) Concentrations of heavy metals in water, sediment and mollusk gastropod, Lanistes carinatus from Lake Manzala, Egypt. Egyptian J Aquatic Res 44: 77-82. doi: 10.1016/j.ejar.2018.05.001
![]() |
[37] | Waykar B, Petare R (2016) Studies on monitoring the heavy metal contents in water, sediment and snail species in Latipada reservoir. J Environ Biol 37: 585-589. |
[38] | Dallinger R (1993) Strategies of metal detoxification in terrestrial invertebrates, in: Dallinger R, Rainbow PS, Ecotoxicology of Metals in Invertebrates, Lewis Publishers, Boca Raton, 245-289. |
[39] |
Alam L, Mohamed CAR, Mokhtar MB (2012) Accumulation pattern of heavy metals in marine organisms collected from a coal burning power plant area of Malacca Strait. Sci Asia, 38: 331-339. doi: 10.2306/scienceasia1513-1874.2012.38.331
![]() |
[40] | Ziyaadini M, Yousefiyanpour Z, Ghasemzadeh J, et al. (2017) Biota-sediment accumulation factor and concentration of heavy metals (Hg, Cd, As, Ni, Pb and Cu) in sediments and tissues of Chiton lamyi (Mollusca: Polyplacophora: Chitonidae) in Chabahar Bay, Iran. Iranian J Fish Sci 16: 1123-1134. |
[41] |
Wang X, Liu X, He Y, et al. (2016) Seasonal variations and health risk of heavy metals in the muscle of Crucian carp (Carassius auratus) cultured in subsidence ponds near Suzhou, Eastcentral China. Expos Hea 8: 79-91. doi: 10.1007/s12403-015-0184-4
![]() |
[42] | Azizi G, Akodad M, Baghour M, et al. (2018) The use of Mytilus spp. mussels as bioindicators of heavy metal pollution in the coastal environment, A review. J Mater Environ Sci 9: 1170-1181. |
[43] |
Casas S, Bacher C (2006) Modelling trace metal (Hg and Pb) bioaccumulation in the Mediterranean mussel, Mytilus galloprovincialis, applied to environmental monitoring. J Sea Res 56: 168-181. doi: 10.1016/j.seares.2006.03.006
![]() |
[44] |
Rivera-Hernández JR, Fernández B, Santos-Echeandia J, et al. (2019) Biodynamics of mercury in mussel tissues as a function of exposure pathway: natural vs microplastic routes. Sci Total Environ 674: 412-423. doi: 10.1016/j.scitotenv.2019.04.175
![]() |
[45] |
González-Fernández C, Albentosa M, Campillo JA, et al. (2015) Influence of mussel biological variability on pollution biomarkers. Environ Res 137: 14-31. doi: 10.1016/j.envres.2014.11.015
![]() |
1. | Bożena Gajdzik, Radoslaw Wolniak, Influence of the COVID-19 Crisis on Steel Production in Poland Compared to the Financial Crisis of 2009 and to Boom Periods in the Market, 2021, 10, 2079-9276, 4, 10.3390/resources10010004 | |
2. | Ahmed S. Elgazzar, Simple mathematical models for controlling COVID-19 transmission through social distancing and community awareness, 2021, 0, 1865-7125, 10.1515/znc-2021-0004 | |
3. | Tahajuddin Sk, Santosh Biswas, Tridip Sardar, The impact of a power law-induced memory effect on the SARS-CoV-2 transmission, 2022, 165, 09600779, 112790, 10.1016/j.chaos.2022.112790 | |
4. | Toshikazu Kuniya, Structure of epidemic models: toward further applications in economics, 2021, 72, 1352-4739, 581, 10.1007/s42973-021-00094-8 | |
5. | Xiaoying Wang, Qing Han, Jude Dzevela Kong, Studying the mixed transmission in a community with age heterogeneity: COVID-19 as a case study, 2022, 7, 24680427, 250, 10.1016/j.idm.2022.05.006 | |
6. | Daisuke Fujii, Taisuke Nakata, Takeshi Ojima, Martial L Ndeffo-Mbah, Heterogeneous risk attitudes and waves of infection, 2024, 19, 1932-6203, e0299813, 10.1371/journal.pone.0299813 |
Author | Country (Period) | Variables | Methodology | Direction of Causality | |||
→ | ← | ↔ | ≠ | ||||
[45] | China (1971–2001) | EC, GDP | Granger Causality |
✔ | ✘ | ✘ | ✘ |
[37] | Pakistan (1972–2009) | EC, RGDP | ARDL, VECM | ✘ | ✔ | ✘ | ✘ |
[46] | Switzerland (1950–2010) | EC, RGDP | ARDL, UECM | ✘ | ✘ | ✔ | ✘ |
[46] | United Kingdom(1975–2010) | EC, RGDP | Panel bootstrap causality | ✔ | ✘ | ✔ | ✘ |
[46] | OECD(1990–2008) | EC, GDP | Panel Granger causality, | ✔ | ✘ | ✔ | ✘ |
[46] | Lebanon (1983–2014) | EC, GDP | Toda-Yamamoto | ✘ | ✘ | ✔ | ✘ |
[46] | Algeria (1971–2010) | EC, GDP | Cointegration, VECM | ✘ | ✘ | ✔ | ✘ |
[3] | Pakistan (1972–2010) | EC, RGDP | Granger Causality | ✔ | ✘ | ✘ | ✘ |
[47] | Poland (2000–2012) | EC, GDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[48] | Turkey (1970–2011) | EC, GDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[19] | Malaysia (2005–2010) | EC, RGDP | Multiplier Approach | ✔ | ✘ | ✘ | ✘ |
[24] | Algeria(1980–2012) | RE, NRE | Granger causality | ✔ | ✘ | ✔ | ✘ |
[49] | Taiwan (1998–2014) | EC, RGDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[50] | U.S.(2005–2015) | EC, RGDP | Wavelet | ✔ | ✘ | ✘ | ✘ |
[51] | Ghana (1971–2014) | EC, Industrial Growth | ARDL, ECM | ✘ | ✔ | ✘ | ✘ |
[52] | Pakistan(1981–2017) | RE, NRE | ARDL, VECM | ✘ | ✘ | ✔ | ✘ |
[53] | Pakistan (1970–2018) | EC, EP, GDP | Granger causality VECM | ✔ | ✘ | ✘ | ✘ |
Note: →, ←, ↔, and ≠ denotes unidirectional bidirectional and no causality, respectively. |
Industrial Variables | Philips-Perron (PP) | Augmented Dickey-Fuller (ADF) | Order of Integration | ||
Levels | First difference | Levels | First difference | ||
EC | (0.32) | (0.00)* | (0.35) | (0.00)* | I(1) |
EP | (0.49) | (0.00)* | (0.51) | (0.00)* | I(1) |
GDP | (0.31) | (0.00)* | (0.29) | (0.00)* | I(1) |
Notes: *, ** asterisk specifies significance at 1% and 5% level, respectively. p-value denotes in () parentheses. |
Hypothesized | r = 0 | r ≤ 1 | r ≤ 2 |
Trace Statistics | 35.40 | 12.83 | 3.29 |
P-Value | (0.01)* | (0.003)* | (0.12) |
Max-Eigen Statistics | 22.57 | 9.54 | 3.33 |
P-Value | (0.03)* | (0.002)* | (0.07) |
Notes: The numbers in the parentheses are p-values, r: indicates the number of co-integration hypothesis relationships, *, ** asterisk implies significance at 5% and 10% level, correspondingly. |
Industrial Sector | VECM (Long-run relationship effects) | ECT(t-Stats) | |||||
Coefficient(p-value) | Standard Error | Coefficient(p-value) | Standard error | Coefficient (p-value) | Standard error | ||
ΔEC | 1 | - | −0.19 | −0.39 | 0.15 | −0.21 | (−2.66)* |
(0.03)* | (0.01)* | ||||||
ΔEP | −0.13 | −0.12 | 1 | - | 0.05 | −0.12 | (−3.97)** |
(0.04)* | (0.12) | ||||||
ΔGDP | 0.09 | −0.14 | −0.16 | −0.24 | 1 | - | (−3.36)* |
(0.001)* | (0.00)** | ||||||
Notes: The number in the parentheses is the p-value. *, ** asterisks indicate significance at 5% and 1% level, respectively. |
Period | Variance Decomposition of EC | |||
S.E. | EC | EP | GDP | |
1 | 0.026154 | 100 | 0 | 0 |
2 | 0.040085 | 96.92866 | 2.657988 | 0.41335 |
3 | 0.055235 | 96.07497 | 3.657699 | 0.267331 |
4 | 0.066524 | 95.30433 | 4.487052 | 0.208615 |
5 | 0.076296 | 93.71978 | 6.059193 | 0.221029 |
6 | 0.084754 | 91.77491 | 7.982465 | 0.242629 |
7 | 0.092929 | 89.84477 | 9.908605 | 0.246629 |
8 | 0.100791 | 87.9735 | 11.77734 | 0.249163 |
9 | 0.108377 | 86.54688 | 13.1834 | 0.269722 |
10 | 0.11552 | 85.69038 | 13.99968 | 0.309941 |
Period | Variance Decomposition of EP: | |||
S.E. | EC | EP | GDP | |
1 | 0.048996 | 56.37411 | 43.62589 | 0 |
2 | 0.080503 | 47.46051 | 52.45839 | 0.081095 |
3 | 0.111427 | 54.85629 | 44.96984 | 0.173865 |
4 | 0.136671 | 58.00428 | 41.5015 | 0.49422 |
5 | 0.158198 | 61.49309 | 36.99908 | 1.507833 |
6 | 0.17542 | 65.03976 | 32.20885 | 2.751394 |
7 | 0.190837 | 68.51778 | 27.79294 | 3.689283 |
8 | 0.20503 | 71.39459 | 24.14619 | 4.459221 |
9 | 0.218298 | 73.62662 | 21.30353 | 5.069856 |
10 | 0.23014 | 75.36009 | 19.18527 | 5.454636 |
Period | Variance Decomposition of GDP: | |||
S.E. | EC | EP | GDP | |
1 | 0.027913 | 7.08547 | 1.074455 | 91.84007 |
2 | 0.047862 | 2.562744 | 1.192117 | 96.24514 |
3 | 0.058702 | 1.780112 | 2.351816 | 95.86807 |
4 | 0.067107 | 2.679764 | 3.801723 | 93.51851 |
5 | 0.076742 | 4.208902 | 5.118915 | 90.67218 |
6 | 0.086717 | 4.775009 | 6.861578 | 88.36341 |
7 | 0.095723 | 5.065867 | 8.767336 | 86.1668 |
8 | 0.103643 | 5.248982 | 10.21453 | 84.53649 |
9 | 0.110786 | 5.330424 | 11.07375 | 83.59582 |
10 | 0.117053 | 5.316104 | 11.42037 | 83.26352 |
Author | Country (Period) | Variables | Methodology | Direction of Causality | |||
→ | ← | ↔ | ≠ | ||||
[45] | China (1971–2001) | EC, GDP | Granger Causality |
✔ | ✘ | ✘ | ✘ |
[37] | Pakistan (1972–2009) | EC, RGDP | ARDL, VECM | ✘ | ✔ | ✘ | ✘ |
[46] | Switzerland (1950–2010) | EC, RGDP | ARDL, UECM | ✘ | ✘ | ✔ | ✘ |
[46] | United Kingdom(1975–2010) | EC, RGDP | Panel bootstrap causality | ✔ | ✘ | ✔ | ✘ |
[46] | OECD(1990–2008) | EC, GDP | Panel Granger causality, | ✔ | ✘ | ✔ | ✘ |
[46] | Lebanon (1983–2014) | EC, GDP | Toda-Yamamoto | ✘ | ✘ | ✔ | ✘ |
[46] | Algeria (1971–2010) | EC, GDP | Cointegration, VECM | ✘ | ✘ | ✔ | ✘ |
[3] | Pakistan (1972–2010) | EC, RGDP | Granger Causality | ✔ | ✘ | ✘ | ✘ |
[47] | Poland (2000–2012) | EC, GDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[48] | Turkey (1970–2011) | EC, GDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[19] | Malaysia (2005–2010) | EC, RGDP | Multiplier Approach | ✔ | ✘ | ✘ | ✘ |
[24] | Algeria(1980–2012) | RE, NRE | Granger causality | ✔ | ✘ | ✔ | ✘ |
[49] | Taiwan (1998–2014) | EC, RGDP | Granger causality | ✘ | ✘ | ✔ | ✘ |
[50] | U.S.(2005–2015) | EC, RGDP | Wavelet | ✔ | ✘ | ✘ | ✘ |
[51] | Ghana (1971–2014) | EC, Industrial Growth | ARDL, ECM | ✘ | ✔ | ✘ | ✘ |
[52] | Pakistan(1981–2017) | RE, NRE | ARDL, VECM | ✘ | ✘ | ✔ | ✘ |
[53] | Pakistan (1970–2018) | EC, EP, GDP | Granger causality VECM | ✔ | ✘ | ✘ | ✘ |
Note: →, ←, ↔, and ≠ denotes unidirectional bidirectional and no causality, respectively. |
Variables | Data Source | Scale Unit |
Industrial electricity consumption | [54] | Gigawatt hours (GWh) |
Electricity prices | [55] | Millions |
Real gross domestic product | [56] | $Millions |
Industrial Variables | Philips-Perron (PP) | Augmented Dickey-Fuller (ADF) | Order of Integration | ||
Levels | First difference | Levels | First difference | ||
EC | (0.32) | (0.00)* | (0.35) | (0.00)* | I(1) |
EP | (0.49) | (0.00)* | (0.51) | (0.00)* | I(1) |
GDP | (0.31) | (0.00)* | (0.29) | (0.00)* | I(1) |
Notes: *, ** asterisk specifies significance at 1% and 5% level, respectively. p-value denotes in () parentheses. |
Hypothesized | r = 0 | r ≤ 1 | r ≤ 2 |
Trace Statistics | 35.40 | 12.83 | 3.29 |
P-Value | (0.01)* | (0.003)* | (0.12) |
Max-Eigen Statistics | 22.57 | 9.54 | 3.33 |
P-Value | (0.03)* | (0.002)* | (0.07) |
Notes: The numbers in the parentheses are p-values, r: indicates the number of co-integration hypothesis relationships, *, ** asterisk implies significance at 5% and 10% level, correspondingly. |
Industrial Sector | VECM (Long-run relationship effects) | ECT(t-Stats) | |||||
Coefficient(p-value) | Standard Error | Coefficient(p-value) | Standard error | Coefficient (p-value) | Standard error | ||
ΔEC | 1 | - | −0.19 | −0.39 | 0.15 | −0.21 | (−2.66)* |
(0.03)* | (0.01)* | ||||||
ΔEP | −0.13 | −0.12 | 1 | - | 0.05 | −0.12 | (−3.97)** |
(0.04)* | (0.12) | ||||||
ΔGDP | 0.09 | −0.14 | −0.16 | −0.24 | 1 | - | (−3.36)* |
(0.001)* | (0.00)** | ||||||
Notes: The number in the parentheses is the p-value. *, ** asterisks indicate significance at 5% and 1% level, respectively. |
Period | Variance Decomposition of EC | |||
S.E. | EC | EP | GDP | |
1 | 0.026154 | 100 | 0 | 0 |
2 | 0.040085 | 96.92866 | 2.657988 | 0.41335 |
3 | 0.055235 | 96.07497 | 3.657699 | 0.267331 |
4 | 0.066524 | 95.30433 | 4.487052 | 0.208615 |
5 | 0.076296 | 93.71978 | 6.059193 | 0.221029 |
6 | 0.084754 | 91.77491 | 7.982465 | 0.242629 |
7 | 0.092929 | 89.84477 | 9.908605 | 0.246629 |
8 | 0.100791 | 87.9735 | 11.77734 | 0.249163 |
9 | 0.108377 | 86.54688 | 13.1834 | 0.269722 |
10 | 0.11552 | 85.69038 | 13.99968 | 0.309941 |
Period | Variance Decomposition of EP: | |||
S.E. | EC | EP | GDP | |
1 | 0.048996 | 56.37411 | 43.62589 | 0 |
2 | 0.080503 | 47.46051 | 52.45839 | 0.081095 |
3 | 0.111427 | 54.85629 | 44.96984 | 0.173865 |
4 | 0.136671 | 58.00428 | 41.5015 | 0.49422 |
5 | 0.158198 | 61.49309 | 36.99908 | 1.507833 |
6 | 0.17542 | 65.03976 | 32.20885 | 2.751394 |
7 | 0.190837 | 68.51778 | 27.79294 | 3.689283 |
8 | 0.20503 | 71.39459 | 24.14619 | 4.459221 |
9 | 0.218298 | 73.62662 | 21.30353 | 5.069856 |
10 | 0.23014 | 75.36009 | 19.18527 | 5.454636 |
Period | Variance Decomposition of GDP: | |||
S.E. | EC | EP | GDP | |
1 | 0.027913 | 7.08547 | 1.074455 | 91.84007 |
2 | 0.047862 | 2.562744 | 1.192117 | 96.24514 |
3 | 0.058702 | 1.780112 | 2.351816 | 95.86807 |
4 | 0.067107 | 2.679764 | 3.801723 | 93.51851 |
5 | 0.076742 | 4.208902 | 5.118915 | 90.67218 |
6 | 0.086717 | 4.775009 | 6.861578 | 88.36341 |
7 | 0.095723 | 5.065867 | 8.767336 | 86.1668 |
8 | 0.103643 | 5.248982 | 10.21453 | 84.53649 |
9 | 0.110786 | 5.330424 | 11.07375 | 83.59582 |
10 | 0.117053 | 5.316104 | 11.42037 | 83.26352 |