Next Article in Journal
Diffusion-Based Causal Representation Learning
Next Article in Special Issue
Optimizing Prognostic Predictions in Liver Cancer with Machine Learning and Survival Analysis
Previous Article in Journal
The Statistics of q-Statistics
Previous Article in Special Issue
Detection of Interaction Effects in a Nonparametric Concurrent Regression Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tail Risk Dynamics under Price-Limited Constraint: A Censored Autoregressive Conditional Fréchet Model

Department of Statistics and Finance, School of Management, University of Science and Technology of China, Hefei 230026, China
*
Authors to whom correspondence should be addressed.
Entropy 2024, 26(7), 555; https://doi.org/10.3390/e26070555
Submission received: 5 May 2024 / Revised: 24 June 2024 / Accepted: 26 June 2024 / Published: 28 June 2024

Abstract

:
This paper proposes a novel censored autoregressive conditional Fréchet (CAcF) model with a flexible evolution scheme for the time-varying parameters, which allows deciphering tail risk dynamics constrained by price limits from the viewpoints of different risk preferences. The proposed model can well accommodate many important empirical characteristics of financial data, such as heavy-tailedness, volatility clustering, extreme event clustering, and price limits. We then investigate tail risk dynamics via the CAcF model in the price-limited stock markets, taking entropic value at risk (EVaR) as a risk measurement. Our findings suggest that tail risk will be seriously underestimated in price-limited stock markets when the censored property of limit prices is ignored. Additionally, the evidence from the Chinese Taiwan stock market shows that widening price limits would lead to a decrease in the incidence of extreme events (hitting limit-down) but a significant increase in tail risk. Moreover, we find that investors with different risk preferences may make opposing decisions about an extreme event. In summary, the empirical results reveal the effectiveness of our model in interpreting and predicting time-varying tail behaviors in price-limited stock markets, providing a new tool for financial risk management.

1. Introduction

Tail risks, which can be quantified by risk measurements such as quantile, expectile, and entropic value at risk [1], highlight the potential for serious losses that could affect investors, financial institutions, and the overall stability of financial markets, making their measurement critical in the financial fields. Understanding and managing tail risk helps mitigate adverse consequences and maintain financial resilience. A voluminous literature provides econometric tools to measure tail risk ([2,3,4,5] and references therein). However, the presence of price-limit trading policies in certain markets complicates the accurate depiction of tail risk, as these constraints may distort the manifestation of tail risk. Price limits, widely endorsed across global stock and futures exchanges, serve as a safeguard procedure for investors and a deterrent against market manipulation [6,7,8]. By imposing restrictions on daily price fluctuations, these policies challenge the applicability of existing tail risk measures, typically designed for unrestricted markets. Consequently, this scenario urgently calls for tailored approaches in modeling tail risk within markets subject to price limits.
The challenges posed by the price-limit policy are twofold. From an econometric point of view, the dispute on the pros and cons of the price-limit trading policy lasts for ages. Some argue that the price-limit policy may lead to ineffectiveness or even destructive market behavior [9,10,11], while others believe that a price-limited policy can reduce market manipulation risk and improve market efficiency [12,13,14]. The disorderly effectiveness of this policy makes it challenging to quantify its impact on tail risk. From a statistical point of view, this policy results in the censoring of observations, and ignoring such censoring could cause substantial bias and size distortion in measuring tail risk, even if the censored probability is tiny. The potential biases also occur in modeling the volatility of returns with price limits. To handle this case, Wei [15] develops a censored-GARCH model to recognize the unobservable feature of price-limited data and Hsieh and Yang [16] subsequently propose a censored stochastic volatility approach based on the censored-GARCH model to further improve the computational efficiency. However, both censored approaches rely on the algebraic relationship between the observed and latent returns (e.g., Equations (2)–(4) in [15]), which breaks down when dealing with the market’s aggregated information, such as tail risk, since it is jointly determined by multiple stocks. In addition, we present a simple analysis of the SSE50 (the SSE50 is a value-weighted price index that represents the performance of the top 50 firms listed on the Shanghai Stock Exchange, selected based on their market capitalization, liquidity, and other criteria) in Figure 1 to further illustrate the damage in measuring the tail risk caused by ignoring the censoring nature. The histogram in the left panel reveals a notable probabilistic stacking of observations triggered by the price limit, whereas the right panel demonstrates that the uncensored fitting method leads to an underestimation of tail risk (quantile).
In the domain of tail risk assessment, Extreme Value Theory (EVT; [17]) stands out as a potent instrument. This theory encompasses two principal methods: one approach is fitting the maximum observations using the generalized extreme value distribution (GEV) and is commonly referred to as the Maxima-GEV or Block Maxima (BM) method (e.g., [18]); another approach involves the Peak-over-Threshold (POT) method [19], which employs the Generalized Pareto Distribution (GPD) to approximate the conditional behavior of random variables that exceed specific high thresholds. Several studies have extensively explored the implications of EVT in assessing tail risk within the context of price limits. Oh et al. [20] assume that the conditional tail distribution of extreme returns obeys a power law and obtain an inferred estimation of tail risk under price limits. Subsequently, Ji et al. [21] introduce a general framework of the self-exciting point process with the truncated generalized Pareto distribution to measure the extreme risks in the price-limited stock markets. Nevertheless, the former ignore the censored nature of extreme returns when estimating the tail index, a critical aspect of accurate tail risk estimation. The latter use a truncated distribution rather than a censored structure to self-adapt to the price limits, which prevents this approach from modeling the latent return and results in their risk measurements continuing to be constrained by price limits, lacking sensitivity to extreme risk events. In addition, both methods ignore a dynamic treatment for the tail index, which has been demonstrated to be necessary by Massacci [22], Zhao et al. [23], and Shen et al. [24], etc. These studies have revealed substantial evidence that the tail risk in financial markets without price limits exhibits significant dynamics over time. Intuitively, these dynamic features of tail risk would also be present in price-limited markets.
For a deeper understanding of tail risk dynamics in price-limited markets, this paper focuses on modeling the time-varying tail features when observations beyond some threshold are censored. We propose a novel censored autoregressive conditional Fréchet model, which accommodates the censoring, heavy-tailed, volatility clustering, and extreme event clustering nature of financial data. The CAcF model incorporates a flexible observation-driven time evolution scheme of the parameters σ t (volatility index) and α t (tail index) of a Fréchet (Type-II GEV) distribution, and the censoring feature into the modeling, allowing for a more explicit exploration of the time-varying tail behavior in price-limited equity markets. Moreover, we employ three typical observation-driven functions to decompose the tail risk from varying risk preference perspectives. (Risk preference is a pivotal factor in economic behavior, directly influencing the choice and behavior of investors in risk investment decisions [25,26,27]).
To empirically illustrate our findings, we utilize stock data from companies included in the SSE50, CSI300 (the CSI300 is a broader index that encompasses the top 300 firms listed on the Shanghai Stock Exchange and Shenzhen Stock Exchange), and TW50 (the TW50 is a market capitalization-weighted stock index developed by the Taiwan Stock Exchange in cooperation with the Financial Times and Stock Exchange (FTSE), which comprises the 50 companies with the largest market capitalization listed on the Taiwan Stock Exchange). In terms of our proposed model, we offer a maximum likelihood estimation (MLE) procedure for model estimation. To quantify tail risk, we adopt the entropic value at risk, which incorporates self-information via entropy and allows for a more flexible and robust representation of risk. We have also derived closed-form expressions for entropic value at risk and censored probability within this framework, providing a convenient approach for out-of-sample prediction. The empirical estimation results demonstrate that the CAcF model can effectively monitor time-varying behaviors of tail risk and provide satisfactory forecasting performance. This suggests its potential value in warning against financial tail risks. In addition, the tail risk of price-limited stock markets is significantly underestimated when censoring is not taken into account. Moreover, our analysis results show that the CAcF-type models with different risk preferences yield varied interpretations of risk. Specifically, risk-preferred investors perceive that hitting limit-down will reduce the potential risk at the next moment, whereas risk-averse investors interpret it conversely. These findings align with the principles of investment psychology and market dynamics. Finally, we study the impact of widening price limits by comparing the performance of the tail risks of TW50 over periods with different price limits. The evidence shows that widening price limits would lead to a decrease in the incidence of extreme events (hitting limit-down), but a significant increase in tail risk.
This paper provides a twofold contribution to the growing literature on tail risk measurement in financial markets. From a statistical modeling point of view, we propose a dynamic tail risk model for price-limited financial markets. The CAcF model incorporates a flexible observation-driven time evolution scheme for the key parameters and accommodates many important empirical characteristics of financial data. We demonstrate that the CAcF model can be derived from a general factor model, which ensures that the dynamic model is theoretically feasible. From an econometric point of view, the CAcF model offers a new perspective to illustrate tail risk dynamics when price limits exist in financial markets. Real applications show that tail risk is seriously underestimated when the price-limited constraint is ignored. Moreover, this study provides valuable insights for policymakers to develop more effective price-limited policies from a risk management perspective.
The rest of this paper is structured as follows. Section 2 introduces the framework of the CAcF model and derives a maximum likelihood estimation procedure. Section 3 presents the empirical results and analysis. Finally, we conclude the paper in Section 4.

