Next Article in Journal
Mapping Risk–Return Linkages and Volatility Spillover in BRICS Stock Markets through the Lens of Linear and Non-Linear GARCH Models
Next Article in Special Issue
Environmental, Social and Governance Awareness and Organisational Risk Perception Amongst Accountants
Previous Article in Journal
Constructing Divisia Monetary Aggregates for the Asian Tigers
Previous Article in Special Issue
Extended Least Squares Making Evident Nonlinear Relationships between Variables: Portfolios of Financial Assets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Lower and Upper Estimates for Ether Option Prices with Conditional Heteroscedasticity and Model Uncertainty

by
Tak Kuen Siu
Department of Actuarial Studies and Business Analytics, Macquarie Business School, Macquarie University, Sydney, NSW 2109, Australia
J. Risk Financial Manag. 2024, 17(10), 436; https://doi.org/10.3390/jrfm17100436
Submission received: 2 August 2024 / Revised: 15 September 2024 / Accepted: 25 September 2024 / Published: 29 September 2024

Abstract

:
This paper aims to leverage Bayesian nonlinear expectations to construct Bayesian lower and upper estimates for prices of Ether options, that is, options written on Ethereum, with conditional heteroscedasticity and model uncertainty. Specifically, a discrete-time generalized conditional autoregressive heteroscedastic (GARCH) model is used to incorporate conditional heteroscedasticity in the logarithmic returns of Ethereum, and Bayesian nonlinear expectations are adopted to introduce model uncertainty, or ambiguity, about the conditional mean and volatility of the logarithmic returns of Ethereum. Extended Girsanov’s principle is employed to change probability measures for introducing a family of alternative GARCH models and their risk-neutral counterparts. The Bayesian credible intervals for “uncertain” drift and volatility parameters obtained from conjugate priors and residuals obtained from the estimated GARCH model are used to construct Bayesian superlinear and sublinear expectations giving the Bayesian lower and upper estimates for the price of an Ether option, respectively. Empirical and simulation studies are provided using real data on Ethereum in AUD. Comparisons with a model incorporating conditional heteroscedasticity only and a model capturing ambiguity only are presented.

1. Introduction

Cryptocurrency options are important instruments for managing risks of trading cryptocurrencies, speculating and leveraging future performances of cryptocurrencies. Bitcoin options, whose underlying assets could be Bitcoin or Bitcoin futures, appear to be the most popular among different cryptocurrency options. Indeed, Bitcoin options are less than ten years old. The first Bitcoin options were launched by Deribit (crypto options and futures exchange for Bitcoin and Ethereum) in 2016, while the first regulated Bitcoin option exchange is LedgerX, which started offering products in 2017. In the aftermath of the enlightening year of 2019 for the cryptocurrency derivatives market, other players started to roll out Bitcoin options. Some of them include Bakkt, CME, Bit.com, OKEx, Binance and Huobi futures. See Yang (2021) for related discussions and more detail about the history of the cryptocurrency derivatives market.
The highly volatile cryptocurrency market poses significant challenges to the valuation of cryptocurrency options. Some researchers have taken up the challenges and proposed quantitative models for pricing cryptocurrency options. Specifically, the uses of stochastic volatility and GARCH-type models for pricing Bitcoin options have received considerable attention in the literature (see, for example, (Chen and Yang 2024; Hou et al. 2020; Jalan et al. 2021; Siu and Elliott 2021; Venter and Maré 2021)). Jump-diffusion models for pricing Bitcoin options have also been considered (see, for example, (Cao and Celik 2021; Hilliard and Ngo 2022)). The uses of machine learning techniques and neural networks for pricing Bitcoin options or cryptocurrency options were considered in (Alexander et al. 2023; Li et al. 2019; Pagnottoni 2019). Alexander et al. (2023) considered the valuation of crypto quanto and inverse options. Lucic (2024) discussed the pricing and hedging of cryptocurrency inverse options. Siu (2024) considered an application of finite-mixture models to pricing Bitcoin options.
While Bitcoin is the largest cryptocurrency in terms of market capitalization, Ethereum (ETH) is the second largest cryptocurrency. The history of the blockchain technology backing Ethereum can be traced back to its creation by Vitalik Buterin in 2013, which was then launched in 2015 (see Copeland (2023)). Ether options, which are also called ETH options, are option contracts written on Ethereum. OKEx and CME launched Ether options in 2020 and 2022, respectively (see Godbole (2020) and Chan (2022)). CME also offered options written on Ether futures. According to McGleenon (2024), the trading volume of Ether options has surged to around USD 20 billion in January 2024. However, there is a relatively small amount of studies on the pricing of Ether options compared with those on the valuation of Bitcoin options. Some works on the pricing of Ether options include, for example, Shegokar (2022) for jump-diffusion stochastic volatility models and Brini and Lenz (2024) for a machine learning regression approach.
Modelling the return dynamics of cryptocurrencies is a challenging task. This is particularly the case in light of the impact of the pandemic. It is imperative to take into account model uncertainty, or ambiguity, about the return dynamics of cryptocurrencies (see, for example, Siu (2023a)). Early works on option pricing with model uncertainty focus on incorporating volatility uncertainty. Specifically, Boyle and Ananthanarayanan (1977) introduced a Bayesian approach to incorporate volatility uncertainty into option pricing. Avellaneda et al. (1995) and Lyons (1995) introduced volatility uncertainty for pricing and hedging options through “uncertain” volatility, where volatility is assumed to take on values in a bounded interval. These works seem to have provided a practical ground for the developments of nonlinear G-expectations, including the work by Peng (1997), which gives a mathematically sound approach to describe volatility uncertainty. Some other works on option pricing with model uncertainty, or volatility uncertainty in particular, include, for example, (Agarwal et al. 2021; Beissner and Riedel 2019; Boyle et al. 2008; Buraschi and Jiltsov 2006; Chau et al. 2022; Cohen and Tegnér 2019; Cont 2006; Elliott et al. 2022; Fouque and Ren 2014; Jaroszkowski and Jensen 2022; Siu 2023a; Wilmott and Oztukel 1998). It does not seem that the impacts of model uncertainty on pricing cryptocurrency options have been well explored or well understood up to now.
This paper aims to leverage Bayesian nonlinear expectations introduced in the work of Siu (2023a) to construct Bayesian lower and upper estimates for the price of an Ether option, that is, an option on Ethereum, in the presence of conditional heteroscedasticity and model uncertainty. Specifically, a discrete-time generalized conditional autoregressive heteroscedastic (GARCH) model introduced by Bollerslev (1986) and Taylor (1986) is adopted to describe conditional heteroscedasticity in the logarithmic returns of Ethereum. To incorporate model uncertainty or ambiguity about the conditional mean and volatility of the logarithmic returns, a pair of Bayesian nonlinear expectations, that is, a Bayesian superlinear expectation and a Bayesian sublinear expectation, are used to construct the Bayesian lower and upper estimates for the price of the Ether option. Firstly, using an extended Girsanov’s principle in discrete time, which is also called a discrete-time Girsanov’s transform (see (Elliott and Madan 1998; Elliott and Siu 2017; Elliott et al. 1995; Siu 2023a)), a family of alternative models capturing model uncertainty about the conditional mean and volatility is introduced via perturbing the conditional mean and volatility of the GARCH model. Then, applying a “risk-neutral” extended Girsanov’s principle to the family of alternative models results in a family of “risk-neutralized” alternative GARCH models. Lastly, using Bayesian credible intervals obtained from conjugate priors and residuals obtained from the estimated GARCH model, Bayesian conditional superlinear and sublinear expectations are constructed, which give the Bayesian lower and upper estimates for the price of the Ether option, respectively. Empirical and simulation studies are provided using real data on Ethereum in AUD (ETH-AUD). The impacts of model uncertainty, conditional heteroscedasticity and their interactions on the Ether option prices are explored by comparing different models.
Boyle and Ananthanarayanan (1977) seem to have been among the first to propose a Bayesian approach to incorporate model uncertainty, particularly volatility uncertainty, into option valuation. Other papers on the use of Bayesian statistics in option pricing include, for example, (Bauwens and Lubrano 2002; Darsinos and Satchell 2007; Foster and Whiteman 1999; Gzyl et al. 2022; Ho et al. 2011; Jang and Lee 2019; Karolyi 1993; Martin et al. 2005; Mozumder et al. 2021; Qin and Almeida 2020; Romboutsa and Stentoft 2014; Wang et al. 2016). Bayesian option pricing under GARCH-type models or stochastic volatility models were considered in (Bauwens and Lubrano 2002; Martin et al. 2005; Mozumder et al. 2021; Romboutsa and Stentoft 2014; Wang et al. 2016). However, it seems that those papers did not use Bayesian nonlinear expectations to describe the impacts of model uncertainty, or ambiguity, about the conditional mean and volatility on estimates for option prices. Furthermore, those papers did not consider the pricing of options on cryptocurrencies such as Ethereum. Siu (2001) considered the use of Bayesian statistics to construct Bayesian credible intervals for coherent risk measures of derivative securities under discrete-time binomial tree models. (Siu et al. 2001, 2004b) proposed the use of Bayesian statistics to construct Bayesian coherent risk measures and Bayesian Value at Risk, respectively, for derivative securities. However, these papers considered risk measures for derivative securities instead of option pricing. They did not consider cryptocurrencies.
The current paper may be thought of as a companion paper of the work of Siu (2023a). In Siu (2023a), the Bayesian nonlinear expectations were applied to provide lower and upper forecasts for Bitcoin returns and to construct lower and upper estimates for tail-based risk metrics, namely Value at Risk and Expected Shortfall, for a Bitcoin position. However, Siu (2023a) did not discuss the use of Bayesian nonlinear expectations for option valuation in general, or pricing Ether options in particular. Now, in the current paper, I endeavour to use the Bayesian nonlinear expectations for constructing Bayesian lower and upper estimates for the price of an Ether option. The current paper also aims to provide novel insights into advancing the understanding for the impacts of model uncertainty, conditional heteroscedasticity and their interactions on Ether option prices through the lens of the Bayesian nonlinear expectations. To this end, the Bayesian lower and upper estimates are interpreted as the “confidence” limits of Bayesian interval estimates for the price of an Ether option along the line of Boyle and Ananthanarayanan (1977). From the perspective of two-price conic finance (Cherny and Madan 2009; Elliott et al. 2022; Madan 2015, 2016; Madan and Cherny 2010), one may explore the possibility of studying a Bayesian version of lower and upper prices of an Ether option using the Bayesian lower and upper estimates. One may also explore the possibility of developing a Bayesian version of nonlinear valuation operators, for example, (Coque et al. 2002; Peng 2004a 2004b, 2004c, 2005). These will be briefly discussed as potential topics for further research in Section 8.
The rest of the paper is organized in the following manner. Section 2 presents a GARCH model without incorporating model uncertainty, namely a reference model, for the dynamics of the logarithmic returns from Ethereum. Then, the pricing of an Ether option under the GARCH model is discussed by reviewing some existing results on GARCH option pricing using an extended Girsanov’s principle. In Section 3, the family of alternative GARCH models and its risk-neutralization were introduced using the extended Girsanov’s principle, and lower and upper estimates for the price of the Ether option are constructed using nonlinear expectations. In Section 4, the Bayesian lower and upper estimates for the price of the Ether option are constructed using the Bayesian conditional superlinear and sublinear expectations, respectively. The numerical and simulation procedures are presented in Section 5. The empirical and simulation results are discussed in Section 6. Model comparisons are also presented in Section 6. The chaotic and random aspects of price fluctuations are briefly commented in Section 7. Concluding remarks are provided, and some limitations of the proposed approach are discussed in Section 8. Some potential topics for further research are briefly discussed in Section 8.

2. The GARCH Model without Model Uncertainty

