Next Article in Journal
Adjacent Markets Influence Over Electricity Trading—Iberian Benchmark Study
Next Article in Special Issue
Do Oil Price Shocks and Other Factors Create Bigger Impacts on Islamic Banks than Conventional Banks?
Previous Article in Journal
Effect of Fins on the Internal Flow Characteristics in the Draft Tube of a Francis Turbine Model
Previous Article in Special Issue
An Efficient Analytical Approach for the Solution of Certain Fractional-Order Dynamical Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of the Risk Quantification in Traditional and Renewable Energy Markets

by
Daniel Velásquez-Gaviria
1,
Andrés Mora-Valencia
2 and
Javier Perote
3,*
1
Departamento de Finanzas, Instituto Tecnológico Metropolitano-ITM, Medellín 050001, Colombia
2
School of Management, Universidad de los Andes, Bogotá 111711, Colombia
3
Department of Economics and IME, University of Salamanca (IME), 37007 Salamanca, Spain
*
Author to whom correspondence should be addressed.
Energies 2020, 13(11), 2805; https://doi.org/10.3390/en13112805
Submission received: 2 April 2020 / Revised: 12 May 2020 / Accepted: 29 May 2020 / Published: 1 June 2020
(This article belongs to the Special Issue Mathematical and Statistical Models for Energy with Applications)

Abstract

:
The transition from traditional energy to cleaner energy sources has raised concerns from companies and investors regarding, among other things, the impact on financial downside risk. This article implements backtesting techniques to estimate and validate the value-at-risk (VaR) and expected shortfall (ES) in order to compare their performance among four renewable energy stocks and four traditional energy stocks from the WilderHill New Energy Global Innovation and the Bloomberg World Energy for the period 2005-2016. The models used to estimate VaR and ES are AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1), and AR(1)-APARCH(1,1), all of them under either normal, skew-normal, Student’s t, skewed-t, Generalized Error or Skew-Generalized Error distributed innovations. Backtesting performance is tested through traditional Kupiec and Christoffersen tests for VaR, but also through recent backtesting ES techniques. The paper extends these tests to the skewed-t, skew-normal and Skew-Generalized Error distributions and applies it for the first time in traditional and renewable energy markets showing that the skewed-t and the Generalized Error distribution are an accurate tool for risk management in those markets. Our findings have important implications for portfolio managers and regulators in terms of capital allocation in renewable and traditional energy stocks, mainly to reduce the impact of possible extreme loss events.

1. Introduction

The demise of fossil fuels is imminent, as can be seen in the recent generation of electricity in the United States. In April 2019, for the first time, renewable sources generated more electricity than coal, according to Bloomberg News based on data from the Energy Information Administration [1]. Moreover, Saudi Arabia is keen to move towards clean energy development, recognizing the changing global agenda in energy resources [2]. This climate shift is leading to concerns among regulators and investors regarding financial downside risk. Consequently, risk quantification and management of financial assets related to traditional and renewable energy have become increasingly important in both academia and industry. This interest is triggered by different issues, such as the high energy price volatility or the demand for mitigation of climate change, which encourages the intensification of energy generation from renewable sources, among others.
In energy markets, it is known that renewable energy stocks are more volatile than traditional energy stocks; the annual volatility between 2005 and 2016 associated with the main traditional energy index, Bloomberg World Energy (BWEI), was 17%, while that of the main renewable energy index, WilderHill New Energy Global Innovation (NEX), was 35%. This difference in volatility can be explained by several factors that more aggressively affect renewable energy stocks; among them, weather conditions, inventory levels, geopolitical events, legislation and exchange rate movements [3]. As a result, investors participating in this market face higher risks in comparison to the traditional energy market.
Regarding investment opportunities, the global green economy represents market capitalization comparable to the fossil fuel sector as of December 2017, according to an FTSE Russell research [4]. Although the market share for each sector represents around 6% (USD 4 trillion), the green economy has been growing since 2013, and it could represent 7% by 2030 based on the current trend. The potential growth could be higher, considering that renewables (excluding hydropower) produce only 8% of the world’s electricity, and financial innovation is helping in developing new ways of investing in renewable energy, as stated in a report of the Economist Intelligence Unit [5]. On the other hand, there has been a decreasing tendency in fossil fuel share since 2011. However, there is a greater role to be played by governments in mitigating the effects of climate change. For instance [6], a Senator in Nebraska argued that renewable energy is very expensive, and replacing traditional energy sources is not foreseeable in the short term.
There are some other risk factors faced by investment in traditional or renewable energy. The performance of traditional energy investment depends greatly on the behavior of oil prices, which is conditional on the supply and demand of the commodity. Oil amounted to 46% of global energy supply during the Arab oil embargo in 1973. In 2014 the stranded assets (fossil fuels) still dominated the energy market, 31% for oil, 29% for coal, and 21% for natural gas, whereas the main renewable sources (wind, solar and geothermal energy) only accounted for around 1% of the market supply.
Additionally, oil companies are under pressure by European and American shareholders who ignore how these companies are going to deal with their business if regulation becomes more oriented to a “cleaner” world energy market. This raises concerns for financial regulators, especially for the current Governor of Bank of England (Mark Carney), about the financial instability potentially caused by the energy transition, since around 80% of fossil fuel reserves could be stranded [7]. Energy markets are becoming more volatile, and shares that were profitable in the past could return huge losses to investors. This is the case of stock prices of E.ON and RWE (the two largest German utilities companies), which declined by two thirds after Japan’s Fukushima nuclear disaster in 2011 [8].
For all the previously stated reasons, there is a need to assess potential losses adequately in traditional and renewable energy financial assets. More specifically, it is necessary to gauge the market risk posed by investments in such assets. To this end, our study implements validation tests on the recently proposed risk measure by the Basel Committee, expected shortfall (ES), which captures the tail risk more adequately than its predecessor, the Value-at-Risk measure. Nonetheless, the accuracy of the risk measure depends on the distribution assumed to fit the returns of the assets.
In this context, the purpose of this study is to provide a technical approach to quantifying market risk measures and compare its performance between traditional and renewable energy markets, considering the financial regulation changes by the BCBS combined with new ES backtesting techniques proposed in the literature. The dataset is composed of four traditional energy stocks, Epistar, Kingspan, Boralex, Nibe, and four traditional energy stocks, Exxon Mobil, Chevron, Total S.A., Schlumberger, the period is from January 2005 to December 2016, for a total of 2901 prices and 2900 log-returns.
Furthermore, we extend the analysis to the skew-normal, skewed-t, Generalized Error and Skew-Generalized Error distributions (GED and SGED henceforth, respectively) under a comprehensive analysis of AR-GARCH/EGARCH/APARCH processes to filter the returns of the analyzed stocks over the period 2005–2016. Our results indicate that the skewed-t and the SGED are the probability distributions that perform best in backtesting when the AR(1)-GARCH(1,1) model is used. In general, there are more VaR violations with renewable energy stocks than with traditional energy stocks; however, renewable energy stocks perform better in ES backtesting than traditional energy stocks. All this reveals the skewed-t and the SGED to be very accurate tools for quantifying ES and implementing risk management policies in energy markets. In terms of ES, the skew-t has better results in the traditional energy sector, and the SGED performs better in the renewable energy sector.
The remainder of this paper is structured as follows. Section 2 provides a literature review in risk quantification. Section 3 reviews the VaR and ES methodologies, as well as the AR-GARCH models. Section 4 presents the data. Section 5 shows an application of VaR and ES backtesting to traditional and renewable energy stocks and the results of the application, whereas Section 6 presents a discussion in light of the new regulation. Lastly, Section 7 summarizes the conclusions.

2. Literature Review: Risk Quantification

This study applies recent ES tests proposed by Acerbi and Székely [9,10] and extends these tests to the skewed Student’s t (skewed-t hereafter) and SGED distributions, which are flexible distributions with good performance for capturing risk in financial markets. We introduce this approach to measuring risk in the renewable energy market and present a comparison in terms of risk assessment with the traditional energy market. The paper aims to provide adequate tools for investors, portfolio and risk managers, and regulators to mitigate the negative impact of several risk factors on investments in traditional and renewable energy stocks.
Risk in energy markets has traditionally been quantified through Value-at-Risk (VaR) [11,12], which can be defined as the maximum expected loss at a given confidence level and time horizon. Statistically, VaR corresponds to a quantile of the asset return distribution. Despite its common use, this measure does not meet the ‘subadditivity’ property or diversification criteria, which means that it is not consistent with portfolio diversification, i.e., the risk of a diversified portfolio may be higher than that incurred independently by its components [13]. Additionally, it is argued that VaR is unable to adequately capture fat-tailed risk and thus performed poorly during the recent financial crisis. Consequently, in 2013, the Basel Committee on Banking Supervision (BCBS) decided to replace VaR at the 1% significance level with ES also known as Tail VaR (TVaR) or Conditional VaR (CVaR), at a 2.5% significance level as the appropriate measure to estimate the regulatory capital under Basel III. ES, defined as the expected loss given that losses have exceeded VaR, is a coherent risk measure [14]. This implies that it not only has better mathematical properties than VaR, but that it also captures the extreme loss events that impact the price of traditional and renewable energy stocks in periods of high instability. However, despite these improvements, ES has not been fully accepted due to criticisms regarding its implementation; these are mainly linked to the lack of ‘elicitability’ [15] and its implications for backtesting. The latter refers to the process of validating the estimates of the risk measure by verifying whether the realized losses, observed ex post, are in line with the ex-ante estimates or forecasts. Research, see, e.g., Acerbi and Szekely [9], Constanzino and Curran [16], Du and Escanciano [17], Fissler and Ziegel [18] and Fissler et al. [19] has demonstrated that elicitability affects the risk measure selection, but not its validation; therefore, backtesting techniques for ES can be performed.
Recent works have assessed expected shortfall in financial markets [20,21], but also in the energy field. For instance, Ref. [22] applied VaR and ES on energy markets, and the author recommended tougher guidelines for the proposal of Basel III; Ref. [23] also employed these risk measures for hedging comparison using future contracts. On the other hand, for systemic risk, marginal and component ES have been estimated to find that Energy and Industrial sectors are the main risk contributors in the Chinese stock markets [24]. Chinese energy stocks have also been examined for diversification purposes in [25]. The authors found that commodity futures could help to reduce the expected shortfall of the diversified portfolios. Expected shortfall or expected shortage was calculated in food, energy, water and environmental security worldwide data by [26]. Other authors [27] have proposed a risk-averse two-stage stochastic mixed-integer optimization model with Conditional Value-at-Risk or ES as a risk measure for power generation capacity expansion, while Ref. [28] applied the financial risk measures but in an energy security context for better understanding of energy accident risks. In related applications, Ref. [29] employed Bayesian CVaR super-quantile regression on the energy price dataset, whereas Ref. [30] used ES and VaR as risk measures in a real options context.

3. Materials and Methods

This section provides the definitions of VaR and ES with their respective closed expressions according to three probability density functions (pdfs): normal, skew-normal, Student’s t, skewed-t, GED and SGED. Subsequently, the complete VaR/ES measure for different mean-variance models is defined. Finally, the three tests proposed by Acerbi and Székely [9,10] for ES backtesting are presented.

3.1. Dynamics of VaR/ES with Mean and Variance Models

It is known that daily financial returns exhibit characteristics of leptokurtosis (fat-tailed risk), skewness, clustering, and persistence in the conditional volatility, in addition to long memory in stochastic shocks [31,32,33,34,35,36], among others. These returns have a predictable component in their conditional mean, which has traditionally been modeled with ARMA structures; in this case, the AR(1) model is used in accordance with most financial literature. The other standard features have traditionally been captured by different conditional volatility models, among them we chose the GARCH [37], EGARCH [38] and APARCH [39] models; the latter two can capture asymmetric effects on conditional volatility (the so-called ‘leverage effect’, which involves a higher severity of negative returns). From an empirical viewpoint, the combination of AR with these different GARCH models was selected according to their performance and accuracy criteria (e.g., AIC or BIC). The formal representation of the models for the conditional mean and the variance are:
A R ( 1 ) μ t = φ + ϕ μ t 1 + ε t ,  
where 1 < ϕ < 1 and φ > 0 .
G A R C H ( 1 , 1 ) σ t 2 = ω + α ε t 1 2 + β σ t 1 2 ,  
where α > 0 ,   β > 0 ,   α + β < 1 and ω > 0 .
E G A R C H ( 1 , 1 ) log ( σ t 2 ) = ω + γ ( | Z t 1 | E | Z t 1 | ) + α Z t 1 + β log ( σ t 1 2 ) ,  
where ω > 0 and γ   .
A P A R C H ( 1 , 1 ) σ t δ = ω + α ( | ε t 1 | γ ε t 1 ) δ + β σ t 1 δ ,  
where δ < 0 and 1 < γ < 1 .
Standardized errors are modeled according to a certain parametric distribution,
Z t = ε t σ t ,   ~   iid   F θ ( Z t ) ,
where Z t represents independent and identically distributed (iid) innovations. The normal ( θ = ( μ = 0 , σ = 1 ) ), the skew-normal ( θ = ( μ = 0 , σ = 1 ,   ξ )), Student’s t ( θ = ν ), skewed-t ( θ = ( ν , λ ) ), GED ( θ = η ) and SGED ( θ = ( η , ζ ) ) distributions are employed to characterize the innovations of GARCH models. Where ξ stands for the skewness parameter in the skew-normal, ν stands for the degrees of freedom in Student’s t, λ stands for skewness parameter in the skewed-t, η stands for the shape parameter in GED, and ζ stands for the skew parameter in SGED.
In the empirical application, we estimate the two risk measures with a time horizon of one day, E S ^ t + 1 ,   1 α and V a R ^ t + 1 ,   1 α , which are subject to the information available up to time t ( t ), which comprises a 500-day rolling window. The parameters of the pdfs associated with the innovations Z t and the AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) model parameters are jointly estimated through maximum likelihood (ML) method and used to estimate the dynamic VaR and ES as follows:
V a R ^ t + 1 , 1 α = μ ^ t + 1 + σ ^ t + 1 V a R 1 α ( Z t + 1 ) ,
E S ^ t + 1 ,   1 α = μ ^ t + 1 + σ ^ t + 1 E S 1 α ( Z t + 1 ) ,
where V a R 1 α ( Z t + 1 )   and E S 1 α ( Z t + 1 )   are computed according to Equations (8) and (9), respectively, depending on the distribution function.

3.2. Risk Measures under Different Distributions