2. Methodology

2.1. Model Specification

Let { Q t } denote the maximal time series representing the cross-sectional maximum of the negative daily return X i t of the stock prices of N companies occurring on day t, i.e., Q t = max 1 i N X i t . Essentially, Q t can offer insight into the potential tail risk of a specific stock market. However, in a price-limited stock market, the econometrician observes that { Q t } is censored by the following mechanism:
Q t = D t M + 1 D t Q t * ,
where D t = I Q t * > M , M is a constant censoring threshold determined by the price limit, Q t * is the latent return but unobserved, and I [ · ] denotes the indicator function.
Remark 1. 
Constrained by the price-limit policy, the observed maximal negative return Q t should be interval censored. However, it is almost impossible for all the stocks in a specific stock market to simultaneously hit the limit-up, which suggests that left-censoring is almost infeasible for Q t . This allows us to relax the assumption on Q t from interval-censored to right-censored as defined in (1).
According to the Fisher–Tippett–Gnedenko theorem [28,29], Q t can typically be modeled, under an independent and identically distributed ( i . i . d . ) assumption, with a truncated GEV distribution. Nevertheless, the classical approach neglects the time dependence among { Q t } and other characteristics, such as censored and heteroskedasticity. To address time-varying behaviors and censored nature of Q t , building upon the dynamic GEV framework (autoregressive conditional Fréchet model, AcF) proposed by Zhao et al. [23], we now introduce the following censored autoregressive conditional Fréchet (CAcF) model:
Q t * = μ + σ t Y t 1 / α t , Q t = D t M + 1 D t Q t * , log σ t = β 0 + β 1 log σ t 1 + β 2 + β 2 * D t 1 G 1 Q t 1 , log α t = γ 0 + γ 1 log α t 1 γ 2 + γ 2 * D t 1 G 2 Q t 1 ,
where 0 β 1 γ 1 < 1 , β 2 * , γ 2 * R , and β 2 , γ 2 > 0 . Further, { Y t } is a sequence of i . i . d . unit Fréchet random variables and ( μ , σ t , α t ) F t 1 = σ ( Q t 1 , D t 1 , Q t 2 , D t 2 , ) represents the location parameter, scale parameter, and shape parameter. The two terms G 1 ( · ) and G 2 ( · ) are the observation-driven functions for { log σ t } and { log α t } , and both are assumed to be monotonically increasing functions of Q t 1 . Combined with the autoregressive scheme for { log σ t } and { log α t } , this setting for the observation-driven functions ensures that the distribution of Q t has a larger scale (larger σ t ) and a heavier tail (smaller α t ) than that of Q t 1 . In other words, this setting offers a joint modeling of both volatility clustering for { σ t } process and extreme event clustering for { α t } process, resulting in a larger tail risk of Q t when an extreme event occurs at time t 1 (large Q t 1 ).
Remark 2. 
In contrast to the conventional AcF model, the threshold M imposes an upper limit on the observable maximum sequence { Q t } . The CAcF model degenerates to the AcF model when M = . Obviously, the CAcF model has a wider range of applications, including financial time series with price limits.
Remark 3. 
Similar to the TGARCH model [30], we insert a threshold structure, β 2 * D t 1 ( γ 2 * D t 1 ), to compensate for the shock of hitting limit-down at time t 1 on the tail risk of Q t . The compensation coefficient β 2 * (as well as γ 2 * ) represents the degree of risk compensation associated with hitting limit-down. A positive value for β 2 * (as well as γ 2 * ) suggests that hitting limit-down at the previous moment intensifies the risk clustering effect. In contrast, a negative value suggests a weakening of the risk clustering effect.
Notice that the boundedness of Q t allows us to relax the boundedness assumption for G 1 ( · ) and G 2 ( · ) in Zhao et al. [23], thus providing a wider range of choices for G 1 ( · ) and G 2 ( · ) compared to AcF model. Drawing on the relationships between risk preference and the convexity–concavity of the utility function in Expected Utility Theory (in this theory, the utility function whose expected value is maximized is convex for a risk-loving agent, linear for a risk-neutral agent, and concave for a risk-averse agent; see more details in Chapter 6 of Mas-Colell et al. [31]), we consider three types of observation-driven functions in the following:
  • Linear:  G 1 z = G 2 z = z .
  • Square:  G 1 z = G 2 z = sgn ( z ) z 2 , where sgn ( · ) is a sign function.
  • Exponential:  G 1 z = exp β 3 z and G 2 z = exp γ 3 z , where β 3 , γ 3 > 0 .
The CAcF models corresponding to the three different functions mentioned above are denoted as CAcF-L (Linear), CAcF-S (Square), and CAcF-E (Exponential), respectively. The main diversity among these three models is that the convexity–concavity of their observation-driven functions differs when faced with loss ( Q t 1 > 0 ). This diversity allows the three models to analyze the impact of historical information on tail risk from different risk preference perspectives. First, the CAcF-L model implies that the impact of the increase in loss is marginally constant. Thus, the CAcF-L model is suitable for those investors or regulators with a risk-neutral preference. Second, the CAcF-S model assumes that an incremental loss can result in a rapid escalation of unease or pessimism about future risk, displaying a high sensitivity to losses. This characteristic is typical of risk-averse individuals. Finally, the CAcF-E model supposes that the impact of the increase in loss is marginally decreasing, which represents the attitude of risk lovers.
Together with the threshold structure in (2), these three types of observation-driven functions allow us to analyze the decisions of investors with different risk perspectives when hitting limit-down. For example, if the compensation coefficients are negative ( β 2 * , γ 2 * < 0 ) in the CAcF-E model, this scenario means that a risk lover may regard hitting limit-down as a signal to buy stocks. A potential explanation is that risk lovers believe that the magnet effect of hitting limit-down leads to an arbitrage opportunity in the stock market. This interpretability of the CAcF model makes it more competitive in analyzing the tail risks of price-limited financial markets.
Remark 4. 
Note that the relationships between risk preference and the convexity–concavity of the observation-driven functions in CAcF models are the reverse of those in Expected Utility Theory. This reversal occurs because the CAcF models focus on losses while Expected Utility Theory concentrates on returns.
Remark 5. 
Although the three CAcF models employ different observation-driven functions, their recovered { σ t } and { α t } processes are highly similar. This is because the recovered { σ t } and { α t } processes by different CAcF models are constrained by the same autoregression structure and the same log-likelihood function. Moreover, we verify this in Section 3. Consequently, the three CAcF models have a consensus on the estimation of tail risk.

2.2. CAcF Model under a Factor Model Framework

In this subsection, we show that the CAcF model can be derived from the limiting form of maximum Q t under a latent general factor model framework. Assume { x i t } i = 1 N is subject to a latent general factor model,
X i t = M , X i t * M , X i t * , M < X i t * < M , M , X i t * M ,
X i t * = f i Z 1 t , Z 2 t , , Z d t + σ i t ε i t ,
where X i t i = 1 p and X i t * i = 1 p are observed and unobserved latent negative returns of company i at time t, respectively, M is a constant censoring threshold determined by the price limit, Z 1 t , Z 2 t , , Z d t consist of observed and unobserved factors, ε i t i = 1 N are i . i . d . noises that are independent of the factors Z i t i = 1 d and σ i t i = 1 N F t 1 are the conditional volatilities of X i t * i = 1 N . The function f i : R d R is a Borel function.
Factor model has been widely used for modeling asset returns [32,33,34]. The general factor model (4) can handle common properties among X i t , such as heterogeneous volatilities and cross-sectional dependence. Together with (3), such a factor structure can well model asset returns in the price-limited market. To incorporate another important characteristic, heavy-tailedness, of many financial time series, we assume that the random noises ε i t i = 1 N are random variables in the Domain of Attraction of Fréchet distribution [35]. Distributions in the Domain of Attraction of Fréchet distribution include a broad class of random variables such as Cauchy, Lévy, Pareto, and t distributions. A specific definition is described in Appendix B.
Further, we consider the following Assumption 1.
Assumption 1. 
(a) 
Under a dynamic model, the tail index α i t of ε i t evolves through time according to certain dynamics and α t F t 1 . In addition, ε i t i = 1 N are i . i . d . random variables in the Domain of Attraction of Fréchet distribution with a tail index α t .
(b) 
sup 1 N < sup 1 i N f i Z 1 t , Z 2 t , , Z d t < , a . s .
(c) 
lim N i = 1 N σ i t α t = and lim N sup 1 i N ( j = 1 N σ j t α t ) 1 σ i t α t = 0 .
Assumption 1(a) commonly addresses the heavy-tailed nature of financial time series within factor models. Assumption 1(b) is a mild assumption about the boundedness of f i . Assumption 1(c) means the magnitudes of conditional volatility σ i t are comparable to each other and there is no single X i t that dominates the total volatility.
Proposition 1. 
Denote Q t = max 1 i N X i t , a N t = 0 , and b N t = i = 1 N σ i t α t 1 / α t . Under Assumptions 1 and given F t 1 , we have, as N ,
Q t a N t b N t d max Ψ α t , ( M a N t ) / b N t ,
where Ψ α t is a Fréchet type random variable with tail index α t and the cumulative distribution function of Ψ α t is exp x α t .
The proof of Proposition 1 is presented in Appendix B. Proposition 1 shows that, under the framework of the general factor model and some mild conditions, the conditional distribution of Q t can be well approximated by a right-censored Fréchet distribution, which provides the rationality of the CAcF model.