A simple discrete-time market consisting of two underlying securities, namely a risk-free (zero-coupon) bond and a risky asset, that is, Ethereum, is considered. The logarithmic return process of Ethereum is modelled by a discrete-time GARCH model. Specifically, we start with a reference model under a real-world reference probability measure for the logarithmic return process of Ethereum in AUD (ETH-AUD). To simplify the discussion and notation, a GARCH(1,1) model for the logarithmic return process is considered. This simple set up can be generalized to a higher-order GARCH model. A discrete-time Girsanov’s transform in, for example, Elliott et al. (1995) and Elliott and Madan (1998), is used to specify a risk-neutral probability measure for valuing an Ether option under the reference GARCH(1,1) model. A discrete-time Girsanov’s transform was employed to specify an equivalent martingale measure for option pricing under discrete-time GARCH-type models (see, for example, Badescu et al. (2008) and Badescu et al. (2011)). Here, we state some relevant results and include the proofs for some of them using the results found in the work of Elliott et al. (1995). This section may serve as a recap on some existing results on GARCH option pricing in the literature. They are included here for the sake of completeness.
Let T denote the (finite) time parameter set { 1 , 2 , , T } of our discrete-time model, where T < . A complete probability space ( Ω , F , P ) is considered on which random processes are defined, where P is a reference (real-world) probability measure on the measurable space ( Ω , F ) . P is interpreted as the real-world probability measure underlying the reference model, which is also called an approximating model. A sequence of independent and identically distributed (i.i.d.) standard normal random variables under P , namely { ϵ t } t T , is considered. In general, one may consider a sequence of i.i.d. non-normal random variables. However, for illustration, a sequence of i.i.d. normal random variables is considered here. Write { Y t } t T for a time series process on ( Ω , F , P ) , which represents the logarithmic return process of Ethereum in AUD. That is, for each t T , Y t denotes the logarithmic return of Ethereum in AUD over the tth period. Assume that the reference model for the logarithmic return process { Y t } t T under P is governed by the GARCH(1,1) model:
Y t = μ Y + h t ϵ t , h t = ω + α h t 1 ϵ t 1 2 + β h t 1 .
Here, μ Y is the mean return; h t is the conditional variance when given the price information up to and including time t; ω , α and β are the parameters in the conditional variance process; and h 0 is a given constant.
In the sequel, a risk-neutral probability measure is specified using a discrete-time Girsanov’s transform. Let ϕ μ , σ 2 ( · ) denote the probability density function (pdf) of a normal distribution with mean μ and variance σ 2 . Write ϕ ( · ) for ϕ 0 , 1 ( · ) . That is, ϕ ( · ) is the pdf of a standard normal distribution with zero mean and unit variance (i.e., N ( 0 , 1 ) ). Let F denote the P -augmentation of the natural filtration { F t } t T generated by the logarithmic return process { Y t } t T . That is, for each t T , the σ -field F t represents the information set generated by the observations { Y 1 , Y 2 , , Y t } about the logarithmic return process up to and including time t.
For each t T , define the notation μ Y ( t ) by
μ Y ( t ) : = r 1 2 h t .
Then, an F -adapted process, denoted by { λ t } t T , is defined by putting
λ t : = ϕ μ Y ( t ) , h t ( Y t ) ϕ μ Y , h t ( Y t ) = ϕ ( Y t μ Y ( t ) h t ) ϕ ( Y t μ Y h t ) .
Furthermore, another F -adapted process, { Λ t } t T , is defined as follows:
Λ t : = k = 1 t λ k .
By definition, { Λ t } t T is a positive-valued process. It is also an ( F , P ) -martingale. Consequently, a probability measure P , which is equivalent to P on F T , is defined by setting
d P d P | F T : = Λ T .
Then, the following lemma provides the conditional distribution of Y t + 1 when given F t under the probability measure P .
Lemma 1. 
Under the probability measure P , the conditional distribution of Y t + 1 , given that F t is a normal distribution with mean r 1 2 h t + 1 and variance h t + 1 . That is, Y t + 1 | F t N ( r 1 2 h t + 1 , h t + 1 ) .
Proof. 
The proof follows some results or arguments in the work of Elliott et al. (1995) (Chapter 3 therein). Using a version of the Bayes’ rule, the martingale property of { Λ t } t T , Equation (4), the law of iterated expectations, and the fact that Λ t is F t -measurable,
P ( Y t + 1 y | F t ) = E [ I { Y t + 1 y } | F t ] = E [ Λ T I { Y t + 1 y } | F t ] E [ Λ T | F t ] = E [ Λ T I { Y t + 1 y } | F t ] Λ t = E { E [ Λ T I { Y t + 1 y } | F t + 1 ] | F t } Λ t = E { I { Y t + 1 y } E [ Λ T | F t + 1 ] | F t } Λ t = E [ ( Λ t ) 1 Λ t + 1 I { Y t + 1 y } | F t ] = E [ λ t + 1 I { Y t + 1 y } | F t ] .
Here, E [ · | F t ] is the conditional expectation given F t under P ; E [ · | F t ] is the conditional expectation given F t under P . Recall that { ϵ t } t T i . i . d . N ( 0 , 1 ) under the real-world reference probability measure P . This, together with Equation (1), imply that Y t + 1 | F t N ( μ Y , h t + 1 ) . Then, by Equations (3) and (6),
P ( Y t + 1 y | F t ) = E [ λ t + 1 I { Y t + 1 y } | F t ] = y ϕ μ Y ( t + 1 ) , h t + 1 ( z ) ϕ μ Y , h t + 1 ( z ) ϕ μ Y , h t + 1 ( z ) d z = y ϕ μ Y ( t + 1 ) , h t + 1 ( z ) d z .
Recall that ϕ μ Y ( t + 1 ) , h t + 1 ( · ) is the pdf of a normal distribution with mean μ Y ( t + 1 ) and variance h t + 1 . Consequently, by Equation (7), Y t + 1 | F t N ( μ Y ( t + 1 ) , h t + 1 ) . Note from Equation (2) that μ Y ( t + 1 ) = r 1 2 h t + 1 . Hence, the result follows. □
By the fundamental theorem of asset pricing, the absence of arbitrage opportunities in a financial market is “essentially” equivalent to the existence of an equivalent martingale measure under which the discounted price process of an asset relative to a numeraire asset is a martingale (see, for example, Delbaen and Schachermayer (1994); Hansen and Richard (1987); Harrison and Kreps (1979); Harrison and Pliska (1981, 1983)). Let { B t } t T denote the price process of a risk-free (zero-coupon) bond with the constant continuously compounded interest rate r > 0 . Then, for each t T ,
B t = e r t .
The bond is the numeraire asset used for discounting.
Let { S t } t T denote the price process of Ethereum in AUD, where S t is the adjusted close price of Ethereum in AUD in the tth period (i.e., the tth trading day). Let S 0 denote the initial price at the beginning of the first period. Then, S t = S t 1 exp ( Y t ) , for each t T . The following theorem shows that the discounted price process { B t 1 S t } t T of Ethereum in AUD by the bond is an ( F , P ) -martingale.
Theorem 1. 
{ B t 1 S t } t T is an ( F , P ) -martingale.
Proof. 
It suffices to prove that
E [ B t + 1 1 S t + 1 | F t ] = B t 1 S t .
The result is then followed by induction.
From Lemma 1, the moment-generating function of Y t + 1 given F t under P evaluated at u is
E [ e u Y t + 1 | F t ] = exp r 1 2 h t + 1 u + 1 2 u 2 h t + 1 .
Setting u = 1 in Equation (10) gives
E [ e Y t + 1 | F t ] = e r .
From Equations (8), (9) and (11),
E [ B t + 1 1 S t + 1 | F t ] = E [ e r ( t + 1 ) S t e Y t + 1 | F t ] = e r ( t + 1 ) S t E [ e Y t + 1 | F t ] = e r t S t = B t 1 S t .
This completes the proof. □
From Theorem 1, P is an equivalent martingale measure, which is a risk-neutral probability measure selected by the discrete-time Giransov’s transform.
Let { ϵ t } t T denote a sequence of i.i.d. standard normal random variables under the probability measure P . Then, from Lemma 1, under P ,
Y t + 1 = r 1 2 h t + 1 + h t + 1 ϵ t + 1 .
Write, for each t T ,
η t : = r μ Y 1 2 h t .
Then, by comparing the real-world and risk-neutral dynamics for { Y t } t T in Equation (1) and Equation (13), respectively,
ϵ t = η t h t + ϵ t .
Using Equation (15) and the conditional variance process in Equation (1),
h t = ω + α h t 1 η t 1 h t 1 + ϵ t 1 2 + β h t 1 .
Putting together Equations (13) and (16), under the risk-neutral probability measure P , the logarithmic return process { Y t } t T has the following GARCH dynamics:
Y t = r 1 2 h t + h t ϵ t , h t = ω + α h t 1 η t 1 h t 1 + ϵ t 1 2 + β h t 1 .
Note that the market described by the discrete-time GARCH model is incomplete. Consequently, there are many equivalent martingale measures under the discrete-time GARCH model. The risk-neutral probability measure P selected by the discrete-time Girsanov’s transform is one of them.
There are other approaches to specify an equivalent martingale measure under a discrete-time GARCH model and their variants. Duan (1995) introduced the locally risk-neutral valuation relationship for option pricing under a discrete-time GARCH model. Siu et al. (2004a) adopted the conditional Esscher transform to select an equivalent martingale measure under discrete-time GARCH models. Elliott et al. (2006) and Siu et al. (2006) discussed the use of the conditional Esscher transform or its variants for pricing European-style options under a discrete-time Markov-switching GARCH model and a discrete-time SETAR-ARCH model, respectively. The self-exciting threshold autoregressive (SETAR) model pioneered by (Tong 1977, 1983, 1990) is one of the major classes of parametric nonlinear time series models. Siu and Elliott (2021) studied the valuation of Bitcoin options under a discrete-time SETAR-GARCH model using the conditional Esscher transform and the variance-dependent pricing kernel.

3. Model Uncertainty via Nonlinear Expectations

In this section, model uncertainty about the conditional mean and volatility of the logarithmic returns of Ethereum is introduced using an extended Girsanov’s principle under the GARCH model in Section 2. Nonlinear expectations are then used to construct lower and upper estimates for the price of an Ether option in light of the model uncertainty. To this end, the following three steps are considered:
(I)
Using an extended Girsanov’s principle, as in (Elliott 2017; Elliott and Siu 2017; Siu 2023a), a family of alternative GARCH models is introduced by perturbing the conditional mean and volatility of the reference GARCH model in Equation (1).
(II)
As in Section 2, a “risk-neutral” extended Girsanov’s principle is adopted to transform the family of alternative GARCH models into a family of “risk-neutralized” alternative GARCH models.
(III)
A pair of conditional nonlinear expectations, namely a conditional superlinear expectation and a conditional sublinear expectation, is constructed based on the family of “risk-neutralized” GARCH alternative models. The conditional superlinear and sublinear expectations then give the lower and upper estimates for the price of an Ether option, respectively.

3.1. A Family of Alternative Models

Step (I) is now considered. Cryptocurrency markets are highly volatile and uncertain. Modelling price fluctuations of cryptocurrencies could be a formidable task. Consequently, model uncertainty is an important and relevant issue as far as the modelling of price fluctuations of cryptocurrencies is concerned. In practice, the mean (or drift) and the volatility of the logarithmic returns are two quantities which interest market practitioners. In the context of a discrete-time GARCH model, the two quantities of interest are the conditional mean and volatility. Accurate modelling or estimation of the conditional mean and volatility of the logarithmic returns of cryptocurrencies could be difficult in practice. Consequently, from a practical perspective, a relevant issue could be how might the model uncertainty, or ambiguity, of the conditional mean and volatility be incorporated.
When faced with model uncertainty, or ambiguity, about the conditional mean and volatility, market practitioners may form interval estimates for misspecifications in the conditional mean and volatility. These interval estimates give the ranges of the “uncertain” parameters, describing the misspecifications in the conditional mean and volatility. Specifically, as in Siu (2023a), a two-dimensional vector of “uncertain” parameters, denoted by θ , describing the misspecifications in the conditional mean and volatility under the GARCH model in Equation (1), is defined as follows:
θ : = ( μ , σ ) × + .
Indeed, μ and σ describe perturbations in the conditional mean and volatility of the logarithmic return process { Y t } t T under the reference GARCH model in Equation (1). μ is the “uncertain” parameter describing the misspecification in the conditional mean (or drift), and σ is the “uncertain” parameter describing the misspecification in the conditional volatility. Specifically, if μ = 0 and σ = 1 , there will be no perturbations/misspecification in the conditional mean and volatility of the logarithmic return process { Y t } t T . y denotes the transpose of a vector, or a matrix, y .
For practical purposes, the interval estimates for misspecifications in the conditional mean and volatility should be finite intervals. Consequently, it may not be unreasonable to assume that the perturbations in the conditional mean and volatility do not take infinite values. This assumption may also be in line with the specifications of the g-expectation for the drift uncertainty in Peng (1997) and the G-expectation for the volatility uncertainty in Peng (2006, 2019). Specifically, the vector of the “uncertain” parameters θ in Equation (18) takes on values in the following rectangular region Θ in × + :
Θ : = [ μ L , μ U ] × [ σ L , σ U ] ,
for μ L , μ U and σ L , σ U + such that μ L < μ U and σ L < σ U . Specifically, μ L and μ U are, respectively, the lower and upper limits of the interval estimate for the “uncertain” parameter μ , while σ L and σ U are, respectively, the lower and upper limits of the interval estimate for the “uncertain” parameter σ .
To capture volatility uncertainty in continuous time, Avellaneda et al. (1995) and Lyons (1995) introduced an “uncertain” volatility parameter and supposed that the “uncertain” volatility parameter lies in an interval whose lower and upper limits could be specified based on the extreme values of the implied volatilities of liquid securities or the high–low peaks of historical volatilities. To incorporate both drift and volatility uncertainties in continuous time, Wilmott and Oztukel (1998) and Wilmott (2007) introduced both the “uncertain” drift and volatility parameters and supposed that they lie in the respective intervals whose lower and upper limits are pre-specified. They called the intervals for the “uncertain” parameters “certainty bands”. Agarwal et al. (2021) incorporated the model uncertainty about the market price of risk in a continuous-time model with multiple risk factors and adopted an ellipsoid-shaped set instead of a rectangular set to specify a statistical confidence region for the “uncertain” market price of risk. The ambiguity specified by an ellipsoid-shaped set is called an elliptical ambiguity, while the ambiguity specified by a rectangular set is called a rectangular ambiguity. The elliptical ambiguity was also considered in Cohen and Tegnér (2019). Siu (2023b) applied some asymptotic results of the maximum likelihood estimates for transition probabilities for the Markov chain to specify the lower and upper limits for the perturbations/misspecifications in the transition rates of the chain.
The aforementioned works considered the incorporation of model uncertainty in continuous time, whereas the incorporation of model uncertainty in discrete time is considered here. Indeed, as noted in Elliott and Siu (2017) and Siu (2023a), besides some advantages from an econometric perspective, in a discrete-time modelling framework, volatility uncertainty can be incorporated via the use of an extended Girsanov’s principle for changing probability measures. Unlike the continuous-time situation where singular probability measures are involved (e.g., Peng (1997) and Fouque and Ren (2014)), the use of an extended Girsanov’s principle involves absolutely continuous changes in probability measures, and it provides a streamlined approach to incorporate volatility uncertainty. Furthermore, in (Agarwal et al. 2021; Avellaneda et al. 1995; Lyons 1995; Wilmott 2007; Wilmott and Oztukel 1998), “uncertain” parameters were introduced to directly describe the ambiguity of the model parameters, such as the drift, volatility, and market price of risk, whereas, as in Siu (2023a), the “uncertain” parameters are introduced here to describe the perturbations/misspecifications in the conditional mean and volatility. Lastly, in some of the existing literature, confidence limits are specified by market information, such as the extreme values of implied volatilities, pre-specified or specified by classical statistical approaches. However, the confidence limits are to be specified here using Bayesian credible intervals along the lines of Boyle and Ananthanarayanan (1977) and Siu (2023a). This will be discussed in more detail in Section 4.
For each θ Θ , let { λ t ( θ ) } t T denote an F -adapted process on ( Ω , F , P ) defined by putting
λ t ( θ ) : = ϕ μ , σ 2 ( ϵ t ) ϕ ( ϵ t ) = ϕ ( ϵ t μ σ ) σ ϕ ( ϵ t ) ,
where { ϵ t } t T i . i . d . N ( 0 , 1 ) under the reference probability measure P as defined in the last section (Section 2).
For each θ Θ , the F -adapted process { Λ t ( θ ) } t T is defined as follows:
Λ t ( θ ) : = k = 1 t λ k ( θ ) .
By definition, for each θ Θ , { Λ t ( θ ) } t T is a positive-valued process. It is also an ( F , P ) -martingale. Consequently, it can be used as a density process for changing probability measures in discrete time. Then, for each θ Θ , a probability measure P θ equivalent to P on F T is defined by putting
d P θ d P | F T : = Λ T ( θ ) .
By a discrete-time Girsanov’s theorem (see, for example, Elliott et al. (1995) and Elliott and Madan (1998)), for each θ Θ , { ϵ t ( θ ) } t T is defined by putting:
ϵ t ( θ ) : = ϵ t μ σ ,
which is a sequence of i.i.d. standard normal random variables under the newly defined probability measure P θ in Equation (22).
The family { ϵ ( θ ) | θ Θ } indexed by Θ , where ϵ ( θ ) : = { ϵ t ( θ ) } t T , defines a family of sequences of Gaussian “uncertain” noises. From Equations (17) and (23), for each θ Θ , under P θ , the logarithmic return process { Y t } t T is governed by
Y t = μ Y + μ h t + σ h t ϵ t ( θ ) , h t = ω + α h t 1 ( μ + σ ϵ t 1 ( θ ) ) 2 + β h t 1 .
When θ : = ( μ , σ ) varies in Θ , a family of alternative GARCH models having the form in Equation (24) is defined under a family of probability measures { P θ | θ Θ } . Borrowing some terminologies from relativity in physics, each of the alternative GARCH models may be considered a reference frame to describe the dynamics of the logarithmic returns. The reference GARCH model in Equation (1) may also be thought of as a reference frame as it is included in the family of alternative GARCH models. In other words, the dynamics of the logarithmic returns are described by a family of reference frames. They may be thought of as “equivalent” reference frames in the sense that they are described by a family of equivalent probability measures.