Given the innovations, Z t   (   t = 1 T ) with a cumulative distribution function (cdf) F Z , the V a R 1 α ( Z t ) , and E S 1 α ( Z t ) can be defined as:
V a R 1 α ( Z ) = i n f { z   :   F Z ( z ) 1 α } ,
E S 1 α ( Z ) = 1 α 0 α V a R 1 q ( Z ) d q ,
where α represents the significance level, which conforms to the regulation is set at 1% and 2.5% for the VaR and ES, respectively. According to this notation, VaR and ES are quantified in the left tail of the innovations (1% and 2.5%), and it is denoted VaR at 99% and ES at 97.5%. Equation (8) can be expressed as:
V a R 1 α ( Z ) = F Z 1 ( α ) ,
where F Z 1 stands for the quantile function, which is represented by Φ 1 ,   Φ ξ 1 , t ν 1 , s t ν , λ 1 , g η 1 and g η , ζ 1 for the cases of the normal, skew-normal, Student’s t, skewed-t, GED and SGED, respectively. Based on Equation (9), and solving the integral, ES under a normal distribution is:
E S 1 α ( Z ) = ϕ ( Φ 1 ( α ) ) α ,  
where ϕ represents the pdf of the standard normal distribution. For the skew-normal distribution, ES is:
E S 1 α ( Z ) = e s 1 α ( ξ ) ,
where e s 1 α ( ξ ) is estimated through numerical integration, using Equation (9) and the pdf of the skew-normal distribution following Fernández and Steel [40] in the general method to transform a unimodal symmetric distribution into a skew-symmetric distribution:
s n ( z , ξ ) = { 2 λ + 1 λ ϕ ( ξ z ) f o r z < 0 , 2 λ + 1 λ ϕ ( z ξ ) f o r z 0 ,
where ξ is the skewness parameter ( ξ   ( 0 , ) ), and ϕ ( z ) denotes the standard normal pdf.
For the Student’s t distribution, the ES [41] is expressed as:
E S 1 α ( Z ) = g ν ( t ν 1 ( α ) ) α ( ν + ( t ν 1 ( α ) ) 2 ν 1 ) ,  
where ν represents the degrees of freedom and t ν   and g ν denote the Student’s t cdf and pdf, respectively.
For the skewed-t distribution, ES is:
E S 1 α ( Z ) = e s 1 α ( v , λ ) ,
where e s 1 α ( v , λ ) is estimated through numerical integration, using Equation (9) and the pdf of the skewed-t distribution proposed by Fernández and Steel [40]:
s t ( z , λ ) = { 2 λ + 1 λ t ( λ z ) f o r z < 0 , 2 λ + 1 λ t ( z λ ) f o r z 0 ,
where λ is the skewness parameter ( λ   ( 0 , ) ), and t ( z ) denotes the Student’s t pdf. For GED, the ES is:
E S 1 α ( Z ) = e s 1 α ( η ) ,
where e s 1 α ( η ) is estimated by numerical integration using Equations (9) and (18). The GED’s pdf as introduced in [38] is:
g ( z ) = η e [ ( 1 2 ) | z κ | η ] κ 2 ( 1 + 1 η ) Γ ( 1 η ) ,
where κ = [ 2 ( 2 / η ) Γ ( 1 / η ) / Γ ( 3 / η ) ] 1 / 2 , and η stands for the shape parameter ( η   ( 0 , ) ). For SGED, the ES is:
E S 1 α ( Z ) = e s 1 α ( η , ζ ) ,
where e s 1 α ( η , ζ ) is estimated throughout numerical integration using Equations (9) and (20), where Equation (20) stands for the SGED’s pdf:
s g ( η , ζ ) = { 2 ζ + 1 ζ g ( ζ z ) f o r z < 0 , 2 ζ + 1 ζ g ( z ζ ) f o r z 0 ,
where ζ is the skewness parameter ( ζ   ( 0 , ) ), and g ( z ) represents the GED’s pdf.

3.3. VaR Backtesting

The literature on VaR backtesting is broad; this article uses the tests by Kupiec [42] and the independence test of Christoffersen [43]. The first test is based on the observation of the VaR estimated with a significance level α , and identifying the occasions when realized losses exceed the VaR during the backtesting period, which are called VaR exceptions. To validate the risk measure, the sequence of exceptions should follow an iid Bernoulli process with probability α ( 1 α ) . This binomial test for the exceptions is known as an unconditional coverage test. The second test examines the independence of the exceptions. The period of both VaR and ES backtesting in the empirical application is T = 2400, and the rolling window size is 500 observations.

3.4. ES Backtesting

Three tests have been proposed for ES backtesting; two of them were presented in [9], and the other was introduced in [10]. To the best of our knowledge, this is the first application of these three tests on traditional and renewable energy markets. To implement the tests, we assume that the model performance derives from a ‘true but unknown’ conditional distribution function F t , which is unknown, but can be estimated by a ‘predictive’ probability distribution P t . In our case, P t represents the normal, skew-normal, Student’s t, and skewed-t distributions. Following the same notation of the authors, the risk measure is denoted by V a R 1 α F and E S 1 α F when the variable X is distributed according to F . It is assumed that the return distribution is continuous and strictly increasing. Therefore, the ES can be expressed as:
E S 1 α ( X t ) = E [ X t | X t + V a R 1 α , t < 0 ] .
For all three tests, 97.5%-VaR backtest is performed first, after which the exceptions are used to build the ES backtests. The null hypothesis assumes that the distribution employed to construct the risk measures is correct, and thus the alternative hypothesis implies risk misspecification.

3.4.1. Z1 Test

The first test proposed by Acerbi and Székely [9] is based on the conditional expectation of the exceptions with VaR at 97.5%. From Equation (17), it follows that
E [ X t E S 1 α , t + 1 |   X t + V a R 1 α , t < 0 ] = 0 ,
Then, the Z 1 statistic is defined as
Z 1 ( X ) = t = 1 T X t I t E S 1 α , t N T + 1 ,
where I t is an indicator function that takes the value 1 if X t + V a R 1 α , t < 0 , and is otherwise equal to 0. N t = t = 1 T I t   is the sum of the exceptions at V a R 1 α . T is the backtesting period, which is 2400 in our empirical application. The null and alternative hypotheses being
H 0 :   P t [ α ] = F t [ α ] ,     t ,  
H 1 :   E S 1 α F E S 1 α , t ,     t ;   and   greater   for   some   t , V a R 1 α F = V a R 1 α , t ,     t .
In test Z 1 , null hypothesis H 0 indicates that the ‘predictive’ probability distribution P t , is correct, since it complies with the characteristics of the F t distribution to characterize financial returns. On the contrary, the alternative hypothesis H 1 indicates an underestimation of the ES, including when the VaR estimation is correct.
The Z 1 statistic reveals that E S 1 α is an appropriate risk measure when the result of the estimate is close to zero, more specifically, E H 0 [ Z 1 |   N T > 0 ] = 0 . Conversely, when E H 1 [ Z 1 |   N T > 0 ]   is negative, it implies that E S 1 α is not a suitable risk measure. Algorithm 1, described in Section 3.4.4, is used to find the critical values of the three backtesting ES tests. These values are tabulated in Table 1, for different values of the parameters of the normal, skew-normal, student’s t, skewed-t, GED and SGED distributions, considering a backtesting period T = 2400 and a 5% confidence level. The critical value at 5% for the normal, skew-normal, GED and SGED distributions is −0.03; for the Student’s t and skewed-t distributions, it is −0.05. The critical values are sensitive to the size T , as shown in [44], and also to the degrees of freedom parameter and shape parameter, but not to the skewness parameter of the skew-normal, skewed-t distributions and SGED. As T increases, the critical value of the test increases until it stabilizes close to zero.

3.4.2. Z2 Test

ES can be expressed as the following unconditional expectation:
E S 1 α , t = E [ X t I t α ] .
The Z 2 statistic is defined as
Z 2 ( X ) = t = 1 T X t I t T α E S 1 α , t + 1 ,
where the null and alternative hypotheses are:
H 0 :   P t [ α ] = F t [ α ] ,   t ,  
H 1 :   E S α F E S α , t ,     t ;   and   greater   for   some   t , V a R α F V a R α , t ,     t .
In test Z 2 , null hypothesis H 0 indicates that the ‘predictive’ probability distribution P t , is correct. On the contrary, the alternative hypothesis, H 1 , indicates an underestimation of the VaR and the ES.
The Z 2 test indicates that E S 1 α is an appropriate risk measure when the result of the estimate is close to zero, i.e., E H 0 [ Z 2 |   N T > 0 ] = 0 . On the contrary, E S 1 α is not an appropriate risk measure when E H 1 [ Z 2 |   N T > 0 ] < 0 , since losses are underestimated. This test penalizes the V a R 1 α due to the number of exceptions; additionally, the alternative hypothesis rejects V a R 1 α and E S 1 α . The critical value at 5% for the normal, skew-normal, GED and SGED distribution is −0.21; for the Student’s t and skewed-t distributions, it is −0.22 (see Table 1).

3.4.3. ZES Test

Tests Z 1 and Z 2 depend not only on the E S 1 α results, but also on the values of V a R 1 α because the latter needs to be backtested previously. Figure 1 illustrates how the critical value of the   Z 2 (ordinate axis) at 5% increases as the value of V a R 1 α (abscissa axis) does. This is crucial, because it increases the possibility of committing a type II error P ( a c c e p t   H 0 |   r e j e c t   H 0 ) . If the type II error occurs, a E S 1 α risk measure would be accepted, which is not suitable, and would leave the investment exposed to market risk.
Given the above feature, Acerbi and Székely [10] proposed the Z E S test, which is not as sensitive to possible VaR values (see dotted dash line in Figure 1); therefore, this test is expected to be more powerful. From Equation (15) and following [14], E S 1 α is equivalent to:
E S 1 α , t = V a R 1 α , t 1 α E [ ( X t + V a R 1 α , t ) I t ] .
Re-organizing the terms and performing the operations, the Z E S test becomes:
Z E S ( X ) = t = 1 T α ( E S 1 α V a R 1 α ) + ( X t + V a R 1 α ) I t T α E S 1 α , t .
The null and alternative hypotheses are:
H 0 :   P t [ α ] = F t [ α ] ,   for   all   t ,
H 1 :   E S 1 α F E S 1 α , t ,     t   and   lesser   for   some   t , V a R 1 α , t ~ V a R 1 α , t F
The Z E S test, just like the other two tests, indicates that E S 1 α is a good risk measure when E H 0 [ Z E S ] = 0 ; and when E H 1 [ Z E S ] < 0 , it indicates that E S 1 α is not a good risk measure. In this test, the alternative hypothesis indicates that E S 1 α is rejected, regardless of the value of V a R 1 α . The critical value at 5% for the normal and skew-normal distributions is −0.20, for Student’s t and skewed-t distributions, it is −0.27, while for GED and SGED, it is −0.22 (see Table 1). The next subsection explains the algorithm used to delimit the rejection area of Z 1 , Z 2 , and Z E S tests, which is based on the computation of the critical values at 5% under the normal, skew-normal, Student’s t, skewed-t, GED and SGED distributions.

3.4.4. Algorithm 1: Critical Values for ES Backtesting under Student’s t Distributions

The critical values of the Z 1 ,   Z 2 , and Z E S tests are determined by generating random numbers from the normal, skew-normal, Student’s t, skewed-t, GED and SGED distributions. In Panel A of Table 1, the degrees of freedom and coefficient of skewness parameters are ( ξ , λ ) = { 0.6 ,   0.8 ,   0.9 ,   1 ,   1.1 ,   1.2 } , and v = { 3 ,   5 ,   7 ,   10 } , respectively. In Panel B, the skewness parameters are ( ζ ) = { 0.6 ,   0.8 ,   0.9 ,   1 ,   1.1 ,   1.2 } , and the shape parameters are η = { 1 ,   1.2 ,   1.4 ,   1.6 , 1.8 } . We choose the parameters according to the model results in Table A2, Table A3 and Table A4 for traditional and renewable energy markets.
Algorithm 1
1. Simulation of X t i under normal, skew-normal, Student’s t, skewed -t, GED and SGED distributions, and M = 10,000 times. The i value indicates the test ( i = 1 , 2 , E S ).
2. Calculation of the Z i = Z ( X i ) statistics, using the data simulated in step 1.
3. Estimation of critical values (cv) at 5%, Z c v = F 5 % 1 ( Z i ) .
4. Resampling the critical values for the different skewness, degrees of freedom and shape parameters.
Table 1 shows the results of critical values.
If Algorithm 1 is applied with T = 2400, it is found that the mean of the Z 1 test with a skewed-t distribution is −0.000318, and the confidence interval at 95% is [−0.05868; 0.05199]. With the Z 2 test, the mean is 0.001611, and the confidence interval at 95% is [−0.25791; 0.24852]. With the Z E S test, the mean is −0.203813, and the confidence interval at 95% is [−0.28459; −0.13449]. For the SGED, the mean of the Z 1 test is −0.000082, and the confidence interval at 95% is [−0.04253; 0.03816]. For the Z 2 test, the mean is −0.001122, and the confidence interval at 95% is [−0.27060; 0.24579]. With the Z E S test, the mean is −0.168741, and the confidence interval at 95% is [−0.23137; −0.11399]. The critical values do not vary consistently over skew-t and SGED distributions. Figure 2 shows these results and the empirical distributions of the three tests under skew-t, and also identifies in blue the critical value at 5%.

4. Data

The dataset consists of the prices of four renewable energy stocks from WilderHill New Energy Global Innovation (NEX), and four traditional energy stocks from the Bloomberg World Energy index (BWEI). Appendix A gathers detailed information on the stocks. The prices are collected on a daily frequency basis from January 2005 to December 2016, for a total of 2900 observations. The logarithmic returns are calculated as r t = 100 log ( V t / V t 1 ) , where V t represents the price of each share at a time t . The units of the variables are logarithmic returns in decimals.
Table 2 displays the descriptive statistics of the data. The typical characteristics of the daily returns are featured, such as mean reversion, positive excess kurtosis, and non-zero skewness. The average annual volatility of renewable energy stocks is 41%, while that of the traditional energy stocks is 19%. Renewable energy stocks have an average excess kurtosis of 10, while that of the traditional energy stocks is 13. The high kurtosis values imply that the empirical distributions of returns exhibit extreme values. The renewable energy stocks Epistar, Kingspan, and Nibe and the traditional energy stock Schlumberger have negative coefficients of skewness; this means that in these stocks, negative returns persist more than positive returns. This finding represents evidence supporting the fact that the analyzed renewable energy stocks are riskier than traditional energy stocks. These stylized facts are illustrated in Figure A1 and Figure A2 (see Appendix B), which depicts prices and log-returns of different energy stocks.

5. Empirical Results

5.1. Conditional Mean and Variance Models in Sample Estimations

The AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) models are used to standardize the returns with the conditional mean and conditional variance. The AR and GARCH orders were selected according to the significance of the coefficients, validation of the assumptions of each model, and information criteria. Each model is estimated jointly by ML method. Table A2, Table A3 and Table A4 (see Appendix B) gather the coefficients associated with each model for the six distributions: normal, skew-normal, Student’s t, skewed-t, GED and SGED. These tables also display the estimates of the conditional mean-variance models, either AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1) or AR(1)-APARCH(1,1). In all cases, the coefficients associated with the AR(1) ( ϕ ) and GARCH(1,1) ( α and β ) processes are significant, both in traditional and renewable energy stocks. Similarly, the results of the coefficients of the skewness ( ξ ) of the skew-normal distribution, the degrees of freedom for the Student’s t distributions ( v ) , the skewness parameter of the skewed-t distributions ( λ ) , and the shape parameter ( η ) from GED, and the skewness parameter from SGED ( ζ ) ,   are consistently significant. For the EGARCH model, the parameter that represents the asymmetry of stochastic shocks in the conditional volatility ( γ ) is positive and significant in all cases. In the APARCH model, the parameter that also represents the asymmetric effects of the stochastic shocks on conditional volatility is significant in all cases, and it is even higher than 1 for the renewable energy stocks Epistar and Boralex, and the traditional energy stocks Exxon, Chevron and Schlumberger.
Table A5 (see Appendix B) shows the results of the two-sample Kolmogorov-Smirnov test, where the null hypothesis indicates that standardized returns and standardized residuals for the models AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1), AR(1)-APARCH(1,1) with the six distributions were drawn from the same continuous distribution. In addition to the estimate of the skewness parameters and excess kurtosis in the standardized residuals of the AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1), AR(1)-APARCH(1,1) models, each fitted with normal, skew-normal, Student’s t, skewed-t, GED and SGED distributions. In the two-sample Kolmogorov-Smirnov test, the null hypothesis, is rejected for all stocks and models, except for Epistar, from the renewable energy market. Nonetheless, the best results were achieved with the skew-t and GED distributions in the AR(1)-EGARCH(1,1) model. In all models, the kurtosis coefficients are greater than three. Moreover, in the AR(1)-GARCH(1,1) model, all the skewness coefficients are negative, except that of the traditional energy stock Nibe, which is positive. In the AR(1)-EGARCH(1,1) model under the Student’s t, skewed-t, GED and SGED distributions, three renewable energy stocks present positive skewness coefficients: Epistar, Kingspan and Nibe. For the remaining stocks, the skewness coefficients are negative. In the AR(1)-APARCH(1,1) model under the Student’s t and skewed-t distributions, two renewable energy stocks have positive skewness coefficients: Kingspan and Nibe. For the remaining stocks, the skewness coefficients are negative. The kurtosis coefficient is consistently lower for the GED and SGED distributions in the AR(1)-APARCH(1,1) model. In summary, in all three models with the three innovation distributions, the standardized residuals still exhibit skewness and kurtosis after filtering the returns. Figure 3 shows the standardized returns and standardized residuals of the AR(1)-GARCH(1,1) model, under the skew-normal, skew-t and SGED distributions.
Table A6 (see Appendix B) shows the Ljung-Box test results to detect autocorrelation of the standardized residuals from the AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1), AR(1)-APARCH(1,1) models, with the different innovation distributions. Considering no seasonal effects in the daily returns included in the sample, it is suggested the number of lags m = ln ( T ) [45] provides better power performance. Accordingly, we choose m = ln ( 2900 ) 8 . In addition, we added the twelfth and twentieth lags. There is no autocorrelation in any of the lags. This indicates that the proposed AR models are suitable for modeling renewable and traditional energy stock returns.
Table A7 (see Appendix B) shows the Ljung-Box test results to detect autocorrelation of the squared standardized residuals from the AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1), AR(1)-APARCH(1,1) models with different innovation distributions. The Kingspan stock shows significant autocorrelation in the first five lags. The Boralex stock, from the renewable sector, evidences autocorrelation until the fifth lag; nonetheless, there is no evidence of autocorrelation for the AR(1)-GARCH(1,1) under the student’s t and skew-t distribution. The Chevron stock displays significant autocorrelation for lags in the orders two, three and four in the AR(1)-GARCH(1,1) model under all distributions; nevertheless, there is no statistical evidence of autocorrelation under all distributions for the AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) models. The remaining stocks work well for all the models and distributions.

5.2. Backtesting of 99%-VaR