2.3. Parameter Estimation

We denote all the parameters in the model by θ = μ , β 0 , β 1 , β 2 , β 2 * , β 3 , γ 0 , γ 1 , γ 2 , γ 2 * , γ 3 and the constraints Θ = θ μ , β 0 , γ 0 , β 2 * , γ 2 * R , 0 β 1 , γ 1 < 1 , β 2 , β 3 , γ 2 , γ 3 > 0 . Using model setting and conditional independence, we can write the log-likelihood function with censored observations { Q t } t = 1 n as
L n θ = 1 n t = 1 n D t log P Q t * > M F t 1 + 1 D t log f Q t * Q t F t 1 = 1 n t = 1 n [ D t log 1 exp ( σ t α t ( M μ ) α t ) + 1 D t [ log α t + α t log σ t α t + 1 log Q t μ σ t α t ( Q t μ ) α t ] ] ,
where f Q t * ( · F t 1 ) is the probability density function of Q t * conditional on the past information F t 1 , and { σ t , α t } t = 1 n can be obtained recursively through (2) with an initial value ( σ 1 , α 1 ) . Subsequently, the maximum likelihood estimator θ ^ can be obtained by iterative search.
Remark 6. 
In practice, we use the estimated ( σ ^ , α ^ ) from the static censored Fréchet distribution with the threshold M as the initial value for ( σ 1 , α 1 ) . Since 0 β 1 , γ 1 < 1 , the influence of ( σ 1 , α 1 ) on recovered ( σ ^ t , α ^ t ) decays exponentially as t increases. Hence, the impact of ( σ 1 , α 1 ) on parameter estimation could be ignored with a sufficiently large sample size.

2.4. Prediction of Maximum Negative Return, Entropic Value at Risk, and Censored Probability

In this subsection, we derive several closed-form expressions for tail risk prediction. First, we give the one-step forward recursive prediction of the maximum negative return Q t as follows:
Q ^ t = E Q t F t 1 = M + σ t Γ 1 1 / α t , σ t α t ( M μ ) α t ( M μ ) exp σ t α t ( M μ ) α t ,
where the upper incomplete gamma function Γ ( a , b ) = b x a 1 exp ( x ) d x .
Second, we employ the entropic value at risk (EVaR) based on Rényi entropy, which is proposed by Pichler and Schlotter [1], as a statistical metric to quantify the level of tail risk associated with an investment portfolio, security, or company. Specifically, we adopt the following definition:
Definition 1 
(Pichler and Schlotter [1]). The entropic value at risk EVaR τ p ( X ) of order p R at confidence level τ [ 0 , 1 ) and X L p based on Rényi entropy is
EVaR τ p ( X ) : = sup E X Z : Z 0 , E Z = 1 and H q ( Z ) log 1 1 τ ,
where 1 p + 1 q = 1 and H q ( Z ) is the Rényi entropy of order q R of a random variable Z, i.e.,
H q ( Z ) : = log P ( Z > 0 ) if q = 0 , Z log Z if q = 1 , log Z if q = , 1 q 1 log E Z q else .
Since the CAcF-type models are based on the Fréchet distribution, whose p-th moment exists if its tail index α > p (i.e., a Fréchet random variable X L p if and only if its tail index α > p ), and the estimated tail index α t ranges from 2 to 12 in our empirical analysis, we mainly focus on the entropic value at risk of order p = 1 at confidence level τ ( EVaR τ 1 ) in the following content. This has another potential benefit in that EVaR τ 1 is a coherent risk measure (see more details in [36]). According to Theorem 12 of Pichler and Schlotter [1], EVaR τ 1 has the dual representation
EVaR τ 1 ( X ) = inf c R c + 1 1 τ · E max ( X c , 0 ) .
Note that EVaR τ 1 is equal to the conditional value at risk (CVaR, [36,37]), also known as the expected shortfall (ES). Then, the one-step forward predictable conditional EVaR τ 1 of Q t * and Q t are
EVaR τ , Q t * F t 1 1 = μ + σ t 1 τ γ 1 1 / α t , log τ ,
where the lower incomplete gamma function γ ( a , b ) = 0 b x a 1 exp ( x ) d x . We call EVaR 1 τ , Q t * F t 1 1  uncensored EVaR in the following content. Both the conditional EVaRs of Q t * and Q t reflect the tail risk of a specific price-limited market. However, the conditional EVaR of Q t , limited by its boundedness, cannot reflect the severity of the risks suffered by a specific price-limited market in a financial crisis. Instead, the uncensored EVaR is a more suitable risk measurement for price-limited markets.
Finally, an alternative indicator that can be used to monitor the tail risk of a price-limited market is the censored probability, whose one-step forward prediction is given as follows:
P ^ t = P Q t * > M F t 1 = 1 exp σ t α t ( M μ ) α t .

3. Empirical Analysis

3.1. Chinese Mainland Stock Market

The implementation of the price-limit policy in the stock markets of mainland China on 16 December 1996 necessitated a ± 10 % limit on the fluctuation of stock prices relative to their previous closing price. To elucidate the dynamics of tail risk within the Chinese equity market, this study applies the CAcF models to examine the cross-sectional maximum of negative daily returns among stocks listed on two pivotal Chinese stock indices: the SSE50 and the CSI300. These indices are not only paramount but also the most frequently cited barometers of Chinese financial market health and are considered reflective of the broader sentiment pervading the Chinese equity landscape. The empirical inquiry spans from 4 January 2005 to 30 December 2022, aggregating 4374 data points. For model calibration, the study delineates the period between 4 January 2005 and 30 December 2020, deploying the subsequent data from 4 January 2021 to 30 December 2022 for out-of-sample predictions.
Note that the price limit was relaxed to ± 20 % on the Science and Technology Innovation Board and Second Board starting 13 June 2019 and 24 August 2020, respectively. For the homogeneity of data, this analysis excludes stocks in SSE50 and CSI300 that are not part of the Main Board, thereby maintaining adherence to the ± 10 % price limit. It is imperative to note that the proportion and temporal span of stocks subjected to the ± 20 % threshold within the SSE50 and CSI300 only account for a small fraction. Consequently, the exclusion of these data subsets does not materially impede the assessment of tail risk associated with the indices.

3.1.1. Fitting Results for SSE50 and CSI300

