Next Article in Journal
Symmetric Collaborative Fault-Tolerant Control of Multi-Intelligence under Long-Range Transmission in Air–Ground Integrated Wireless High-Mobility Self-Organizing Networks
Next Article in Special Issue
Refinement of a Lyapunov-Type Inequality for a Fractional Differential Equation
Previous Article in Journal
On Uniformly Starlike Functions with Respect to Symmetrical Points Involving the Mittag-Leffler Function and the Lambert Series
Previous Article in Special Issue
A Functional Inequality and a New Class of Probabilities in the N-Person Red-and-Black Game
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Bilinear Time Series Models with Time-Varying and Symmetric GARCH Coefficients: Estimation and Simulation

1
Department of Mathematics, Al-Zaytoonah University of Jordan, Amman 11733, Jordan
2
Department of Basic Sciences, Al-Ahliyya Amman University, Amman 19328, Jordan
3
Department of Mathematics, Seville University, 41004 Sevilla, Spain
4
Department of Mathematics, University of Oum El Bouaghi, Oum El Bouaghi 04000, Algeria
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(5), 581; https://doi.org/10.3390/sym16050581
Submission received: 26 January 2024 / Revised: 21 April 2024 / Accepted: 26 April 2024 / Published: 8 May 2024
(This article belongs to the Special Issue Advance in Functional Equations, Second Edition)

Abstract

:
This paper makes a significant contribution by focusing on estimating the coefficients of a sample of non-linear time series, a subject well-established in the statistical literature, using bilinear time series. Specifically, this study delves into a subset of bilinear models where Generalized Autoregressive Conditional Heteroscedastic (GARCH) models serve as the white noise component. The methodology involves applying the Klimko–Nilsen theorem, which plays a crucial role in extracting the asymptotic behavior of the estimators. In this context, the Generalized Autoregressive Conditional Heteroscedastic model of order (1,1) noted that the GARCH (1,1) model is defined as the white noise for the coefficients of the example models. Notably, this GARCH model satisfies the condition of having time-varying coefficients. This study meticulously outlines the essential stationarity conditions required for these models. The estimation of coefficients is accomplished by applying the least squares method. One of the key contributions lies in utilizing the fundamental theorem of Klimko and Nilsen, to prove the asymptotic behavior of the estimators, particularly how they vary with changes in the sample size. This paper illuminates the impact of estimators and their approximations based on varying sample sizes. Extending our study to include the estimation of bilinear models alongside GARCH and GARCH symmetric coefficients adds depth to our analysis and provides valuable insights into modeling financial time series data. Furthermore, this study sheds light on the influence of the GARCH white noise trace on the estimation of model coefficients. The results establish a clear connection between the model characteristics and the nature of the white noise, contributing to a more profound understanding of the relationship between these elements.

1. Introduction

In recent decades, a diverse range of time series models has emerged, thanks to the contributions of an expanding community of scholars. This interest has attracted not only mathematicians but also economists and sociologists, who recognize the broad applicability of these models in finance and economics. Nonlinear time series models, in particular, have gained prominence due to their ability to capture the nonlinear characteristics often observed in empirical finance models describing volatility and returns. These models find applications across various disciplines, including marketing, insurance, virology [1], and chemistry.
One of the two models introduced by Ganger and Andersen [2] in 1978 is a bilinear time series model. Applications of bilinear models can be found in the domains of engineering, medicine, and biology. Later research, like that carried out by Subba [2], concentrated on these models’ statistical and probabilistic aspects.
Numerous approaches, such as moments and least squares (LS), have been used to address the estimating problem for bilinear models. Meanwhile, both Autoregressive Conditional Heteroskedasticity (ARCH) and Generalized Autoregressive Conditional Heteroskedasticity (GARCH) models are important in financial time series analysis. They are used to estimate economic activities, analyze holding risk, and assess option pricing.
As a subclass of nonlinear models, bilinear time series models have attracted interest because of their capacity to capture complex relationships in data. These models find use in a variety of industries, including engineering, medicine, and biology. Bilinear models have been studied since Ganger and Andersen’s 1978 paper [1], when curiosity about their statistical and probabilistic characteristics first arose.
In financial time series analysis, Generalized Autoregressive Conditional Heteroskedasticity (GARCH) models have become essential tools. GARCH models were developed to capture time-varying variances and volatility clustering. They offer insights into economic forecasting, asset pricing, risk management, and option valuation. Prominent scholarly contributions, including those from the authors, highlighted in the introduction (references [3,4,5,6], for example), have demonstrated GARCH models’ applicability in various financial contexts.
When it comes to bilinear time series models, the least squares method is a basic estimation technique. The least squares approach, which is popular due to its ease of use and effectiveness, attempts to reduce the sum of squared residuals to provide estimates for the model parameters. The least squares method is a useful technique for estimating parameters in the domain of bilinear models, where interactions between variables may display nonlinearity. Its use in bilinear models makes it easier to understand how variables relate to one another and makes it possible to simulate intricate, time-varying events.
This paper is organized as follows. Section 2 provides notifications and significant results on G A R C H with bilinear time series as preliminaries. Section 3, using the nature of the proposed model, presents the estimation approach for all coefficients of bilinear mixed models using white noise G A R C H . Section 4 illustrates the asymptotic behavior of the estimators and their mechanisms through numerical examples and simulations, employing the Klimko–Nilsen theorem. Section 5 discusses the importance of the paper and its differences from other published works, defining the content of the new work related to this research. To enrich our study, a comparison was made between two bilinear models, the first followed by GARCH and the second by Symmetric GARCH.

2. Preliminaries