The results of the backtesting for 99%-VaR are presented in Table 3, including, for all the models and stocks, the observed and expected number of violations as well as the Kupiec and independence tests. Model performance is considered accurate if it simultaneously rejects the null hypothesis in both tests.
In the AR(1)-GARCH(1,1) model, the normal distribution does not perform accurately in any case; however, the renewable energy stock Nibe exceeds the Kupiec test. The model under the skew-normal only works well for Chevron. As for the Student’s t distribution, model performance is found to be accurate for four stocks: Epistar and Boralex for renewable energy and Total S.A. and Schlumberger for traditional energy. The model under skewed-t distribution performs well in all cases except with two traditional energy stocks—Chevron and Total S.A.—which do not reject the null hypothesis in the independence test. The null hypothesis is not rejected for any of the stocks in either test, except for Exxon Mobil and Chevron, from the traditional energy sector for the GED distribution. The model for the SGED distribution works well in all cases, except for Exxon Mobil.
For the AR(1)-EGARCH(1,1) model, the normal and skew-normal distributions are not appropriate in any case; however, as in the previous model for the normal, the renewable energy stock Nibe does not reject the null hypothesis in the Kupiec test. Regarding the Student’s t distribution, the model only works well in two cases, with the renewable energy stock Boralex and the traditional energy stock Schlumberger. The model under the skewed-t distribution performs accurately in four cases: renewable energy stock Boralex, and traditional energy stocks Exxon Mobil, Total S.A., and Schlumberger. For the GED distribution case, the null hypothesis is not rejected for all of the stocks in either of the tests, except for Exxon Mobil, Chevron and Total S.A. from the traditional energy industry. On the other hand, the null hypothesis is rejected for Epistar, Exxon Mobil and Chevron for the SGED distribution.
As far as the AR(1)-APARCH(1,1) model is concerned, under normal and skew-normal distributions, only the renewable energy stock Nibe performs well. Besides, the Student’s t distribution behaves adequately in three cases: Kingspan and Nibe for renewable energy and Schlumberger for traditional energy. The skewed-t model, however, presents adequate performance in six cases, two from renewable energy stocks (Kingspan and Boralex) and for all the four traditional energy stocks. For the GED and SGED distributions, both tests perform well in all cases, except for Exxon Mobil, Chevron and Total S.A. from the traditional energy sector.

5.3. Backtesting of 97.5%-VaR

Backtesting analyses for VaR are extended in Table 4 by considering 97.5%-VaR.
The AR(1)-GARCH(1,1) model under normal distribution performs well in three stocks: Kingspan and Nibe from the renewable energy sector, and Schlumberger from the traditional energy sector. The model under the skew-normal distribution performs well in Nibe, Exxon Mobil, and Schlumberger. All stocks pass the test with the Student’s t distribution, except two from the traditional power sector: Exxon Mobil and Chevron. The model performs accurately with the skewed-t distribution in five cases, but the renewable energy stocks Epistar and Nibe and the traditional energy stock Chevron reject the null hypothesis in the Kupiec and independence test. The model under the GED distribution does not reject both tests in any of the cases, except for Exxon Mobil and Chevron, from the traditional energy sector. Whereas for the SGED distribution, the model performs well in all cases except for Epistar, from the renewable energy sector.
The AR(1)-EGARCH(1,1) model under normal distribution is only satisfactory for the traditional energy stock Schlumberger. Under the skew-normal distribution, four stocks pass the tests: Nibe from the renewable sector, and Exxon Mobil, Chevron, and Schlumberger from the traditional energy sector. As for the Student’s t distribution, the model seems adequate in two cases, the renewable energy stock Nibe and the traditional energy stock Schlumberger. The renewable energy stocks do not pass the test for the skewed-t distribution. On the contrary, all the traditional energy stocks except for Chevron reject the null hypothesis in the Kupiec and independence test. The model under the GED distribution accurately exhibits results in all cases, except for Exxon Mobil, Chevron and Total S.A. In addition, for the SGED distribution, the model works well in all cases except for Epistar, from the renewable energy sector.
The AR(1)-APARCH(1,1) model performs well under normal and skew-normal distributions in two cases: with the renewable energy stock Nibe, and with the traditional energy stock Schlumberger. The model behaves adequately under the Student’s t distribution with the renewable energy stocks Kingspan, Boralex and Nibe, and the traditional energy stock Schlumberger. Under the skewed-t distribution, the model performs accurately with the renewable energy stock Boralex and with all the traditional energy stocks, with the exception of Chevron. The model under the GED distribution works well in all cases, except for Exxon Mobil, Chevron, and Total S.A. from the traditional energy sector. For the SGED distribution, the model works precisely in all cases except for Epistar from the renewable energy sector.

5.4. Backtesting of 97.5%-ES

Table 5 displays the results for the 97.5%-ES backtesting for the three tests (Z1, Z2 y ZES).

5.4.1. Z1 Test

With the AR(1)-GARCH(1,1) model under normal and skew-normal distributions, all of the stocks reject the null hypothesis of the test. Under the Student’s t distribution, the test works well with the renewable energy stocks Epistar, Kingspan and Boralex, and the traditional energy stock Total S.A. The test presents satisfactory results under the skewed-t distribution; in all cases, the null hypothesis of the test is not rejected, except with the traditional energy stock Schlumberger. For the GED and SGED, the test works well for Epistar, Kingspan and Boralex from the renewable energy sector, and for Total S.A. from the traditional energy sector.
In the AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) model under normal and skew-normal distributions, all of the stocks rejects the null hypothesis of the test. The four renewable energy stocks pass the test with the Student’s t and skewed-t distributions, and the traditional energy stocks Total S.A. and Chevron also pass, but the latter only passes the test in the AR(1)-APARCH(1,1) model. Under the GED and SGED, the test works well for Epistar, Kingspan and Boralex, and for Total S.A.

5.4.2. Z2 Test

In the AR(1)-GARCH(1,1) model under normal and skew-normal distributions, all of the stocks reject the null hypothesis of the test. For the Student’s t distribution, the test performs well with all four renewable energy stocks and with the traditional energy stock Schlumberger. The results improve noticeably with the skewed-t distribution; none of the stocks reject the null hypothesis of the test except the renewable energy stocks Epistar and Nibe. Under the GED, the test performs well for all four stocks from the renewable sector, and for Total S.A. and Schlumberger from the traditional energy sector. In the case of SGED, the test presents satisfactory results for all the stocks, except for Epistar, from the renewable energy sector.
In the AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) models under normal and skew-normal distributions, all of the stocks rejects the null hypothesis of the test. All the renewable energy stocks, except for Epistar in AR(1)-EGARCH(1,1), pass the test with the Student’s t distribution; Schlumberger from the renewable energy stocks also passes the test. None of the stocks reject the null hypothesis of the test under the skewed-t distribution, except the renewable energy stocks Epistar and Nibe. Under the GED, the test performs well for all the four stocks from the renewable sector, and for Schlumberger from the traditional energy sector. For the SGED, the test performs well for all the stocks, except for Epistar and Kingspan, from the renewable energy sector.

5.4.3. ZES Test

In the AR(1)-GARCH(1,1) model under normal and skew-normal distributions, all of the stocks reject the null hypothesis of the test. For the Student’s t distribution, the test works well with three renewable energy stocks (Epistar, Kingspan and Boralex) and with the traditional energy stock Schlumberger. With the skewed-t distribution, the null hypothesis of the test cannot be rejected for all the stocks, with the exception of the renewable energy stocks Epistar and Nibe and the traditional energy stock Exxon Mobil. Under the GED and the SGED, the test exhibits satisfactory results for Epistar, Kingspan, Boralex, Total S.A., and Schlumberger.
None of the stocks pass the test with the AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) models under normal and skew-normal distributions. With the Student’s t distribution, however, none of the renewable energy stocks reject the null hypothesis of the test; nor does the traditional energy stock Schlumberger. Under the skewed-t distribution, all stocks pass the test except the renewable energy stocks Epistar and Nibe. Under the GED and SGED, the test works well in all cases, except for Exxon Mobil and Chevron from the traditional energy market.
In summary, the higher flexibility of the skewed-t and the SGED provides the best results for VaR and ES testing. The skewed-t has better results in the traditional energy sector; meanwhile, the SGED performs better in the renewable energy sector. Nevertheless, we recommend filtering the returns by employing AR(1)-GARCH(1,1), since the results are similar to other GARCH(1,1) family models tested in our application. Furthermore, Z E S statistic is proposed to validate ES quantities, since this test is more potent than Z1 and Z2 tests.

6. Discussion

The empirical analyses in the previous section shed light on the best techniques for risk management in energy markets, revealing the skewed-t distribution to be an accurate tool, but also the adequacy of ES as a new paradigm for regulatory capital. In this section, we present a discussion on the implications in terms of capital reserves of the shift in the risk measure from VaR99% to ES97.5% in renewable and traditional energy stocks considering skewed-t distribution for the innovations. Table 6 shows the daily VaR99%, ES97.5% and the difference between these two risk measures diff = ES97.5%VaR99% calculated for the renewable and traditional energy stock returns with the AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) models in the analyzed sample. This difference (diff) represents the percentage of extra reserve capital that must be established in adopting ES.
In Table 6, we can observe that VaR, ES, and diff values do not vary considerably between AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1), and AR(1)-APARCH(1,1) models. On average, VaR and ES quantities are higher in renewable (−4.6338 and −4.8391, respectively) than in traditional energy stocks (−3.2264 and −3.3087), this can be attributed to the higher volatility faced in the renewable energy market. In addition, ES estimates in all cases are considerably higher than VaR. This is a result of extreme values that fatten the tail of the return distribution, which leads ES to assign higher regulatory capital when the skewed-t distribution is employed to fit the standardized residuals of stock returns. These findings have important implications for portfolio managers and regulators in terms of capital allocation in renewable and traditional energy stocks mainly to reduce the impact of possible extreme losses events.
To have a better understanding of the effects of switching the risk measures, we compute AR(1)-GARCH(1,1)-VaR and -ES for two stocks, Epistar (renewable energy) and Exxon Mobil (traditional energy), by employing a 500-day rolling window, and then we calculate the daily difference between the two risk measures. The results are illustrated in Figure 4. In this graph, it can be observed that, on average, the difference between both risk measures is much higher in Epistar (0.2402%) than in Exxon Mobil (0.1119%). Nevertheless, in 2008, within the context of the subprime crisis, differences for Exxon Mobil sharply rose to 1%, while those for Epistar remained at 0.25%. On the contrary, in relative calm periods such as in 2012 to 2013, the difference between the risk measures is around 0.2%. On the other hand, the difference between the two risk measures has increased consistently for Epistar since 2012, with a difference between 0.3% and 0.7% in 2014 and 0.6% in 2016.
Figure 5 shows the difference between the 97.5%-ES and 99%-VaR, estimated with the AR(1)-GARCH(1,1) and the skewed-t distribution for innovations during the 2008 subprime crisis. Throughout this period, the ES was capable of adapting the extreme values and of allocating higher capital reserves than the VaR. This evidence emphasizes the need to search for the most accurate risk measures, especially in highly volatile scenarios. In these contexts, the use of ES with skewed-t distribution is advisable as a protection mechanism against extreme loss events in renewable energy stocks due to the capacity of skewed-t to adapt to potential volatility shocks and the asymmetric impact of market news.
Our results support the ES as a suitable risk measure for renewable and traditional energy stocks in the context of periods of high volatility and turmoil when investors face extreme market risk. Consequently, this study helps energy market participants to have a better understanding of the market risk posed in renewable and traditional energy stock investments. Moreover, the efficient allocation of regulatory capital through accurate risk management techniques also helps regulators to fight against market instability and policymakers to decide on real investment decisions that might have a significant impact on the future of traditional and renewable energies.

7. Conclusions

The increasing demand for renewable energy and the high volatility of the prices of these energy sources have increased the need to accurately measure the risk of these markets. This is particularly important for the decision-making of firms, investors, and regulators and represents a cornerstone for the stability and sustainability of these markets with profound implications for the whole economy. This paper fits into this framework, implementing a rigorous methodology for risk management. This section summarizes the main contributions and results of the paper, which are enumerated below.
(1) The article studies the risk in traditional and renewable energy markets by examining four renewable energy stocks and four traditional energy stocks. The prices of these stocks exhibit the same features as other financial assets, particularly leptokurtosis and skewness, calling for the use of conditional mean-variance models and flexible non-normal distributions. The performance of three models—(AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1))—with either normal, skew-normal, Student’s t, skewed-t, GED or SGED innovations are compared in terms of both VaR and ES backtesting techniques.
(2) The paper implements three tests proposed by Acerbi and Székely [9,10] to backtest 97.5%-ES. Backtesting of 99%-VaR and 97.5%-VaR is also carried out with Kupiec and independence tests. Our study contributes to the literature by extending the proposed ES tests to the skewed-t, GED and SGED distributions. A straightforward algorithm is proposed to compute critical values, and we tabulate them for different values of the skewed-t parameters. Furthermore, to the best of our knowledge, this is the first application of these tests to renewable and traditional energy stocks.
(3) The renewable energy stocks generally have higher price volatility than the traditional energy stocks, and thus the renewable energy stocks on average have more VaR violations than do the traditional energy stocks. The results indicate that the fitting of the normal distribution is weak in the three GARCH models for both traditional and renewable energy stocks. The best-fitting distribution in the backtests is the skewed-t, and AR(1)-GARCH(1,1) is the model that best fits the sample data.
(4) In the ES backtest, renewable energy stocks perform better than traditional energy stocks. We show that the tests proposed by Acerbi and Székely [9,10] are adequate for solving the drawback of the backtestability of ES, especially the Z E S   test, which is more powerful. For this test, the skewed-t has better results in the traditional energy sector; meanwhile, the SGED performs better in the renewable energy sector. Consequently, we provide evidence on the appropriateness of ES, the new risk measure proposed by the last Basel Accords, for quantifying the risk of traditional and renewable energy stocks.
(5) Our results have important implications for the efficient allocation of regulatory capital, particularly in periods of high instability. The use of flexible distributions as the skewed-t and the SGED may help to provide accurate capital provisions in case of potential losses provoked by extreme events but also to free resources in more calm periods that can be used for undertaking real investments in traditional or renewable energies.
All in all, our results constitute a methodological framework with empirical evidence that can be used by regulators and investors to mitigate the market risk emerging from renewable and traditional energy markets.

Author Contributions

The empirical analysis was mainly developed by D.V.-G. For the remaining issues related to the research, all authors contributed equally. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministry of Economics and Competitiveness, project ECO2016-75631-P, Castilla and León Government, project SA049G19, FAPA-Uniandes, project PR.3.2016.2807, and Instituto Tecnológico Metropolitano de Medellín.

Acknowledgments

The authors acknowledge the institutions above for funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Description of the stocks.
Table A1. Description of the stocks.
StockIndexDescription
EPISTARNEXEpistar Corporation manufactures and markets light-emitting diode (LED) chips and epitaxial wafers. The Company sells its products in Taiwan and exports worldwide.
Kingspan Group PLCNEXKingspan Group PLC is a global market player in high-performance insulation and building envelope technologies.
Boralex Inc.NEXBoralex Inc. is an electricity producer whose core business is the development and operation of renewable energy power stations. The Corporation operates assets in the following power generation types—wind, hydroelectric, thermal and solar-based in Canada, the Northeastern United States, and France.
NIBE IndustrierNEXNIBE Industrier AB is an international heating technology company. The Company is organized around three business areas, all united under a shared vision to create world-class solutions in sustainable energy. NIBE produces and sells heat pumps, boiler and water heaters, electrical heating elements as well as freestanding fireplaces.
Exxon Mobil CorporationBWEIExxon Mobil Corporation operates petroleum and petrochemical businesses on a worldwide basis. The Company operations include exploration and production of oil and gas, electric power generation, and coal and mineral operations. Exxon Mobil also manufactures and markets fuels, lubricants and chemicals.
Chevron CorporationBWEIChevron Corporation is an integrated energy company with operations in countries located around the world. The Company produces and transports crude oil and natural gas. Chevron also refines, markets, and distributes fuels, as well as is involved in chemical and mining operations, power generation and energy services.
TOTAL S.A.BWEITOTAL S.A. explores for, produces, refines, transports, and markets oil and natural gas. The Company also operates a chemical division that produces polypropylene, polyethylene, polystyrene, rubber, paint, ink, adhesives, and resins. TOTAL operates gasoline filling stations in Europe, the United States and Africa.
Schlumberger LimitedBWEISchlumberger Limited is an oil services company. The Company, through its subsidiaries, provides a wide range of services, including technology, project management, and information solutions to the international petroleum industry as well as advanced acquisition and data processing surveys.
NEX stands for WilderHill New Energy Global Innovation index, and BWEI stands for Bloomberg World Energy index. Source: Bloomberg LP.

Appendix B. Tables and Figures