We employed the CAcF models to fit the in-sample data and conducted a comparative analysis with the AcF model, which ignores the censored nature. The results of parameter estimation are presented in Table 1. This comparison yielded several noteworthy insights. Specifically, analysis of the estimated location parameter, μ , reveals that, for both the SSE50 and CSI300 datasets, all three variants of the CAcF model exhibit lower absolute values in the estimated location parameters compared to the AcF model, with correspondingly reduced standard deviations. Furthermore, it is observed that the location parameter μ estimated by the AcF model significantly deviates from the observed value range. This evidence suggests that neglecting the censored characteristics of the data introduces substantial bias and size distortion to the location parameter estimation, potentially leading to inaccurate tail risk assessment when censored observations are modeled using the AcF approach.
Subsequently, we focus on analyzing the estimated autoregressive coefficients β 1 , γ 1 and the estimated compensation coefficients β 2 * , γ 2 * for the volatility index process { σ t } and the tail index process { α t } . The remaining coefficients only determine the movement range of { σ t } and { α t } without other particular meanings. For the SSE50 dataset, all the estimated autoregressive coefficients β 1 and γ 1 range from 0.83 to 0.88 , indicating a strong persistence in both the { σ t } process and the { α t } process while, for the CSI300 dataset, all the estimated autoregressive coefficients are uniformly smaller than those for the SSE50 dataset, which suggests that there is a stronger observation-driven latent factor. This is understandable because the CSI300 index includes more constituent stocks than the SSE50 index, making the observation process { Q t } contain more latent information. Recall model (2); the compensation coefficients β 2 * and γ 2 * represent the degree of risk compensation associated with hitting limit-down, and reflect the decisions of investors with different risk perspectives when hitting limit-down. For both the SSE50 dataset and the CSI300 dataset, the compensation coefficients of the CAcF-E model are negative, while the compensation coefficients of the CAcF-L and CAcF-S models are positive. This result indicates that, in the Chinese mainland stock market, risk lovers regard hitting limit-down as a signal to buy stocks, while risk-neutral or risk-averse agents regard hitting limit-down as a signal that triggers further deterioration of stock markets. We can also see from Table 1 that the significance of the CAcF model coefficients is significantly better than that of the AcF model coefficients. This result is also strong evidence that the CAcF models are more suitable for price-limited markets.
Furthermore, we present the CAcF-E modeling results for the SSE50 in detail, where Figure 2 and Figure 3 show the estimation results of the tail index { α ^ t } and the volatility index { σ ^ t } , respectively. More results for other models and the CSI300 dataset are given in Appendix D. In Figure 2, the red line represents the estimated { α ^ t } and the green line represents the one-step prediction of the tail index base on historical information F t 1 . A smaller α ^ t indicates greater tail risk. Overall, the estimated tail index of the CAcF-E model is roughly in the range of 2 to 12. The two periods in which the tail index is consistently smaller ( α t < 4 ) coincide with the Subprime Crisis and the Chinese stock market crash in 2015. There is a significant correlation between the upward and downward trends in the closing price series of SSE50 and the estimated tail index, making α ^ t a valid indicator of the tail risk of the underlying market, i.e., the market stability index. In addition, before the SSE50 index declined significantly during the 2015 Chinese stock market crash, the estimated tail index based on the CAcF-E model had already begun to fall suddenly and sharply, providing an early warning of significant risks in the market. The performance of the tail index estimated via the CAcF-E model is consistent with those empirical findings of Massacci [22] and Zhao et al. [23] in markets without a price-limit policy. Moreover, the top-right subfigure presents the point-wise 90% confidence interval for the tail index based on historical information F t 2 . (The closed-form expressions for confidence intervals of α ^ t and σ ^ t are shown in Appendix C.) Instead of using a one-step point forecast, investors can use the lower bound of the confidence interval of the tail index as an earlier risk warning. Analogously, volatility risk in the market is also of great concern to investors. Figure 3 represents the estimated { σ ^ t } and the one-step prediction of the volatility together with the point-wise 90% confidence interval in the top-right subfigure. For comparison, we also fit a GARCH(1,1) model for the cross-sectional maxima of the negative daily returns of the stocks in SSE50 and plot the standardized daily volatility estimated by the GARCH model in Figure 3 (blue dashed line). The two series are very close, suggesting that CAcF’s dynamic scaling parameter σ t (including its forecasts) accurately measures market volatility.
To verify the discussion in Remark 5, we employ two measures of sequence similarity, namely Cosine Distance (CosD) and mean Euclidean Distance (MED), defined as follows
CosD ( a t , b t ) = 1 a t b t / a t 2 b t 2 and MED ( a t , b t ) = 1 n ( a t b t ) 2 ,
where { a t , 1 t n } and { b t , 1 t n } denote two different estimated dynamic parameter sequences. A pairwise comparison of the estimates of the dynamic scale parameter and the shape parameter for the three different CAcF models is presented in Table 2. It is observed that the two index sequences obtained from the different CAcF-type models are highly similar, regardless of whether it is an angular measure (CosD) or a distance measure (MED). This result suggests that the three CAcF models only decompose the { σ t } process and the { α t } process from different risk preference perspectives through different observation-driven functions, but these three models are consistent in their estimates of the { σ t } process and the { α t } process. In addition, combining the close estimated location parameters of the three CAcF models, this result further indicates that the three CAcF models have a consensus on the estimation of tail risk.

3.1.2. Tail Risk Analysis and Forecasting Performance

In this subsection, we apply the CAcF models to assess time-varying tail risk. Figure 4 shows both in-sample and out-of-sample forecasts of CAcF-based entropic value at risk with order p = 1 of the Q t * series in SSE50 at the 90%, 95%, and 99% confidence levels (solid lines denote in-sample and dashed lines denote out-of-sample). Notably, CAcF-based EVaR shows its sensitivity to market turbulence and reaches high levels around 2007–2008 and 2015, which is consistent with the actual events of the Subprime Crisis and the Chinese stock market crash in 2015. The strength of its fluctuation is also consistent with the brown line below indicating the negative daily return of the SSE50 index. Figure 5 elucidates the probability of actual observations being censored—specifically, the likelihood of encountering the limit-down scenario, as inferred from the latent maximum series Q t * distribution, which is derived from the parameters estimated by both the CAcF-E and AcF models. A significant takeaway from this analysis is the discernibly higher probability of breaching price limits when acknowledging the censored nature of the data. This observation intimates a more dire market scenario than might otherwise be perceived. The rationale behind this finding is intuitively straightforward: when the censored data characteristic is overlooked, a 10 % decline is interpreted as a genuine market movement. However, this interpretation often overlooks the critical role of price-limit constraints, which, if absent, would likely result in declines exceeding 10 % , thereby indicating a heightened risk scenario.
To illustrate the forecasting performance of the CAcF models, we employ three evaluation measures to evaluate the out-of-sample performance of recursive prediction from 4 January 2021 to 30 December 2022, and the results are detailed in Table 3. The three metrics include mean absolute error (MAE), mean absolute percentage error (MAPE), and mean censored probability (MCP), defined as follows,
MAE = 1 n t | Q ^ t Q t | , MAPE = 1 n t | ( Q ^ t Q t ) / Q t | , MCP = 1 n t | P ^ t D t | ,
where Q ^ t denotes the prediction of Q t and P ^ t = E ^ [ D t F t 1 ] denotes the probability of hitting limit-down calculated based on the distribution of Q t * in model (2). In both the SSE50 or CSI300 datasets, all three CAcF models obtained smaller values of MAE, MAPE, and MCP, i.e., better predictive performance, than the benchmark AcF model. For example, for the out-of-sample prediction of Q t , the MAPE of the CAcF-L model for CSI300 is 0.2706 , which is 9 % less than 0.2973 in the AcF model. Further, for the out-of-sample prediction of D t , there is a 96 % reduction. The comparison results demonstrate that the censored structure of CAcF models effectively improves the prediction performance.

3.2. Chinese Taiwan Stock Market

For the Taiwan stock market, the price-limit policy has been revised several times. The price limits were ± 5 % until January 1989 , ± 7 % from January 1989 to 30 May 2015, and ± 10 % thereafter. We collect the historical stock data from 4 January 2005 to 30 December 2022, with a total of 4431 observations. Since the sample period covers two different price-limit policies, we divide the data into two periods and fit them with the CAcF model. Period I is from 4 January 2005 to 30 May 2015, while period II is from 1 June 2015 to 30 December 2022. The price limits of these two periods are ± 7 % and ± 10 % , respectively.
In this part of our study, we are mainly concerned with two questions. First, we want to verify whether the patterns we observed in the stock market of mainland China are also present in other markets. Second, we are interested in how the relaxation of the price limits affects the tail risk of the stock market.

3.2.1. Fitting Results for TW50

To accurately capture the variations in tail risk after the implementation of a new price limit in the Taiwan stock market, we do not split the sample data for out-of-sample validation in this study. The parameter estimation results for the TW50 index are described in Table 4. These findings echo the insights garnered from the analysis of the SSE50 and CSI300 indices, underscoring that employing the AcF model for censored observations could lead to an inaccurate assessment of tail risk. Notably, for both periods under review, the estimated autoregressive coefficients β 1 and γ 1 exhibit values ranging between 0.82 and 0.90 , with minimal variance between the two periods. This observation suggests a notable consistency in the persistence of both the σ t and the α t processes across the different periods. Regarding the compensation coefficients β 2 * , γ 2 * , the estimated outcomes suggest that the investment behavior in the Taiwan stock market aligns with that observed within the Chinese mainland stock market, indicating a similar risk disposition among investors. Additionally, we compute the pairwise similarity of the dynamic scale parameter and the shape parameter across the three distinct CAcF models applied to the TW50 dataset. The outcomes of this analysis are presented in Table 5, offering yet another strong confirmation of Remark 5.

3.2.2. Empirical Analysis on Widening Price Limits

The adjustment of price limits in the Taiwan stock market offers an opportunity to examine the impact of broader price limits on tail risk. To this end, we reconstruct the tail index { α t } , volatility index { σ t } , and tail risk processes for each of the two periods using the CAcF models. Figure 6 showcases the estimated tail index α ^ t processes for both periods as derived through the CAcF-E model. Notably, there is no abrupt transition between the two tail index processes at the juncture of price-limit adjustment, with the estimated tail index for Period II exhibiting a broader range of movement compared to Period I. This suggests that the immediate probability of extreme events does not significantly alter following the introduction of new price limits. However, over the longer term, the likelihood of extreme events during tranquil market periods decreases, yet remains significant during turbulent times. This implies that, while broader price limits may reduce the frequency of hitting limit-down events, they do not mitigate tail risks during periods of market distress. This observation is partially supported by the fact that the incidence of limit-down events was 3.73% in Period I and reduced to 2.95% in Period II.
The estimated volatility index processes σ ^ t are depicted in Figure 7, alongside the daily returns of the TW50 index. For both periods, the estimated volatility index process σ ^ t accurately captures the fluctuations in the TW50’s daily returns. However, unlike the tail index estimate α ^ t , there is a noticeable jump between the volatility indices for the two periods. This discrepancy indicates that the volatility of the cross-sectional maximum of negative daily returns for stocks listed on the TW50 index significantly improved after the expansion of price limits.
Figure 8 and Figure 9 depict the evolution of time-varying tail risk for the TW50 index, demonstrating the behavior of CAcF-based EVaR with order p = 1 and the estimated probability of encountering censored events. These illustrations reveal pronounced increases in both the CAcF-based EVaR and the estimated censored probability during the 2007–2008 and 2020–2021 intervals, corresponding to the periods of the Subprime Mortgage Crisis and the COVID-19 pandemic. This correlation underscores the responsiveness of these risk metrics to market turmoil. Furthermore, the apex observed during the 2020–2021 period surpasses that of the 2007–2008 phase, suggesting that the widening of price limits not only enhances market volatility, as depicted in Figure 7, but also exacerbates tail risks. After relaxing price limits, the CAcF-based EVaR initially rises, then subsequently falls, reflecting a sophisticated market’s reaction to regulatory changes through the lens of price fluctuations. This pattern suggests that, while broader price limits can lead to initial instability, they may ultimately contribute to a more resilient market environment by allowing more freedom for price adjustments in response to market conditions.

4. Conclusions

This paper proposes a novel censored autoregressive conditional Fréchet model for analyzing the dynamic tail risk of financial markets constrained with price limits. Through empirically investigating the entropic value at risk of stock data from the Chinese mainland and Taiwan, we find that tail risk is seriously underestimated in price-limited stock markets when ignoring the censoring nature of data. We simultaneously decompose the tail risk from different risk preference perspectives via CAcF-type models with different observation-driven functions. The results offer guidelines and explanations for the decision-making for investors with diverse risk preferences dealing with risky events (hitting limit-down). Moreover, we explore the similarities and differences in the tail risk of the TW50 index over two periods with different price limits to study the impact of widening the price limit. The results suggest that widening the price limit would reduce the incidence of events of hitting limit-down; however, it also results in significant increases in the level of tail risk. The conclusions have important implications for policymakers in reconsidering market stabilization mechanisms. Finally, the out-of-sample performance of the CAcF model demonstrates its efficacy in forecasting tail risk and thus proves its potential for widespread application in tail risk monitoring.
However, the current study has some limitations, such as the lack of a way to analyze the risk relationship between multiple markets and the inability to include more lagged information. Therefore, in future applications, we would extend the CAcF model to the multivariate setting via a flexible extreme value copula to study the risk spillover effect among different stock markets. Another extended solution is to incorporate an ARMA structure for the time-varying parameters to utilize more historical information to enhance the forecasting performance of tail risk.

Author Contributions

Methodology, T.X. and L.S.; software, T.X.; validation, T.X.; formal analysis, T.X., L.S. and Y.C.; investigation, Y.C.; writing—original draft, T.X. and L.S.; writing—review and editing, Y.C.; visualization, T.X.; supervision, L.S. and Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Social Science Fund of China (No. 22BTJ027).

Data Availability Statement

Data available from Yahoo Finance: https://finance.yahoo.com/ (accessed on 4 May 2024).

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. The Significance of the Autoregression Coefficients

This section demonstrates the significance of the autoregression coefficients in model (2) concerning the μ , β 0 , β 1 , β 2 , β 2 * , β 3 , γ 0 , γ 1 , γ 2 , γ 2 * , and γ 3 parameters of the three types of CAcF models.
Table A1. The significance of the autoregression coefficients in model (2).
Table A1. The significance of the autoregression coefficients in model (2).
CoefficientSignificance
μ The location parameter in the conditional Fréchet distribution.
β 0 The constant term, determining the range of movement of { log σ t } process.
β 1 The autoregression coefficient for { log σ t } process.
β 2 The coefficient related to G 1 Q t 1 , representing the impact of Q t 1 through some observation-driven function G 1 ( · ) on log σ t .
β 2 * The compensation coefficient related to G 1 Q t 1 and the state D t 1 , indicating the impact of hitting limit-down on { log σ t } .
β 3 The coefficient embedded in the observation-driven function G 1 ( · ) of exponential type for { log σ t } process, which controls the shape of the function.
γ 0 The constant term, determining the range of movement of { log α t } process.
γ 1 The autoregression coefficient for log α t 1 .
γ 2 The coefficient related to G 2 Q t 1 , representing the impact of Q t 1 through some observation-driven function G 2 ( · ) on log α t .
γ 2 * The compensation coefficient related to G 1 Q t 1 and the state D t 1 , indicating the impact of hitting limit-down on { log α t } .
γ 3 The coefficient embedded in the observation-driven function G 2 ( · ) of exponential type for { log α t } process, which controls the shape of the function.

Appendix B. Proof of Proposition 1

We first present the following definition of a random variable in the Domain of Attraction of Fréchet distribution in Leadbetter et al. [35]:
Definition A1. 
A random variable ε is in the Domain of Attraction of Fréchet distribution with tail index α if and only if x F = and 1 F ε ( x ) l ( x ) x α , α > 0 , where F ε is the cdf of ε , l ( x ) is a slowly varying function, and x F = sup x : F ε ( x ) < 1 . Here and after, for two positive functions m 1 ( x ) and m 2 ( x ) , m 1 ( x ) m 2 ( x ) means m 1 ( x ) m 2 ( x ) 1 , as x .
Proof of Proposition 1. 
In the following, the conditioning on F t 1 is implicit and we drop the time index t for notation simplicity. Namely, Q N = max 1 i N X i , a N = 0 , and b N = i = 1 N σ i t α t 1 / α t . Denote Q N * = max 1 i N X i * and set A N = { exist i such that X i * M , 1 i N } , for any δ > 0 , we have
P lim N | Q N max Q N * , M | < δ = P A N P lim N | M M | < δ A N + P A N c P lim N | Q N * Q N * | < δ A N c = P A N + P A N c = 1 .
Namely, Q N = max Q N * , M , a . s . , as N .
By Proposition 1 of Zhao et al. [23], we have
Q N * a N b N d Ψ α .
Then, we can conclude that
Q N a N b N d max Ψ α , ( M a N ) / b N ,
by the continuous mapping theorem. □

Appendix C. The Confidence Intervals of α ^ t and σ ^ t

Here, we construct confidence intervals of α ^ t and σ ^ t based on historical information F t 2 . Recall Equation (2); it is obvious that Y t 1 , Q t 1 F t 1 and α t 1 , σ t 1 F t 2 . Combining the monotonic relationship among α t , σ t , and Q t 1 , given the historical information F t 2 , we can obtain the τ % confidence intervals for α ^ t and σ ^ t , respectively, as follows:
γ ^ 0 + γ ^ 1 log α ^ t 1 + γ ^ 2 + γ ^ 2 * D t 1 ( 1 + τ ) / 2 G ^ 2 Q t 1 ( 1 + τ ) / 2 , γ ^ 0 + γ ^ 1 log α ^ t 1 + γ ^ 2 + γ ^ 2 * D t 1 ( 1 τ ) / 2 G ^ 2 Q t 1 ( 1 τ ) / 2 , and β ^ 0 + β ^ 1 log σ ^ t 1 + β ^ 2 + β ^ 2 * D t 1 ( 1 τ ) / 2 G ^ 1 Q t 1 ( 1 τ ) / 2 , β ^ 0 + β ^ 1 log σ ^ t 1 + β ^ 2 + β ^ 2 * D t 1 ( 1 + τ ) / 2 G ^ 1 Q t 1 ( 1 + τ ) / 2 ,
where D t τ = I [ Q t τ , * > M ] , Q t τ = min ( M , Q t τ , * ) , and Q t τ , * = μ ^ + α ^ t ( log τ ) 1 / σ ^ t is τ -quantile of Q t * .