3.2. Risk-Neutralization

In this subsection, Step (II) is considered. For each θ Θ and each t T , let
μ Y θ ( t ) : = μ Y + μ h t .
Then, an F -adapted process { λ t ( θ ) } t T is defined by putting
λ t ( θ ) : = ϕ μ Y ( t ) , σ 2 h t ( Y t ) ϕ μ Y θ ( t ) , σ 2 h t ( Y t ) = ϕ ( Y t μ Y ( t ) σ h t ) ϕ ( Y t μ Y θ ( t ) σ h t ) ,
where μ Y ( t ) and μ Y θ ( t ) are given by Equations (2) and (25), respectively.
From Equation (26), for each θ Θ , an F -adapted process { Λ t ( θ ) } t T is defined by putting
Λ t ( θ ) : = k = 1 t λ k ( θ ) .
Note that for each θ Θ , { Λ t ( θ ) } t T is a positive-valued ( F , P θ ) -martingale. Consequently, for each θ Θ , a new probability measure P ( θ ) , which is equivalent to P θ on F T , is defined by setting
d P ( θ ) d P θ | F T : = Λ T ( θ ) .
Using arguments similar to the proof of Lemma 1, it can be shown that for each θ Θ , under the probability measure P ( θ ) , the conditional distribution of Y t + 1 with F t given is a normal distribution with mean μ Y ( t + 1 ) : = r 1 2 h t + 1 and variance σ h t + 1 . Similarly to Equation (13), for each θ Θ , under P ( θ ) ,
Y t + 1 = r 1 2 h t + 1 + σ h t + 1 ϵ t + 1 ( θ ) ,
where { ϵ t ( θ ) } t T i . i . d . N ( 0 , 1 ) under P ( θ ) .
For each θ Θ and each t T , let
η t ( θ ) : = r μ Y 1 2 h t μ h t .
Comparing the dynamics for { Y t } t T in Equation (24) and those in Equation (29) gives
ϵ t ( θ ) = η t ( θ ) σ h t + ϵ t ( θ ) .
Substituting Equation (31) into the conditional variance process in Equation (24) gives
h t = ω + α h t 1 μ + η t 1 ( θ ) h t 1 + σ ϵ t 1 ( θ ) 2 + β h t 1 .
Putting together Equations (29) and (32), for each θ Θ , under P ( θ ) , the logarithmic return process { Y t } t T has the following dynamics:
Y t = r 1 2 h t + σ h t ϵ t ( θ ) , h t = ω + α h t 1 μ + η t 1 ( θ ) h t 1 + σ ϵ t 1 ( θ ) 2 + β h t 1 .
Using arguments similar to the proof of Theorem 1 and Equation (33), it can be shown that for each θ Θ , the discounted price process of ETH-AUD, i.e., { B t 1 S t } t T , is an ( F , P ( θ ) ) -martingale. Consequently, for each θ Θ , P ( θ ) is a “risk-neutral” probability measure. It defines a “risk-neutralized” alternative GARCH model in Equation (33), which incorporates the impact of model uncertainty about the conditional mean and volatility via the vector of the “uncertain” parameters θ : = ( μ , σ ) Θ . Then, the family of “risk-neutral” probability measures { P ( θ ) | θ Θ } defines a family of “risk-neutralized” alternative GARCH models, each having the form in Equation (33).

3.3. Lower and Upper Estimates for Ether Option Prices

Step (III) is considered in this subsection. A European-style Ether option with maturity at time T and payoff described by a random variable X on ( Ω , F , P ) is considered. Assume that the random variable X is F T Y -measurable and integrable under the probability measure P . That is, E [ | X | ] < .
For each θ Θ , let E θ , [ · | F t ] denote the conditional expectation given F t under a “risk-neutral” probability measure P ( θ ) . Then, an F -adapted process { X t θ } t T is defined by putting
X t θ B t : = E θ , X B T | F t .
From Equation (34), for each θ Θ , { X t θ B t } t T is an ( F , P ( θ ) ) -martingale. For a given θ Θ , if { X t θ } t T is taken as a price process of the Ether option, and then the discounted price process of the Ether option { X t θ B t } t T is an ( F , P ( θ ) ) -martingale, where the numeraire asset for the discounting is chosen to be the bond.
The conditional superlinear and sublinear expectations for the discounted payoff of the Ether option X B T with respect to the family of “risk-neutral” probability measures { P ( θ ) | θ Θ } are now defined. Specifically, for each t T , given F t , the conditional superlinear expectation for X B T with respect to { P ( θ ) | θ Θ } , denoted by E Θ , L [ B T 1 X | F t ] , is defined as follows:
E Θ , L [ B T 1 X | F t ] : = ess inf θ Θ E θ , [ B T 1 X | F t ] ,
where the essential infimum ess inf is taken over θ Θ .
Similarly, for each t T , given F t , the conditional sublinear expectation for X B T with respect to { P ( θ ) | θ Θ } , denoted by E Θ , U [ B T 1 X | F t ] , is defined as follows:
E Θ , U [ B T 1 X | F t ] : = ess sup θ Θ E θ , [ B T 1 X | F t ] ,
where the essential supremum ess sup is taken over θ Θ .
Note that the conditional sublinear expectation in Equation (36) satisfies the following five properties:
(i)
Measurability: For any integrable random variable X F t (i.e., X is an integrable, F t -measurable random variable),
E Θ , U [ X | F t ] = X .
(ii)
Monotonicity: For any two integrable random variables X , Z F T , satisfying X : = X ( ω ) Z : = Z ( ω ) for all ω Ω ,
E Θ , U [ X | F t ] E Θ , U [ Z | F t ] .
(iii)
Translation Invariance: For any two integrable random variables X F T and C F t ,
E Θ , U [ X + C | F t ] = E Θ , U [ X | F t ] + C .
(iv)
Positive Homogeneity: For any two integrable random variables X F T and Λ F t such that Λ : = Λ ( ω ) 0 , for all ω Ω ,
E Θ , U [ Λ X | F t ] = Λ E Θ , U [ X | F t ] .
(v)
Sub-additivity: For any two integrable random variables X , Z F T ,
E Θ , U [ X + Z | F t ] E Θ , U [ X | F t ] + E Θ , U [ Z | F t ] .
The conditional superlinear expectation in Equation (35) satisfies Properties (i)–(iv) above. However, instead of satisfying the sub-additivity property in Equation (41), the conditional superlinear expectation in Equation (35) satisfies the following super-additivity property:
(v*)
Super-additivity: For any two integrable random variables X , Z F T ,
E Θ , L [ X + Z | F t ] E Θ , L [ X | F t ] + E Θ , L [ Z | F t ] .
For each t T , let L 1 ( Ω , F t , P ) be the space of integrable random variables defined on ( Ω , F t , P ) . Since { E θ , [ B T 1 X | F t ] | θ Θ } L 1 ( Ω , F t , P ) , and L 1 ( Ω , F t , P ) is a complete lattice, both E Θ , U [ B T 1 X | F t ] and E Θ , L [ B T 1 X | F t ] exist, and they are F t -measurable random variables (see Dunford and Schwartz (1958), p. 302 therein, and Elliott (1982), p. 225 and Lemma 16.A.2 therein).
Using the conditional superlinear expectation in Equation (35), the lower estimate for the price of the Ether option at time t, denoted by X t L ( Θ ) , with respect to the family of “risk-neutral” probability measures { P ( θ ) | θ Θ } is defined by
X t L ( Θ ) : = B t E Θ , L X B T | F t = e r ( T t ) E Θ , L [ X | F t ] .
Similarly, using the conditional sublinear expectation in Equation (36), the upper estimate for the price of the Ether option at time t, denoted by X t U ( Θ ) , with respect to the family of “risk-neutral” probability measures { P ( θ ) | θ Θ } is defined by
X t U ( Θ ) : = B t E Θ , U X B T | F t = e r ( T t ) E Θ , U [ X | F t ] .

4. Bayesian Lower and Upper Estimates

In this section, the Bayesian lower and upper estimates for the price of an Ether option under the GARCH model are constructed using Bayesian nonlinear expectations in Siu (2023a). The basic idea is to adopt the “confidence” limits of Bayesian credible intervals to specify the lower and upper limits for each of the “uncertain” parameters μ and σ describing the perturbations/misspecifications in the conditional mean and volatility, respectively. Specifically, the “confidence” limits of the Bayesian credible intervals are used to specify the lower and upper limits of the rectangular region Θ in Equation (19).
Bayesian statistics provides a data analytical approach, which is supplemented by expert opinion or subjective judgement, to specify the conditional nonlinear expectations in Equations (35) and (36). The two sources of information are combined “rationally” using the Bayes’ theorem. Indeed, Bayesian statistics contributes to many developments of machine learning, artificial intelligence and data science. Its relevance to model uncertainty, or ambiguity, was featured in the smooth ambiguity approach to model uncertainty proposed by Klibanoff et al. (2005). Its connection with decision theory was highlighted in the classic monograph by Bernardo and Smith (2000). See also Siu (2023a) for some related discussions. Bayesian statistics and computation can also be used to handle missing data, which may be particularly relevant as far as big data are concerned. See, for example, (Kong et al. 1994; Ramoni and Sebastiani 2001; Rubin 1976; Tanner and Wong 1987). The constructions and notation below follow those in Siu (2023a), Section 4 therein, which were based on results in Bernardo and Smith (2000).
From Equation (23),
ϵ t = μ + σ ϵ t ( θ ) .
Recall that for each θ Θ , under P θ , ϵ 1 ( θ ) , ϵ 2 ( θ ) , , ϵ T ( θ ) i . i . d . N ( 0 , 1 ) . Consequently, by Equation (45), ϵ 1 , ϵ 2 , , ϵ T | θ i . i . d . N ( μ , σ 2 ) under P θ . The above statements still hold if the two-dimensional vector of the “uncertain” parameters θ : = ( μ , σ ) ( , ) × ( 0 , ) . In Bayesian statistics, the vector of the “uncertain” parameters ( μ , σ ) is regarded as a random vector. When evaluating the Bayesian credible intervals for the random vector ( μ , σ ) , it is assumed that the state space of the random vector ( μ , σ ) is ( , ) × ( 0 , ) . Following Bernardo and Smith (2000), Siu (2023a) and some relevant references therein, instead of considering σ , the precision λ is considered, where λ : = 1 σ 2 . In other words, the random vector ( μ , λ ) , with the same state space ( , ) × ( 0 , ) , is considered. Since σ is an “uncertain” parameter, λ may also be interpreted as an “uncertain” parameter.
A Normal–Gamma conjugate prior distribution is assigned for the random vector ( μ , λ ) . That is,
μ | λ N μ 0 , 1 t 0 λ , λ G a ( α 0 , β 0 ) ,
where μ 0 and t 0 λ are the prior mean and precision of the normal prior distribution for μ | λ , respectively. Note that μ 0 ( , ) and t 0 ( 0 , ) . The prior distribution of λ is a Gamma distribution with the shape parameter α 0 and the rate parameter β 0 . Note also that α 0 ( 0 , ) and β 0 ( 0 , ) . In practice, those prior parameters, which are also called hyper-parameters, can be chosen based on expert opinion or subjective judgement. See Siu (2023a), Section 4, for discussions about the intuitions behind and interpretations for those prior parameters.
Suppose that the reference model under the reference (real-world) probability measure P (e.g., the GARCH model in Equation (1)) is estimated, given F n . That is, the GARCH model is estimated using observations about the logarithmic returns { Y 1 , Y 2 , , Y n } . The estimation for the GARCH model based on the conditional maximum likelihood estimation method is performed using the R package “rugarch” Galanos and Kley (2023). Let I : = { 1 , 2 , , n } , which is a set of indices of the observations. Then, from Equation (1), the standardized residuals { e i } i I from the model fitting are computed as follows:
e i = Y i μ ^ Y h ^ i ,
where μ ^ Y is the estimate for the mean parameter μ Y , and h ^ i is the estimate for the conditional variance h i such that
h ^ i = ω ^ + α ^ h ^ i 1 e i 1 2 + β ^ h ^ i 1 .
Note that ω ^ , α ^ and β ^ are the conditional maximum likelihood estimates for the GARCH parameters ω , α and β , respectively.
Recall that under P θ , ϵ 1 , ϵ 2 , , ϵ n | θ i . i . d . N ( μ , σ 2 ) . It may not be unreasonable to assume that under P θ , e 1 , e 2 , , e n | θ i . i . d . N ( μ , σ 2 ) . Let e ( n ) : = { e 1 , e 2 , , e n } . Write γ 1 and γ 2 for two probability levels such that γ 1 , γ 2 ( 0 , 1 ) . Let t 2 α 0 + n ( γ 1 2 ) denote the upper-tail critical value of a t-distribution with the degree of freedom 2 α 0 + n corresponding to the probability level γ 1 2 . Then, using results in Bernardo and Smith (2000), p. 440 (see also Siu (2023a), Section 4), the 100 ( 1 γ 1 ) % Bayesian credible interval for the “uncertain” parameter μ given e ( n ) is ( L γ 1 ( μ | e ( n ) ) , U γ 1 ( μ | e ( n ) ) ) such that
L γ 1 ( μ | e ( n ) ) = μ n t 2 α 0 + n γ 1 2 σ n n + t 0 , U γ 1 ( μ | e ( n ) ) = μ n + t 2 α 0 + n γ 1 2 σ n n + t 0 ,
where
μ n = t 0 μ 0 + n e ¯ n t 0 + n , σ n = β n α 0 + n 2 , β n = β 0 + 1 2 n s n 2 + t 0 n ( μ 0 e ¯ n ) 2 2 ( t 0 + n ) , e ¯ n = 1 n i = 1 n e i , s n 2 = 1 n i = 1 n ( e i e ¯ n ) 2 .
Using some terminologies in machine learning, the sample mean e ¯ n may be related to the training error, while the unknown mean μ may be related to the test error that is to be learned or estimated from the training error (see, for example, Zhang (2023), page 9 therein).
Furthermore, the 100 ( 1 γ 2 ) % Bayesian credible interval for the “uncertain” parameter λ given e ( n ) is ( L γ 2 ( λ | e ( n ) ) , U γ 2 ( λ | e ( n ) ) ) such that
L γ 2 ( λ | e ( n ) ) = Γ α 0 + n 2 , β n 1 γ 2 2 , U γ 2 ( λ | e ( n ) ) = Γ α 0 + n 2 , β n γ 2 2 ,
where Γ α 0 + n 2 , β n ( γ ) is the upper-tail critical value of a Gamma distribution with the shape parameter α 0 + n 2 and the rate parameter β n corresponding to the probability level γ .
Take
L γ 2 ( σ | e ( n ) ) : = 1 U γ 2 ( λ | e ( n ) ) , U γ 2 ( σ | e ( n ) ) : = 1 L γ 2 ( λ | e ( n ) ) .
Then, the following rectangular region is defined:
Θ γ 1 , γ 2 ( e ( n ) ) : = [ L γ 1 ( μ | e ( n ) ) , U γ 1 ( μ | e ( n ) ) ] × [ L γ 2 ( σ | e ( n ) ) , U γ 2 ( σ | e ( n ) ) ] .
Note that Θ γ 1 , γ 2 ( e ( n ) ) may be called a Bayesian (posterior) rectangular region. It may be thought of as a Bayesian (posterior) estimate for the rectangular region in Equation (19). To simplify the notation, we write Θ ( e ( n ) ) for Θ γ 1 , γ 2 ( e ( n ) ) .
Let N T denote the set of non-negative integers truncated at the terminal time T < . That is, N T : = { 0 , 1 , 2 , , T } . Define the product space Θ ( e ) by
Θ ( e ) : = n = 1 N T Θ ( e ( n ) ) .
We consider an { F n } n N T -adapted stochastic process, denoted by θ ( e ) : = { θ n ( e ) } n N T , such that for each n N T and each ω Ω , θ n ( e ) : = θ n ( e , ω ) Θ ( e ( n ) ) .
Define, for each n N T , θ n ( e ) : = ( μ n ( e ) , σ n ( e ) ) Θ ( e ( n ) ) . Note that for each n N T , μ n ( e ) is an F n -measurable random variable such that for each ω Ω , μ n ( e ) : = μ n ( e , ω ) [ L γ 1 ( μ | e ( n ) ) , U γ 1 ( μ | e ( n ) ) ] . Similarly, for each n N T , σ n ( e ) is an F n -measurable random variable such that for each ω Ω , σ n ( e ) : = σ n ( e , ω ) [ L γ 2 ( σ | e ( n ) ) , U γ 2 ( σ | e ( n ) ) ] .
For each process θ ( e ) : = { θ n ( e ) } n N T Θ ( e ) , each n N T and each k = n + 1 , n + 2 , , T , let
μ Y ( k , θ n ( e ) ) : = μ Y + μ n ( e ) h k .
Then, for each θ ( e ) Θ ( e ) , each n N T and each k = n + 1 , n + 2 , , T , we define
λ k ( θ n ( e ) ) : = ϕ μ Y ( k ) , ( σ n ( e ) ) 2 h k ( Y k ) ϕ μ Y ( k , θ n ( e ) ) , ( σ n ( e ) ) 2 h k ( Y k ) × ϕ μ n ( e ) , ( σ n ( e ) ) 2 ( ϵ k ) ϕ ( ϵ k ) = ϕ ( Y k μ Y ( k ) σ n ( e ) h k ) ϕ ( Y k μ Y ( k , θ n ( e ) ) σ n ( e ) h k ) × ϕ ( ϵ k μ n ( e ) σ n ( e ) ) σ n ( e ) ϕ ( ϵ k ) ,
where μ Y ( k ) and μ Y ( k , θ n ( e ) ) are given by Equation (2) and Equation (55), respectively. Furthermore, we take
λ n ( θ n ( e ) ) = 1 .
For each n N T , let F n denote the “tail” filtration { F n , k } k = n , n + 1 , , T generated by the logarithmic returns { Y n , Y n + 1 , , Y T } such that for each k = n , n + 1 , , T ,
F n , k : = σ { Y m | m = n , n + 1 , , k } N ,
where A 1 A 2 : = σ ( A 1 A 2 ) for any two σ -algebras A 1 and A 2 on Ω ; N is a collection of P -null subsets in the σ -field F .
By the definitions in Equations (56)–(58), for each n N T , the “tail” process { λ k ( θ n ( e ) ) } k = n , n + 1 , , T is an F n -adapted stochastic process.
For each θ ( e ) Θ ( e ) and each n N T , using Equation (56), the F n -adapted process { Λ k ( θ n ( e ) ) } k = n , n + 1 , , T is defined such that for each k = n , n + 1 , , T ,
Λ k ( θ n ( e ) ) : = m = n k λ m ( θ n ( e ) ) .
Then, by the definitions in Equations (56), (57) and (59), for each n N T , { Λ k ( θ n ( e ) ) } k = n , n + 1 , , T is an ( F n , P ) -martingale.
Consequently, for each θ ( e ) Θ ( e ) and each n N T , a new probability measure P ( θ n ( e ) ) , which is equivalent to P on F n , T , is defined by setting
d P ( θ n ( e ) ) d P | F n , T : = Λ T ( θ n ( e ) ) .
For each n N T , a family of Bayesian (posterior) “risk-neutral” probability measures, given F n , which is indexed by the space of processes Θ ( e ) , is defined by
P n ( Θ ( e ) ) : = { P ( θ n ( e ) ) | θ ( e ) Θ ( e ) } .
The family P n ( Θ ( e ) ) defines a family of probability measures underlying a family of Bayesian (posterior) “risk-neutral” alternative GARCH models given F n . Note that for each n N T , only the n t h coordinate, i.e., θ n ( e ) , of the process θ ( e ) Θ ( e ) is involved in defining the probability measure P ( θ n ( e ) ) corresponding to θ ( e ) within the family P n ( Θ ( e ) ) .
For each n N T , a Bayesian (posterior) superlinear expectation for the discounted payoff of the Ether option X B T given F n with respect to P n ( Θ ( e ) ) in Equation (61), denoted by E Θ ( e ) , L [ B T 1 X | F n ] , is defined as follows:
E Θ ( e ) , L [ B T 1 X | F n ] : = ess inf θ ( e ) Θ ( e ) E θ n ( e ) , [ B T 1 X | F n ] ,
where E θ n ( e ) , [ · | F n ] is the conditional expectation given F n under the probability measure P ( θ n ( e ) ) P n ( Θ ( e ) ) . The essential infimum in Equation (62) is taken over all the probability measures in the family P n ( Θ ( e ) ) . Note that the Bayesian (posterior) superlinear expectation in Equation (62) satisfies the five properties (Properties (i)–(iv) in Equations (37)–(40) and Property (v*) in Equation (42)).
Similarly, for each n N T , a Bayesian (posterior) sublinear expectation for the discounted payoff of the Ether option X B T given F n with respect to P n ( Θ ( e ) ) in Equation (61), denoted by E Θ ( e ) , U [ B T 1 X | F n ] , is defined as follows:
E Θ ( e ) , U [ B T 1 X | F n ] : = ess sup θ ( e ) Θ ( e ) E θ n ( e ) , [ B T 1 X | F n ] .
The Bayesian (posterior) sublinear expectation in Equation (63) satisfies the five properties (Properties (i)–(v) in Equations (37)–(41)).
Using the Bayesian (posterior) superlinear expectation in Equation (62), the Bayesian (posterior) lower estimate for the price of the Ether option at time n with respect to P n ( Θ ( e ) ) , denoted by X n L ( Θ ( e ) ) , is defined by
X n L ( Θ ( e ) ) : = B n E Θ ( e ) , L X B T | F n = e r ( T n ) E Θ ( e ) , L [ X | F n ] .
Similarly, using the Bayesian (posterior) sublinear expectation in Equation (63), the Bayesian (posterior) upper estimate for the price of the Ether option at time n with respect to P n ( Θ ( e ) ) , denoted by X n U ( Θ ( e ) ) , is defined by
X n U ( Θ ( e ) ) : = B n E Θ ( e ) , U X B T | F n = e r ( T n ) E Θ ( e ) , U [ X | F n ] .
The Bayesian (posterior) lower and upper estimates for the price of the Ether option are interpreted as the Bayesian (posterior) lower and upper confidence limits on the option price, respectively, in the line with the interpretations adopted in Boyle and Ananthanarayanan (1977).

5. Numerical/Simulation Procedure

In this section, the numerical/simulation procedure for computing the Bayesian lower and upper estimates for the price of the Ether option is described. Given e ( n ) (i.e., the standardized residuals { e i } i I ), the Bayesian (posterior) rectangular region Θ ( e ( n ) ) : = Θ γ 1 , γ 2 ( e ( n ) ) can be computed using Equations (49)–(53). Then, the Bayesian lower and upper estimates for the price of the Ether option in Equations (64) and (65) are computed using the procedures described in Section 4. However, it does not seem that closed-form expressions for the Bayesian lower and upper estimates in Equations (64) and (65) are available. Consequently, they are computed by the numerical/simulation procedures to be described in the sequel.
To compute the Bayesian lower and upper estimates, a numerical/simulation procedure based on the Monte Carlo simulation is adopted. A European call Ether option with strike price K and maturity at time T is considered. In this case, the terminal payoff of the Ether option at time T is
X = max ( S T K , 0 ) ,
where S T is the price of Ethereum in AUD at the maturity T.
To approximate the Bayesian (posterior) superlinear expectation E Θ ( e ) , L [ X | F n ] , it is conjectured that the infirmum of E θ , [ X | F n ] over θ Θ ( e ( n ) ) is attained at the extremal point ( μ , σ ) = ( L γ 1 ( μ | e ( n ) ) , L γ 2 ( σ | e ( n ) ) ) . Similarly, to approximate the Bayesian (posterior) sublinear expectation E Θ ( e ) , U [ X | F n ] , it is conjectured that the supremum of E θ , [ X | F n ] over θ Θ ( e ( n ) ) is attained at the extremal point ( μ , σ ) = ( U γ 1 ( μ | e ( n ) ) , U γ 2 ( σ | e ( n ) ) ) .
It is possible that the supremum and infimum may be attained at some interior points of the rectangular region Θ ( e ( n ) ) . However, the approximations to the Bayesian (posterior) superlinear and sublinear expectations based on the extremal points look simpler than using grid search. The latter could be computationally intensive when the number of grid points is large and it is used coupled with the Monte Carlo simulation method. Furthermore, the approximation results based on extremal, or boundary, points to be presented look reasonable at an intuitive level.
To sum up, when approximating the Bayesian (posterior) superlinear expectation E Θ ( e ) , L [ X | F n ] , the following extremal, or boundary, point is adopted:
θ ¯ ( L ) : = ( μ ¯ ( L ) , σ ¯ ( L ) ) = ( L γ 1 ( μ | e ( n ) ) , L γ 2 ( σ | e ( n ) ) ) .
Similarly, when approximating E Θ ( e ) , U [ X | F n ] , the following extremal, or boundary, point is used:
θ ¯ ( U ) : = ( μ ¯ ( U ) , σ ¯ ( U ) ) = ( U γ 1 ( μ | e ( n ) ) , U γ 2 ( σ | e ( n ) ) ) .
Consequently, the Bayesian (posterior) superlinear and sublinear expectations may be, respectively, approximated as follows:
E Θ ( e ) , L [ X | F n ] E θ ¯ ( L ) , [ X | F n ] ,
and
E Θ ( e ) , U [ X | F n ] E θ ¯ ( U ) , [ X | F n ] .
The Monte Carlo method is used to estimate the conditional expectations E θ ¯ ( L ) , [ X | F n ] and E θ ¯ ( U ) , [ X | F n ] in Equations (69) and (70), respectively. The use of the Monte Carlo method for option pricing was pioneered by Boyle (1977). See also Boyle et al. (1997) for an excellent account. Here, the Monte Carlo method with control variates is used to estimate E θ ¯ ( L ) , [ X | F n ] and E θ ¯ ( U ) , [ X | F n ] . To illustrate this, the estimation of E θ ¯ ( U ) , [ X | F n ] based on the Monte Carlo method with control variates is presented here. The estimation of E θ ¯ ( L ) , [ X | F n ] follows similarly. For more detail about the Monte Carlo method with control variates and other variance reduction techniques, such as antithetic variates, please refer to Boyle (1977) and Boyle et al. (1997).
Under the probability measure P ( θ ¯ ( U ) ) , the estimated dynamics for the “risk-neutralized” alternative GARCH in Equation (33) are given by
Y t = r 1 2 h ^ t + σ ¯ ( U ) h ^ t ϵ t ( θ ¯ ( U ) ) , h ^ t = ω ^ + α ^ h ^ t 1 μ ¯ ( U ) + η ^ t 1 ( θ ¯ ( U ) ) h ^ t 1 + σ ¯ ( U ) ϵ t 1 ( θ ¯ ( U ) ) 2 + β ^ h ^ t 1 ,
whereby, from Equation (30),
η ^ t 1 ( θ ¯ ( U ) ) : = r μ ^ Y 1 2 h ^ t 1 μ ¯ ( U ) h ^ t 1 .
Recall that μ ^ Y , ω ^ , α ^ and β ^ are the estimates for the unknown parameters in the GARCH model with drift in Equation (1) under the reference probability measure P . Under the probability measure P ( θ ¯ ( U ) ) , { ϵ t ( θ ¯ ( U ) ) } t T i . i . d . N ( 0 , 1 ) .
Suppose now that a set of simulated values for the terminal payoff X in Equation (66), that is, { X ( m ) ( θ ¯ ( U ) ) } m = 1 , 2 , , M , is generated by simulating from the estimated dynamics in Equation (71). Then, the Monte Carlo estimate for the conditional expectation E θ ¯ ( U ) , [ X | F n ] is given by
E θ ¯ ( U ) , [ X | F n ] 1 M m = 1 M X ( m ) ( θ ¯ ( U ) ) .
Using Equations (65), (70) and (73), without control variates, the Monte Carlo estimate for the Bayesian (posterior) upper price is given by
X n U ( Θ ( e ) ) e r ( T n ) 1 M m = 1 M X ( m ) ( θ ¯ ( U ) ) : = V ^ n ( θ ¯ ( U ) ) , say .
Let V n G denote the known price of the European call Ether option at time n, which is taken as the Black–Scholes–Merton call price at time n. Write V ^ n G for the Monte Carlo estimate for the known price V n G . Then, with control variates, the Monte Carlo estimate, denoted by V ^ n c v ( θ ¯ ( U ) ) , for the Bayesian (posterior) upper estimate for the price of the Ether option is given by
V ^ n c v ( θ ¯ ( U ) ) = V ^ n ( θ ¯ ( U ) ) + β ( V n G V ^ n G ) ,
where
β = Cov ( V ^ n ( θ ¯ ( U ) ) , V ^ n G ) Var ( V ^ n G ) .
The standard error of the Monte Carlo estimate in Equation (75) is given by
SD ( V ^ n c v ( θ ¯ ( U ) ) ) = Var ( V ^ n c v ( θ ¯ ( U ) ) ) + β 2 Var ( V ^ n G ) 2 β Cov ( V ^ n ( θ ¯ ( U ) ) , V ^ n G ) .
In line with Boyle and Ananthanarayanan (1977), a Bayesian model for estimating the unknown volatility parameter in the Black–Scholes–Merton (BSM) model is adopted to compute Bayesian confidence limits and a Bayesian point estimate for the price of an Ether option in the current paper. For convenience, this model is called the Bayesian BSM model here. The Bayesian BSM model corresponds to a model with volatility uncertainty (but without conditional heteroscedasticity), and it is briefly discussed in the rest of this section for the sake of completeness. In Section 6, the Bayesian confidence limits obtained from the Bayesian BSM model will be compared with those obtained from the proposed model, which incorporates both model uncertainty and conditional heteroscedasticity.
In the Bayesian BSM model, the logarithmic returns of Ethereum are directly used to compute the Bayesian credible interval for the unknown volatility parameter. Let Y ( n ) : = { Y 1 , Y 2 , , Y n } (i.e., the observations about the logarithmic returns of Ethereum up to and including time n). Define
Y ¯ n : = i = 1 n Y i n , s Y 2 ( n ) : = 1 n i = 1 n ( Y i Y ¯ n ) 2 .
As in Equation (50), define β n Y and σ n Y as follows:
β n Y = β 0 + 1 2 n s Y 2 ( n ) + t 0 n ( μ 0 Y ¯ n ) 2 2 ( t 0 + n ) , σ n Y = β n Y α 0 + n 2 .
Similarly to Equations (51) and (52), the 100 ( 1 γ 2 ) % Bayesian credible interval for the unknown annualized volatility of the logarithmic returns of Ethereum, i.e., σ Y , given Y ( n ) , is ( L γ 2 ( σ Y | Y ( n ) ) , U γ 2 ( σ Y | Y ( n ) ) ) such that
L γ 2 ( σ Y | Y ( n ) ) = Γ α 0 + n 2 , β n Y γ 2 2 1 2 × 252 , U γ 2 ( σ Y | Y ( n ) ) = Γ α 0 + n 2 , β n Y 1 γ 2 2 1 2 × 252 ,
where the scaling factor 252 is applied to annualize the lower and upper limits of the Bayesian credible interval for σ Y under the assumption that there are 252 trading days per annum.
Let r Y denote the annualized risk-free interest rate. Then, given the volatility σ Y , the Black–Scholes–Merton model gives the following closed-form solution to the price of the European call Ether option with strike price K and time-to-maturity τ :
C ( s , r Y , σ Y , K , τ ) = s Φ ( d 1 ) K e r Y τ Φ ( d 2 ) ,
where
d 1 : = ln ( s K ) + ( r Y + 1 2 σ Y 2 ) τ σ Y τ , d 2 : = d 1 σ Y τ .
Note that s is the current price of Ethereum in AUD; τ is measured in year(s).
The Bayesian point estimate for σ Y given Y ( n ) , denoted by M ( σ Y | Y ( n ) ) , is taken as follows:
M ( σ Y | Y ( n ) ) = σ n Y × 252 .
Then, the (closed-form) Bayesian lower and upper estimates for the price of the European call Ether option given Y ( n ) are, respectively, given by
C n L ( s , r Y , K , τ ) = C ( s , r Y , L γ 2 ( σ Y | Y ( n ) ) , K , τ ) ,
and
C n U ( s , r Y , K , τ ) = C ( s , r Y , U γ 2 ( σ Y | Y ( n ) ) , K , τ ) .
Furthermore, the (closed-form) Bayesian point estimate for the price of the European call Ether option given Y ( n ) is given by
C n P ( s , r Y , K , τ ) = C ( s , r Y , M ( σ Y | Y ( n ) ) , K , τ ) .

6. Empirical and Simulation Results

Real data on the daily adjusted close prices of Ethereum in AUD (ETH-AUD) are considered. The period of the data is from 1 January 2018 to 21 May 2024 (a total of 2333 observations) covering some periods of the COVID-19 pandemic. The data were downloaded from Ethereum AUD (ETH-AUD) price history and historical data—Yahoo Finance (https://au.finance.yahoo.com/quote/ETH-AUD/, (accessed on 15 September 2024)). Data on the daily logarithmic returns from ETH-AUD (a total of 2332 observations) are used to estimate the reference model, i.e., the GARCH(1,1) model with drift, in Equation (1). The estimation of the GARCH(1,1) model was performed using the R package “rugarch” Galanos and Kley (2023). The Monte Carlo estimates and their standard errors for the Bayesian lower and upper estimates for prices of Ether options with different strike prices and maturities are computed. The strike prices are 80%, 90%, 100%, 110% and 120% of the current adjusted close price of ETH-AUD. Consequently, in-the-money (ITM), at-the-money (ATM) and out-of-the-money (OTM) Ether options are considered. The maturities are three months, six months, nine months and one year. It is supposed that there are 252 trading days per year. Besides the Bayesian lower and upper estimates, the prices of Ether options without incorporating model uncertainty, which are called the standard prices, are also considered. Specifically, Monte Carlo estimates and their standard errors for the standard prices with different strike prices and maturities are computed. The standard prices refer to the prices of the Ether options from a model incorporating conditional heteroscedasticity only, that is, the reference model in Equation (1). The Bayesian lower and upper estimates, as well as the Bayesian point estimates, for the prices of the Ether options are computed from the model incorporating volatility uncertainty only in line with Boyle and Ananthanarayanan (1977), as described in Section 5. All the statistical analyses, simulations and computations were performed using R.
Table 1 presents some summary statistics of the daily logarithmic returns data.
From Table 1, the daily logarithmic returns from ETH-AUD have a negatively skewed and heavy-tailed distribution. Furthermore, the range of the returns is large. Similar findings were discussed in Siu (2023a) for Bitcoin returns data covering some periods of the COVID-19 pandemic.
Figure 1 depicts the time series plots for the daily adjusted close prices (a) and logarithmic returns (b) from ETH-AUD. [subfigure]justification=centering
From Figure 1a, there was an upward trend in the adjusted close prices of ETH-AUD starting around the 1000th day and ending at around the 1500th day. This corresponds to the first half year of 2020 during which the impact of the COVID-19 pandemic started emerging. Furthermore, from Figure 1b, the daily logarithmic returns from ETH-AUD exhibit volatility clustering, that is, periods of high (low) volatility followed by periods of high (low) volatility. Similar observations were made in Siu (2023a) for Bitcoin returns data covering some periods of the COVID-19 pandemic.
Figure 2 depicts the sample ACF plot (a) and the sample PACF plot (b) for the daily logarithmic returns from ETH-AUD.
From Figure 2, the daily logarithmic returns from ETH-AUD appear to be stationary. To provide further evidence for the stationarity, three types of the augmented Dickey–Fuller (ADF) test were performed using the R package “aTSA” Qiu (2024). They are “Type 1: no drift no trend”, “Type 2: with drift no trend”, and “Type 3: with drift and trend”. The alternative hypothesis of the test is “stationary”. The p values of the three types of the ADF test are all less than 0.01 , indicating the stationarity of the logarithmic returns data.
Table 2 provides the estimation results for the GARCH(1,1) model with drift. Note that the standard errors reported are robust standard errors.
From the estimation results in Table 2, the standard error of the estimate for the drift parameter μ Y is large compared with the estimate for μ Y . This indicates that estimating the drift parameter could be noisy, which has been documented in Merton (1980) for market returns data. Furthermore, the (G)ARCH coefficients α and β are significant at the 5% significance level. This provides evidence for the presence of conditional heteroscedasticity. To provide further evidence, a weighted Ljung–Box test on the standardized squared residuals and a weighted ARCH LM test were performed. Table 3 presents the results of the weighted Ljung–Box test on standardized squared residuals and the weighted ARCH LM test. The results of the two tests support the presence of conditional heteroscedasticity or ARCH effect in the daily logarithmic returns from ETH-AUD.
Table 4 provides the lower and upper limits of the 90% Bayesian credible intervals for the “uncertain” parameters μ and σ , which were computed using the Normal–Gamma conjugate prior and the standardized residuals from the estimated reference model. As an illustration, the assumed values for the prior parameters are μ 0 = 0.05 , t 0 = 100 , α 0 = 2 and β 0 = 1 . These assumed values for the prior parameters were considered in Siu (2023a) for Bitcoin returns data.
From Table 4, the Bayesian credible interval for μ covers the point zero. This indicates that the Bayesian credible interval includes the situation where there are no perturbations/misspecifications in the conditional mean. Likewise, the Bayesian credible interval for σ covers the point one. This indicates that the Bayesian credible interval includes the situation where there are no perturbations/misspecifications in the conditional volatility. The Bayesian credible interval for μ appears to be relatively wide compared with that for σ . This indicates that model uncertainty about the conditional mean could be more significant than model uncertainty about the conditional volatility for the daily logarithmic returns from ETH-AUD. This appears to be in line with the finding in Merton (1980) that the drift estimate could be more noisy than the volatility estimate for market returns data. The lower and upper limits of the Bayesian credible intervals for μ and σ in Table 4 are used to compute the Monte Carlo estimates and their standard errors for the Bayesian lower and upper estimates for the prices of Ether options following the procedures described in Section 5.
Table 5, Table 6 and Table 7 give the Monte Carlo estimates for the Bayesian lower and upper estimates for the prices of the Ether options as well as their standard prices with different strike prices K and maturities T. The respective standard errors are presented in parentheses below the Monte Carlo estimates. Each Monte Carlo estimate is computed by 10,000 simulation runs. The risk-free interest rate r is taken as 5% per annum. The current adjusted close price of ETH-AUD s 0 = 5693.319 .
From Table 5, Table 6 and Table 7, the Bayesian lower and upper estimates as well as the standard prices increase as the maturity T increases. Furthermore, they decrease as the strike price K increases. These results are intuitively reasonable. Specifically, the results can be explained by the time value of the call option and the payoff of the call option as a decreasing function of the strike price. Furthermore, the Bayesian lower estimates are lower than the respective Bayesian upper estimates, and the standard prices lie between the respective Bayesian lower and upper estimates. This may lend support to the conjecture that the supremum and infimum of E θ , [ X | F n ] over θ Θ ( e ( n ) ) may be attained at the extremal, or boundary, points of the rectangular region Θ ( e ( n ) ) . The standard errors are relatively small compared with the respective Monte Carlo estimates. This indicates that the Monte Carlo estimates for the Bayesian lower and upper prices as well as the standard prices are quite precise.
Table 8 presents the differences in the Monte Carlo estimates for the Bayesian upper and lower estimates for the prices of the Ether options with different strikes K and maturities T.
From Table 8, for a fixed strike price K, the differences in the Monte Carlo estimates for the Bayesian upper and lower estimates for the prices of the Ether options increase as the maturity T increases. This indicates that as the maturity T of an Ether option becomes longer, the effect of model uncertainty about the conditional mean and volatility on the option prices could be more significant.
Table 9 gives the Bayesian upper, point and lower estimates for the prices of the Ether options with different strikes K and maturities T under the Bayesian BSM model. Note that the prior parameters used here are the same as those values assumed before. The 90% Bayesian credible interval for the annualized volatility parameter σ Y is given by ( 0.862603 , 0.8955688 ) . The Bayesian point estimate for σ Y is given by 0.878703 .
From Table 9, the Bayesian lower, point and upper estimates for the prices of the Ether options increase as the maturity T increases. They decrease as the moneyness increases. Again, these results are intuitively reasonable. The Bayesian lower, point and upper estimates for the prices of the Ether options are in correct order. That is, the Bayesian lower estimates are below the respective Bayesian upper estimates, and the Bayesian point estimates lie between the respective Bayesian lower and upper estimates. This is explained by the fact that the Black–Scholes–Merton formula for the call price is an increasing function of the annualized volatility parameter σ Y , and the numerical results that the Bayesian point estimate for σ Y lies within the Bayesian credible interval for σ Y .
Table 10 presents the differences in the upper and lower estimates for the prices of the Ether options with different strikes K and maturities T under the Bayesian BSM model.
From Table 10, for a fixed strike price, the differences in the upper and lower estimates for the prices of the Ether options increase as the maturity T increases. This result is consistent with that from Table 8. This reflects that as the maturity T of an Ether option becomes longer, the impact of volatility uncertainty on the option prices could be more pronounced. Comparing Table 8 and Table 10, the differences in the upper and lower estimates for the prices of the Ether options under the proposed model with both model uncertainty and conditional heteroscedasticity are substantially higher than those under the Bayesian BSM model. This indicates that the combined effect of model uncertainty about the conditional mean and conditional heteroscedasticity on the option prices are more significant than the impact of volatility uncertainty.

7. Price Fluctuations: Chaotic or Random?

Are price fluctuations chaotic or random? This question has captured considerable attention from academic researchers and market practitioners. If price fluctuations are chaotic, they may be considered a deterministic (nonlinear) dynamical system. If price fluctuations are random, they are totally unpredictable and may be described by the random walk theory. The latter is related to a (weak) form of market efficiency, although rejecting random walk theory does not imply market efficiency when trading costs are taken into account (see (Taylor 1986, 2005)). Henri Poincaré realized more than a century ago that a deterministic (nonlinear) dynamical system, which is sensitive to its initial state, can generate behaviour which is reminiscent of random behaviour. In a seminal article by May (1976), a logistic map was used to illustrate chaotic dynamics in the study of animal population dynamics. Although it may not be possible to define the notion of chaos precisely, the sensitivity to the initial state is an important property of a chaotic dynamical system. See (Tong 1990, 1995, 2002) for related discussions and more detail. See also the monographs by (Gleick 1987; Palmer 2023) for an excellent account to the chaos theory and its applications to diverse fields.
Motivated by the market crash in October 1987, an early study on deterministic chaotic and nonlinear dynamical systems in financial markets was conducted in a seminal article by Hsieh (1991). Through performing a statistical test, namely the BDS test Broock et al. (1996), on weekly stock returns from the Center for Research in Securities Prices (CRSP) at the University of Chicago, the empirical results in Hsieh (1991) revealed that logarithmic stock returns are not independent and identically distributed (i.i.d.) and that conditional heteroscedasticity is the cause of the rejection of the i.i.d. logarithmic stock returns. Hsieh (1991) also pointed out that to model nonlinear behaviour in the logarithmic stock returns, instead of focusing on the conditional mean returns including chaotic dynamics, attention should be given to modelling conditional heteroscedasticity, i.e., time-varying conditional volatility. (Hsieh 1993, 1995) further explored implications of nonlinear dynamical systems for financial risk management and financial markets, respectively. Using data on four currency futures markets (the British Pound, Deutsche Mark, Japanese Yen and Swiss Franc), Hsieh (1993) found empirically that while conditional mean returns may not be predictable, the conditional variances are predictable. Using data on stock index futures prices for the S&P, Nikkei, DAX, FTSE and CAC, the empirical findings in Hsieh (1995) indicated that conditional variance models such as the ARCH-type models can account for a great deal of non-random behaviour in asset returns. For excellent accounts on chaos and nonlinear dynamical systems in financial markets and the economy, one may refer to the monographs by (Broock et al. 1991; Farmer 2024; Peters 1996). Additionally, Schefczyk (2012) may provide some discussions from a philosophical perspective.
Deterministic nonlinear dynamical systems and chaos may not provide realistic descriptions for real-world dynamics. To address some limitations of deterministic nonlinear dynamical systems and chaos, stochastic dynamical systems and noisy chaos were introduced by researchers in the 1990s. See, for example, (Bartlett 1990; Berliner 1992; Casdagli 1992; Cheng and Tong 1992, 1994; Fan et al. 1993; Tong 1990, 1992, 1995; Yao and Tong 1994a, 1994b; Chan and Tong 1994). Specifically, stochastic dynamical systems and noisy chaos are described by stochastic difference equations or nonlinear autoregressive models which generalize deterministic difference equations by introducing stochastic noise. Two important approaches have been introduced to quantify the sensitivity of stochastic dynamical systems with respect to their initial conditions. They are the conditional distribution approach and the conditional mean approach. The conditional distribution approach was considered in (Yao and Tong 1994a, 1994b). The key idea of this approach is to describe the sensitivity of the conditional distribution with respect to the initial state through the divergence of the conditional distribution relating to the Kullback–Leibler information. The conditional mean approach describes the sensitivity of the conditional mean with respect to the initial state. It aims to generalize the notion of the Lyapunov exponent in a stochastic environment. See (Chan and Tong 2001; Tong 1995, 2002; Yao and Tong 1994a) for more detailed discussions on the two approaches. Indeed, the sensitivity of stochastic dynamical systems with respect to their initial states has an important implication for nonlinear prediction. Specifically, how well you can predict depends on where you are now (Chan and Tong 2001; Fan and Yao 2003; Tong 1995; Yao and Tong 1994a). This highlights the importance of state-dependency for nonlinear prediction. Stochastic (nonlinear) dynamical systems have been explored for modelling, estimating and predicting volatility in (Elliott et al. 2012; Tong 2007). Yao and Tong (1998) introduced a bootstrap test for deterministic chaos and stochastic randomness. As a potential topic for further research, it might be interesting to explore the possibility of using this bootstrap test to examine whether price fluctuations in financial markets, or particularly in crypto markets, are chaotic or random. Indeed, Garnier and Solna (2019) and Pietrych (2021) studied the nonlinearity, chaos and multifractality in the cryptocurrency markets as well as their implications for the market efficiency of cryptocurrencies.

8. Conclusions

Using Bayesian nonlinear expectations introduced in Siu (2023a), Bayesian lower and upper estimates for the price of an Ether option are constructed with the objective of incorporating the impacts of model uncertainty and conditional heteroscedasticity. Specifically, model uncertainty about the conditional mean and volatility of the return dynamics of Ethereum modelled by a reference GARCH model is incorporated by introducing the “uncertain” parameters describing the perturbations/misspecifications in the conditional mean and volatility. Then, a family of alternative GARCH models are introduced via changing probability measures using an extended Girsanov’s principle. The risk-neutralization of the family of alternative GARCH models is also performed via the extended Girsanov’s principle. The Bayesian nonlinear expectations are constructed using closed-form Bayesian credible intervals with conjugate priors and the residuals from the estimated reference GARCH model. The numerical/simulation procedures for the computation of the Bayesian lower and upper estimates for the prices of Ether options are introduced based on the Monte Carlo simulation method with control variates. Empirical results from real data on the daily logarithmic returns from Ethereum in AUD are presented. Monte Carlo estimates and their standard errors for the Bayesian lower and upper estimates for the prices of the Ether options as well as their standard prices are provided. The Monte Carlo estimates are precise and consistent with intuition. It is also found that the impact of model uncertainty about the conditional mean and volatility on the Ether option prices is more significant as the maturity increases. Comparisons between the Bayesian lower and upper estimates for the prices of the Ether options under the Bayesian BSM model, which incorporates volatility uncertainty only, are provided. It is found that the combined effect of model uncertainty about the conditional mean and conditional heteroscedasticity on the option prices is more pronounced than the impact of volatility uncertainty. Some discussions about the chaotic and random aspects of price fluctuations in financial markets in general, or crypto markets in particular, are provided. The chaotic and random aspects of price fluctuations may have implications for market efficiency in cryptocurrencies (see, for example, Garnier and Solna (2019) and Pietrych (2021)). The latter is relevant to crypto option pricing as discussed in Siu and Elliott (2021). It is hoped that those discussions might provide insights and inspiration, which may be helpful for future endeavours to explore the chaotic and random aspects of price fluctuations, as well as their implications for financial markets in general, or crypto markets in particular.
Some limitations of the proposed model are briefly discussed here. Firstly, the conjugate priors for the “uncertain” parameters describing the perturbations/misspecifications in the conditional mean and volatility are adopted here. Although the conjugate priors lead to closed-form Bayesian credible intervals, it is possible that the “uncertain” parameters may be better modelled by other distributions. Indeed, one may explore the possibility of using a Bayesian nonparametric approach to construct Bayesian credible intervals for the “uncertain” parameters as noted in Siu (2023a), Online Appendix F. Secondly, to simplify the discussion, it is assumed here that the random shocks of the GARCH model are normally distributed. It is possible that the random shocks of the GARCH model follow non-normal distributions. A similar remark was made in Boyle and Ananthanarayanan (1977) for the Bayesian BSM model. Thirdly, for illustration, the GARCH model is adopted here as a reference model. It is possible that some variants of the GARCH model, such as the GJR-GARCH model and the SETAR-GARCH model, may be better candidates for a reference model. Here, it is assumed that the risk-free interest rate is a constant. However, in practice, it may be realistic to consider the situation of a stochastic interest rate. It is conjectured here that the infimum and supremum of a conditional expectation are attained at the extremal points of a rectangular region when approximating the Bayesian (posterior) superlinear and sublinear expectations. Although the numerical results presented here may lend support to this conjecture, it is possible that the the infimum and supremum may be attained at some interior points of the rectangular region. Although the chaotic and random aspects of price fluctuations are briefly touched on in Section 7, empirical studies on the chaotic and random aspects of price fluctuations in Ether and other cryptocurrencies are not provided in the current paper. The papers by Garnier and Solna (2019) and Pietrych (2021) may provide an excellent account on this topic.
As noted in Godbole (2023), bid and ask prices prevail in Bitcoin and Ether option markets. However, it seems that the literature on modelling the bid and ask prices in cryptocurrency option markets may be rather limited. In a recent paper by Madan et al. (2019), they calibrated the Bitcoin price dynamics under several Markov models using vanilla Bitcoin options and studied the implied liquidity of Bitcoin call options using conic finance theory. How might the prevalence of bid and ask prices in Bitcoin and Ether options be explained? This may be a potential topic for further research. From the perspective of a two-price economy (Cherny and Madan 2009; Madan 2015, 2016; Madan and Cherny 2010), a recent paper by Elliott et al. (2022) considered a two-price system based on a lower price process and an upper price process attributed to model uncertainty, or uncertainty in probability, for financial valuation. Suppose the Bayesian lower and upper estimates for the price of an Ether option are interpreted as a Bayesian version of the lower and upper prices of the Ether option. Then, one may consider a Bayesian two-price system for the Ether option, in which the prevalence of bid and ask prices of the Ether option is explained by the model uncertainty about the conditional mean and volatility of the logarithmic returns of Ethereum. In this case, some further questions may arise. For example, besides model uncertainty, can the prevalence of bid and ask prices of the Ether option be explained by other economic factors? If yes, what other economic factors might be included? Another potential topic for further research may be to explore the possibility of developing a Bayesian version of nonlinear valuation operators in (Coque et al. 2002; Peng 2004a, 2004b, 2004c, 2005). For example, can a Bayesian version of nonlinear valuation operators be developed in continuous time?

Funding

The author wishes to acknowledge the Discovery Grant from the Australian Research Council (ARC) (project number: DP190102674).

Data Availability Statement

The data supporting reported results can be found from Yahoo Finance (https://au.finance.yahoo.com/quote/ETH-AUD/ (accessed on 15 September 2024)).

Acknowledgments

The author sincerely thanks the three reviewers for helpful and insightful comments. The author would like to thank the comment from a referee which inspired the discussions in Section 7.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Agarwal, Ankush, Christian Oliver Ewald, and Yihan Zou. 2021. An RBSDE Based Monte Carlo Method for Solving Optimal Stopping Problems under Parameter Uncertainty. Available online: https://ssrn.com/abstract=3549891 (accessed on 21 June 2021).
  2. Alexander, Carol, Ding Chen, and Arben Imeraj. 2023. Crypto quanto and inverse options. Mathematical Finance 33: 1005–43. [Google Scholar] [CrossRef]
  3. Avellaneda, Marco, Arnon Levy, and Antonio Paras. 1995. Pricing and hedging derivative securities in markets with uncertain volatilities. Applied Mathematical Finance 2: 73–88. [Google Scholar] [CrossRef]
  4. Badescu, Alexandru M., and Reg J. Kulperger. 2008. GARCH option pricing: A semiparametric approach. Insurance: Mathematics and Economics 43: 69–84. [Google Scholar] [CrossRef]
  5. Badescu, Alexandru M., Robert J. Elliott, Reg J. Kulperger, Jarkko Miettinen, and Tak Kuen Siu. 2011. A comparison of pricing kernels for GARCH option pricing with generalized hyperbolic distributions. International Journal of Theoretical and Applied Finance 14: 669–708. [Google Scholar] [CrossRef]
  6. Bartlett, Maurice S. 1990. Chance or chaos (with discussion). Journal of Royal Statistical Society A 153: 321–47. [Google Scholar] [CrossRef]
  7. Bauwens, Luc, and Michel Lubrano. 2002. Bayesian option pricing using asymmetric GARCH models. Journal of Empirical Finance 9: 321–42. [Google Scholar] [CrossRef]
  8. Beissner, Patrick Beissner, and Frank Riedel. 2019. Equilibria under Knightian price uncertainty. Econometrica 87: 37–64. [Google Scholar] [CrossRef]
  9. Berliner, L. Mark. 1992. Statistics, probability and chaos. Statistical Science 7: 69–122. [Google Scholar] [CrossRef]
  10. Bernardo, José M., and Adrian F. M. Smith. 2000. Bayesian Theory. New York: John Wiley & Sons, Inc. [Google Scholar]
  11. Bollerslev, Tim. 1986. Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics 31: 307–27. [Google Scholar] [CrossRef]
  12. Boyle, Phelim P. 1977. Options: A Monte Carlo approach. Journal of Financial Economics 4: 323–38. [Google Scholar] [CrossRef]
  13. Boyle, Phelim P., and A. L. Ananthanarayanan. 1977. The impact of variance estimation in option valuation models. Journal of Financial Economics 5: 375–87. [Google Scholar] [CrossRef]
  14. Boyle, Phelim P., Mark Broadie, and Paul Glasserman. 1997. Monte Carlo methods for security pricing. Journal of Economic Dynamics and Control 21: 1267–321. [Google Scholar] [CrossRef]
  15. Boyle, Phelim P., Shui Feng, Weidong Tian, and Tan Wang. 2008. Robust stochastic discount factors. Review of Financial Studies 21: 1077–122. [Google Scholar] [CrossRef]
  16. Brini, Alessio, and Jimmie Lenz. 2024. Pricing cryptocurrency options with machine learning regression for handling market volatility. Economic Modelling 136: 106752. [Google Scholar] [CrossRef]
  17. Broock, William A., David A. Hsieh, and Blake LeBaron. 1991. Nonlinear Dynamics, Chaos and Instability: Statistical Theory and Economic Evidence. Cambridge: MIT Press. [Google Scholar]
  18. Broock, William A., José A. Scheinkman, W. Davis Dechert, and Blake Lebaron. 1996. A test for independence based on the correlation dimension. Econometric Reviews 15: 197–235. [Google Scholar] [CrossRef]
  19. Buraschi, Andrea, and Alexei Jiltsov. 2006. Model uncertainty and option markets with heterogeneous beliefs. Journal of Finance 61: 2841–97. [Google Scholar] [CrossRef]
  20. Cao, Melanie, and Batur Celik. 2021. Valuation of bitcoin options. The Journal of Futures Markets 41: 1007–26. [Google Scholar] [CrossRef]
  21. Casdagli, Martin. 1992. Chaos and deterministic versus stochastic nonlinear modelling. Journal of Royal Statistical Society Series B 54: 303–28. [Google Scholar] [CrossRef]
  22. Chan, Jeremy. 2022. CME to Launch Ether Options Contract. Financial News. Dow Jones & Company, Inc. September. Available online: https://www.fnlondon.com/articles/cme-ether-bitcoin-options-futures-crypto-merge-derivatives-contract-20220913 (accessed on 15 September 2024).
  23. Chan, Kung-Sik, and Howell Tong. 1994. A note on noisy chaos. Journal of Royal Statistical Society Series B 56: 301–11. [Google Scholar] [CrossRef]
  24. Chan, Kung-Sik, and Howell Tong. 2001. Chaos: A Statistical Approach. New York: Springer. [Google Scholar]
  25. Chau, Huy N., Masaaki Fukasawa, and Miklós Rásonyi. 2022. Super-replication with transaction costs under model uncertainty for continuous processes. Mathematical Finance 32: 1066–85. [Google Scholar] [CrossRef]
  26. Chen, Kuo Shing, and J. Jimmy Yang. 2024. Price dynamics and volatility jumps in bitcoin options. Financial Innovation 10: 132. [Google Scholar] [CrossRef]
  27. Cheng, Bing, and Howell Tong. 1992. Consistent nonparametric order determination and chaos (with discussion). Journal of Royal Statistical Society Series B 54: 427–449, 451–474. [Google Scholar] [CrossRef]
  28. Cheng, Bing, and Howell Tong. 1994. Orthogonal projection, embedding dimension and sample size in chaotic time series from a statistical perspective. Philosophical Transactions of the Royal Society (London) A 348: 325–41. [Google Scholar]
  29. Cherny, Alexander, and Dilip B. Madan. 2009. New measures for performance evaluation. The Review of Financial Studies 22: 2571–606. [Google Scholar] [CrossRef]
  30. Cohen, Samuel N., and Martin Tegnér. 2019. European option pricing with stochastic volatility models under parameter uncertainty. In Frontiers in Stochastic Analysis–BSDEs, SPDEs and Their Applications. Edited by Samuel N. Cohen, Istvan Gyongy, Goncalo dos Reis, David Siska and Lukasz Szpruch. BSDE-SPDE 2017. Springer Proceedings in Mathematics & Statistics. Cham: Springer, vol. 289. [Google Scholar] [CrossRef]
  31. Cont, Rama. 2006. Model uncertainty and its impact on the pricing of derivative instruments. Mathematical Finance 16: 519–47. [Google Scholar] [CrossRef]
  32. Copeland, Tim. 2023. A Brief History of Ethereum. The Block. September. Available online: https://www.theblock.co/learn/245716/a-brief-history-of-ethereum (accessed on 15 September 2024).
  33. Coque, François, Ying Hu, Jean Mémin, and Shige Peng. 2002. Filtration-consistent nonlinear expectations and related g-expectations. Probability Theory and Related Fields 123: 1–27. [Google Scholar] [CrossRef]
  34. Darsinos, Theo, and Stephen Satchell. 2007. Bayesian analysis of the Black-Scholes option price. In Forecasting Expected Returns in the Financial Markets. Edited by Stephen Satchell. Cambridge: Academic Press, pp. 117–50. [Google Scholar]
  35. Delbaen, Freddy, and Walter Schachermayer. 1994. A general version of the fundamental theorem of asset pricing. Mathematische Annalen 300: 463–520. [Google Scholar] [CrossRef]
  36. Duan, Jin-Chuan. 1995. The GARCH option pricing model. Mathematical Finance 5: 13–32. [Google Scholar] [CrossRef]
  37. Dunford, Nelson, and Jacob T. Schwartz. 1958. Linear Operators: Part I. New York: Interscience. [Google Scholar]
  38. Elliott, Robert J. 1982. Stochastic Calculus and Applications. Berlin/Heidelberg and New York: Springer. [Google Scholar]
  39. Elliott, Robert J. 2017. Filtering with uncertain noise. IEEE Transactions on Automatic Control 62: 876–81. [Google Scholar] [CrossRef]
  40. Elliott, Robert J., and Dilip B. Madan. 1998. A discrete time equivalent martingale measure. Mathematical Finance 8: 127–52. [Google Scholar] [CrossRef]
  41. Elliott, Robert J., and Tak Kuen Siu. 2017. Stochastic volatility with regime switching and uncertain noise: Filtering with sub-linear expectations. Discrete and Continuous Dynamical Systems Series B 22: 59–81. [Google Scholar] [CrossRef]
  42. Elliott, Robert J., Dilip B. Madan, and Tak Kuen Siu. 2022. Lower and upper pricing of financial assets. Probability, Uncertainty and Quantitative Risk 7: 45–66. [Google Scholar] [CrossRef]
  43. Elliott, Robert J., Lakhdar Aggoun, and John B. Moore. 1995. Hidden Markov Models: Estimation and Control, 1st ed. New York: Springer. [Google Scholar]
  44. Elliott, Robert J., Tak Kuen Siu, and Eric S. Fung. 2012. Filtering a nonlinear stochastic volatility model. Nonlinear Dynamics 67: 1295–313. [Google Scholar] [CrossRef]
  45. Elliott, Robert J., Tak Kuen Siu, and Leunglung Chan. 2006. Option pricing for GARCH models with Markov switching. International Journal of Theoretical and Applied Finance 9: 825–41. [Google Scholar] [CrossRef]
  46. Fan, Jianqing, and Qiwei Yao. 2003. Nonlinear Time Series: Nonparametric and Parametric Methods. New York: Springer. [Google Scholar]
  47. Fan, Jianqing, Qiwei Yao, and Howell Tong. 1993. Estimation of conditional densities and sensitivity measures in nonlinear dynamical systems. Biometrika 83: 189–206. [Google Scholar] [CrossRef]
  48. Farmer, J. Doyne. 2024. Making Sense of Chaos: A Better Economics for a Better World. London: Penguin. [Google Scholar]
  49. Foster, F. Douglas, and Charles H. Whiteman. 1999. An application of Bayesian option pricing to the soybean market. American Journal of Agricultural Economics 81: 722–27. [Google Scholar] [CrossRef]
  50. Fouque, Jean-Pierre, and Bin Ren. 2014. Approximation for option prices under uncertain volatility. SIAM Journal on Financial Mathematics 5: 360–83. [Google Scholar] [CrossRef]
  51. Galanos, Alexios, and Tobias Kley. 2023. Univariate GARCH Models. R Package “Rugarch”. Version 1.5-1. September 20. Available online: https://cran.r-project.org/web/packages/rugarch/rugarch.pdf (accessed on 15 September 2024).
  52. Garnier, Josselin, and Knut Solna. 2019. Chaos and order in the bitcoin market. Physica A: Statistical Mechanics and its Applications 524: 708–21. [Google Scholar] [CrossRef]
  53. Gleick, James. 1987. Chaos: Making a New Science. New York: Viking Books. [Google Scholar]
  54. Godbole, Omkar. 2020. Crypto Derivatives Exchange OKEx Launches Options on Ether. CoinDesk. June. Available online: https://www.coindesk.com/markets/2020/06/05/crypto-derivatives-exchange-okex-launches-options-on-ether/ (accessed on 15 September 2024).
  55. Godbole, Omkar. 2023. Bitcoin, Ether Options Order Books Signal Calm Amid Mounting Risks. CoinDesk. October. Available online: https://www.coindesk.com/markets/2023/10/10/bitcoin-ether-options-order-books-signal-calm-amid-mounting-risks/ (accessed on 15 September 2024).
  56. Gzyl, Henryk, Enrique ter Horst, and Samuel Malone. 2022. Towards a Bayesian Framework for Option Pricing. Available online: https://arxiv.org/pdf/cs/0610053 (accessed on 15 September 2024).
  57. Hansen, Lars Peter, and Scott F. Richard. 1987. The role of conditioning information in deducing testable restrictions implied by dynamic asset pricing models. Econometrica 55: 587–613. [Google Scholar] [CrossRef]
  58. Harrison, J. Michael, and David M. Kreps. 1979. Martingales and arbitrage in multiperiod securities markets. Journal of Economic Theory 20: 381–408. [Google Scholar] [CrossRef]
  59. Harrison, J. Michael, and Stanley R. Pliska. 1981. Martingales and stochastic integrals in the theory of continuous trading. Stochastic Processes and Their Applications 11: 215–80. [Google Scholar] [CrossRef]
  60. Harrison, J. Michael, and Stanley R. Pliska. 1983. A stochastic calculus model of continuous trading: Complete markets. Stochastic Processes and Their Applications 15: 313–16. [Google Scholar] [CrossRef]
  61. Hilliard, Jimmy E., and Julie T. D. Ngo. 2022. Bitcoin: Jumps, convenience yields, and option prices. Quantitative Finance 22: 2079–91. [Google Scholar] [CrossRef]
  62. Ho, Shu Wing, Alan Lee, and Alastair Marsden. 2011. Use of Bayesian estimates to determine the volatility parameter input in the Black-Scholes and binomial option pricing models. Journal of Risk and Financial Management 3: 74–96. [Google Scholar] [CrossRef]
  63. Hou, Ai Jun, Weining Wang, Cathy Y. H. Chen, and Wolfgang Karl Härdle. 2020. Pricing cryptocurrency options. Journal of Financial Econometrics 18: 250–79. [Google Scholar]
  64. Hsieh, David A. 1991. Chaos and nonlinear dynamics: Application to financial markets. Journal of Finance 46: 1839–77. [Google Scholar] [CrossRef]
  65. Hsieh, David A. 1993. Implications of nonlinear dynamics for financial risk management. Journal of Financial and Quantitative Analysis 28: 41–64. [Google Scholar] [CrossRef]
  66. Hsieh, David A. 1995. Nonlinear dynamics in financial markets: Evidence and implication. Financial Analyst Journal 51: 55–62. [Google Scholar] [CrossRef]
  67. Jalan, Akanksha, Roman Matkovskyy, and Saqib Aziz. 2021. The Bitcoin options market: A first look at pricing and risk. Applied Economics 53: 2026–41. [Google Scholar] [CrossRef]
  68. Jang, Huisu, and Jaewook Lee. 2019. Generative Bayesian neural network model for risk-neutral pricing of American index options. Quantitative Finance 19: 587–603. [Google Scholar] [CrossRef]
  69. Jaroszkowski, Bartosz, and Max Jensen. 2022. Valuation of European options under an uncertain market price of volatility risk. Applied Mathematical Finance 29: 213–26. [Google Scholar] [CrossRef]
  70. Karolyi, G. Andrew. 1993. A Bayesian approach to modeling stock return volatility for option valuation. Journal of Financial and Quantitative Analysis 28: 579–94. [Google Scholar] [CrossRef]
  71. Klibanoff, Peter, Massimo Marinacci, and Sujoy Mukerji. 2005. A smooth model of decision making under ambiguity. Econometrica 73: 1849–92. [Google Scholar] [CrossRef]
  72. Kong, Augustine, Jun S. Liu, and Wing Hung Wong. 1994. Sequential imputations and Bayesian missing data problems. Journal of the American Statistical Association 89: 278–88. [Google Scholar] [CrossRef]
  73. Li, Lun, Ali Arab, Jiqiang Liu, Jingxian Liu, and Zhu Han. 2019. Bitcoin options pricing using LSTM-based prediction model and blockchain statistics. Paper presented at IEEE International Conference on Blockchain (Blockchain), Atlanta, GA, USA, July 14–17; pp. 67–74. [Google Scholar]
  74. Lucic, Vladimir, and Artur Sepp. 2024. Valuation and hedging of cryptocurrency inverse options. Quantitative Finance 24: 851–69. [Google Scholar] [CrossRef]
  75. Lyons, Terry J. 1995. Uncertain volatility and the risk-free synthesis of derivatives. Applied Mathematical Finance 2: 117–33. [Google Scholar] [CrossRef]
  76. Madan, Dilip B. 2015. Asset pricing theory for two price economies. Annals of Finance 11: 1–35. [Google Scholar] [CrossRef]
  77. Madan, Dilip B. 2016. Benchmarking in two price financial markets. Annals of Finance 12: 201–19. [Google Scholar] [CrossRef]
  78. Madan, Dilip B., and Alexander Cherny. 2010. Markets as a counterparty: An introduction to conic finance. International Journal of Theoretical and Applied Finance 13: 1149–77. [Google Scholar] [CrossRef]
  79. Madan, Dilip B., Sofie Reyners, and Wim Schoutens. 2019. Advanced model calibration on bitcoin options. Digital Finance 1: 117–37. [Google Scholar] [CrossRef]
  80. Martin, Gael M., Catherine S. Forbes, and Vance L. Martin. 2005. Implicit Bayesian inference using option prices. Journal of Time Series Analysis 26: 437–62. [Google Scholar] [CrossRef]
  81. May, Robert M. 1976. Simple mathematical models with very complicated dynamics. Nature 261: 459–67. [Google Scholar] [CrossRef]
  82. McGleenon, Brian. 2024. Ether Options Trading Volume Hits All-Time High as Price Sits around $2300. The Block. January. Available online: https://www.theblock.co/post/275300/ether-options-trading-volume-hits-all-time-high-as-price-sits-around-2300 (accessed on 15 September 2024).
  83. Merton, Robert C. 1980. On estimating the expected return on the market: An exploratory investigation. Journal of Financial Economics 8: 323–61. [Google Scholar] [CrossRef]
  84. Mozumder, Sharif, Taufiq Choudhry, and Michael Dempsey. 2021. Option pricing model biases: Bayesian and Markov Chain Monte Carlo regression analysis. Computational Economics 57: 1287–305. [Google Scholar] [CrossRef]
  85. Pagnottoni, Paolo. 2019. Neural network models for Bitcoin option pricing. Frontiers in Artificial Intelligence 2: 5. [Google Scholar] [CrossRef]
  86. Palmer, Tim. 2023. The Primacy of Doubt. Oxford: Oxford University Press. [Google Scholar]
  87. Peng, Shige. 1997. Backward SDE and related g-expectation. In Backward Stochastic Differential Equations. Pitman Research Notes in Mathematics Series 364; London: Longman Scientific & Technical, pp. 141–59. [Google Scholar]
  88. Peng, Shige. 2004a. Dynamical evaluations. C. R. Acad. Sci. Paris Ser. I 339: 585–89. [Google Scholar] [CrossRef]
  89. Peng, Shige. 2004b. Filtration consistent nonlinear expectations and evaluations of contingent claims. Acta Mathematicae Applicatae Sinica English Series 20: 191–214. [Google Scholar] [CrossRef]
  90. Peng, Shige. 2004c. Nonlinear expectations, nonlinear evaluations and risk measures. In Stochastic Methods in Finance. Lecture notes in mathematics 1856. New York: Springer, pp. 138–65. [Google Scholar]
  91. Peng, Shige. 2005. Dynamically Consistent Nonlinear Evaluations and Expectations. Available online: https://arxiv.org/abs/math/0501415 (accessed on 15 September 2024).
  92. Peng, Shige. 2006. G-expectation, G-Brownian motion and related stochastic calculus of Itô’s type. In The Abel Symposium 2005, Abel Symposia 2. Berlin: Springer, pp. 541–67. [Google Scholar]
  93. Peng, Shige. 2019. Nonlinear Expectations and Stochastic Calculus under Uncertainty with Robust CLT and G-Brownian Motion. New York: Springer. [Google Scholar]
  94. Peters, Edgar E. 1996. Chaos and Order in the Capital Markets: A New View of Cycles, Prices, and Market Volatility, 2nd ed. New York: John Wiley & Sons, Inc. [Google Scholar]
  95. Pietrych, Lukasz, Julio E. Sandubete, and Lorenzo Escot. 2021. Solving the chaos model-data paradox in the cryptocurrency market. Communications in Nonlinear Science and Numerical Simulation 102: 105901. [Google Scholar] [CrossRef]
  96. Qin, Zhang, and Caio Almeida. 2020. A Bayesian nonparametric approach to option pricing. Brazilian Review of Finance 18: 115–37. [Google Scholar] [CrossRef]
  97. Qiu, Debin. 2024. Alternative Time Series Analysis. R Package “aTSA”. Version 3.1.2.1. February 21. Available online: https://cran.r-project.org/web/packages/aTSA/aTSA.pdf (accessed on 15 September 2024).
  98. Ramoni, Marco, and Paola Sebastiani. 2001. Robust learning with missing data. Machine Learning 45: 147–70. [Google Scholar] [CrossRef]
  99. Romboutsa, Jeroen V. K., and Lars Stentoft. 2014. Bayesian option pricing using mixed normal heteroskedasticity models. Computational Statistics and Data Analysis 76: 588–605. [Google Scholar] [CrossRef]
  100. Rubin, Donald B. 1976. Inference and missing data. Biometrika 63: 581–92. [Google Scholar] [CrossRef]
  101. Schefczyk, Michael. 2012. The financial crisis, the exemption view and the problem of the harmless torturer. Philosophy of Management 11: 25–38. [Google Scholar] [CrossRef]
  102. Shegokar, Nikhilesh. 2022. Option Pricing on Cryptocurrency Underlying Using Diffusion and Jump-Diffusion Stochastic Volatility Models. May. Available online: https://cryptoquant.uk/assets/documentation/main.pdf (accessed on 15 September 2024).
  103. Siu, Tak Kuen. 2001. Risk Measures in Finance and Insurance. Ph.D. thesis, The University of Hong Kong, Hong Kong SAR, China. [Google Scholar]
  104. Siu, Tak Kuen. 2023a. Bayesian nonlinear expectation for time series modelling and its application to Bitcoin. Empirical Economics 64: 505–37. [Google Scholar] [CrossRef]
  105. Siu, Tak Kuen. 2023b. European option pricing with market frictions, regime switches and model uncertainty. Insurance: Mathematics and Economics 113: 233–50. [Google Scholar] [CrossRef]
  106. Siu, Tak Kuen. 2024. Finite mixture models for option pricing: An application to Bitcoin options. Submitted. [Google Scholar]
  107. Siu, Tak Kuen, and Robert J. Elliott. 2021. Bitcoin option pricing with a SETAR-GARCH model. European Journal of Finance 27: 564–95. [Google Scholar] [CrossRef]
  108. Siu, Tak Kuen, Howell Tong, and Hailiang Yang. 2001. Bayesian risk measures for derivatives via random Esscher transform. North American Actuarial Journal 5: 78–91. [Google Scholar] [CrossRef]
  109. Siu, Tak Kuen, Howell Tong, and Hailiang Yang. 2004a. On pricing derivatives under GARCH models: A dynamic Gerber-Shiu approach. North American Actuarial Journal 8: 17–31. [Google Scholar] [CrossRef]
  110. Siu, Tak Kuen, Howell Tong, and Hailiang Yang. 2004b. On Bayesian Value at Risk: From linear to non-linear portfolios. Asia-Pacific Financial Markets 11: 161–84. [Google Scholar] [CrossRef]
  111. Siu, Tak Kuen, Howell Tong, and Hailiang Yang. 2006. Option pricing under threshold autoregressive models by threshold Esscher transform. Journal of Industrial and Management Optimization 2: 177–97. [Google Scholar]
  112. Tanner, Martin A., and Wing Hung Wong. 1987. The calculation of posterior distributions by data augmentation. Journal of the American Statistical Association 82: 528–50. [Google Scholar] [CrossRef]
  113. Taylor, Stephen J. 1986. Modelling Financial Time Series. Chichester: Wiley. [Google Scholar]
  114. Taylor, Stephen J. 2005. Asset Price Dynamics, Volatility, and Prediction. Princeton: Princeton University Press. [Google Scholar]
  115. Tong, Howell. 1977. Discussion of a paper by A. J. Lawrance and N. T. Kottegoda. Journal of Royal Statistical Society: Series A 140: 34–35. [Google Scholar]
  116. Tong, Howell. 1983. Threshold Models in Nonlinear Time Series Analysis. Berlin, Heidelberg and New York: Springer Verlag. [Google Scholar]
  117. Tong, Howell. 1990. Non-linear Time Series: A Dynamical System Approach. Oxford: Oxford University Press. [Google Scholar]
  118. Tong, Howell. 1992. Some comments on a bridge between nonlinear dynamicists and statisticians. Physica D: Nonlinear Phenomena 58: 299–303. [Google Scholar] [CrossRef]
  119. Tong, Howell. 1995. A personal overview of non-linear time series analysis from a chaos perspective (with discussion and rejoinder). Scandinavian Journal of Statistics 22: 399–445. [Google Scholar]
  120. Tong, Howell. 2002. Nonlinear time series analysis since 1990: Some personal reflections. Acta Mathematicae Applicatae Sinica, English Series 18: 177–84. [Google Scholar] [CrossRef]
  121. Tong, Howell. 2007. Exploring volatility from a dynamical system perspective. In 56th Session of the International Statistical Institute. Lisbon: PRT. [Google Scholar]
  122. Venter, Pierre J., and Eben Maré. 2021. Univariate and multivariate GARCH models applied to Bitcoin futures option pricing. Journal of Risk and Financial Management 14: 261. [Google Scholar] [CrossRef]
  123. Wang, Ying, Sai Tsang Boris Choy, and Hoi Ying Wong. 2016. Bayesian option pricing framework with stochastic volatility for FX data. Risks 4: 51. [Google Scholar] [CrossRef]
  124. Wilmott, Paul. 2007. Paul Wilmott Introduces Quantitative Finance, 2nd ed. Hoboken: Wiley. [Google Scholar]
  125. Wilmott, Paul, and Asli Oztukel. 1998. Uncertain parameters, an empirical stochastic volatility model and confidence limits. International Journal of Theoretical and Applied Finance 1: 175–89. [Google Scholar] [CrossRef]
  126. Yang, K. 2021. History of Crypto Derivatives. Crypto.com. February. Available online: https://crypto.com/research/history-of-crypto-derivatives (accessed on 15 September 2024).
  127. Yao, Qiwei, and Howell Tong. 1994a. Quantifying the influence of initial values on non-linear prediction. Journal of Royal Statistical Society Series B 56: 701–25. [Google Scholar] [CrossRef]
  128. Yao, Qiwei, and Howell Tong. 1994b. On prediction and chaos in stochastic systems. Philosophical Transactions of the Royal Society (London) A 348: 357–69. [Google Scholar]
  129. Yao, Qiwei, and Howell Tong. 1998. A bootstrap detection for operational determinism. Physica D: Nonlinear Phenomena 115: 49–55. [Google Scholar] [CrossRef]
  130. Zhang, Tong. 2023. Mathematical Analysis of Machine Learning Algorithms. Cambridge: Cambridge University Press. [Google Scholar]
Figure 1. Time series plots for adjusted close prices and log returns.
Figure 1. Time series plots for adjusted close prices and log returns.
Jrfm 17 00436 g001
Figure 2. SACF and SPACF plots for log returns.
Figure 2. SACF and SPACF plots for log returns.
Jrfm 17 00436 g002
Table 1. Summary statistics of daily log returns from ETH-AUD.
Table 1. Summary statistics of daily log returns from ETH-AUD.
Min.1st Qu.MedianMean
−0.5189823−0.01855880.00088160.0007498
3rd Qu.Max.SkewKurtosis
0.02251480.2283313−0.8810.32
Table 2. The estimation results of the GARCH(1,1) model with drift.
Table 2. The estimation results of the GARCH(1,1) model with drift.
ParameterEstimateStd. Errort ValuePr(> | t | )
μ Y 0.0010970.0008841.24070.214723
ω 0.0000410.0000351.16000.246037
α 0.0717530.0304672.35510.018517
β 0.9135440.04007122.79840.000000
Table 3. Weighted Ljung–Box test on standardized squared residuals and ARCH LM test.
Table 3. Weighted Ljung–Box test on standardized squared residuals and ARCH LM test.
Ljung–Box ARCH LM
Lag Statistic p -Value Lag Statistic p -Value
10.051120.821130.23060.6311
52.517410.502254.64610.1240
94.341780.531975.00070.2245
Table 4. 90% Bayesian credible intervals for μ and σ .
Table 4. 90% Bayesian credible intervals for μ and σ .
Lower Limit μ Upper Limit μ Lower Limit σ Upper Limit σ
−0.017951620.034037960.98149261.01901
Table 5. Monte Carlo estimates for lower limits.
Table 5. Monte Carlo estimates for lower limits.
T3 Months6 Months9 Months1 Year
K
0.8 s 0 1614.87651924.2072144.8912291.315
(8.543328)(14.27263)(20.09029)(20.35247)
0.9 s 0 1296.40031646.1551889.4792052.171
(7.794544)(13.55024)(19.41157)(19.73576)
s 0 1038.99991406.4441669.1711845.062
(6.909680)(12.77218)(18.69219)(19.10111)
1.1 s 0 832.30411201.5811474.6491662.490
(6.186688)(12.03232)(17.98552)(18.49396)
1.2 s 0 664.57701028.4631305.1081500.547
(5.841916)(11.40278)(17.32325)(17.94511)
Table 6. Monte Carlo estimates for standard prices.
Table 6. Monte Carlo estimates for standard prices.
T3 Months6 Months9 Months1 Year
K
0.8 s 0 1664.50922018.7302273.2482435.673
(9.616927)(17.23003)(25.58117)(26.54255)
0.9 s 0 1349.99411745.8842024.4932203.923
(8.793997)(16.40912)(24.78752)(25.82538)
s 0 1095.04761508.9481807.9922002.367
(7.814137)(15.51467)(23.94530)(25.08685)
1.1 s 0 888.87221304.8801615.3891823.577
(7.001017)(14.64918)(23.10734)(24.37146)
1.2 s 0 719.66151130.8751446.2861663.437
(6.609144)(13.90233)(22.31443)(23.71967)
Table 7. Monte Carlo estimates for upper limits.
Table 7. Monte Carlo estimates for upper limits.
T3 Months6 Months9 Months1 Year
K
0.8 s 0 1734.35392160.1232473.2062663.540
(11.278233)(22.56787)(36.63348)(40.43543)
0.9 s 0 1425.08321894.0892233.0742442.569
(10.355051)(21.61634)(35.69513)(39.60193)
s 0 1173.39861661.0702022.3592249.001
(9.248022)(20.56020)(34.68488)(38.73560)
1.1 s 0 968.18311458.1231832.5482076.176
(8.306346)(19.51948)(33.65647)(37.88148)
1.2 s 0 797.25681283.1321664.8711919.336
(7.834245)(18.59461)(32.66048)(37.08164)
Table 8. Differences in the Monte Carlo estimates for upper and lower limits.
Table 8. Differences in the Monte Carlo estimates for upper and lower limits.
T3 Months6 Months9 Months1 Year
K
0.8 s 0 119.4774235.9164328.3146372.2243
0.9 s 0 128.6830247.9340343.5949390.3978
s 0 134.3987254.6257353.1871403.9392
1.1 s 0 135.8790256.5423357.8990413.6860
1.2 s 0 132.6797254.6698359.7624418.7889
Table 9. Upper (U), point (P) and lower (L) estimates under the Bayesian BSM model.
Table 9. Upper (U), point (P) and lower (L) estimates under the Bayesian BSM model.
T3 Months6 Months9 Months1 Year
K
0.8 s 0 (U)1619.84481976.5002254.7822487.521
(M)1605.41981955.4362229.2452458.637
(L)1591.71161935.3522204.8542431.015
0.9 s 0 (U)1301.25941703.1992008.7442261.713
(M)1284.25161679.5551980.5442230.078
(L)1268.02341656.9631953.5702199.788
s 0 (U)1038.10281468.5301793.4412061.766
(M)1019.53821443.0831763.2342027.959
(L)1001.80091418.7491734.3211995.572
1.1 s 0 (U)824.02431267.5521604.8071884.203
(M)804.88231241.0091573.1771848.736
(L)786.60531215.6321542.9051814.758
1.2 s 0 (U)651.84961095.6351439.2561726.047
(M)632.94281068.6011406.7101689.364
(L)614.92881042.7791375.5781654.235
Table 10. Differences in the upper and lower estimates under the Bayesian BSM model.
Table 10. Differences in the upper and lower estimates under the Bayesian BSM model.
T3 Months6 Months9 Months1 Year
K
0.8 s 0 28.1331841.1477249.9278456.50651
0.9 s 0 33.2359746.2364755.1736161.92441
s 0 36.3019549.7812559.1197266.19463
1.1 s 0 37.4189551.9199561.9017969.44560
1.2 s 0 36.9207752.8559263.6779971.81156
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

Siu, T.K. Bayesian Lower and Upper Estimates for Ether Option Prices with Conditional Heteroscedasticity and Model Uncertainty. J. Risk Financial Manag. 2024, 17, 436. https://doi.org/10.3390/jrfm17100436

AMA Style

Siu TK. Bayesian Lower and Upper Estimates for Ether Option Prices with Conditional Heteroscedasticity and Model Uncertainty. Journal of Risk and Financial Management. 2024; 17(10):436. https://doi.org/10.3390/jrfm17100436

Chicago/Turabian Style

Siu, Tak Kuen. 2024. "Bayesian Lower and Upper Estimates for Ether Option Prices with Conditional Heteroscedasticity and Model Uncertainty" Journal of Risk and Financial Management 17, no. 10: 436. https://doi.org/10.3390/jrfm17100436

APA Style

Siu, T. K. (2024). Bayesian Lower and Upper Estimates for Ether Option Prices with Conditional Heteroscedasticity and Model Uncertainty. Journal of Risk and Financial Management, 17(10), 436. https://doi.org/10.3390/jrfm17100436

Article Metrics

Back to TopTop