In 2004, Bibi studied, in reference [6], a sample of bilinear models with time-varying coefficients that we define in a probability space ( Ω , Γ , P ) according to the following general stochastic expression:
X t = i = 1 p u i , t ( a ) X t i + j = 1 q v t , j ( a ) ε t j + i = 1 p 0 j = 1 q 0 ϕ t , i j ( a ) X t i ε t j + ε t .
The process ( X t ) t Z denoted by B L ( p , q , p 0 , q 0 ) , where
{ u i , t ( a ) t } 1 i p , { v j ( t ) } 1 j q , { ϕ i j ( t ) } , 1 i p 0 , 1 j q 0
are usually constants, but there will be time-varying coefficients that take their values in R ; where a is a vector, the white noise sequence ( ε t ) t Z takes many forms and is generally distributed according to Gaussian law with zero mean and variance σ 2 < . But we will aim at a specified sample of bilinear coefficient models provided by G A R C H ( 1 , 1 ) white noise.
ε t = h t η t , h t 2 = γ 0 + α t ( α ) ε t 1 2 + β t ( β ) h t 1 2 .
where the general GARCH ( p 0 , q 0 ) model is defined as follows:
ε 2 ( t ) = h 2 ( t ) η 2 ( t ) h 2 ( t ) = γ 0 + i = 1 p 0 α i ε t i 2 + j = 1 q 0 β i h t j 2
Note that with the distribution η t N ( 0 , 1 ) , the sequences α t ( α ) and β t ( β ) are time-varying coefficients while maintaining γ 0 > 0 as a constant. h t is independent of the σ −field generated by η t + k , k 0 . And ε t is a measurable function of the variables η t , 0 . It is well known that the strong GARCH(1,1) model is strictly stationary if and only if
E log α t ( α ) η t 2 + β t ( β ) < 0 .
We define the function φ as following φ ( η t 2 ) = α t ( α ) η t 2 + β t ( β ) , and through this function, we can write h t 2 as the following form:
h t 2 = γ 0 + α t ( α ) h t 1 2 η t 1 2 + β t ( β ) h t 1 2 = γ 0 + α t ( α ) η t 1 2 + β t ( β ) h t 1 2 = γ 0 + φ ( η t 1 2 ) h t 1 2 .
The last expression can be written recurrently using the following formula:
h t 2 = γ 0 + γ 0 i = 1 q j = 1 i φ ( η t j 2 ) + j = 1 q + 1 φ ( η t j 2 ) h t q 1 2 .
According to the stability conditions, we can demonstrate that E ( h t 2 ) < ; the proof will be checked, where it suffices to show that the following product is less than 1. Then, according to Jensen inequality, q tends to infinity, so
h t 2 = γ 0 1 + i = 1 j = 1 i φ ( η t j 2 ) .
We can prove that lim q j = 1 q + 1 φ ( η t j 2 ) = 0 ; it is clear that
lim q E j = 1 q + 1 φ ( η t j 2 ) j = 1 E φ ( η t j 2 ) = j = 1 α t ( α ) + β t ( β ) .
By considering the following inequality α t ( α ) + β t ( β ) = δ < 1 , we will ensure that
j = 1 δ j , δ j = δ = c ,
where c is a constant, (7) tends towards zero, which demonstrates the following:
lim q j = 1 q + 1 φ ( η t j 2 ) = 0 .
In another way, we can prove that the series is almost surely convergent, h t 2 , if we take the general term of the series and apply the Cauchy criterion of convergence:
lim n j = 1 n φ ( η t j 2 ) 1 n = lim n e 1 n j = 1 n ln φ ( η t j 2 ) .
Also, if v t is an i.i.d sequence of random variables admitting an expectation which can be infinite, then 1 n k = 1 n v k is approximate towards E ( v 1 ) . So, through this construction and using Jensen’s inequality we obtain the following:
lim n e 1 n j = 1 n ln φ ( η t j 2 ) = lim n e E ( ln φ ( η t 1 2 ) e ln E φ ( η t 1 2 ) = e ln ( α t ( α ) + β t ( β ) ) < 1 ,
which will show that the series is convergent.

2.1. Stationarity Study

In the following subsection, we will deal with the following bilinear model:
X t = ϕ t ( a ) X t s ε t 1 + ε t ,
where s 1 , and ε t is the white noise that follows the GARCH (1,1) model with time-varying coefficients, which was defined above using Expression (1); ϕ t ( a ) , t Z is a sequence of time-varying coefficients, and a is defined as a vector a = ( a 1 , a 2 , , a m ) included in subset Θ of R m . It is well known that stationary solutions exist for this model if ϕ t 2 ( a ) E ( ε t 2 ) < 1 , where we can give some extensions of stability according to the white noise and, recurrently, ε t is written as follows:
ε t = X t + k = 1 t 1 ( 1 ) k i = 0 k 1 ϕ t k ( a ) × i = 0 k 1 X t i s X t k .
In the case where
E ( ε t 2 ) = E ( η t 2 ) E ( h t 2 ) = E ( h t 2 ) = γ 0 + α t ( α ) E ( h t 1 2 ) E ( η t 1 2 ) + β t ( β ) E ( h t 1 2 ) = γ 0 + α t ( α ) E ( ε t 1 2 ) + β t ( β ) E ( ε t 1 2 ) .
And, where E ( ε t 2 ) = E ( ε t 1 2 ) , we find that E ( ε t 2 ) = γ 0 1 α t ( α ) β t ( β ) , such as α t ( α ) + β t ( β ) < 1 ; then, the necessary conditions for model stability will be
ϕ t 2 ( a ) γ 0 1 α t ( α ) β t ( β ) < 1 .
And the recurring expression of the model in case, where s 1 , will be
X t = ε t + j = 1 t / s 1 i = 0 j 1 ϕ t i s ( a ) ε t i s 1 ε t s j .
where y denotes the integer part of value y , Model (3) under a condition of stability and the unique solution with recursive form when s = 1 :
X t = ε t + j = 1 i = 0 j 1 ϕ t i ( a ) ε t i 1 ε t j .
Theorem 1.
The model solution series (12) under a condition of stability is almost surely convergent.
Proof. 
Our method of proof will be based on showing that E ( X t ) < . So, by using the Schwartz inequality, we have
S = E ε t j i = 0 j 1 ϕ t i ( a ) ε t i 1 { E ( ε t j 2 ) } 0.5 i = 0 j 1 ϕ t i ( a ) E ( ε t i 1 2 ) 0.5 { E ( h t j 2 ) } 0.5 i = 0 j 1 ϕ t i ( a ) E ( ε t i 1 2 ) 0.5 .
and as ϕ t i ( a ) E ( ε t i 1 2 ) 0.5 = ρ < 1 , { E ( h t j 2 ) } 0.5 = M is bounded:
j = 1 E ε t j i = 0 j 1 ϕ t i ( a ) ε t i 1 M j = 1 ρ j < M 1 ρ .
By showing that the series is a converging geometric series, our proof is complete here. □
Theorem 2.
Consider the stochastic process defined by Model (3), driven by the noise ARCH(1). If the following conditions hold,
1.
For all t, α t ( α ) [ 0 , 1 ) ,
2.
ϕ t ( a ) γ 0 1 α t ( α ) < 1 ,
then the series
X t = ε t + j = 1 i = 0 j 1 ϕ t ( a ) ε t i 1 ε t j = f
almost surely converges, and the defined X t is the unique, strictly stationary solution.
Proof. 
Firstly, we show that E i = 0 j 1 ϕ t ( a ) ε t i 1 < , where j is any integer, by using the Schwartz inequality and the condition of stability:
E i = 0 j 1 ϕ t ( a ) ε t i k 1 i = 0 j 1 E ϕ t ( a ) ε t i k 1 2 0.5 = E ϕ t ( a ) ε t i k 1 2 0.5 = ϕ t ( a ) γ 0 1 α t ( α ) ζ < 1 ,
then
j = 1 i = 0 j 1 E ϕ t ( a ) ε t i k 1 j = 1 ξ i < ,
where s 1 , because of the recurrent model
X t = ε t + j = 1 t / s 1 i = 0 j 1 ϕ t i s ( a ) ε t i s 1 ε t s j = g .
Here, y denotes the greatest integer less than or equal to y. It is worth noticing that g f implies E g E f . □
Theorem 3.
Let s = 1 in Model (3) be driven by the noise ARCH(1). If the condition E ln ϕ t ( a ) ε t , 0 holds, then it implies that, for α t ( α ) [ 0 , 1 ] ,
ϕ t ( a ) γ 0 1 α t ( α ) < 1 ,
where α t ( α ) represents the parameter in the range [ 0 , 1 ] .
Proof. 
Using Jensen’s inequality, we have
E ln ϕ t ( a ) ε t ln E ϕ t ( a ) ε t .
Now, by applying the Schwartz inequality, we obtain
ln E ϕ t ( a ) ε t ln E ϕ t ( a ) ε t 2 0.5 = ln ϕ t ( a ) γ 0 1 α t ( α ) 0.5 < 0 ,
where we observe that ϕ t ( a ) γ 0 1 α t ( α ) 0.5 < 1 , ensuring the logarithm is negative.
Thus, we have shown that the condition E ln ϕ t ( a ) ε t , 0 implies that ϕ t ( a ) γ 0 1 α t ( α ) < 1 , for α t ( α ) [ 0 , 1 ] , which completes the proof. □
Theorem 4.
For the ARCH(1) model, under the condition α t ( α ) < 1 and γ 0 = 0 , the process is bounded.
Proof. 
Since ε t = η t h t and h t = α t ( α ) ε t 1 , we find that
h t i = k = 1 i α t ( α ) ε t i 1 ,
then, lim i h t i = 0 , demonstrating that there exists M such that ε t M . □

2.2. Klimko–Nilsen Theorem and Estimation Approach

This theorem played a fundamental role in proving the existence and uniqueness of the estimators, as well as their asymptotic behaviors according to the approach of square worlds. Here is its text in the form of the hypotheses.
Firstly, let H N ( X t ) be the σ -field generated by the set of observations { X t , t = 1 , , N } and the vectors a = ( a 1 , a 2 , , a n 1 ) ; as for the white noise coefficients, we consider here the vectors α = ( α 1 , α 2 , , α n 2 ) and β = ( β 1 , β 2 , , β n 3 ) , where the image of these vectors with their functions exists in R . Then, the parameter that we are going to estimate will be ω = ( a , γ 0 , α , β ) , and we assume ω 0 is a true value of ω included in an open Ω of R n 1 + n 2 + n 3 + 1 . We introduce the predictor or orthogonal projection ρ t t 1 on the observations up to time t 1 by the following difference: ρ t t 1 ( ω ) = X t ε t ( ω ) .
So, the principle of the least squares method is based on looking for the parameter ω , such that the penalty function will be defined with the following standard algorithm of estimation:
w 1 = a 1 , , ω n 1 = a n 1 , ω n 1 + 1 = γ 0 , ω n 1 + 2 = α 1 , , ω n 1 + n 2 + 1 = α n 2 , ω n 1 + n 2 + 2 = β 1 , , ω n 1 + n 2 + n 3 + 1 = β n 3 .
Theorem 5.
Under stable conditions for ( X t ) t Z generated by Equation (8), and by adding ϕ ( a ) < 1 , where its white noise follows GARCH(1,1), we suppose that ρ t t 1 ( ω ) is almost surely doubly continuously differentiable in an open subset containing the true value ω 0 of the vector ω. We assume two constants, K 0 and K 1 , such that
(a) 
E ω 0 𝜕 ε t 2 ( ω ) 𝜕 ω i 4 K 0 , i = n 1 + n 2 + n 3 + 1 .
(b) 
E ω 0 𝜕 2 ε t 2 ( ω ) 𝜕 ω i 𝜕 ω j E ω 0 𝜕 2 ε t 2 ( ω ) 𝜕 ω i 𝜕 ω j H t 1 2 K 1 , i , j = n 1 + n 2 + n 3 + 1 .
(c) 
1 2 N t = 1 N E ω 0 𝜕 2 ε t 2 ( ω ) 𝜕 ω i 𝜕 ω j H N 1 converges almost surely to the matrix M ( ω ) which is strictly positive.
(d) 
lim N sup δ 0 1 δ N t = 1 N 𝜕 2 ε t 2 ( ω ) 𝜕 ω i 𝜕 ω j ω = ω ˜ 𝜕 2 ε t 2 ( ω ) 𝜕 ω i 𝜕 ω j ω = ω 0 < , where ω ω 0 < δ , δ > 0 , where ω ˜ represents the intermediate value between ω and ω 0 .
Then, there exists an estimator ω ^ N such that ω ^ N ω 0 as N , if these conditions are satisfied as well as the following assumption:
(e) 
1 N t = 1 N E ω 0 𝜕 ε t 2 ( ω ) 𝜕 ω 𝜕 ε t 2 ( ω ) 𝜕 ω T H t 1 E ω 0 𝜕 ε t 2 ( ω ) 𝜕 ω 𝜕 ε t 2 ( ω ) 𝜕 ω T 0 as N .
Proof. 
See [7]. □
We will prove some elements of this theorem according to the model proposed in this paper. First, we will give some techniques to calculate the derivatives. Let
𝜕 ε t 2 ( ω ) 𝜕 ω i = 2 ε t ( ω ) 𝜕 ε t ( ω ) 𝜕 ω i = 2 ε t ( ω ) 𝜕 ρ t t 1 ( ω ) 𝜕 ω i .
Then, 8 E ω 0 ε t 4 ( ω ) = E ω 0 𝜕 ρ t t 1 ( ω ) 𝜕 ω i 4 .
So, in case s = 1 , we obtain
ρ t t 1 ( ω ) = j = 1 i = 0 j 1 ϕ t i ( a ) ε t i 1 ε t j .
In the case where ω p = a p , p = 1 , , n , we will use
𝜕 ρ t t 1 ( ω ) 𝜕 ω i = k = 1 t 1 ( 1 ) k 𝜕 𝜕 a i i = 0 k 1 ϕ t k ( a ) × i = 0 k 1 X t i s X t k .
And
Z 1 = 𝜕 𝜕 a p i = 0 k 1 ϕ t k ( a ) = p = 0 k 1 i = 0 i p k 1 ϕ t i ( a ) 𝜕 ϕ t p ( a ) 𝜕 a p ,
where we put q = max 𝜕 ϕ t i ( a ) 𝜕 a p , and ϕ t i ( a ) θ < 1 ; then, it will be
Z 1 q ( k ) θ k 1 .
We will obtain
𝜕 ρ t t 1 ( ω ) 𝜕 ω i k = 1 t 1 ( 1 ) k q ( k ) θ k 1 i = 0 k 1 X t i s X t k < .
In addition, we establish E 𝜕 ρ t t 1 ( ω ) 𝜕 ω i < . According to previous analysis, given E ε t 2 ( ω ) = E η 2 E ( h t 2 ) , where E ( h t 2 ) is bounded, it follows that E ε t 4 ( ω ) < . Consequently, we deduce E ω 0 𝜕 ε t 2 ( ω ) 𝜕 ω i 4 K 0 , thus validating the hypothesis. Further details are provided in the subsequent sections (see [7]).

3. Least Squares Approach

The least squares estimators for ARCH models are asymptotically normal, though less efficient than methods such as the generalized method of moments (GMM) and maximum likelihood estimation (MLE) [7,8,9]. However, in this section, we will focus on the least squares method, particularly in situations where the coefficients are time-varying and driven by white noise with time-varying coefficients.

3.1. Algorithm

The primary algorithm of the least squares method aims to find the estimator ω ^ N by minimizing the penalty function q N ( ω ) , defined as follows:
arg min ω Ω q N ( ω ) ,
where
q N ( ω ) = 1 N t = 1 N ε t 2 ( ω ) .

3.2. Derivation Techniques

The least squares approach relies on Taylor’s second-degree formula. For any i , j { 1 , , n 1 + n 2 + n 3 + 1 } , we have the following:
q N ( ω ) = q N ( ω 0 ) + ( ω ω 0 ) 𝜕 q N 𝜕 ω T ( ω 0 ) + 1 2 ( ω ω 0 ) 𝜕 2 q N 𝜕 ω T 𝜕 ω ( ω ˜ ) ( ω ω 0 ) T ,
where ω ˜ is an intermediate point between ω and ω 0 with ω ω ˜ δ , δ > 0 . The derivatives 𝜕 q N ( ω ) 𝜕 ω i form a vector with n 1 + n 2 + n 3 + 1 coordinates, and 𝜕 2 q N ( ω ) 𝜕 ω i 𝜕 ω j presents a matrix of size ( n 1 + n 2 + n 3 + 1 ) 2 .

3.3. Derivatives Calculation

We derive expressions for the first and second derivatives concerning the parameters ω i and ω j . Firstly, for ω = a , we have
𝜕 ε t 2 ( ω ) 𝜕 ω = 2 ε t ( ω ) 𝜕 ε t ( ω ) 𝜕 ω = 2 ε t ( ω ) 𝜕 ρ t t 1 ( ω ) 𝜕 ω ,
where
𝜕 ρ t t 1 ( ω ) 𝜕 ω = k = 1 t 1 ( 1 ) k 𝜕 𝜕 a i = 0 k 1 ϕ t k ( a ) i = 0 k 1 X t i s X t k .
Similarly, for ω = α , the expression is given by
𝜕 ε t 2 ( ω ) 𝜕 ω = 4 η t 2 α t ( α ) 1 β t ( β ) ε t 1 𝜕 ε t 1 𝜕 ω ,
with
𝜕 ε t 1 𝜕 ω i = k = 1 t 2 ( 1 ) k 𝜕 𝜕 ω i i = 0 k 1 ϕ t 1 k ( a ) i = 0 k 1 X t 1 i s X t 1 k .
For the second derivative, considering ω i = a i and ω j = α j , we have
𝜕 ε t 2 ( ω ) 𝜕 ω i 𝜕 ω j = 4 η t 3 𝜕 ε t 1 𝜕 ω i × 𝜕 𝜕 ω j α t ( α ) 1 β t ( β ) h t + 𝜕 h t 𝜕 ω j α t ( α ) 1 β t ( β ) .
Additionally, when ω i = α i and ω j = α j , where i , j { 1 , , n 2 } , the second derivative is given by
𝜕 ε t 2 ( ω ) 𝜕 ω i 𝜕 ω j = 4 η t 3 𝜕 ε t 1 𝜕 ω i × 𝜕 𝜕 ω j α t ( α ) 1 β t ( β ) h t + 𝜕 h t 𝜕 ω j α t ( α ) 1 β t ( β ) .
These derivative expressions facilitate the computation of second-order terms in the least squares optimization process.
The derived equations lay the foundation for implementing the least squares approach, enabling the estimation of time-varying coefficients within the framework of ARCH models. This method offers a practical means of modeling volatility dynamics, particularly in financial time series data, where coefficients may exhibit temporal variations.
In summary, the least squares approach outlined here provides a systematic methodology for estimating parameters in ARCH models, leveraging optimization principles and derivative calculus to refine parameter estimates iteratively. Through careful consideration of the model structure and the properties of the underlying data, this approach offers valuable insights into the dynamics of volatility and enables effective forecasting in various domains.

4. Simulation and Graphic Illustrations

This section presents a simulation study and graphical illustrations for the model under consideration. We begin by defining a vector of parameters a = ( a 1 , a 2 , a 3 ) belonging to an open subset I 1 × I 2 × I 3 within R 3 . We impose the constraint a 2 = 1 a 1 and consider a 3 = max ( a 1 , a 2 ) . The time-varying coefficients are expressed as follows:
ϕ t ( a ) = a 1 , t 1 3 , a 2 , t 2 3 , a 3 , t 0 3 .
The objective of this section is to simulate the model X t = ϕ t ( a ) X t s ε t 1 + ε t with GARCH(1,1) white noise, where γ 0 = 0.005 , α = ( α 1 , α 2 ) , and β = ( β 1 , β 2 ) , such that
α t ( α 1 , α 2 ) = α 1 , t 2 N + 1 , α 2 , t 2 N . β t ( β 1 , β 2 ) = β 1 , t 2 N + 1 , β 2 , t 2 N .
We employ the notation of N s for the number of simulations, N for the sample size, the true parameter values ( a 1 , a 2 , a 3 ) , and their estimated counterparts ( a ^ 1 , a ^ 2 , a ^ 3 ) . In our simulations, we set s = 2 and adopt the GARCH white noise formulation:
ε t = η t h t h t 2 = 0.005 + α t ( 0.3 , 0.4 ) ε t 2 + β t ( 0.003 , 0.001 ) h t 2 .
where ω = ( 0.005 , 0.3 , 0.4 , 0.003 , 0.001 ) is the real value proposed for the model and ω ^ is its estimated value. In a situation where β t ( β ) = 0 , we have the following simulation.
By considering Table 1, Table 2 and Table 3, we can take the same model proposed but with constant coefficients and with the same GARCH proposed in our paper to extend our study of this model on the side of simulations and get the following Table 4.
X t = a X t s ε t 1 + ε t .
Notably, models featuring estimated coefficients and bilinear models guided by the GARCH(1,1) model with true coefficients exhibit compatibility as shown in Figure 1, Figure 2 and Figure 3. This compatibility effectively demonstrates that the estimation of these proposed models delivers efficient outcomes. Furthermore, these models validate the asymptotic behaviors of the estimators.

4.1. Asymmetric and Symmetric GARCH Models

The significance and impact of symmetric and asymmetric GARCH models in financial econometrics cannot be overstated. These models serve as indispensable tools for capturing the intricate dynamics of volatility in financial time series data.
Symmetric and asymmetric GARCH models offer distinct methodologies for volatility modeling. Symmetric GARCH models assume equal effects of positive and negative shocks on volatility, while asymmetric models such as EGARCH or GJR-GARCH allow for varying responses to such shocks. Depending on the characteristics of the data, one model type may yield more precise estimates of volatility dynamics than the other.
The choice of a GARCH model profoundly influences the accuracy of results in simulating future scenarios or conducting Monte Carlo simulations in financial modeling. Asymmetric GARCH models, adept at capturing the asymmetry and skewness commonly observed in financial returns data, often provide more realistic simulations of future volatility and asset prices. Moreover, the selection between symmetric and asymmetric GARCH models carries implications for risk management practices. Asymmetric GARCH models excel in capturing volatility clustering, a critical aspect for assessing and managing financial risk. However, it is essential to note that asymmetric GARCH models may require more computational resources for estimation and simulation, particularly when dealing with large datasets or extensive simulations. Researchers and practitioners must carefully consider computational efficiency alongside model accuracy when choosing between symmetric and asymmetric GARCH models.
Extending the analysis, GARCH models are recognized for their superiority in capturing large shocks and volatility clustering compared to ARCH models, as highlighted by Patton and Sheppard (2015). Nevertheless, GARCH models have notable limitations, particularly their assumption of symmetric volatility regardless of the nature of shocks, as emphasized by Lubrano (2001) [10]. In contrast, ARCH and GARCH models are considered more suitable for risk measurement than descriptive statistics, especially in scenarios where stock returns undergo frequent fluctuations and the distribution of returns is markedly non-symmetric, according to Ang and Bekaert (2007) [11]. Despite this, variance remains a trusted measure for positive volatility, an aspect often overlooked by stockholders, underscoring the issue of instability and the uneven distribution of stock returns. Also, see [12,13,14].
To address the challenge of asymmetric instability, the present study explores various models for measuring GARCH effects, with a specific focus on GARCH symmetry.
ε ( t ) = z ( t ) η ( t ) z 2 ( t ) = θ 0 + θ 1 ε 2 ( t 1 ) + θ 2 z 2 ( t 1 ) .
And in another way, the following model explains the Asymmetric GARCH model:
ε ( t ) = z ( t ) η ( t ) z ( t ) = θ 0 + θ 1 ε ( t 1 ) + θ 2 z ( t 1 ) .
To construct the asymmetry in a model driven by GARCH, we can write our model using the following expression:
X ( t ) = θ 0 X ( t s ) ε ( t 1 ) + ε ( t ) ε ( t ) G A R C H ( 1 , 1 ) .
The asymmetric or symmetric nature of the white noise GARCH component imparts its characteristic properties to the bilinear models of the time series. This influence is demonstrated through the resulting curves, which depict the intricate interplay between volatility dynamics and the underlying characteristics of the GARCH model. By incorporating asymmetric or symmetric GARCH into bilinear models, we can capture the nuanced features of volatility, including asymmetry, clustering, and the impact of shocks, thereby enhancing our ability to model and forecast financial time series data accurately. The curves serve as visual representations of how the properties of GARCH reverberate through the bilinear models, providing valuable insights into the complex interactions within the data and informing our understanding of financial market behavior. See [15,16,17,18].
The estimation of model coefficients oriented by a GARCH model gives us the estimated data with the following Table 5, where a is the true value of model (32), a ^ is its estimators, and ns represents the number of simulations
It can be observed that estimators of variable values in the bilinear model, as defined by its expression, yield efficient results when driven by Symmetric GARCH compared to Asymmetric GARCH, indicating that symmetry plays a crucial role in enhancing the accuracy of estimation for the proposed model [19,20,21,22,23,24,25].

4.2. Graphic Illustration

It is evident from the graphs that the impact of GARCH white noise on the bilinear model is pronounced. The asymmetry and symmetry inherent in GARCH manifest clearly in the original bilinear models depicted in the following Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9. Additionally, some deviations from symmetry are observable in the models, indicating that the time-varying coefficients, specifically the constant coefficients of the models, contribute to maintaining the properties of white noise. The final graph illuminates the overall symmetry of graphs driven by symmetric GARC.

5. Concluding Comments

The simulations conducted on the bilinear model outlined in our paper illustrate a notable trend: as the sample size increases, the estimators converge toward their true values. This underscores the critical role of sample size in achieving accurate approximations. Additionally, we observe that increasing the number of simulations enhances the proximity between the estimated values and their true counterparts.
Interestingly, the white noise ARCH(1) model demonstrates commendable accuracy in approximating the estimators compared to the GARCH model for the coefficients. Moreover, we consistently observe a close resemblance between the two graphs, indicating robustness in the estimation process.
It is noteworthy that there exists a critical sample size, beyond which slight perturbations may arise in our simulations. However, the asymptotic behaviors of the estimators remain unaffected, suggesting the efficiency of our estimation approach in capturing temporal variations in coefficients.
In summary, our estimation approach, incorporating time-varying coefficients in the bilinear model along with a GARCH white noise component, yields significant theoretical and numerical insights, as evidenced by the simulations. Future research endeavors will delve into exploring various models in economics that integrate a combination of GARCH models to better capture and model realistic data patterns.

Author Contributions

Conceptualization, M.A.H. and S.A.; Data curation, S.A.; Formal analysis, N.L. and A.O.; Funding acquisition, S.A.; Investigation, S.A.; Methodology, R.A.; Project administration, S.A.; Resources, N.L., A.O. and M.A.H.; Software, M.A.H.; Supervision, S.A.; Validation, S.A.; Visualization, A.O.; Writing—original draft, S.A.; Writing—review & editing, S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in thearticle, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declared that they have no conflict of interest.

References

  1. Granger, C.W.J.; Andersen, A.P. An introduction to bilinear time series models. Angew. Stat. Okonom. 1978, 8, 7–94. [Google Scholar]
  2. Rao, T.S.; Gabr, M.M. An Introduction to Bispectral Analysis and Bilinear Time Series Models; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; Volume 24. [Google Scholar]
  3. Pan, J.; Li, G.; Xie, Z. Stationary solution and parametric estimation for Bilinear model driven by ARCH noises. Sci. China Ser. A Math. 2002, 45, 1523–1537. [Google Scholar]
  4. Terdik, G. Bilinear Stochastic Models and Related Problems of Nonlinear Time Series Analysis: A Frequency Domain Approach; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1999; Volume 142. [Google Scholar]
  5. Tong, H. Non-Linear Time Series: A Dynamical System Approach; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  6. Bibi, A.; Oyet, A.J. Estimation of some bilinear time series models with time varying coefficients. Stoch. Anal. Appl. 2004, 22, 355–376. [Google Scholar] [CrossRef]
  7. Klimko, L.A.; Nelson, P.I. On conditional least squares estimation for stochastic processes. Ann. Stat. 1978, 6, 629–642. [Google Scholar] [CrossRef]
  8. Weiss, A.A. ARCH and bilinear time series models: Comparison and combination. J. Bus. Econ. 1986, 4, 59–70. [Google Scholar] [CrossRef]
  9. Pantula, S.G. Estimation of autoregressive models with ARCH errors. Sankhyā Indian J. Stat. Ser. B 1988, 50, 119–138. [Google Scholar]
  10. Lubrano, M. Smooth Transition Garch Models: A Bayesian Perspective. Louvain Econ. Rev. 2001, 67, 257–287. [Google Scholar] [CrossRef]
  11. Patton, A.J.; Sheppard, K. Good Volatility, Bad Volatility: Signed Jumps and The Persistence of Volatility. Rev. Econ. Stat. 2015, 97, 683–697. [Google Scholar] [CrossRef]
  12. Fiorucci, J.A.; Silva, G.N.; Barboza, F. Reaction trend system with GARCH quantiles as action points. Expert Syst. Appl. 2022, 198, 116750. [Google Scholar] [CrossRef]
  13. Hammad, M.A.; Alsharif, S.; Shmasnah, A.; Khalil, R. Fractional Bessel differential equation and fractional Bessel functions. Ital. J. Pure Appl. Math. 2022, 47, 521–531. [Google Scholar]
  14. Hammad, M.A.; Rasool, S.; Alotaibi, B.M.; Alotiby, M.; Tiofack, C.L.; Alrowaily, A.W.; El-Tantawy, S.A. On the modified versions of G’G-expansion technique for analyzing the fractional coupled Higgs system. AIP Adv. 2023, 13, 105131. [Google Scholar] [CrossRef]
  15. Gherghina, Ș.C.; Armeanu, D.Ș.; Joldeș, C.C. COVID-19 pandemic and Romanian stock market volatility: A GARCH approach. J. Risk Financ. Manag. 2021, 14, 341. [Google Scholar] [CrossRef]
  16. Zraiqat, A.; Al-Hwawcha, L.K. On exact solutions of second order nonlinear ordinary differential equations. Appl. Math. 2015, 6, 953–957. [Google Scholar] [CrossRef]
  17. Khan, M.S.; Khan, K.I.; Mahmood, S.; Sheeraz, M. Symmetric and asymmetric volatility clustering via GARCH family models: An evidence from religion dominant countries. Paradigms 2019, 13, 20–25. [Google Scholar]
  18. Warren, R.; Kaseke, F. Modeling South African Stock Market Volatility Using Univariate Symmetric and Asymmetric GARCH Models. Indian J. Financ. Bank. 2021, 6, 1–16. [Google Scholar]
  19. Carol, A.; Lazar, E.; Stanescu, S. Analytic moments for GJR-GARCH (1, 1) processes. Int. J. Forecast. 2021, 37, 105–124. [Google Scholar]
  20. Wang, Z.; Xing, T.; Wang, X. Economic uncertainty and stock market asymmetric volatility: Analysis based on the asymmetric GARCH-MIDAS model. Int. J. Emerg. Mark. 2024. [Google Scholar] [CrossRef]
  21. Xiong, L.; Zhu, F. Robust estimation for the one-parameter exponential family integer-valued GARCH (1,1) models based on a modified Tukey’s biweight function. Comput. Stat. 2024, 39, 495–522. [Google Scholar] [CrossRef]
  22. Zuo, H.; Guo, H. Structural nonlinear damage identification based on the information distance of GNPAX/GARCH model and its experimental study. Struct. Health Monit. 2024, 23, 991–1012. [Google Scholar] [CrossRef]
  23. Zhang, R. Asymmetric beta-binomial GARCH models for time series with bounded support. Appl. Math. Comput. 2024, 470, 128556. [Google Scholar] [CrossRef]
  24. Liu, H.; Yang, P.; Song, H.; Wu, D.C. Global and domestic economic policy uncertainties and tourism stock market: Evidence from China. Tour. Econ. 2024, 30, 567–591. [Google Scholar] [CrossRef]
  25. Amirshahi, B.; Lahmiri, S. Hybrid deep learning and GARCH-family models for forecasting volatility of cryptocurrencies. Mach. Learn. Appl. 2023, 12, 100465. [Google Scholar] [CrossRef]
Figure 1. Simulation for bilinear time series model with its true values N = 900 , N s = 250 .
Figure 1. Simulation for bilinear time series model with its true values N = 900 , N s = 250 .
Symmetry 16 00581 g001
Figure 2. Simulation for bilinear time series model replaced by its estimated values N = 900 , N s = 250 .
Figure 2. Simulation for bilinear time series model replaced by its estimated values N = 900 , N s = 250 .
Symmetry 16 00581 g002
Figure 3. Comparisonbetween the two models.
Figure 3. Comparisonbetween the two models.
Symmetry 16 00581 g003
Figure 4. The time plot of daily Iran stock returns.
Figure 4. The time plot of daily Iran stock returns.
Symmetry 16 00581 g004
Figure 5. The number of simulations.
Figure 5. The number of simulations.
Symmetry 16 00581 g005
Figure 6. Bilinear time series model with constant coefficients driven by Asymmetric GARCH, where N s = 1200 and N = 900 .
Figure 6. Bilinear time series model with constant coefficients driven by Asymmetric GARCH, where N s = 1200 and N = 900 .
Symmetry 16 00581 g006
Figure 7. The number of simulations.
Figure 7. The number of simulations.
Symmetry 16 00581 g007
Figure 8. Bilinear time series model with constant coefficients driven by Symmetric GARCH, where N s = 5000 and N = 1000 .
Figure 8. Bilinear time series model with constant coefficients driven by Symmetric GARCH, where N s = 5000 and N = 1000 .
Symmetry 16 00581 g008
Figure 9. Curve illustrates the symmetry of the bilinear model driven by symmetric GRACH compared with several lines.
Figure 9. Curve illustrates the symmetry of the bilinear model driven by symmetric GRACH compared with several lines.
Symmetry 16 00581 g009
Table 1. Estimation of model coefficients according to different sizes, where s = 1 .
Table 1. Estimation of model coefficients according to different sizes, where s = 1 .
N s NReal Values
( a 1 , a 2 , a 3 )
Estimates
( a ^ 1 , a ^ 2 , a ^ 3 )
250300 (0.0097, 0.3230, 0.3706)
600 (0.0182, 0.2872, 0.3797)
900(0.02, 0.25, 0.45)(0.0129, 0.2879, 0.3802)
500300 (0.0166, 0.2704, 0.4119)
600 (0.0225, 0.2787, 0.4152)
900 (0.0220, 0.2795, 0.4456)
Table 2. Estimation of model coefficients and ARCH white noise coefficients, where s = 1 .
Table 2. Estimation of model coefficients and ARCH white noise coefficients, where s = 1 .
N s NReal ValueEstimates ω ^
250300 0.015 , 0.098 , 0.235 , 0.082 , 0.033
600 0.009 , 0.165 , 0.385 , 0.007 , 0.014
900 ω 0.003 , 0.263 , 0.388 , 0.007 , 0.014
500300 0.003 , 0.369 , 0.387 , 0.006 , 0.003
600 0.003 , 0.369 , 0.387 , 0.006 , 0.002
900 0.004 , 0.318 , 0.407 , 0.010 , 0.002
Table 3. Estimation of the model according to other coordinates, for ( a 1 , a 2 , a 3 , γ 0 , α 1 , α 2 ) .
Table 3. Estimation of the model according to other coordinates, for ( a 1 , a 2 , a 3 , γ 0 , α 1 , α 2 ) .
N s N ( a ^ 1 , a ^ 2 , a ^ 3 , γ ^ 0 , α ^ 1 , α ^ 2 )
250300 ( 0.022 , 0.234 , 0.453 , 0.005 , 0.289 , 0.406 )
500600 ( 0.0198 , 0.252 , 0.398 , 0.005 , 0.301 , 0.411 )
where ( a 1 , a 2 , a 3 , γ 0 , α 1 , α 2 ) = ( 0.02 , 0.25 , 0.45 , 0.005 , 0.3 , 0.4 ) .
Table 4. Estimation of GARCH model coefficients for a bilinear model with constant coefficients.
Table 4. Estimation of GARCH model coefficients for a bilinear model with constant coefficients.
True Values: θ 0 = ( a , γ 0 , α 1 , α 2 , β 1 , β 2 ) = ( 0.07 , 0.01 , 0.1 , 0.3 , 0.6 , 0.7 )
N N s θ ^ = ( a ^ , γ ^ 0 , α ^ 1 , α ^ 2 , β ^ 1 , β ^ 2 )
300250 ( 0.0694 , 0.0124 , 0.0999 , 0.4179 , 0.5825 , 0.5125 )
600500 ( 0.0642 , 0.0115 , 0.1004 , 0.4145 , 0.5997 , 0.5122 )
Table 5. The estimation of model coefficients oriented by a GARCH model.
Table 5. The estimation of model coefficients oriented by a GARCH model.
ns N θ 0 θ ^ 0 Using Symmetric GARCH θ ^ 0 Using Asymmetric GARCH
1001200.090.08790.1672
1002400.090.09020.1677
1003000.090.08980.1399
3001200.20.19440.1944
3002400.20.19670.1967
3003000.20.20130.2013
9001200.70.59720.5972
9002400.70.68970.6897
9003000.70.70110.7011
120012000.90.90020.9038
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

Abu Hammad, M.; Alkhateeb, R.; Laiche, N.; Ouannas, A.; Alshorm, S. Comparative Analysis of Bilinear Time Series Models with Time-Varying and Symmetric GARCH Coefficients: Estimation and Simulation. Symmetry 2024, 16, 581. https://doi.org/10.3390/sym16050581

AMA Style

Abu Hammad M, Alkhateeb R, Laiche N, Ouannas A, Alshorm S. Comparative Analysis of Bilinear Time Series Models with Time-Varying and Symmetric GARCH Coefficients: Estimation and Simulation. Symmetry. 2024; 16(5):581. https://doi.org/10.3390/sym16050581

Chicago/Turabian Style

Abu Hammad, Ma’mon, Rami Alkhateeb, Nabil Laiche, Adel Ouannas, and Shameseddin Alshorm. 2024. "Comparative Analysis of Bilinear Time Series Models with Time-Varying and Symmetric GARCH Coefficients: Estimation and Simulation" Symmetry 16, no. 5: 581. https://doi.org/10.3390/sym16050581

APA Style

Abu Hammad, M., Alkhateeb, R., Laiche, N., Ouannas, A., & Alshorm, S. (2024). Comparative Analysis of Bilinear Time Series Models with Time-Varying and Symmetric GARCH Coefficients: Estimation and Simulation. Symmetry, 16(5), 581. https://doi.org/10.3390/sym16050581

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