Appendix D. Appendix of Empirical Results

Appendix D.1. Complementary Empirical Results for SSE50 and CSI300

We present here complementary empirical results for the SSE50 and CSI300 datasets. The two graphs on the left side in Figure A1 show the in-sample fitting and out-of-sample forecasts of the tail index { α ^ t } in the SSE50 based on the CAcF-L and CAcF-S models. Correspondingly, the two graphs on the right show the in-sample fitting and out-of-sample predicted volatility { σ ^ t } . These results are highly similar to Figure 2 and Figure 3, and further verify that the different CAcF models have a consensus on the estimated { α ^ t } and { σ ^ t } processes, but different CAcF models decompose the two processes from different risk preference perspectives. Thus, although investors have different risk preferences, the risk factors embedded in the market are clear, only with different interpretations. The two graphs on the left in Figure A2 show the in-sample fitting and out-of-sample predicted results of the CAcF-L and CAcF-S models for the entropic value at risk of Q t * in SSE50. The brown line indicates the negative daily returns of the SSE50 index (the positive daily returns of the SSE50 index are cut off to reflect the downward risk better). Correspondingly, the two graphs on the right show the in-sample fitting and out-of-sample predicted censored probability. Similarly, this is consistent with the results presented in Figure 4 and Figure 5. The straight fact is that, if the censored features of the data are not accounted for, the tail risks in the price-limited market will be underestimated.
Figure A1. The estimated tail index { α ^ t } and volatility { σ ^ t } by CAcF-L and CAcF-S (red line) and the corresponding predicted tail index (blue line) for the SSE50 dataset.
Figure A1. The estimated tail index { α ^ t } and volatility { σ ^ t } by CAcF-L and CAcF-S (red line) and the corresponding predicted tail index (blue line) for the SSE50 dataset.
Entropy 26 00555 g0a1
Figure A2. The uncensored entropic VaR and censored probability estimated via CAcF-L and CAcF-S models for the SSE50 dataset.
Figure A2. The uncensored entropic VaR and censored probability estimated via CAcF-L and CAcF-S models for the SSE50 dataset.
Entropy 26 00555 g0a2
In addition, we also present the results for the CSI300 dataset, as shown in Figure A3 and Figure A4. These results are consistent with those of the SSE50 dataset, which demonstrates that our conclusions are also applicable to the CSI300 dataset.
Figure A3. The estimated tail index { α ^ t } and volatility { σ ^ t } by CAcF-E, CAcF-L, and CAcF-S models for the CSI300 dataset.
Figure A3. The estimated tail index { α ^ t } and volatility { σ ^ t } by CAcF-E, CAcF-L, and CAcF-S models for the CSI300 dataset.
Entropy 26 00555 g0a3
Figure A4. The uncensored entropic VaR and censored probability estimated via CAcF-E, CAcF-L, and CAcF-S models for the CSI300 dataset.
Figure A4. The uncensored entropic VaR and censored probability estimated via CAcF-E, CAcF-L, and CAcF-S models for the CSI300 dataset.
Entropy 26 00555 g0a4

Appendix D.2. Complementary Empirical Results for TW50

This section presents complementary empirical results for the TW50 index based on the CAcF-L and CAcF-S models. As shown in Figure A5 and Figure A6, these results are highly similar to the results based on the CAcF-E model (see more details in Section 3.2).
Figure A5. The estimated tail index { α ^ t } and volatility index { σ ^ t } by CAcF-L and CAcF-S for the TW50 dataset. The white area represents the result of Period I and the yellow area represents the result of Period II.
Figure A5. The estimated tail index { α ^ t } and volatility index { σ ^ t } by CAcF-L and CAcF-S for the TW50 dataset. The white area represents the result of Period I and the yellow area represents the result of Period II.
Entropy 26 00555 g0a5
Figure A6. The uncensored entropic VaR and censored probability estimated via CAcF-L and CAcF-S models for the TW50 dataset. The white area represents the result of Period I and the yellow area represents the result of Period II.
Figure A6. The uncensored entropic VaR and censored probability estimated via CAcF-L and CAcF-S models for the TW50 dataset. The white area represents the result of Period I and the yellow area represents the result of Period II.
Entropy 26 00555 g0a6

References

  1. Pichler, A.; Schlotter, R. Entropy based risk measures. Eur. J. Oper. Res. 2020, 285, 223–236. [Google Scholar] [CrossRef]
  2. Kou, S.; Peng, X. On the measurement of economic tail risk. Oper. Res. 2016, 64, 1056–1072. [Google Scholar] [CrossRef]
  3. Liu, F.; Wang, R. A theory for measures of tail risk. Math. Oper. Res. 2021, 46, 1109–1128. [Google Scholar] [CrossRef]
  4. Tan, K.; Chen, Y.; Chen, D. A New Risk Measure MMVaR: Properties and Empirical Research. J. Syst. Sci. Complex. 2023, 36, 2026–2045. [Google Scholar] [CrossRef]
  5. Zhu, D.; Zhang, C.; Pan, D. Extreme Risk Measurement of Carbon Market Considering Multifractal Characteristics. J. Syst. Sci. Complex. 2023, 36, 2497–2514. [Google Scholar] [CrossRef]
  6. Brennan, M.J. A theory of price limits in futures markets. J. Financ. Econ. 1986, 16, 213–233. [Google Scholar] [CrossRef]
  7. Veld-Merkoulova, Y.V. Price limits in futures markets: Effects on the price discovery process and volatility. Int. Rev. Financ. Anal. 2003, 12, 311–328. [Google Scholar] [CrossRef]
  8. Tang, S. Price limit performance: New evidence from a quasi-natural experiment in China’s ChiNext market. Int. Rev. Financ. Anal. 2023, 89, 102747. [Google Scholar] [CrossRef]
  9. Chan, S.H.; Kim, K.A.; Rhee, S.G. Price limit performance: Evidence from transactions data and the limit order book. J. Empir. Financ. 2005, 12, 269–290. [Google Scholar] [CrossRef]
  10. Chen, T.; Gao, Z.; He, J.; Jiang, W.; Xiong, W. Daily price limits and destructive market behavior. J. Econom. 2019, 208, 249–264. [Google Scholar] [CrossRef]
  11. Jin, S.; Zhou, C.; Peng, H. Does price limit reduce stock price volatility on the limit up and down day? Financ. Res. Lett. 2023, 58, 104302. [Google Scholar] [CrossRef]
  12. Kodres, L.E.; O’Brien, D.P. The existence of Pareto-superior price limits. Am. Econ. Rev. 1994, 84, 919–932. [Google Scholar]
  13. Fernandes, M.; Rocha, M. Are price limits on futures markets that cool? Evidence from the Brazilian Mercantile and Futures Exchange. J. Financ. Econom. 2007, 5, 219–242. [Google Scholar] [CrossRef]
  14. Deb, S.S.; Kalev, P.S.; Marisetty, V.B. Are price limits really bad for equity markets? J. Bank. Financ. 2010, 34, 2462–2471. [Google Scholar] [CrossRef]
  15. Wei, S.X. A censored–GARCH model of asset returns with price limits. J. Empir. Financ. 2002, 9, 197–223. [Google Scholar] [CrossRef]
  16. Hsieh, P.H.; Yang, J.J. A censored stochastic volatility approach to the estimation of price limit moves. J. Empir. Financ. 2009, 16, 337–351. [Google Scholar] [CrossRef]
  17. Embrechts, P.; Klüppelberg, C.; Mikosch, T. Modelling Extremal Events: For Insurance and Finance; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; Volume 33. [Google Scholar]
  18. McNeil, A.J.; Frey, R. Estimation of tail-related risk measures for heteroscedastic financial time series: An extreme value approach. J. Empir. Financ. 2000, 7, 271–300. [Google Scholar] [CrossRef]
  19. Pickands, J. Statistical inference using extreme order statistics. Ann. Stat. 1975, 3, 119–131. [Google Scholar]
  20. Oh, S.; Kee, H.; Park, K. Tail risk under price limits. Econ. Model. 2019, 77, 113–123. [Google Scholar] [CrossRef]
  21. Ji, J.; Wang, D.; Xu, D.; Xu, C. Combining a self-exciting point process with the truncated generalized Pareto distribution: An extreme risk analysis under price limits. J. Empir. Financ. 2020, 57, 52–70. [Google Scholar] [CrossRef]
  22. Massacci, D. Tail risk dynamics in stock returns: Links to the macroeconomy and global markets connectedness. Manag. Sci. 2017, 63, 3072–3089. [Google Scholar] [CrossRef]
  23. Zhao, Z.; Zhang, Z.; Chen, R. Modeling maxima with autoregressive conditional Fréchet model. J. Econom. 2018, 207, 325–351. [Google Scholar] [CrossRef]
  24. Shen, Z.; Chen, Y.; Shi, R. Modeling tail index with autoregressive conditional Pareto model. J. Bus. Econ. Stat. 2022, 40, 458–466. [Google Scholar] [CrossRef]
  25. Andreoni, J.; Sprenger, C. Risk preferences are not time preferences. Am. Econ. Rev. 2012, 102, 3357–3376. [Google Scholar] [CrossRef]
  26. Rieger, M.O.; Wang, M.; Hens, T. Risk preferences around the world. Manag. Sci. 2015, 61, 637–648. [Google Scholar] [CrossRef]
  27. Schildberg-Hörisch, H. Are risk preferences stable? J. Econ. Perspect. 2018, 32, 135–154. [Google Scholar] [CrossRef] [PubMed]
  28. Fisher, R.A.; Tippett, L.H.C. Limiting forms of the frequency distribution of the largest or smallest member of a sample. In Mathematical Proceedings of the Cambridge Philosophical Society; Cambridge University Press: Cambridge, UK, 1928; Volume 24, pp. 180–190. [Google Scholar]
  29. Gnedenko, B. Sur la distribution limite du terme maximum d’une serie aleatoire. Ann. Math. 1943, 44, 423–453. [Google Scholar] [CrossRef]
  30. Zakoian, J.M. Threshold heteroskedastic models. J. Econ. Dyn. Control. 1994, 18, 931–955. [Google Scholar] [CrossRef]
  31. Mas-Colell, A.; Whinston, M.D.; Green, J.R. Microeconomic Theory; Oxford University Press: New York, NY, USA, 1995; Volume 1. [Google Scholar]
  32. Fama, E.F.; French, K.R. Common risk factors in the returns on stocks and bonds. J. Financ. Econ. 1993, 33, 3–56. [Google Scholar] [CrossRef]
  33. Bai, J.; Ng, S. Evaluating latent and observed factors in macroeconomics and finance. J. Econom. 2006, 131, 507–537. [Google Scholar] [CrossRef]
  34. Fama, E.F.; French, K.R. A five-factor asset pricing model. J. Financ. Econ. 2015, 116, 1–22. [Google Scholar] [CrossRef]
  35. Leadbetter, M.R.; Lindgren, G.; Rootzén, H. Extremes and Related Properties of Random Sequences and Processes; Springer Series in Statistics; Springer: Berlin/Heidelberg, Germany, 1983. [Google Scholar]
  36. Rockafellar, R.T.; Uryasev, S. Conditional value-at-risk for general loss distributions. J. Bank. Financ. 2002, 26, 1443–1471. [Google Scholar] [CrossRef]
  37. Arici, G.; Campi, M.C.; Carè, A.; Dalai, M.; Ramponi, F.A. A theory of the risk for empirical CVaR with application to portfolio selection. J. Syst. Sci. Complex. 2021, 34, 1879–1894. [Google Scholar] [CrossRef]
Figure 1. The left panel shows the histogram of the cross-sectional maxima of the negative daily returns of stocks in the SSE50 index for the period from 4 January 2005 to 30 December 2022; and the right panel shows the tail quantiles estimated by Fréchet distribution and censored Fréchet distribution, respectively.
Figure 1. The left panel shows the histogram of the cross-sectional maxima of the negative daily returns of stocks in the SSE50 index for the period from 4 January 2005 to 30 December 2022; and the right panel shows the tail quantiles estimated by Fréchet distribution and censored Fréchet distribution, respectively.
Entropy 26 00555 g001
Figure 2. The figure shows the estimated tail index { α ^ t } by CAcF-E (red line) and the corresponding predicted tail index (green line). The blue line represents the closing price of SSE50. The top-right subfigure presents the 90% confidence interval for the predicted tail index.
Figure 2. The figure shows the estimated tail index { α ^ t } by CAcF-E (red line) and the corresponding predicted tail index (green line). The blue line represents the closing price of SSE50. The top-right subfigure presents the 90% confidence interval for the predicted tail index.
Entropy 26 00555 g002
Figure 3. Estimated volatility: CAcF v.s. GARCH. The red solid line is the estimated { σ ^ t } by CAcF-E, the green solid line is the corresponding predicted value, and the blue dashed line is the estimated volatility by GARCH(1,1) model. All series are standardized to be zero mean and unit variance for comparison. The top-right subfigure presents the 90% confidence interval for the predicted volatility.
Figure 3. Estimated volatility: CAcF v.s. GARCH. The red solid line is the estimated { σ ^ t } by CAcF-E, the green solid line is the corresponding predicted value, and the blue dashed line is the estimated volatility by GARCH(1,1) model. All series are standardized to be zero mean and unit variance for comparison. The top-right subfigure presents the 90% confidence interval for the predicted volatility.
Entropy 26 00555 g003
Figure 4. The figure shows the in-sample and out-of-sample predicted results for CAcF-based entropic value at risk (VaR) of the Q t * in SSE50. The brown line indicates the negative daily returns of the SSE50 index (the positive daily returns of the SSE50 index are cut off to reflect the downward risk better).
Figure 4. The figure shows the in-sample and out-of-sample predicted results for CAcF-based entropic value at risk (VaR) of the Q t * in SSE50. The brown line indicates the negative daily returns of the SSE50 index (the positive daily returns of the SSE50 index are cut off to reflect the downward risk better).
Entropy 26 00555 g004
Figure 5. The figure shows the in-sample and out-of-sample predicted results for CAcF-based censored probability of the cross-sectional maxima of the negative daily returns of the stocks in SSE50.
Figure 5. The figure shows the in-sample and out-of-sample predicted results for CAcF-based censored probability of the cross-sectional maxima of the negative daily returns of the stocks in SSE50.
Entropy 26 00555 g005
Figure 6. The estimated tail index { α ^ t } of TW50 by CAcF-E, where the white area indicates results for Period I and the yellow area for Period II. The blue line illustrates the TW50’s closing price.
Figure 6. The estimated tail index { α ^ t } of TW50 by CAcF-E, where the white area indicates results for Period I and the yellow area for Period II. The blue line illustrates the TW50’s closing price.
Entropy 26 00555 g006
Figure 7. The estimated volatility index { σ ^ t } of TW50 by CAcF-E, where the white area indicates results for Period I and the yellow area for Period II. The blue line illustrates the TW50’s closing price.
Figure 7. The estimated volatility index { σ ^ t } of TW50 by CAcF-E, where the white area indicates results for Period I and the yellow area for Period II. The blue line illustrates the TW50’s closing price.
Entropy 26 00555 g007
Figure 8. The CAcF-based entropic value at risk (EVaR) for the cross-sectional maximum, Q t * , of negative daily returns for TW50 stocks is depicted, with the white area illustrating results from Period I and the yellow area from Period II. The brown line represents the negative daily returns of the TW50 index, with positive daily returns excluded to better emphasize downward risk.
Figure 8. The CAcF-based entropic value at risk (EVaR) for the cross-sectional maximum, Q t * , of negative daily returns for TW50 stocks is depicted, with the white area illustrating results from Period I and the yellow area from Period II. The brown line represents the negative daily returns of the TW50 index, with positive daily returns excluded to better emphasize downward risk.
Entropy 26 00555 g008
Figure 9. The censored probability derived from the CAcF-E model for the cross-sectional maxima of negative daily returns of TW50 stocks, with the white and yellow areas indicating the outcomes for Period I and Period II, respectively.
Figure 9. The censored probability derived from the CAcF-E model for the cross-sectional maxima of negative daily returns of TW50 stocks, with the white and yellow areas indicating the outcomes for Period I and Period II, respectively.
Entropy 26 00555 g009
Table 1. In-sample estimation results for SSE50 and CSI300.
Table 1. In-sample estimation results for SSE50 and CSI300.
SSE50CSI300
CAcF-L CAcF-S CAcF-E AcF CAcF-L CAcF-S CAcF-E AcF
μ 0.124 *** 0.123 *** 0.122 *** 3.538 0.182 *** 0.180 *** 0.179 *** 5.864
0.004 0.005 0.004 7.240 0.006 0.005 0.006 ( 4.938 )
β 0 0.297 *** 0.309 *** 0.206 *** 0.167 0.514 *** 0.497 *** 0.415 *** 1.901 ***
0.048 0.051 0.011 0.260 0.056 0.056 0.018 0.066
β 1 0.852 *** 0.842 *** 0.884 ***0.869 *** 0.680 *** 0.683 *** 0.703 *** 0.636 ***
0.023 0.026 0.006 0.018 0.035 0.036 0.013 0.031
β 2 0.420 *** 4.638 *** 0.035 *** 0.002 0.723 *** 5.596 *** 0.087 *** 1.257
0.066 0.775 0.004 0.005 0.088 0.799 0.009 0.758
β 2 * 0.077 0.011 0.194 *** 0.305 *** 2.139 ** 0.331 ***
0.105 1.731 0.012 0.069 0.851 0.019
β 3 24.867 *** 13.577 *** 20.839 *** 0.025
( 2.556 ) ( 5.236 ) ( 1.476 ) ( 0.091 )
γ 0 0.447 *** 0.358 *** 6.421 *** 4.829 *** 0.839 *** 0.733 *** 2.847 *** 0.932
0.056 0.045 0.376 0.392 0.120 0.107 0.573 ( 2.7044 )
γ 1 0.834 *** 0.852 *** 0.848 *** 0.845 *** 0.700 *** 0.713 *** 0.712 *** 0.646 ***
0.021 0.019 0.005 0.004 0.045 0.044 0.027 0.034
γ 2 3.007 *** 31.229 *** 6.831 *** 5.750 *** 3.108 *** 27.584 *** 3.652 *** 1.206
0.336 3.301 0.389 0.238 0.474 4.068 0.637 2.689
γ 2 * 1.117 ** 1.453 0.118 *** 1.699 *** 11.849 ** 0.184 ***
0.624 8.694 0.023 0.435 4.799 0.027
γ 3 0.411 *** 0.564 *** 0.854 *** 3.047
( 0.018 ) ( 0.025 ) ( 0.133 ) ( 2.133 )
Note: The numbers in the brackets are the corresponding standard deviations. The symbol − indicates that the model excludes the corresponding parameter and the symbols *, **, and *** represent statistical significance at 10%, 5%, and 1% levels, respectively.
Table 2. Sequence similarity for the scale and shape parameters as estimated across the three CAcF models.
Table 2. Sequence similarity for the scale and shape parameters as estimated across the three CAcF models.
SSE50 ( σ ^ t ( L ) , σ ^ t ( E ) ) ( σ ^ t ( L ) , σ ^ t ( S ) ) ( σ ^ t ( E ) , σ ^ t ( S ) ) ( α ^ t ( L ) , α ^ t ( E ) ) ( α ^ t ( L ) , α ^ t ( S ) ) ( α ^ t ( E ) , α ^ t ( S ) )
CosD 3.84 × 10 5 0.0001 2.24 × 10 5 4.99 × 10 5 0.0008 0.0010
MED 0.0021 0.0008 0.0022 0.0974 0.2791 0.2711
CSI300 ( σ ^ t ( L ) , σ ^ t ( E ) ) ( σ ^ t ( L ) , σ ^ t ( S ) ) ( σ ^ t ( E ) , σ ^ t ( S ) ) ( α ^ t ( L ) , α ^ t ( E ) ) ( α ^ t ( L ) , α ^ t ( S ) ) ( α ^ t ( E ) , α ^ t ( S ) )
CosD 1.57 × 10 5 4.66 × 10 5 1.12 × 10 5 1.73 × 10 5 0.0003 0.0003
MED 0.0034 0.0024 0.0018 0.1034 0.1487 0.1360
Table 3. Out-of-sample forecasting performance of the three CAcF models and the AcF model.
Table 3. Out-of-sample forecasting performance of the three CAcF models and the AcF model.
SSE50CSI300
CAcF-L CAcF-S CAcF-E AcF CAcF-L CAcF-S CAcF-E AcF
MAE 0.0187 0.0188 0.0186 0.0195 0.0186 0.0188 0.0187 0.0195
MAPE 0.3160 0.3102 0.3143 0.3312 0.2706 0.2713 0.2709 0.2973
MCP 0.0051 0.0042 0.0052 0.0113 0.0055 0.0056 0.0068 0.1483
Note: Bold characters indicate the best prediction result.
Table 4. Parameter estimation results for TW50.
Table 4. Parameter estimation results for TW50.
Period IPeriod II
CAcF-L CAcF-S CAcF-E AcF CAcF-L CAcF-S CAcF-E AcF
μ 0.112 *** 0.105 *** 0.108 *** 2.769 0.136 *** 0.133 *** 0.136 *** 0.376
0.007 0.006 0.007 5.428 0.009 0.005 0.010 ( 0.242 )
β 0 0.338 *** 0.373 *** 0.250 *** 1.957 0.264 *** 0.289 *** 0.180 *** 0.077 ***
0.069 0.082 0.038 2.162 0.058 0.061 0.025 0.056
β 1 0.837 *** 0.821 *** 0.853 *** 0.872 *** 0.861 *** 0.846 *** 0.870 *** 0.897 ***
0.033 0.039 0.024 0.112 0.029 0.032 0.020 0.022
β 2 0.472 *** 6.438 *** 0.065 *** 1.827 0.377 *** 4.479 *** 0.071 *** 0.020
0.097 1.391 0.017 4.427 0.082 0.934 0.017 0.020
β 2 * 0.257 * 3.362 0.039 *** 0.022 0.101 0.011
0.148 2.513 0.009 0.171 0.924 0.008
β 3 10.391 *** 0.012 6.738 *** 8.833
( 3.268 ) ( 0.336 ) ( 1.630 ) ( 9.927 )
γ 0 0.438 *** 0.394 *** 5.312 *** 3.276 * 0.454 *** 0.400 *** 4.936 *** 4.700 ***
0.089 0.087 0.204 1.955 0.115 0.103 0.253 ( 0.184 )
γ 1 0.850 *** 0.844 *** 0.854 *** 0.858 *** 0.844 *** 0.847 *** 0.847 *** 0.878 ***
0.032 0.035 0.003 0.121 0.041 0.040 0.004 0.002
γ 2 2.939 *** 35.548 *** 5.734 ***4.120 ** 2.470 *** 25.857 *** 5.383 *** 5.183 ***
0.585 7.388 0.212 2.034 0.601 6.032 0.263 0.176
γ 2 * 1.118 12.960 0.079 ** 1.090 4.186 0.112 **
0.847 10.510 0.037 1.248 8.918 0.051
γ 3 0.510 *** 0.747 0.462 *** 0.472 ***
( 0.062 ) ( 0.818 ) ( 0.066 ) ( 0.053 )
Note: The numbers in the brackets are the corresponding standard deviation. The symbol − denotes the exclusion of a parameter from the model. The symbols *, **, and *** indicate statistical significance at 10%, 5%, and 1% levels, respectively.
Table 5. Sequence similarity for the scale and shape parameters estimated by the three CAcF models.
Table 5. Sequence similarity for the scale and shape parameters estimated by the three CAcF models.
Period I ( σ ^ t ( L ) , σ ^ t ( E ) ) ( σ ^ t ( L ) , σ ^ t ( S ) ) ( σ ^ t ( E ) , σ ^ t ( S ) ) ( α ^ t ( L ) , α ^ t ( E ) ) ( α ^ t ( L ) , α ^ t ( S ) ) ( α ^ t ( E ) , α ^ t ( S ) )
CosD 2.11 × 10 6 1.71 × 10 5 8.82 × 10 6 3.41 × 10 6 0.0003 0.0003
MED 7.97 × 10 5 0.0001 6.17 × 10 5 0.0052 0.0102 0.0059
Period II ( σ ^ t ( L ) , σ ^ t ( E ) ) ( σ ^ t ( L ) , σ ^ t ( S ) ) ( σ ^ t ( E ) , σ ^ t ( S ) ) ( α ^ t ( L ) , α ^ t ( E ) ) ( α ^ t ( L ) , α ^ t ( S ) ) ( α ^ t ( E ) , α ^ t ( S ) )
CosD 1.36 × 10 6 2.32 × 10 5 1.35 × 10 5 1.75 × 10 6 0.0004 0.0004
MED 6.38 × 10 6 5.67 × 10 5 5.78 × 10 5 0.0004 0.0077 0.0078
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, T.; Shu, L.; Chen, Y. Tail Risk Dynamics under Price-Limited Constraint: A Censored Autoregressive Conditional Fréchet Model. Entropy 2024, 26, 555. https://doi.org/10.3390/e26070555

AMA Style

Xu T, Shu L, Chen Y. Tail Risk Dynamics under Price-Limited Constraint: A Censored Autoregressive Conditional Fréchet Model. Entropy. 2024; 26(7):555. https://doi.org/10.3390/e26070555

Chicago/Turabian Style

Xu, Tao, Lei Shu, and Yu Chen. 2024. "Tail Risk Dynamics under Price-Limited Constraint: A Censored Autoregressive Conditional Fréchet Model" Entropy 26, no. 7: 555. https://doi.org/10.3390/e26070555

APA Style

Xu, T., Shu, L., & Chen, Y. (2024). Tail Risk Dynamics under Price-Limited Constraint: A Censored Autoregressive Conditional Fréchet Model. Entropy, 26(7), 555. https://doi.org/10.3390/e26070555

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