Table A2. Estimation of AR(1)-GARCH(1,1) models.
Table A2. Estimation of AR(1)-GARCH(1,1) models.
EpistarKingspanBoralexNibeExxonChevronTotalSchlumberger
CoefficientNormal Distribution
μ0.0000
(0.001)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
ϕ0.0499
(0.019)
0.0283
(0.020)
−0.0524
(0.019)
0.0017
(0.020)
−0.0517
(0.019)
−0.0342
(0.019)
−0.0172
(0.019)
−0.0267
(0.019)
ω0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
α0.0414
(0.007)
0.0952
(0.017)
0.0319
(0.004)
0.0931
(0.014)
0.0812
(0.010)
0.0837
(0.010)
0.0758
(0.010)
0.0496
(0.006)
β0.9433
(0.010)
0.8894
(0.020)
0.9659
(0.004)
0.8652
(0.020)
0.8989
(0.012)
0.9012
(0.011)
0.9074
(0.012)
0.9428
(0.007)
Skew-Normal Distribution
μ0.0001
(0.0005)
0.0012
(0.0004)
0.0007
(0.0003)
0.0013
(0.0004)
0.0004
(0.0002)
0.0005
(0.0002)
0.0004
(0.0003)
0.0007
(0.0003)
ϕ0.0531
(0.0197)
0.0298
(0.0198)
−0.0593
(0.0192)
0.0001
(0.0207)
−0.0611
(0.0198)
−0.0409
(0.0195)
−0.0219
(0.0196)
−0.0235
(0.0198)
ω0.0000
(0.0000)
0.0000
(0.0000)
0.0000
(0.0000)
0.0000
(0.0000)
0.0000
(0.0000)
0.0000
(0.0000)
0.0000
(0.0000)
0.0000
(0.0000)
α0.0405
(0.0028)
0.0485
(0.0182)
0.0317
(0.0022)
0.0862
(0.0142)
0.0857
(0.0093)
0.0866
(0.0111)
0.0682
(0.0095)
0.0507
(0.0076)
β0.9435
(0.0050)
0.9456
(0.0217)
0.9663
(0.0020)
0.8740
(0.0207)
0.8945
(0.0123)
0.8997
(0.0125)
0.9184
(0.0119)
0.9412
0.0089
ξ1.0119
(0.0235)
1.0162
(0.0221)
1.0032
(0.0194)
1.0340
(0.0212)
0.9062
(0.0211)
0.8597
(0.0225)
0.9295
(0.0216)
0.9729
(0.0226)
Student’s t Distribution
μ0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
ϕ0.0274
(0.019)
0.0216
(0.018)
−0.0653
(0.018)
−0.0065
(0.0019)
−0.0524
(0.018)
−0.0352
(0.019)
−0.0196
(0.018)
−0.0285
(0.018)
ω0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
α0.0500
(0.008)
0.0446
(0.009)
0.1131
(0.028)
0.0883
(0.019)
0.0831
(0.013)
0.0805
(0.011)
0.0689
(0.011)
0.0478
(0.007)
β0.9421
(0.010)
0.9497
(0.010)
0.8494
(0.039)
0.8762
(0.026)
0.9026
(0.014)
0.9089
(0.011)
0.9151
(0.014)
0.9484
(0.007)
8.7578
(1.529)
5.2540
(0.498)
5.8807
(0.298)
5.2573
(0.488)
6.6061
(0.766)
10.000
(1.452)
7.6678
(1.035)
7.8742
(1.028)
Skewed-t Distribution
μ0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
ϕ0.0286
(0.019)
0.0216
(0.018)
−0.0653
(0.018)
−0.0030
(0.019)
−0.0572
(0.019)
−0.0427
(0.019)
−0.0268
(0.018)
−0.0285
(0.018)
ω0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
α0.0501
(0.008)
0.0445
(0.009)
0.1132
(0.027)
0.0887
(0.018)
0.0833
(0.012)
0.0809
(0.011)
0.0670
(0.011)
0.0478
(0.007)
β0.9421
(0.010)
0.9497
(0.010)
0.8491
(0.038)
0.8762
(0.025)
0.9018
(0.014)
0.9093
(0.011)
0.9172
(0.014)
0.9484
(0.007)
ν8.6184
(1.478)
5.2529
(0.498)
5.8799
(0.298)
5.2686
(0.489)
6.7975
(0.812)
10.000
(1.422)
7.6179
(1.031)
7.8740
(1.028)
λ1.0288
(0.025)
0.9987
(0.025)
1.0167
(0.024)
1.0483
(0.027)
0.9385
(0.024)
0.882
(0.024)
0.9253
(0.023)
1.0003
(0.026)
GED
μ−0.0003
(0.000)
0.0009
(0.000)
0.0003
(0.000)
0.0007
(0.000)
0.0005
(0.000)
0.0007
(0.000)
0.0006
(0.000)
0.0007
(0.000)
ϕ0.0220
(0.019)
0.0103
(0.019)
−0.0676
(0.015)
−0.0165
(0.019)
−0.0490
(0.019)
−0.0356
(0.019)
−0.023
(0.020)
−0.0283
(0.019)
ω0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
α0.0444
(0.004)
0.0431
(0.008)
0.0394
(0.006)
0.0760
(0.017)
0.0851
(0.013)
0.0852
(0.012)
0.0676
(0.012)
0.0501
(0.018)
β0.9443
(0.006)
0.9516
(0.009)
0.9521
(0.007)
0.8873
(0.027)
0.8977
(0.015)
0.9014
(0.014)
0.9188
(0.015)
0.9439
(0.021)
η1.4626
(0.053)
1.2959
(0.046)
1.1214
(0.038)
1.2557
(0.041)
1.4437
(0.051)
1.6629
(0.061)
1.4599
(0.054)
1.4819
(0.056)
SGED
μ−0.0001
(0.001)
0.0010
(0.000)
0.0005
(0.000)
0.0010
(0.000)
0.0003
(0.000)
0.0005
(0.000)
0.0004
(0.000)
0.0006
(0.000)
ϕ0.0224
(0.019)
0.0111
(0.016)
−0.0695
(0.013)
−0.0124
(0.022)
−0.0601
(0.020)
−0.0429
(0.019)
−0.0287
(0.019)
−0.0290
(0.019)
ω0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
α0.0446
(0.004)
0.0436
(0.009)
0.0394
(0.006)
0.0777
(0.017)
0.0850
(0.013)
0.0846
(0.013)
0.0649
(0.012)
0.0500
(0.018)
β0.9439
(0.006)
0.9511
(0.010)
0.9520
(0.007)
0.8865
(0.026)
0.8972
(0.015)
0.9025
(0.015)
0.9217
(0.016)
0.9440
(0.021)
η1.4597
(0.023)
1.2970
(0.020)
1.1220
(0.019)
1.2533
(0.028)
1.4618
(0.023)
1.7051
(0.023)
1.4639
(0.023)
1.4834
(0.025)
ζ1.0266
(0.053)
1.0128
(0.046)
1.019
(0.038)0
1.0576
(0.040)
0.9339
(0.052)
0.8713
(0.064)
0.9331
(0.054)
0.9920
(0.057)
Standard error in parentheses.
Table A3. Estimation of the AR(1)-EGARCH(1,1) model.
Table A3. Estimation of the AR(1)-EGARCH(1,1) model.
EpistarKingspanBoralexNibeExxonChevronTotalSchlumberger
CoefficientNormal Distribution
μ0.0000 (0.0009)0.0000 (0.0004)0.0000 (0.0004)0.0000
(0.0003)
0.0000 (0.0002)0.0000 (0.0002)0.0000 (0.0002)0.0000 (0.0003)
ϕ0.055
(0.077)
0.019
(0.020)
−0.052
(0.020)
−0.002
(0.003)
−0.051
(0.020)
−0.022
(0.011)
−0.007
(0.019)
−0.027
(0.019)
ω−2.025
(0.491)
−0.307
(0.061)
−0.390
(0.075)
−0.258
(0.027)
−0.213
(0.021)
−0.157
(0.035)
−0.335
(0.043)
−0.329
(0.045)
α−0.095
(0.022)
−0.042
(0.013)
−0.069
(0.011)
−0.060
(0.011)
−0.069
(0.011)
−0.078
(0.012)
−0.104
(0.005)
−0.060
(0.012)
β0.712
(0.070)
0.957
(0.008)
0.948
(0.010)
0.965
(0.004)
0.975
(0.002)
0.981
(0.004)
0.959
(0.005)
0.957
(0.006)
γ0.278
(0.019)
0.232
(0.026)
0.188
(0.005)
0.154
(0.038)
0.177
(0.016)
0.158
(0.041)
0.160
(0.005)
0.178
(0.020)
Skew-Normal Distribution
μ−0.0001
(0.0006)
0.0007
(0.0003)
0.0002
(0.0002)
0.0011
(0.0004)
0.0001
(0.0002)
0.0002
(0.0003)
0.0001
(0.0003)
0.0003
(0.0003)
ϕ0.0488
(0.0192)
0.0219
(0.0172)
−0.0591
(0.0152)
−0.0085
(0.0196)
−0.0570
(0.0188)
−0.0232
(0.0189)
−0.0155
(0.0188)
−0.0225
(0.0191)
ω−0.1515
(0.0086)
−0.0308
(0.0005)
−0.0195
(0.0004)
−0.2447
(0.0089)
−0.2277
(0.0235)
−0.1516
(0.0328)
−0.1556
(0.0089)
−0.0828
(0.0011)
α−0.0324
(0.0079)
−0.0327
(0.0065)
−0.0456
(0.0047)
−0.0585
(0.0060)
−0.0700
(0.0111)
−0.0784
(0.0112)
−0.0803
(0.0107)
−0.0413
(0.0078)
β0.9782
(0.0012)
0.9955
(0.0001)
0.9968
(0.0001)
0.9667
(0.0013)
0.9731
(0.0027)
0.9820
(0.0039)
0.9808
(0.0011)
0.9890
(0.0002)
γ0.0935
(0.0118)
0.0770
(0.00275
0.0690
(0.0027)
0.1481
(0.0156)
0.1822
(0.0165)
0.1561
(0.0423)
0.1107
(0.0202)
0.1008
(0.0036)
ξ1.0226
(0.0240)
1.0219
(0.0219)
1.0123
(0.0193)
1.0390
(0.0211)
0.9006
(0.0212)
0.8569
(0.0230)
0.9351
(0.0219)
0.9757
(0.0226)
Student’s t Distribution
μ0.000
(0.001)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
ϕ0.036
(0.020)
0.023
(0.020)
−0.069
(0.019)
−0.007
(0.018)
−0.052
(0.018)
−0.027
(0.019)
−0.014
(0.020)
−0.025
(0.019)
ω−1.491
(0.341)
−0.562
(0.605)
−1.550
(0.534)
−0.311
(0.026)
−0.150
(0.013)
−0.138
(0.029)
−0.176
(0.002)
−0.068
(0.003)
α−0.095
(0.021)
−0.069
(0.029)
−0.038
(0.027)
−0.066
(0.015)
−0.063
(0.012)
−0.072
(0.013)
−0.083
(0.011)
−0.051
(0.009)
β0.788
(0.049)
0.925
(0.081)
0.798
(0.070)
0.959
(0.003)
0.983
(0.002)
0.984
(0.004)
0.979
(0.000)
0.991
(0.000)
γ0.289
(0.019)
0.290
(0.148)
0.414
(0.035)
0.174
(0.024)
0.172
(0.020)
0.150
(0.055)
0.118
(0.009)
0.104
(0.011)
ν9.972
(2.372)
5.947
(0.574)
5.621
(0.242)
5.344
(0.505)
7.021
(0.898)
13.361
(5.56)
8.410
(1.176)
8.356
(0.801)
Skewed-t Distribution
μ0.000
(0.001)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
0.000
(0.000)
ϕ−0.008
(0.018)
0.011
(0.018)
−0.070
(0.019)
−0.005
(0.023)
−0.057
(0.018)
−0.041
(0.019)
−0.027
(0.018)
−0.025
(0.015)
ω−0.117
(0.072)
−0.132
(0.075)
−1.574
(0.530)
−0.313
(0.026)
−0.150
(0.013)
−0.350
(0.093)
−0.365
(0.110)
−0.068
(0.002)
α−0.698
(0.116)
−0.558
(0.175)
−0.038
(0.027)
−0.066
(0.015)
−0.063
(0.012)
−1.083
(0.278)
−1.417
(0.200)
−0.051
(0.009)
β0.935
(0.032)
0.942
(0.028)
0.795
(0.070)
0.959
(0.003)
0.983
(0.001)
0.889
(0.027)
0.888
(0.032)
0.991
(0.000)
γ2.410
(0.330)
2.553
(0.517)
0.416
(0.034)
0.176
(0.024)
0.171
(0.020)
3.587
(0.145)
2.464
(0.155)
0.104
(0.011)
ν4.010
(0.001)
5.456
(0.000)
5.613
(0.242)
5.354
(0.506
7.229
(0.957)
5.812
(0.001)
4.513
(0.000)
8.357
(0.808)
λ1.034
(0.029)
0.988
(0.023)
1.011
(0.025)
1.045
(0.027)
0.935
(0.024)
0.926
(0.028)
0.945
(0.025)
0.995
(0.026)
GED
μ−0.0005
(0.001)
0.0007
(0.000)
0.0002
(0.000)
0.0004
(0.000)
0.0003
(0.000)
0.0003
(0.000)
0.0003
(0.000)
0.0004
(0.000)
ϕ0.0191
(0.014)
0.0073
(0.015)
−0.0685
(0.021)
−0.0180
(0.020)
−0.0477
(0.021)
−0.0232
(0.019)
−0.0215
(0.019)
−0.0245
(0.019)
ω−0.1301
(0.004)
−0.0371
(0.001)
−0.0612
(0.001)
−0.2399
(0.037)
−0.2105
(0.024)
−0.1526
(0.019)
−0.1609
(0.003)
−0.0761
(0.002)
α−0.0328
(0.009)
−0.0337
(0.008)
−0.0395
(0.009)
−0.0613
(0.014)
−0.0688
(0.014)
−0.0772
(0.012)
−0.0800
(0.010)
−0.0471
(0.009)
β0.9816
(0.001)
0.9952
(0.000)
0.9922
(0.000)
0.9684
(0.005)
0.9759
(0.003)
0.9823
(0.002)
0.9807
(0.001)
0.9903
(0.000)
γ0.1035
(0.008)
0.0796
(0.007)
0.0924
(0.009)
0.1409
(0.040)
0.1804
(0.019)
0.1570
(0.023)
0.1144
(0.009)
0.1016
(0.001)
η1.4724
(0.057)
1.3183
(0.046)
1.1360
(0.037)
1.2649
(0.041)
1.4612
(0.054)
1.7093
(0.068)
1.5078
(0.055)
1.4797
(0.051)
SGED
μ−0.0003
(0.001)
0.0008
(0.000)
0.0003
(0.000)
0.0007
(0.000)
0.0001
(0.000)
0.0002
(0.000)
0.0001
(0.000)
0.0003
(0.000)
ϕ0.0181
(0.020)
0.0075
(0.016)
−0.0698
(0.013)
−0.0159
(0.019)
−0.0579
(0.019)
−0.0275
(0.019)
−0.0239
(0.019)
−0.0254
(0.019)
ω−0.1310
(0.004)
−0.0369
(0.001)
−0.0608
(0.001)
−0.2357
(0.030)
−0.2234
(0.042)
−0.1475
(0.028)
−0.1566
(0.001)
−0.0766
(0.002)
α−0.0332
(0.009)
−0.0335
(0.008)
−0.0390
(0.009)
−0.0593
(0.014)
−0.0712
(0.014)
−0.0764
(0.012)
−0.0794
(0.010)
−0.0470
(0.009)
β0.9815
(0.001)
0.9952
(0.000)
0.9923
(0.000)
0.9690
(0.004)
0.9743
(0.005)
0.9827
(0.003)
0.9811
(0.000)
0.9902
(0.000)
γ0.1039
(0.008)
0.0800
(0.007)
0.0923
(0.009)
0.1446
(0.033)
0.1815
(0.021)
0.1534
(0.038)
0.1092
(0.007)
0.1014
(0.001)
η1.4707
(0.057)
1.3192
(0.046)
1.1362
(0.037)
1.2652
(0.041)
1.4805
(0.055)
1.7464
(0.082)
1.5056
(0.055)
1.4806
(0.051)
ζ1.0310
(0.024)
1.0063
(0.020)
1.0131
(0.019)
1.0514
(0.026)
0.9247
(0.022)
0.8641
(0.024)
0.9294
(0.022)
0.9879
(0.025)
The standard error in parentheses.
Table A4. Estimation of the AR(1)-APARCH(1,1) model.
Table A4. Estimation of the AR(1)-APARCH(1,1) model.
EpistarKingspanBoralexNibeExxonChevronTotalSchlumberger
CoefficientNormal Distribution
μ0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
ϕ0.0463
(0.018)
0.0145
(0.014)
−0.0511
(0.019)
0.0019
(0.019)
−0.0529
(0.020)
−0.0246
(0.020)
−0.0151
(0.016)
−0.0218
(0.19)
ω0.0005
(0.000)
0.0003
(0.000)
0.0000
(0.000)
0.0018
(0.000)
0.0001
(0.000)
0.0001
(0.000)
0.0009
(0.000)
0.0000
(0.000)
α0.0498
(0.007)
0.0361
(0.05)
0.0307
(0.005)
0.0844
(0.010)
0.0881
(0.011)
0.0802
(0.010)
0.0653
(0.009)
0.0500
(0.007)
β0.9396
(0.010)
0.9690
(0.003)
0.9697
(0.003)
0.9020
(0.013)
0.8988
(0.011)
0.9124
(0.009)
0.9279
(0.011)
0.9471
(0.006)
γ0.3283
(0.094)
0.4062
(0.150)
0.4534
(0.108)
0.4670
(0.106)
0.3935
(0.076)
0.5248
(0.097)
0.8203
(0.112)
0.3714
(0.098)
δ1.0739
(0.263)
0.7834
(0.253)
1.6283
(0.232)
0.7959
(0.159)
1.3440
(0.214)
1.2293
(0.202)
0.8168
(0.167)
1.4029
(0.226)
Skew-Normal Distribution
μ−0.0001
(0.0005)
0.0007
(0.0003)
0.0004
(0.0003)
0.0011
(0.0002)
0.0001
(0.0002)
0.0002
(0.0002)
0.0000
(0.0001)
0.0004
(0.0003)
ϕ0.0498
(0.0215)
0.0206
(0.0130)
−0.0590
(0.0190)
−0.0029
(0.0038)
−0.0553
(0.0196)
−0.0262
(0.0195)
−0.0167
(0.0042)
−0.0188
(0.0194)
ω0.0005
(0.0005)
0.0002
(0.0000)
0.0000
(0.0000)
0.0018
(0.0012)
0.0001
(0.0001)
0.0001
(0.0001)
0.0010
(0.0007)
0.0001
(0.0001)
α0.0495
(0.0103)
0.0373
(0.0013)
0.0305
(0.0020)
0.0795
(0.0108)
0.0947
(0.0114)
0.0807
(0.0103)
0.0544
(0.0155)
0.0499
(0.0066)
β0.9378
(0.0156)
0.9675
(0.0003)
0.9682
(0.0015)
0.9062
(0.0132)
0.8916
(0.0117)
0.9118
(0.0100)
0.9422
(0.0216)
0.9474
(0.0061)
γ0.3413
(0.1006)
0.4763
(0.0991)
0.4041
(0.0759)
0.4901
(0.1130)
0.3914
(0.0754)
0.4775
(0.0945)
0.9037
(0.1513)
0.3765
(0.1012)
δ1.0937
(0.2763)
0.8959
(0.0482)
1.7392
(0.0380)
0.7835
(0.1661)
1.3313
(0.2018)
1.2970
(0.1934)
0.7297
(0.2298)
1.3378
(0.2035)
ξ1.0235
(0.0239)
1.0258
(0.0220)
1.0135
(0.0198)
1.0402
(0.0215)
0.8990
(0.0212)
0.8575
(0.0225)
0.9351
(0.0219)
0.9732
(0.0226)
Student’s t Distribution
μ0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
ϕ0.0240
(0.019)
0.0169
(0.019)
−0.0660
(0.018)
−0.0056
(0.017)
−0.0516
(0.018)
−0.0295
(0.019)
−0.0196
(0.020)
−0.0228
(0.018)
ω0.0003
(0.000)
0.0002
(0.000)
0.0001
(0.000)
0.0011
(0.000)
0.0001
(0.000)
0.0001
(0.000)
0.0010
(0.000)
0.0001
(0.000)
α0.0599
(0.009)
0.0440
(0.008)
0.1036
(0.029)
0.0842
(0.014)
0.0940
(0.014)
0.0764
(0.011)
0.0629
(0.009)
0.0514
(0.008)
β0.9374
(0.010)
0.9611
(0.006)
0.8826
(0.037)
0.9087
(0.016)
0.8995
(0.013)
0.9190
(0.010)
0.9309
(0.012)
0.9494
(0.007)
γ0.2804
(0.090)
0.4878
(0.140)
0.1524
(0.073)
0.4344
(0.112)
0.3965
(0.089)
0.4985
(0.108)
0.8705
(0.115)
0.4807
(0.115)
δ1.1112
(0.253)
0.9558
(0.173)
1.6047
(0.266)
0.8777
(0.182)
1.2454
(0.235)
1.2586
(0.227)
0.7786
(0.180)
1.2967
(0.230)
ν8.9644
(0.1609)
5.4858
(0.534)
5.9731
(0.315)
5.3611
(0.500)
7.0192
(0.854)
10.000
1.381
8.5843
(1.272)
8.1418
(1.068)
Skewed-t Distribution
μ0.0000
(0.001)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
0.0000
(0.000)
ϕ0.0247
(0.020)
0.0170
(0.019)
−0.0663
(0.018)
−0.0042
(0.018)
−0.0547
(0.020)
−0.0335
(0.019)
−0.0226
(0.014)
−0.0227
(0.018)
ω0.0003
(0.000)
0.0002
(0.000)
0.0001
(0.000)
0.0010
(0.000)
0.0001
(0.000)
0.0001
(0.000)
0.0010
(0.000)
0.0001
(0.000)
α0.0604
(0.009)
0.0440
(0.008)
0.1040
(0.029)
0.0850
(0.014)
0.0939
(0.014)
0.0759
(0.011)
0.0613
(0.009)
0.0514
(0.008)
β0.9369
(0.011)
0.9611
(0.006)
0.8819
(0.037)
0.9084
(0.016)
0.8985
(0.013)
0.9200
(0.010)
0.9324
(0.012)
0.9494
(0.007)
γ0.2915
(0.091)
0.4885
(0.141)
0.1508
(0.074)
0.4317
(0.112)
0.3997
(0.089)
0.4903
(0.108)
0.8812
(0.116)
0.4808
(0.116)
δ1.0959
(0.251)
0.9560
(0.173)
1.6053
(0.266)
0.8782
(0.183)
1.2520
(0.233)
1.2570
(0.219)
0.7915
(0.183)
1.2972
(0.231)
ν8.8193
(1.553)
5.4835
(0.534)
5.9717
(0.315)
5.3631
(0.500)
7.2550
(0.914)
10.0000
(1.352)
8.4353
(1.239)
8.1419
(1.068)
λ1.0373
(0.026)
0.9979
(0.025)
1.0150
(0.024)
1.0496
(0.028)
0.9347
(0.024)
0.8833
(0.025)
0.9256
(0.023)
1.0023
(0.026)
GED
μ−0.0005
(0.000)
0.0007
(0.000)
0.0002
(0.000)
0.0004
(0.000)
0.0003
(0.000)
0.0004
(0.000)
0.0002
(0.000)
0.0004
(0.000)
ϕ0.0192
(0.015)
0.0069
(0.018)
−0.0684
(0.017)
−0.0151
(0.020)
−0.0483
(0.019)
−0.0270
(0.020)
−0.0245
(0.005)
−0.0232
(0.019)
ω0.0005
(0.000)
0.0001
(0.000)
0.0000
(0.000)
0.0010
(0.001)
0.0001
(0.000)
0.0001
(0.000)
0.0011
(0.002)
0.0001
(0.000)
α0.0550
(0.013)
0.0386
(0.002)
0.0412
(0.010)
0.0770
(0.014)
0.0933
(0.013)
0.0807
(0.011)
0.0551
(0.043)
0.0505
(0.007)
β0.9378
(0.018)
0.9657
(0.000)
0.9536
(0.003)
0.9094
(0.018)
0.8958
(0.013)
0.9116
(0.011)
0.9416
(0.060)
0.9488
(0.006)
γ0.3166
(0.104)
0.4740
(0.124)
0.2562
(0.067)
0.4785
(0.132)
0.3864
(0.088)
0.4722
(0.102)
0.8899
(0.355)
0.4355
(0.116)
δ1.0750
(0.289)
0.9619
(0.060)
1.7863
(0.045)
0.9253
(0.210)
1.3296
(0.232)
1.2952
(0.218)
0.7131
(0.607)
1.2793
(0.223)
η1.4721
(0.059)
1.3237
(0.046)
1.1369
(0.038)
1.2715
(0.041)
1.4668
(0.052)
1.7166
(0.064)
1.5221
(0.056)
1.4899
(0.052)
SGED
μ−0.0003
(0.001)
0.0008
(0.000)
0.0003
(0.000)
0.0007
(0.000)
0.0001
(0.000)
0.0002
(0.000)
0.0001
(0.000)
0.0004
(0.000)
ϕ0.0184
(0.020)
0.0075
(0.018)
−0.0704
(0.017)
−0.0130
(0.010)
−0.0565
(0.020)
−0.0303
(0.020)
−0.0258(0.006)−0.0241
(0.019)
ω0.0005
(0.000)
0.0001
(0.000)
0.0000
(0.000)
0.0010
(0.001)
0.0001
(0.000)
0.0001
(0.000)
0.0011
(0.001)
0.0001
(0.000)
α0.0552
(0.014)
0.0388
(0.002)
0.0409
(0.018)
0.0785
(0.014)
0.0946
(0.014)
0.0794
(0.011)
0.0532
(0.003)
0.0504
(0.007)
β0.9373
(0.019)
0.9656
(0.000)
0.9537
(0.002)
0.9088
(0.017)
0.8936
(0.013)
0.9133
(0.011)
0.9438
(0.003)
0.9488
(0.006)
γ0.3225
(0.105)
0.4689
(0.124)
0.2533
(0.011)
0.4546
(0.124)
0.3995
(0.088)
0.4713
(0.101)
0.9123
(0.024)
0.4366
(0.116)
δ1.0626
(0.289)
0.9603
(0.060)
1.7912
(0.109)
0.9259
(0.178)
1.3185
(0.234)
1.3070
(0.209)
0.7026
(0.129)
1.2759
(0.223)
η1.4705
(0.059)
1.3252
(0.046)
1.1374
(0.038)
1.2718
(0.041)
1.4876
(0.053)
1.7570
(0.067)
1.5176
(0.056)
1.4911
(0.052)
ζ1.0316
(0.024)
1.0082
(0.021)
1.0160
(0.023)
1.0526
(0.024)
0.9239
(0.023)
0.8652
(0.023)
0.9278
(0.022)
0.9878
(0.024)
Standard error in parentheses.
Table A5. Two-sample Kolmogorov-Smirnov tests, skewness and kurtosis estimation of standardized residuals from AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) models.
Table A5. Two-sample Kolmogorov-Smirnov tests, skewness and kurtosis estimation of standardized residuals from AR(1)-GARCH(1,1), AR(1)-EGARCH(1,1) and AR(1)-APARCH(1,1) models.
Normal Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
KS testSkewnessKurtosisKS testSkewnessKurtosisKS testSkewnessKurtosis
Epistar0.34974−0.025093.560910.585780.010603.513670.564050.013693.52055
Kingspan0.00001−0.026225.976710.000130.046305.458960.000110.058065.42608
Boralex0.00644−0.100206.772460.011800.013806.062450.007680.011856.06490
Nibe0.001280.029818.257700.001570.094678.134980.000850.097497.76289
Exxon Mobil0.00001−0.324294.453710.00017−0.338254.471720.00015−0.339224.43026
Chevron0.00001−0.352443.870620.00003−0.348353.769990.00003−0.345223.73576
Total S.A0.00095−0.195594.097680.00914−0.150773.878620.01180−0.146973.83074
Schlumberger0.00077−0.146724.844470.00281−0.130055.284810.00157−0.153775.20671
Skew-Normal Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
KS testSkewnessKurtosisKS testSkewnessKurtosisKS testSkewnessKurtosis
Epistar0.34974−0.024733.559320.542570.011673.512330.521380.015073.51885
Kingspan0.00001−0.023685.964160.000090.048785.456970.000090.062105.42004
Boralex0.00492−0.099726.763910.010850.015326.058710.007040.012976.05604
Nibe0.000690.034018.262140.000950.102238.140030.000450.104697.75844
Exxon Mobil0.00002−0.326794.457980.00026−0.341304.476210.00029−0.342874.43726
Chevron0.00001−0.357263.881380.00005−0.351963.782880.00004−0.348443.74654
Total S.A0.00128−0.198844.104540.01515−0.152593.880380.01180−0.148123.83148
Schlumberger0.00095−0.147354.841260.00309−0.131185.274910.00191−0.154415.19833
Student’s t Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
KS testSkewnessKurtosisKS testSkewnessKurtosisKS testSkewnessKurtosis
Epistar0.56405−0.039863.701860.781830.004933.639430.739660.008403.64400
Kingspan0.00007−0.037566.038030.000330.063065.461810.000370.071295.44417
Boralex0.07290−0.197739.380790.05533−0.116417.516640.07797−0.102448.02954
Nibe0.012840.042348.330480.035830.122758.246160.026450.124797.96779
Exxon Mobil0.00000−0.325804.480210.00003−0.337004.482180.00002−0.339004.45428
Chevron0.00000−0.356613.885720.00001−0.352613.788920.00001−0.349273.75188
Total S.A0.00077−0.198714.104700.00492−0.152343.880050.00538−0.147873.83225
Schlumberger0.00041−0.137534.891480.00173−0.115635.508720.00142−0.143055.42437
Skewed-t Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
KS testSkewnessKurtosisKS testSkewnessKurtosisKS testSkewnessKurtosis
Epistar0.34974−0.038873.699550.607720.006663.638090.585780.010313.64268
Kingspan0.00003−0.035786.028670.000210.063305.461920.000230.071615.44377
Boralex0.02264−0.193649.330150.03325−0.119927.561090.04153−0.104098.04632
Nibe0.001730.045928.330990.005380.125638.235260.003400.127007.95628
Exxon Mobil0.00002−0.326524.482990.00011−0.337884.482420.00012−0.340414.45491
Chevron0.00001−0.359873.896530.00004−0.355303.798120.00002−0.351183.75814
Total S.A0.00173−0.201794.115350.01932−0.153293.880200.01515−0.148293.83103
Schlumberger0.00041−0.137684.891710.00173−0.115665.507290.00142−0.143025.42275
GED
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
KS testSkewnessKurtosisKS testSkewnessKurtosisKS testSkewnessKurtosis
Epistar0.71801−0.039083.659050.629780.004833.603800.629780.009273.60931
Kingspan0.00002−0.039736.043440.000110.052765.469520.000110.061415.44694
Boralex0.01782−0.128917.322840.01395−0.048016.605000.00838−0.028646.57144
Nibe0.013950.034818.269400.033250.109898.180710.028570.111327.89423
Exxon Mobil0.00000−0.323784.464750.00006−0.337624.475450.00004−0.338924.44059
Chevron0.00000−0.355073.874630.00001−0.350643.778250.00001−0.347563.74339
Total S.A0.00037−0.198354.099000.00309−0.152313.879940.00492−0.148173.83287
Schlumberger0.00062−0.142454.866410.00173−0.121855.399270.00128−0.147575.31753
SGED
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
KS testSkewnessKurtosisKS testSkewnessKurtosisKS testSkewnessKurtosis
Epistar0.50054−0.038163.655810.607720.006243.603130.564050.011033.60817
Kingspan0.00001−0.037596.032330.000080.053355.469250.000070.062305.44584
Boralex0.00838−0.126637.309070.00644−0.048046.608180.00538−0.027656.56856
Nibe0.001420.039078.274740.006440.114838.177570.005380.115617.89235
Exxon Mobil0.00003−0.326394.469750.00023−0.339644.477820.00017−0.341404.44407
Chevron0.00001−0.358793.887910.00005−0.353973.790380.00003−0.350023.75210
Total S.A0.00142−0.201694.110090.01395−0.153563.880530.01515−0.148363.83142
Schlumberger0.00085−0.142574.864690.00173−0.122185.394940.00157−0.147735.31342
KS test column shows p-values in two-sample Kolmogorov-Smirnov test. The null hypothesis indicates that standardized returns and standardized residuals were drawn from the same continuous distribution.
Table A6. Ljung-Box test of standardized residuals from AR(1)-GARCH(1,1), AR(1)-GARCH(1,1) and AR(1)-APARCH(1,1) models.
Table A6. Ljung-Box test of standardized residuals from AR(1)-GARCH(1,1), AR(1)-GARCH(1,1) and AR(1)-APARCH(1,1) models.
Normal Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.750.570.660.600.640.970.660.660.630.650.850.600.900.520.960.680.640.610.930.750.730.650.980.74
20.940.680.900.790.860.400.850.500.890.610.970.810.990.500.970.550.900.591.000.880.920.520.970.54
30.930.780.970.630.590.460.850.200.890.730.990.650.710.630.890.210.900.701.000.740.610.620.870.18
40.580.310.990.270.710.580.920.290.570.351.000.270.810.730.960.280.580.361.000.310.710.710.950.26
50.710.100.480.310.440.670.970.300.700.100.530.330.640.840.980.350.710.100.520.370.590.820.980.31
60.770.080.590.430.560.470.940.190.770.060.660.450.750.590.970.200.780.070.640.490.700.560.970.17
70.640.110.700.370.590.570.970.260.620.090.760.360.750.690.990.250.630.100.750.380.730.670.980.23
80.580.150.120.470.490.450.770.270.540.120.100.470.570.590.820.290.550.140.110.480.590.550.840.25
120.840.090.110.270.580.610.440.290.820.060.080.270.690.690.500.330.840.060.090.250.690.650.490.28
200.450.080.230.230.490.450.520.400.440.040.210.220.590.580.670.430.460.050.240.230.610.550.680.41
Skew-Normal Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.780.590.670.630.390.540.460.580.660.660.830.610.630.790.950.640.670.610.920.760.610.920.990.70
20.950.690.910.810.670.340.710.480.910.610.970.810.890.590.960.540.910.590.990.880.860.570.960.53
30.930.790.970.640.500.410.760.190.900.730.990.650.680.710.890.210.900.701.000.740.600.660.870.18
40.590.320.990.280.630.530.870.280.580.351.000.270.780.800.960.280.590.361.000.310.700.750.950.26
50.710.110.480.320.390.640.940.300.700.100.530.330.630.890.980.350.710.100.520.370.590.850.980.31
60.770.080.590.440.510.440.910.190.770.060.650.450.730.640.970.200.780.070.640.500.700.580.970.17
70.640.120.700.370.540.550.950.260.620.090.760.370.740.730.990.250.630.100.740.380.730.690.980.23
80.580.150.120.480.450.430.740.260.550.120.100.470.560.620.820.290.550.140.110.480.590.570.840.25
120.840.100.110.280.540.600.420.280.830.060.080.270.680.730.510.330.840.060.090.250.690.670.500.28
200.450.080.230.230.460.440.500.390.440.040.210.220.570.600.670.430.460.050.230.240.610.560.690.41
Student’s t Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.170.230.200.390.610.820.530.550.130.270.160.560.880.730.790.650.130.270.230.600.770.880.810.63
20.380.380.380.610.840.400.760.450.320.380.350.750.980.590.940.520.320.370.430.770.930.590.950.49
30.540.520.560.540.580.460.790.180.470.520.540.640.690.690.870.180.470.500.610.690.600.660.850.15
40.360.190.720.230.700.580.890.260.320.240.710.290.790.790.950.240.330.250.770.310.710.750.940.21
50.490.050.360.280.450.680.950.290.450.070.350.360.630.880.980.320.450.080.400.390.590.850.980.28
60.580.040.480.390.570.470.920.190.540.040.470.480.740.640.960.190.550.050.530.510.710.600.960.16
70.470.060.600.330.600.580.960.250.440.070.590.400.750.730.980.240.440.080.640.410.730.700.980.21
80.430.080.180.430.490.450.750.270.390.100.130.500.570.630.800.280.390.110.170.520.590.590.830.23
120.730.060.110.270.570.620.430.300.700.050.090.300.690.730.470.330.720.050.110.290.690.690.470.27
200.390.050.220.220.500.450.510.400.380.030.190.240.590.590.650.430.390.030.240.250.620.560.670.40
Skewed-t Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.170.230.200.390.610.820.530.550.130.270.160.560.880.730.790.650.130.270.230.600.770.880.810.63
20.380.380.380.610.840.400.760.450.320.380.350.750.980.590.940.520.320.370.430.770.930.590.950.49
30.540.520.560.540.580.460.790.180.470.520.540.640.690.690.870.180.470.500.610.690.600.660.850.15
40.360.190.720.230.700.580.890.260.320.240.710.290.790.790.950.240.330.250.770.310.710.750.940.21
50.490.050.360.280.450.680.950.290.450.070.350.360.630.880.980.320.450.080.400.390.590.850.980.28
60.580.040.480.390.570.470.920.190.540.040.470.480.740.640.960.190.550.050.530.510.710.600.960.16
70.470.060.600.330.600.580.960.250.440.070.590.400.750.730.980.240.440.080.640.410.730.700.980.21
80.430.080.180.430.490.450.750.270.390.100.130.500.570.630.800.280.390.110.170.520.590.590.830.23
120.730.060.110.270.570.620.430.300.700.050.090.300.690.730.470.330.720.050.110.290.690.690.470.27
200.390.050.220.220.500.450.510.400.380.030.190.240.590.590.650.430.390.030.240.250.620.560.670.40
GED
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.070.150.380.210.790.730.430.430.050.230.370.330.990.770.740.540.050.220.470.360.890.930.710.53
20.180.270.680.410.920.390.690.410.160.330.670.561.000.600.920.490.150.320.760.590.960.600.920.46
30.300.400.840.390.620.450.740.160.260.460.840.510.710.710.850.180.250.440.900.570.620.670.830.15
40.200.150.940.160.730.570.850.240.180.210.930.220.800.800.940.240.180.210.960.250.720.760.930.22
50.300.040.450.200.460.670.930.270.280.060.480.290.640.890.970.320.270.060.500.320.600.860.970.28
60.380.030.570.290.580.470.900.170.360.030.610.400.740.650.960.180.360.040.630.430.710.600.960.15
70.310.050.680.250.610.570.950.240.290.050.720.330.750.740.980.230.280.060.740.350.740.710.980.20
80.280.060.140.340.500.460.720.250.250.080.130.430.570.640.790.270.250.080.140.450.590.600.820.23
120.570.040.110.210.580.620.400.280.550.040.090.260.690.740.460.320.550.040.100.260.690.690.460.27
200.280.040.230.180.500.450.490.380.280.020.210.210.590.600.640.430.280.030.240.230.620.570.660.40
SGED
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.070.160.330.300.420.470.290.410.050.230.340.390.590.960.630.510.050.230.410.430.570.920.660.49
20.190.290.620.520.690.320.540.400.140.330.630.630.870.620.870.480.140.330.710.660.830.580.880.45
30.310.420.800.470.510.390.630.160.240.460.810.550.660.730.830.170.230.450.860.600.580.670.830.15
40.200.150.910.190.630.520.770.240.170.210.920.230.760.810.930.240.170.220.940.260.690.760.920.21
50.300.040.430.240.400.620.870.260.260.060.460.300.620.900.970.310.250.060.480.330.580.860.970.27
60.390.030.550.340.520.430.850.170.340.030.590.410.730.650.950.180.340.040.610.450.690.590.960.15
70.310.050.660.290.550.540.910.230.270.050.700.340.740.740.980.230.270.060.720.360.730.700.980.20
80.280.070.130.380.450.420.680.240.240.080.120.440.560.630.800.270.240.090.130.460.580.580.820.23
120.580.040.100.240.550.590.380.280.530.040.080.260.680.740.480.320.540.040.100.260.680.680.480.27
200.290.040.220.200.470.430.470.380.270.020.200.220.580.600.650.430.270.030.230.230.610.570.680.39
Figures correspond to p-values for the Ljung-Box test, which null hypothesis is the absence of autocorrelation, m stands for lag order.
Table A7. Ljung-Box test of standardized squared residuals from AR(1)-GARCH(1,1), AR(1)-GARCH(1,1) and AR(1)-APARCH(1,1) models.
Table A7. Ljung-Box test of standardized squared residuals from AR(1)-GARCH(1,1), AR(1)-GARCH(1,1) and AR(1)-APARCH(1,1) models.
Normal Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.150.000.000.720.280.170.240.600.190.000.000.090.430.220.380.570.230.000.000.030.300.160.580.98
20.250.010.000.930.080.010.360.300.360.010.000.220.250.150.400.340.430.010.000.090.340.120.260.45
30.430.030.000.750.170.020.260.480.560.020.000.260.430.090.480.530.640.020.000.120.540.070.260.63
40.380.040.000.780.280.030.230.620.390.040.000.390.580.140.360.680.410.030.000.210.710.100.210.77
50.190.060.000.830.400.050.340.500.220.040.000.530.720.180.500.620.220.040.000.320.820.160.320.64
60.260.080.000.900.510.070.300.620.310.080.000.650.780.210.280.740.310.070.000.430.860.180.150.76
70.180.120.000.950.480.080.110.540.250.120.000.730.790.280.150.530.250.110.000.520.840.240.100.52
80.230.170.000.970.410.060.160.650.270.180.000.820.410.120.220.620.260.160.000.620.490.100.150.63
120.350.370.000.990.440.030.290.820.360.390.000.950.260.030.390.740.360.360.000.850.400.020.320.79
200.120.410.001.000.650.150.610.850.200.330.000.990.700.140.690.840.180.280.000.970.770.120.630.89
Skew-Normal Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.150.000.000.730.270.160.280.600.190.000.000.100.420.250.430.570.230.000.000.030.290.170.610.98
20.250.010.000.940.100.020.380.310.360.010.000.250.300.130.390.340.430.010.000.100.370.110.240.45
30.420.030.000.740.200.020.240.490.560.020.000.280.490.090.450.530.640.020.000.120.580.070.230.63
40.380.040.000.780.320.040.210.620.390.040.000.410.640.140.320.680.410.040.000.210.740.100.180.77
50.190.060.000.820.450.060.320.500.220.050.000.540.770.180.450.630.220.040.000.320.840.150.280.64
60.260.080.000.900.570.070.260.630.320.080.000.670.820.210.240.750.320.070.000.440.880.180.120.76
70.190.130.000.940.520.080.100.540.260.130.000.750.830.280.130.530.260.120.000.530.860.230.080.52
80.230.180.000.970.450.070.150.650.280.190.000.830.420.120.190.620.270.170.000.630.480.100.130.63
120.350.380.000.990.490.040.270.820.360.400.000.950.260.030.350.740.360.370.000.850.380.020.280.79
200.120.410.001.000.690.180.590.850.200.330.000.990.700.130.650.840.180.280.000.970.760.120.580.89
Student’s t Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.370.000.140.570.230.210.340.700.400.010.000.070.370.290.390.480.470.010.020.050.270.200.620.79
20.620.010.280.830.100.010.360.310.690.020.000.190.280.100.380.400.770.020.070.130.340.100.220.50
30.780.020.340.720.190.020.190.490.840.040.000.220.460.080.450.570.870.040.090.150.530.070.220.66
40.430.030.430.770.320.030.150.640.360.060.000.340.620.130.340.730.340.060.140.250.700.100.180.80
50.140.00.580.830.440.050.240.490.150.080.000.470.750.160.480.640.130.070.230.370.810.150.280.66
60.180.060.670.910.560.060.180.620.200.130.010.600.810.190.270.760.180.120.290.490.850.170.130.77
70.090.100.770.950.490.070.070.560.130.190.010.690.810.260.140.540.110.180.390.590.820.230.090.52
80.140.150.770.980.470.060.110.660.160.270.010.780.490.110.200.630.140.250.390.690.530.090.130.63
120.230.360.870.980.540.030.210.840.230.480.030.930.370.030.370.800.210.450.510.880.470.020.280.82
200.100.420.340.990.690.150.510.850.170.360.010.980.780.130.680.870.140.320.150.970.810.120.590.91
Skewed-t Distribution
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.370.000.140.570.230.210.340.700.400.010.000.070.370.290.390.480.470.010.020.050.270.200.620.79
20.620.010.280.830.100.010.360.310.690.020.000.190.280.100.380.400.770.020.070.130.340.100.220.50
30.780.020.340.720.190.020.190.490.840.040.000.220.460.080.450.570.870.040.090.150.530.070.220.66
40.430.030.430.770.320.030.150.640.360.060.000.340.620.130.340.730.340.060.140.250.700.100.180.80
50.140.050.580.830.440.050.240.490.150.080.000.470.750.160.480.640.130.070.230.370.810.150.280.66
60.180.060.670.910.560.060.180.620.200.130.010.600.810.190.270.760.180.120.290.490.850.170.130.77
70.090.100.770.950.490.070.070.560.130.190.010.690.810.260.140.540.110.180.390.590.820.230.090.52
80.140.150.770.980.470.060.110.660.160.270.010.780.490.110.200.630.140.250.390.690.530.090.130.63
120.230.360.870.980.540.030.210.840.230.480.030.930.370.030.370.800.210.450.510.880.470.020.280.82
200.100.420.340.990.690.150.510.850.170.360.010.980.780.130.680.870.140.320.150.970.810.120.590.91
GED
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.240.000.000.630.250.190.290.670.300.010.000.090.400.250.390.500.340.000.000.050.280.180.600.83
20.430.010.000.880.080.010.360.310.550.010.000.220.260.130.370.370.620.010.000.140.330.110.230.47
30.630.020.000.730.170.020.220.500.740.030.000.240.440.090.450.550.790.030.000.160.530.070.230.65
40.420.030.000.780.280.030.180.640.380.050.000.370.600.140.340.710.370.050.000.270.700.100.190.78
50.150.050.000.840.410.050.280.500.160.060.000.510.730.170.480.630.150.060.000.400.810.150.290.65
60.200.070.000.910.520.060.230.620.230.110.000.640.790.210.270.750.210.100.000.520.850.180.130.76
70.120.100.000.950.470.070.090.550.160.160.000.720.800.280.150.540.150.150.000.610.830.240.090.52
80.160.150.000.980.430.060.130.660.190.230.000.810.430.110.210.630.170.220.000.710.500.090.140.63
120.280.360.000.990.470.030.240.840.280.450.000.950.300.030.370.780.270.420.010.900.420.020.290.81
200.110.420.001.000.670.150.560.850.180.350.000.990.730.130.690.860.160.310.000.980.780.110.600.90
SGED
AR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
mEpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKingspanBoralexNibeExxon MobilChevronTotal S.ASchlum
berger
EpistarKing
span
BoralexNibeExxon MobilChevronTotal S.ASchlum
berger
10.240.000.000.640.260.180.340.670.300.010.000.100.420.270.440.490.340.000.000.050.290.180.650.82
20.430.010.000.880.100.010.380.310.550.010.000.230.310.110.380.370.610.010.000.150.380.100.220.47
30.630.020.000.730.200.020.210.500.750.030.000.260.510.090.430.550.790.030.000.170.580.070.200.64
40.430.030.000.780.320.030.160.640.380.050.000.390.660.140.310.710.370.050.000.280.740.100.160.78
50.160.050.000.830.450.050.260.490.160.060.000.520.790.170.440.630.150.060.000.400.850.150.250.65
60.210.070.000.910.570.070.200.620.230.110.000.650.840.200.230.750.220.100.000.530.880.180.110.76
70.120.110.000.950.510.070.080.550.160.160.000.740.840.270.120.540.150.150.000.620.860.230.070.52
80.170.160.000.970.460.060.120.660.190.230.000.820.460.110.180.630.180.220.000.720.510.090.110.63
120.290.370.000.990.520.040.230.840.280.460.000.950.300.030.330.780.270.430.010.900.420.020.250.80
200.120.420.001.000.700.170.530.850.190.360.000.990.740.120.640.860.160.310.000.980.790.120.550.90
Figures correspond to p-values for the Ljung-Box test, which null hypothesis is the absence of autocorrelation, m stands for lag order.
Figure A1. Prices for different energy stocks.
Figure A1. Prices for different energy stocks.
Energies 13 02805 g0a1
Figure A2. Log-returns for different energy stocks.
Figure A2. Log-returns for different energy stocks.
Energies 13 02805 g0a2

References

  1. Bloomberg News. Available online: https://www.bloomberg.com/news/articles/2019-06-25/for-first-time-ever-renewables-surpass-coal-in-u-s-power-mix. (accessed on 25 December 2019).
  2. World Economic Forum. Available online: https://www.weforum.org/agenda/2019/09/what-the-saudis-can-teach-australia-about-the-end-of-fossil-fuels/. (accessed on 25 December 2019).
  3. Askari, H.; Krichene, N. Oil price dynamics (2002–2006). Energy Econ. 2008, 30, 2134–2153. [Google Scholar] [CrossRef]
  4. Russell, F.T.S.E. Available online: http://www.ftserussell.com/sites/default/files/ftse_russell_investing_in_the_global_green_economy_busting_common_myths_may_2018.pdf. (accessed on 25 December 2019).
  5. The economist. Available online: https://www.economist.com/special-report/2018/03/15/switching-to-renewables-will-not-be-as-rapid-as-many-hope. (accessed on 25 December 2019).
  6. The economist. Available online: https://www.economist.com/united-states/2017/11/25/the-keystone-xl-pipeline-has-won-approval-in-nebraska. (accessed on 25 December 2019).
  7. The economist. Available online: https://www.economist.com/special-report/2016/11/24/the-future-of-oil. (accessed on 25 December 2019).
  8. The economist. Available online: https://www.economist.com/business/2016/09/08/breaking-bad. (accessed on 25 December 2019).
  9. Acerbi, C.; Székely, B. Back-testing expected shortfall. Risk Mag. 2014, 27, 76–81. [Google Scholar]
  10. Acerbi, C.; Székely, B. General Properties of Backtestable Statistics. 2017. Available online: https://ssrn.com/abstract=2905109 (accessed on 25 December 2019).
  11. Aloui, C.; Mabrouk, S. Value-at-risk estimations of energy commodities via long-memory, asymmetry, and fat-tailed GARCH models. Energy Policy 2010, 38, 2326–2339. [Google Scholar] [CrossRef]
  12. Feng, Z.H.; Wei, Y.M.; Wang, K. Estimating risk for the carbon market via extreme value theory: An empirical analysis of the EU ETS. Appl. Energy 2012, 99, 97–108. [Google Scholar] [CrossRef]
  13. Artzner, P.; Delbaen, F.; Eber, J.M.; Heath, D. Coherent measures of risk. Math. Financ. 1999, 9, 203–228. [Google Scholar] [CrossRef]
  14. Acerbi, C.; Tasche, D. On the coherence of expected shortfall. J. Bank. Financ. 2002, 26, 1487–1503. [Google Scholar] [CrossRef] [Green Version]
  15. Gneiting, T. Making and evaluating point forecasts. J. Am. Stat. Assoc. 2011, 106, 746–762. [Google Scholar] [CrossRef] [Green Version]
  16. Costanzino, N.; Curran, M. Backtesting General Spectral Risk Measures with Application to Expected Shortfall. SSRN 2514403. 2015. Available online: https://ssrn.com/abstract=2514403 (accessed on 25 December 2019).
  17. Du, Z.; Escanciano, J.C. Backtesting expected shortfall: Accounting for tail risk. Manag. Sci. 2016, 63, 940–958. [Google Scholar] [CrossRef] [Green Version]
  18. Fissler, T.; Ziegel, J.F. Higher order elicitability and Osband’s principle. Ann. Stat. 2016, 44, 1680–1707. [Google Scholar] [CrossRef]
  19. Fissler, T.; Ziegel, J.F.; Gneiting, T. Expected Shortfall is jointly elicitable with Value at Risk-Implications for backtesting. Risk 2016, 29, 58–61. [Google Scholar]
  20. Del Brio, E.B.; Mora-Valencia, A.; Perote, J. Expected shortfall assessment in commodity (L)ETF portfolios with semi-nonparametric specifications. Eur. J. Financ. 2019, 25, 1746–1764. [Google Scholar] [CrossRef]
  21. Del Brío, E.B.; Mora-Valencia, A.; Perote, J. Risk quantification for commodity ETFs: Backtesting value-at-risk and expected shortfall. Int. Rev. Financ. Anal. 2020. [Google Scholar] [CrossRef]
  22. Tarrant, W. The utility of Basel III rules on excessive violations of internal risk models. J. Risk Model Valid. 2019, 13, 25–37. [Google Scholar] [CrossRef]
  23. Hanly, J. Managing Energy Price Risk using Futures Contracts: A Comparative Analysis. Energy J. 2017, 38, 93–112. [Google Scholar] [CrossRef]
  24. Wu, F. Sectoral contributions to systemic risk in the Chinese stock market. Financ. Res. Lett. 2019, 31, 386–390. [Google Scholar] [CrossRef]
  25. Wen, X.; Nguyen, D.K. Can investors of Chinese energy stocks benefit from diversification into commodity futures? Econ. Model. 2017, 66, 184–200. [Google Scholar] [CrossRef]
  26. Ermolieva, T.; Havlík, P.; Ermoliev, Y.; Mosnier, A.; Obersteiner, M.; Leclère, D.; Khabarav, N.; Valin, H.; Reuter, W. Integrated Management of Land Use Systems under Systemic Risks and Security Targets: A Stochastic Global Biosphere Management Model. J. Agric. Econ. 2016, 67, 584–601. [Google Scholar] [CrossRef] [Green Version]
  27. Vespucci, M.T.; Bertocchi, M.; Pisciella, P.; Zigrino, S. Two-stage stochastic mixed integer optimization models for power generation capacity expansion with risk measures. Optim. Methods Softw. 2016, 31, 305–327. [Google Scholar] [CrossRef]
  28. Spada, M.; Parashiv, F.; Burgherr, P. A comparison of risk measures for accidents in the energy sector and their implications on decision-making strategies. Energy 2018, 154, 277–288. [Google Scholar] [CrossRef] [Green Version]
  29. Tsionas, M.G.; Izzeldin, M. Bayesian CV@R/super-quantile regression. J. Appl. Stat. 2018, 45, 2943–2957. [Google Scholar] [CrossRef] [Green Version]
  30. Abadie, L.M.; Goicochea, N.; Galarraga, I. Carbon risk and optimal retrofitting in cement plants: An application of stochastic modelling, MonteCarlo simulation and Real Options Analysis. J. Clean. Prod. 2017, 142, 3117–3130. [Google Scholar] [CrossRef]
  31. Mandelbrot, B. The variation of certain speculative prices. J. Bus. 1963, 36, 394–419. [Google Scholar] [CrossRef]
  32. Fama, E.F. The behavior of stock-market prices. J. Bus. 1965, 38, 34–105. [Google Scholar] [CrossRef]
  33. Bera, A.K.; Higgins, M.L. ARCH models: Properties, estimation and testing. J. Econ. Surv. 1993, 7, 305–366. [Google Scholar] [CrossRef]
  34. Granger, C.W.J.; Ding, Z. Some properties of absolute returns. An alternative measure of risk. Ann. Deconomie Stat. 1995, 40, 67–91. [Google Scholar] [CrossRef]
  35. Pagan, A. The econometrics of financial markets. J. Empir. Financ. 1996, 3, 15–102. [Google Scholar] [CrossRef]
  36. Cont, R. Empirical properties of asset returns: Stylized facts and statistical issues. Quant. Financ. 2001, 1, 223–236. [Google Scholar] [CrossRef]
  37. Bollerslev, T. Generalized autoregressive conditional heteroskedasticity. J. Econom. 1986, 31, 307–327. [Google Scholar] [CrossRef] [Green Version]
  38. Nelson, D.B. Conditional Heteroskedasticity in asset returns: A new approach. Econometrica 1991, 59, 347. [Google Scholar] [CrossRef]
  39. Ding, Z.; Granger, C.W.; Engle, R.F. A long memory property of stock market returns and a new model. J. Empir. Financ. 1993, 1, 83–106. [Google Scholar] [CrossRef]
  40. McNeil, A.J.; Frey, R.; Embrechts, P. Quantitative Risk Management: Concepts, Techniques and Tools, revised ed.; Princeton University Press: Princeton, NJ, USA, 2015. [Google Scholar]
  41. Fernández, C.; Steel, M.F. On Bayesian modeling of fat tails and skewness. J. Am. Stat. Assoc. 1998, 93, 359–371. [Google Scholar]
  42. Kupiec, P. Techniques for verifying the accuracy of risk measurement models. J. Deriv. 1995, 3, 73–84. [Google Scholar] [CrossRef]
  43. Christoffersen, P. Evaluating interval forecasts. Int. Econ. Rev. 1998, 39, 841–862. [Google Scholar] [CrossRef]
  44. Cardona, E.; Mora-Valencia, A.; Velásquez-Gaviria, D. Testing expected shortfall: An application to emerging market stock indices. Risk Manag. 2019, 21, 153–182. [Google Scholar] [CrossRef]
  45. Tsay, R.S. Analysis of Financial Time Series, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2010. [Google Scholar]
Figure 1. Sensitivity of the test Z2 and ZES to VaR.
Figure 1. Sensitivity of the test Z2 and ZES to VaR.
Energies 13 02805 g001
Figure 2. Critical values tests for skewed-t distribution.
Figure 2. Critical values tests for skewed-t distribution.
Energies 13 02805 g002
Figure 3. Epistar’s standardized returns and standardized residuals of the AR(1)-GARCH(1,1) model under the skew-normal, skew-t and SGED distributions.
Figure 3. Epistar’s standardized returns and standardized residuals of the AR(1)-GARCH(1,1) model under the skew-normal, skew-t and SGED distributions.
Energies 13 02805 g003
Figure 4. The difference in risk measures in Epistar and Exxon Mobil stock returns.
Figure 4. The difference in risk measures in Epistar and Exxon Mobil stock returns.
Energies 13 02805 g004
Figure 5. Estimated VaR-99% and ES-97.5% with the skewed-t distribution.
Figure 5. Estimated VaR-99% and ES-97.5% with the skewed-t distribution.
Energies 13 02805 g005
Table 1. Selection of critical values for Z1, Z2 y ZES tests, and different distributions. T = 2400 and significance level of 5%.
Table 1. Selection of critical values for Z1, Z2 y ZES tests, and different distributions. T = 2400 and significance level of 5%.
Panel A: Normal, Skew-Normal, t-Student and Skew tPanel B: GED and SGED
Test Z1Skew Parameter (λ,ξ) Skew Parameter (ζ)
Degrees of freedom (ν)0.60.80.911.11.2Shape parameter (η)0.60.80.911.11.2
3−0.13−0.13−0.14−0.14−0.14−0.141−0.06−0.06−0.06−0.06−0.05−0.05
5−0.07−0.07−0.07−0.08−0.07−0.081.2−0.05−0.05−0.05−0.05−0.05−0.05
7−0.05−0.06−0.06−0.06−0.06−0.061.4−0.05−0.05−0.04−0.04−0.04−0.04
10−0.05−0.05−0.05−0.05−0.05−0.051.6−0.04−0.04−0.04−0.04−0.04−0.04
Skew−Normal−0.03−0.03−0.03−0.03−0.03−0.031.8−0.04−0.04−0.04−0.03−0.03−0.03
Test Z2Skew Parameter (λ,ξ) Skew Parameter (ζ)
Degrees of freedom (ν)0.60.80.911.11.2Shape parameter (η)0.60.80.911.11.2
3−0.25−0.26−0.26−0.26−0.26−0.261−0.22−0.22−0.22−0.22−0.22−0.22
5−0.22−0.23−0.22−0.23−0.23−0.221.2−0.22−0.22−0.22−0.22−0.21−0.22
7−0.22−0.21−0.22−0.22−0.23−0.221.4−0.22−0.22−0.22−0.22−0.22−0.21
10−0.22−0.22−0.22−0.22−0.22−0.221.6−0.22−0.22−0.21−0.22−0.22−0.22
Skew−Normal−0.21−0.21−0.21−0.21−0.21−0.211.8−0.22−0.22−0.22−0.21−0.22−0.22
Test ZESSkew Parameter (λ,ξ) Skew Parameter (ζ)
Degrees of freedom (ν)0.60.80.911.11.2Shape parameter (η)0.60.80.911.11.2
3−0.50−0.51−0.52−0.52−0.52−0.521−0.34−0.33−0.33−0.32−0.31−0.30
5−0.34−0.35−0.35−0.36−0.36−0.371.2−0.31−0.30−0.29−0.29−0.28−0.27
7−0.29−0.30−0.30−0.30−0.31−0.311.4−0.28−0.28−0.27−0.26−0.25−0.24
10−0.25−0.26−0.27−0.27−0.27−0.281.6−0.25−0.25−0.25−0.24−0.23−0.22
Skew−Normal−0.19−0.20−0.20−0.20−0.21−0.211.8−0.25−0.23−0.23−0.22−0.21−0.21
Z1 and Z2 tests are proposed by Acerbi and Székely [9], while ZES is proposed in Acerbi and Székely [10]. Skew-normal, Skewed-t and SGED distributions collapse to standard normal, Student’s t and GED distribution, respectively, when the skew parameter is equal to 1.
Table 2. Descriptive statistics of renewable and traditional energy stocks.
Table 2. Descriptive statistics of renewable and traditional energy stocks.
Total Period (January 2005–December 2016)
StockTypeMeanMedianStd. Dev.Excess KurtosisSkewnessMinMaxObs
EpistarRenewable−0.00010.00070.030 (47.62%)3.494−0.014−0.1090.0982900
KingspanRenewable0.00030.00080.027 (42.86%)11.997−0.546−0.2540.1792900
BoralexRenewable0.00040.00040.023 (36.51%)13.6770.062−0.1950.2202900
NibeRenewable0.00050.00020.024 (38.10%)12.231−0.119−0.2550.2242900
Exxon MobilTraditional0.00020.00030.016 (25.40%)16.2830.028−0.1500.1592900
ChevronTraditional0.00030.00080.017 (26.99%)15.4950.086−0.1330.1892900
Total S.A.Traditional0.00000.00040.019 (30.16%)10.1300.184−0.1130.1622900
SchlumbergerTraditional0.00030.00020.023 (36.51%)10.806−0.533−0.2030.1392900
Annual volatility is presented in parentheses.
Table 3. 99%-VaR backtesting. Expected violations (EV) = 24.
Table 3. 99%-VaR backtesting. Expected violations (EV) = 24.
AR(1)-GARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
Epistar364.57
(0.03)
4.92
(0.09)
365.25
(0.02)
5.58
(0.06)
240.02
(0.89)
1.45
(0.48)
270.21
(0.65)
1.28
(0.53)
250.04
(0.84)
1.30
(0.52)
280.64
(0.42)
1.56
(0.46)
Kingspan397.10
(0.01)
9.14
(0.01)
4413.51
(0.00)
14.82
(0.00)
200.97
(0.32)
7.70
(0.02)
270.21
(0.65)
4.68
(0.10)
250.04
(0.84)
4.97
(0.08)
240.00
(1.00)
5.23
(0.07)
Boralex4311.20
(0.00)
12.70
(0.00)
4818.79
(0.00)
19.69
(0.00)
230.12
(0.73)
1.65
(0.48)
260.07
(0.80)
1.25
(0.54)
240.00
(1.00)
1.38
(0.50)
270.36
(0.55)
1.39
(0.50)
Nibe311.50
(0.22)
5.00
(0.08)
387.01
(0.01)
9.12
(0.01)
172.73
(0.10)
5.35
(0.07)
270.21
(0.65)
4.68
(0.10)
210.40
(0.53)
6.64
(0.05)
270.36
(0.55)
4.73
(0.09)
Exxon Mobil4918.75
(0.00)
14.48
(0.00)
4514.76
(0.00)
16.48
(0.00)
397.10
(0.01)
9.14
(0.01)
280.42
(0.51)
4.63
(0.10)
408.97
(0.00)
10.33
(0.01)
365.25
(0.02)
6.35
(0.04)
Chevron4614.78
(0.00)
15.92
(0.00)
343.73
(0.05)
4.70
(0.10)
311.50
(0.22)
5.00
(0.08)
172.73
(0.10)
5.35
(0.07)
408.97
(0.00)
10.33
(0.01)
270.36
(0.55)
0.98
(0.61)
Total S.A.4513.54
(0.00)
13.58
(0.00)
4110.03
(0.00)
10.15
(0.01)
270.21
(0.65)
4.68
(0.10)
191.45
(0.23)
4.36
(0.08)
333.05
(0.08)
3.56
(0.17)
280.64
(0.42)
1.56
(0.46)
Schlumberger397.10
(0.01)
7.31
(0.03)
397.96
(0.00)
9.25
(0.01)
230.12
(0.73)
1.65
(0.48)
260.07
(0.80)
1.25
(0.54)
311.89
(0.17)
2.70
(0.26)
301.40
(0.24)
2.16
(0.34)
AR(1)-EGARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
Epistar4513.54
(0.00)
14.81
(0.00)
4211.14
(0.00)
12.70
(0.00)
321.99
(0.16)
5.27
(0.07)
386.20
(0.01)
8.40
(0.02)
333.05
(0.08)
6.04
(0.05)
365.25
(0.02)
7.69
(0.02)
Kingspan4311.20
(0.00)
15.27
(0.00)
4413.51
(0.00)
14.82
(0.00)
220.31
(0.58)
6.30
(0.04)
290.71
(0.40)
8.96
(0.01)
210.40
(0.53)
2.22
(0.33)
270.36
(0.55)
4.73
(0.09)
Boralex5223.11
(0.00)
23.74
(0.00)
5021.68
(0.00)
21.68
(0.00)
270.42
(0.51)
4.63
(0.10)
260.07
(0.80)
1.25
(0.54)
280.64
(0.42)
1.56
(0.46)
311.89
(0.17)
2.54
(0.28)
Nibe301.07
(0.30)
8.93
(0.01)
365.25
(0.02)
7.69
(0.02)
210.59
(0.44)
6.94
(0.03)
270.21
(0.65)
9.28
(0.01)
210.40
(0.53)
2.22
(0.33)
260.16
(0.69)
4.80
(0.09)
Exxon Mobil5832.85
(0.00)
36.45
(0.00)
4211.14
(0.00)
12.64
(0.00)
4614.78
(0.00)
15.94
(0.00)
321.99
(0.16)
2.60
(0.27)
4211.14
(0.00)
12.64
(0.00)
376.10
(0.01)
7.26
(0.03)
Chevron5426.20
(0.00)
28.29
(0.00)
408.97
(0.00)
10.33
(0.01)
4918.75
(0.00)
21.64
(0.00)
220.31
(0.58)
6.30
(0.04)
4717.40
(0.00)
19.28
(0.00)
376.10
(0.01)
7.26
(0.03)
Total S.A.4918.75
(0.88)
18.75
(0.00)
408.97
(0.00)
9.12
(0.01)
386.20
(0.01)
6.45
(0.04)
280.42
(0.51)
4.63
(0.10)
397.96
(0.00)
8.15
(0.02)
333.05
(0.08)
3.97
(0.14)
Schlumberger386.20
(0.01)
8.40
(0.00)
387.01
(0.01)
8.23
(0.01)
290.71
(0.40)
1.58
(0.45)
280.42
(0.51)
4.63
(0.10)
290.99
(0.32)
1.70
(0.43)
270.36
(0.55)
0.98
(0.61)
AR(1)-APARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
Epistar5527.81
(0.00)
29.76
(0.00)
4412.20
(0.00)
12.44
(0.00)
386.20
(0.01)
8.40
(0.02)
4210.10
(0.00)
14.40
(0.00)
333.05
(0.08)
6.04
(0.05)
280.64
(0.42)
1.56
(0.46)
Kingspan4311.20
(0.00)
11.29
(0.00)
419.05
(0.00)
9.31
(0.00)
210.59
(0.44)
2.46
(0.29)
280.42
(0.51)
4.63
(0.10)
270.36
(0.55)
4.73
(0.09)
240.02
(0.89)
0.00
(0.99)
Boralex408.05
(0.00)
8.22
(0.02)
4511.20
(0.00)
13.13
(0.00)
220.31
(0.58)
6.30
(0.04)
240.02
(0.89)
0.00
(0.99)
240.02
(0.89)
0.00
(0.99)
240.02
(0.89)
0.00
(0.99)
Nibe290.71
(0.40)
4.67
(0.10)
321.99
(0.16)
5.51
(0.08)
250.00
(0.95)
1.30
(0.52)
397.10
(0.01)
15.92
(0.00)
220.31
(0.58)
4.63
(0.10)
280.64
(0.42)
4.73
(0.09)
Exxon Mobil5629.45
(0.00)
30.23
(0.00)
4311.20
(0.00)
12.41
(0.00)
4614.78
(0.00)
16.82
(0.00)
280.42
(0.51)
4.63
(0.10)
4211.14
(0.00)
12.64
(0.00)
4413.51
(0.00)
14.82
(0.00)
Chevron5020.16
(0.00)
20.16
(0.00)
4614.78
(0.00)
16.01
(0.00)
386.20
(0.01)
6.45
(0.04)
200.97
(0.32)
3.56
(0.08)
4413.51
(0.00)
14.82
(0.00)
4717.40
(0.00)
19.28
(0.00)
Total S.A.5121.61
(0.00)
22.32
(0.00)
408.97
(0.00)
11.46
(0.01)
332.54
(0.11)
5.52
(0.05)
250.00
(0.95)
1.23
(0.59)
386.20
(0.01)
8.40
(0.02)
376.10
(0.01)
7.26
(0.03)
Schlumberger419.05
(0.00)
9.31
(0.00)
397.10
(0.01)
8.35
(0.00)
250.00
(0.95)
1.23
(0.59)
260.07
(0.80)
1.25
(0.54)
280.42
(0.51)
4.63
(0.10)
270.21
(0.65)
2.60
(0.27)
v stands for violations depending on the distribution (normal, Skew-normal,Student’s t, Skewed-t, GED and SGED). p-values for each test in parentheses. Bold figures represent adequate performance of the model for each test.
Table 4. 97.5%-VaR backtesting. Expected violations (EV) = 60.
Table 4. 97.5%-VaR backtesting. Expected violations (EV) = 60.
AR(1)-GARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
Epistar815.60
(0.02)
8.98
(0.01)
795.62
(0.02)
9.16
(0.01)
660.29
(0.59)
1.05
(0.58)
773.57
(0.06)
7.73
(0.02)
743.12
(0.08)
4.26
(0.12)
785.07
(0.02)
7.09
(0.03)
Kingspan731.98
(0.16)
3.33
(0.19)
848.77
(0.00)
10.07
(0.01)
610.01
(0.92)
1.23
(0.54)
670.44
(0.51)
1.13
(0.57)
610.02
(0.90)
1.15
(0.56)
640.27
(0.60)
2.67
(0.26)
Boralex742.34
(0.13)
7.13
(0.03)
743.12
(0.08)
10.02
(0.01)
620.00
(0.98)
1.12
(0.57)
630.02
(0.87)
1.05
(0.59)
530.87
(0.35)
1.38
(0.50)
550.44
(0.51)
0.83
(0.66)
Nibe610.01
(0.92)
0.16
(0.92)
670.81
(0.37)
1.42
(0.49)
570.39
(0.53)
0.71
(0.70)
784.04
(0.05)
6.22
(0.04)
511.46
(0.23)
2.11
(0.35)
630.15
(0.70)
0.22
(0.89)
Exxon Mobil9112.41
(0.00)
13.17
(0.00)
795.62
(0.02)
6.32
(0.04)
815.60
(0.02)
5.79
(0.06)
640.08
(0.78)
0.15
(0.93)
795.62
(0.02)
5.68
(0.06)
722.32
(0.13)
2.33
(0.31)
Chevron879.39
(0.00)
9.84
(0.01)
701.62
(0.20)
2.31
(0.31)
794.53
(0.03)
4.66
(0.10)
436.54
(0.01)
6.62
(0.04)
859.48
(0.00)
9.89
(0.01)
640.27
(0.60)
0.63
(0.73)
Total S.A.805.05
(0.02)
6.40
(0.04)
785.07
(0.02)
6.34
(0.04)
721.65
(0.20)
2.40
(0.30)
521.68
(0.20)
1.68
(0.43)
722.32
(0.13)
3.14
(0.21)
660.60
(0.44)
1.06
(0.59)
Schlumberger660.29
(0.59)
0.32
(0.85)
670.81
(0.37)
0.82
(0.66)
531.34
(0.25)
1.90
(0.39)
531.34
(0.25)
1.90
(0.39)
610.02
(0.90)
0.14
(0.93)
600.00
(1.00)
0.16
(0.92)
AR(1)-EGARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
Epistar9616.68
(0.00)
20.76
(0.00)
859.48
(0.00)
12.00
(0.00)
8810.11
(0.00)
13.99
(0.00)
9414.90
(0.00)
23.57
(0.00)
743.12
(0.08)
5.74
(0.06)
764.04
(0.04)
6.35
(0.04)
Kingspan8810.11
(0.00)
13.99
(0.00)
8912.54
(0.00)
12.58
(0.00)
742.34
(0.13)
9.53
(0.01)
763.13
(0.08)
9.77
(0.01)
732.71
(0.10)
2.73
(0.26)
753.57
(0.06)
3.62
(0.16)
Boralex879.39
(0.00)
11.80
(0.00)
806.20
(0.01)
6.40
(0.04)
742.34
(0.13)
7.13
(0.03)
731.98
(0.16)
7.00
(0.03)
650.42
(0.52)
0.45
(0.80)
670.81
(0.37)
0.82
(0.66)
Nibe620.00
(0.98)
5.19
(0.07)
691.32
(0.25)
3.02
(0.22)
600.05
(0.82)
1.38
(0.50)
868.70
(0.00)
15.04
(0.00)
521.14
(0.28)
1.72
(0.42)
590.02
(0.90)
0.21
(0.90)
Exxon Mobil9616.68
(0.00)
19.19
(0.00)
743.12
(0.08)
3.16
(0.21)
9112.41
(0.00)
14.26
(0.00)
701.08
(0.30)
4.56
(0.10)
806.20
(0.01)
6.24
(0.04)
691.32
(0.25)
1.32
(0.52)
Chevron9011.62
(0.00)
15.12
(0.00)
660.60
(0.44)
1.06
(0.59)
879.39
(0.00)
11.80
(0.00)
492.91
(0.09)
11.94
(0.00)
8710.96
(0.00)
10.97
(0.00)
650.42
(0.52)
0.83
(0.66)
Total S.A.9818.55
(0.00)
18.86
(0.00)
795.62
(0.02)
5.98
(0.04)
868.70
(0.00)
9.92
(0.01)
660.29
(0.59)
1.05
(0.59)
785.07
(0.02)
5.15
(0.08)
722.32
(0.13)
2.63
(0.27)
Schlumberger752.72
(0.59)
3.87
(0.14)
722.32
(0.13)
3.14
(0.21)
660.29
(0.59)
2.53
(0.28)
630.02
(0.87)
2.73
(0.26)
610.02
(0.90)
0.25
(0.88)
600.00
(1.00)
0.20
(0.91)
AR(1)-APARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
vKupiec TestInd.
Test
Epistar10222.53
(0.00)
27.10
(0.00)
10222.53
(0.00)
24.10
(0.00)
9011.62
(0.00)
11.78
(0.00)
8910.85
(0.00)
11.79
(0.00)
764.04
(0.04)
6.35
(0.04)
806.20
(0.01)
6.40
(0.04)
Kingspan9112.41
(0.00)
12.55
(0.00)
8810.11
(0.00)
11.42
(0.00)
711.35
(0.25)
1.77
(0.41)
763.13
(0.08)
5.61
(0.06)
743.12
(0.08)
5.74
(0.06)
701.08
(0.30)
4.56
(0.10)
Boralex815.60
(0.02)
5.79
(0.06)
805.05
(0.02)
6.51
(0.02)
680.62
(0.43)
0.63
(0.73)
690.84
(0.36)
0.84
(0.66)
650.42
(0.52)
0.83
(0.66)
691.32
(0.25)
1.32
(0.52)
Nibe630.02
(0.87)
1.05
(0.59)
630.02
(0.87)
2.73
(0.26)
620.00
(0.98)
0.24
(0.89)
9213.22
(0.00)
18.12
(0.00)
521.14
(0.28)
1.72
(0.42)
570.39
(0.53)
0.71
(0.70)
Exxon Mobil10121.50
(0.00)
21.87
(0.00)
806.20
(0.01)
6.40
(0.04)
9213.22
(0.00)
13.28
(0.00)
711.35
(0.25)
1.35
(0.51)
806.20
(0.01)
6.40
(0.04)
731.98
(0.16)
1.05
(0.59)
Chevron9414.90
(0.00)
15.01
(0.00)
9011.62
(0.00)
10.08
(0.00)
794.53
(0.03)
5.79
(0.06)
436.54
(0.01)
6.62
(0.04)
8810.11
(0.00)
11.42
(0.00)
711.35
(0.25)
1.77
(0.41)
Total S.A.9011.62
(0.01)
11.78
(0.00)
9011.62
(0.00)
13.22
(0.00)
826.17
(0.01)
6.20
(0.05)
640.08
(0.78)
4.80
(0.09)
806.20
(0.01)
6.40
(0.04)
743.12
(0.08)
3.16
(0.21)
Schlumberger763.13
(0.82)
4.16
(0.12)
731.98
(0.16)
1.05
(0.59)
600.05
(0.82)
0.22
(0.90)
560.57
(0.45)
0.64
(0.73)
620.00
(0.98)
0.24
(0.89)
690.84
(0.36)
0.84
(0.66)
10222.53
(0.00)
27.10
(0.00)
10222.53
(0.00)
24.10
(0.00)
9011.62
(0.00)
11.78
(0.00)
8910.85
(0.00)
11.79
(0.00)
630.02
(0.87)
2.73
(0.26)
600.05
(0.82)
1.38
(0.50)
v stands for violations depending on the distribution (normal, Skew-normal,Student’s t, Skewed-t, GED and SGED). p-values for each test in parentheses. Bold figures represent the adequate performance of the model for each test.
Table 5. 97.5%-ES backtesting. Values for Z1, Z2 y ZES statistics.
Table 5. 97.5%-ES backtesting. Values for Z1, Z2 y ZES statistics.
AR(1)-GARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
Z1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZES
Epistar−0.09−0.45−0.32−0.10−0.44−0.340.02−0.17−0.250.00−0.37−0.33−0.02−0.20−0.20−0.02−0.32−0.20
Kingspan−0.10−0.52−0.36−0.12−0.57−0.390.02−0.05−0.240.01−0.22−0.24−0.03−0.06−0.10−0.03−0.11−0.16
Boralex−0.21−0.35−0.41−0.19−0.46−0.430.030.08−0.260.020.01−0.23−0.030.05−0.22−0.020.00−0.19
Nibe−0.18−0.27−0.37−0.21−0.35−0.41−0.040.03−0.28−0.01−0.49−0.42−0.130.04−0.28−0.12−0.18−0.34
Exxon Mobil−0.16−0.72−0.48−0.14−0.50−0.40−0.05−0.48−0.39−0.05−0.09−0.29−0.10−0.45−0.39−0.09−0.20−0.34
Chevron−0.11−0.70−0.41−0.10−0.28−0.32−0.05−0.46−0.37−0.010.25−0.16−0.09−0.54−0.38−0.08−0.15−0.28
Total S.A.−0.09−0.57−0.36−0.09−0.41−0.32−0.02−0.31−0.300.000.05−0.20−0.02−0.21−0.22−0.03−0.15−0.16
Schlumberger−0.15−0.36−0.36−0.15−0.29−0.35−0.05−0.05−0.27−0.06−0.06−0.25−0.09−0.11−0.18−0.09−0.09−0.18
AR(1)-EGARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
Z1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZES
Epistar−0.08−0.52−0.34−0.09−0.54−0.350.00−0.28−0.27−0.05−0.41−0.36−0.02−0.21−0.22−0.01−0.33−0.19
Kingspan−0.08−0.60−0.36−0.10−0.63−0.380.03−0.20−0.260.01−0.21−0.26−0.00−0.190.21−0.01−0.26−0.17
Boralex−0.14−0.52−0.40−0.16−0.54−0.420.03−0.04−0.280.02−0.05−0.27−0.03−0.13−0.19−0.03−0.13−0.10
Nibe−0.19−0.23−0.37−0.20−0.38−0.41−0.050.02−0.280.00−0.42−0.38−0.120.03−0.16−0.14−0.12−0.18
Exxon Mobil−0.18−0.70−0.49−0.17−0.44−0.41−0.10−0.50−0.43−0.07−0.11−0.26−0.13−0.51−0.43−0.12−0.19−0.40
Chevron−0.14−0.63−0.43−0.12−0.23−0.32−0.09−0.62−0.42−0.070.22−0.17−0.11−0.61−0.41−0.10−0.19−0.36
Total S.A.−0.12−0.58−0.39−0.10−0.45−0.34−0.01−0.42−0.35−0.01−0.03−0.22−0.02−0.40−0.22−0.03−0.21−0.20
−0.15−0.30−0.35−0.12−0.35−0.34−0.07−0.08−0.28−0.07−0.01−0.27−0.09−0.11−0.19−0.09−0.09−0.20
AR(1)-APARCH(1,1)
NormalSkew-Normalt-StudentSkew-tGEDSGED
Z1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZESZ1Z2ZES
Epistar−0.11−0.63−0.29−0.08−0.52−0.320.04−0.21−0.27−0.02−0.41−0.360.03−0.21−0.21−0.03−0.31−0.17
Kingspan−0.09−0.49−0.35−0.11−0.58−0.360.01−0.20−0.260.01−0.19−0.22−0.02−0.180.15−0.00−0.28−0.21
Boralex−0.12−0.72−0.44−0.12−0.61−0.400.05−0.04−0.250.02−0.05−0.24−0.02−0.11−0.16−0.03−0.12−0.14
Nibe−0.21−0.35−0.36−0.23−0.42−0.37−0.030.02−0.220.00−0.42−0.38−0.120.03−0.22−0.12−0.12−0.17
Exxon Mobil−0.11−0.59−0.47−0.26−0.51−0.47−0.10−0.50−0.43−0.09−0.11−0.21−0.15−0.55−0.39−0.11−0.17−0.41
Chevron−0.09−0.45−0.45−0.14−0.28−0.36−0.03−0.62−0.42−0.020.22−0.23−0.10−0.60−0.41−0.10−0.18−0.36
Total S.A.−0.11−0.69−0.42−0.11−0.45−0.38−0.02−0.42−0.35−0.01−0.03−0.25−0.03−0.41−0.21−0.01−0.21−0.22
Schlumberger−0.18−0.22−0.35−0.17−0.35−0.33−0.12−0.08−0.25−0.08−0.01−0.19−0.09−0.11−0.17−0.11−0.11−0.16
The critical value for Z1 test is -0.03 for normal, skew-normal, GED and SGED, -0.05 for Student’s t and skewed-t. The critical value for Z2 test is -0.21 for normal, skew-normal, GED and SGED, -0.22 for Student’s t and skewed-t. The critical value for ZES test is -0.20 for normal and skew-normal, -0.27 for Student’s t and skewed-t, and -0.22 for GED and SGED. Bold figures represent the adequate performance of the model for each test.
Table 6. Capital requirements percentage for energy stocks under skewed-t innovations.
Table 6. Capital requirements percentage for energy stocks under skewed-t innovations.
Total Period (January 2005–December 2016)
ModelAR(1)-GARCH(1,1)AR(1)-EGARCH(1,1)AR(1)-APARCH(1,1)
StockVaR99%ES97.5%diffVaR99%ES97.5%diffVaR99%ES97.5%diff
Renewable
Epistar−4.52166−4.627060.10540−4.50284−4.604160.10131−4.58420−4.687960.10376
Kingspan−5.68544−5.939380.25394−5.52422−5.757600.23338−5.29004−5.512060.22202
Boralex−3.74001−3.998520.25851−3.48382−3.709960.22614−3.85014−4.104330.25420
Nibe−4.58810−4.791420.20332−4.91927−5.130770.21150−4.71179−4.913500.20171
Traditional
Exxon Mobil−2.75996−2.843870.08391−2.43863−2.506950.06832−2.44536−2.513130.06777
Chevron−3.03334−3.094450.06111−3.06976−3.111850.04209−3.15134−3.214260.06291
Total S.A.−3.87745−3.980490.10305−3.93551−4.031670.09616−3.83333−3.924980.09166
Schlumberger−3.23484−3.315880.08104−3.08337−3.158760.07540−3.06615−3.139850.07370
The table presents the differences in capital requirements from replacing VaR99% with ES97.5%, computed as diff = ES97.5%VaR99%.

Share and Cite

MDPI and ACS Style

Velásquez-Gaviria, D.; Mora-Valencia, A.; Perote, J. A Comparison of the Risk Quantification in Traditional and Renewable Energy Markets. Energies 2020, 13, 2805. https://doi.org/10.3390/en13112805

AMA Style

Velásquez-Gaviria D, Mora-Valencia A, Perote J. A Comparison of the Risk Quantification in Traditional and Renewable Energy Markets. Energies. 2020; 13(11):2805. https://doi.org/10.3390/en13112805

Chicago/Turabian Style

Velásquez-Gaviria, Daniel, Andrés Mora-Valencia, and Javier Perote. 2020. "A Comparison of the Risk Quantification in Traditional and Renewable Energy Markets" Energies 13, no. 11: 2805. https://doi.org/10.3390/en13112805

APA Style

Velásquez-Gaviria, D., Mora-Valencia, A., & Perote, J. (2020). A Comparison of the Risk Quantification in Traditional and Renewable Energy Markets. Energies, 13(11), 2805. https://doi.org/10.3390/en13112805

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop