The aim of this thesis is a development of a new class of financial models with random parameters, which are computationally efficient and have the same level of performance as existing ones. In particular, this research is threefold. I have studied the evolution of storable commodity and commodity futures prices in time using a new random parameter model coupled with a Kalman filter. Such a combination allows one to forecast arbitrage-free futures prices and commodity spot prices one step ahead. Another direction of my research is a new volatility model, where the volatility is a random variable. The main advantage of this model is high calibration speed compared to the existing stochastic volatility models such as the Bates model or the Heston model. However, the performance of the new model is comparable to the latter. Comprehensive numerical studies demonstrate that the new model is a very competitive alternative to the Heston or the Bates model in terms of accuracy of matching option prices or computing hedging parameters. Finally, a new futures pricing model for electricity futures prices was developed. The new model has a random volatility parameter in its underlying process. The new model has less parameters, as compared to two-factor models for electricity commodity pricing with and without jumps. Numerical experiments with real data illustrate that it is quite competitive with the existing two-factor models in terms of pricing one step ahead futures prices, while being far simpler to calibrate. Further, a new heuristic for calibrating two-factor models was proposed. The new calibration procedure has two stages, offline and online. The offline stage calibrates parameters under a physical measure, while the online stage is used to calibrate the risk-neutrality parameters on each iteration of the particle filter. A particle filter was used to estimate the values of the underlying stochastic processes and to forecast futures prices one step ahead. The contributory material from two chapters of this thesis have been submitted to peer reviewed journals in terms of two papers: • Chapter 4: “A fast calibrating volatility model” has been submitted to the European Journal of Operational Research. • Chapter 5: “Electricity futures price models : calibration and forecasting” has been submitted to the European Journal of Operational Research.

The famous Black-Scholes formula provided the first mathematically sound mechanism to price financial options. It is based on the assumption, that daily random stock returns are identically normally distributed and hence stock prices follow a stochastic process with a constant volatility. Observed prices, at which options trade on the markets, don¡¯t fully support this hypothesis. Options corresponding to different strike prices trade as if they were driven by different volatilities. To capture this so-called volatility smile, we need a more sophisticated option-pricing model assuming that the volatility itself is a random process. The price we have to pay for this stochastic volatility model is that such models are computationally extremely intensive to simulate and hence difficult to fit to observed market prices. This difficulty has severely limited the use of stochastic volatility models in the practice. In this project we propose to overcome the obstacle of computational complexity by executing the simulations in a massively parallel fashion on the graphics processing unit (GPU) of the computer, utilizing its hundreds of parallel processors. We succeed in generating the trillions of random numbers needed to fit a monthly options contract in 3 hours on a desktop computer with a Tesla GPU. This enables us to accurately price any derivative security based on the same underlying stock. In addition, our method also allows extracting quantitative measures of the riskiness of the underlying stock that are implied by the views of the forward-looking traders on the option markets.

This study includes three issues about the dynamic of 30-days Taiwan Commercial Paper rate (CP2).The first issue focuses on the estimation of continuous-time short-term interest rate models. We discretize the continuous-time models by using two different approaches, and then use weekly and monthly data to estimate the parameters. The models are evaluated by data fit. We find that the estimated parameters are similar for different discretization approaches and would be more stable and efficient under quasi-maximum likelihood (QML) with weekly data. There exists mean reversion for Taiwan CP rate and the relationship between the volatility and the level of interest rates are less than 1 and smaller than that of American T-Bill rates reported by CKLS (1992) and Nowman (1997). We also find that CIR-SR model performs best for Taiwan CP rate. The second issue compares the continuous-time short-term interest rate models empirically both by predictive accuracy test and encompassing test. Having the estimated parameters of the models by discretization of Nowman(1997) and QML, we produce the forecasts on conditional mean and volatility for the interest rate over multiple-step-ahead horizons. The results indicate that the sophisticated models outperform the simpler models in the in-sample data fit, but have a distinct performance in the out-of-sample forecasting. The models equipped with mean reversion can produce better forecasts on conditional means during some period, and the heteroskedasticity variance model with outperform counterparts in volatility forecasting in some periods. The third issue concerns the persistent and massive volatility of short-term interest rates. This part inquires how the realizations on Taiwan short-term interest rates can be best described empirically. Various popular volatility specifications are estimated and tested. The empirical findings reveal that the mean reversion is an important characteristic for the Taiwan interest rates, and the level effect exists. Overall, the GARCH-L model fits well to Taiwan interest rates.

Ce travail de thèse poursuit une perspective double dans l'usage conjoint des méthodes de Monte Carlo séquentielles (MMS) et de l'algorithme Espérance-Maximisation (EM) dans le cadre des modèles de Markov cachés présentant une structure de dépendance markovienne d'ordre supérieur à 1 au niveau de la composante inobservée. Tout d'abord, nous commençons par un exposé succinct de l'assise théorique des deux concepts statistiques à travers les chapitres 1 et 2 qui leurs sont consacrés. Dans un second temps, nous nous intéressons à la mise en pratique simultanée des deux concepts au chapitre 3 et ce dans le cadre usuel où la structure de dépendance est d'ordre 1. L'apport des méthodes MMS dans ce travail réside dans leur capacité à approximer efficacement des fonctionnelles conditionnelles bornées, notamment des quantités de filtrage et de lissage dans un cadre non linéaire et non gaussien. Quant à l'algorithme EM, il est motivé par la présence à la fois de variables observables et inobservables (ou partiellement observées) dans les modèles de Markov Cachés et singulièrement les mdèles de volatilité stochastique étudié. Après avoir présenté aussi bien l'algorithme EM que les méthodes MCs ainsi que quelques unes de leurs propriétés dans les chapitres 1 et 2 respectivement, nous illustrons ces deux outils statistiques au travers de la calibration d'un modèle de volatilité stochastique. Cette application est effectuée pour des taux change ainsi que pour quelques indices boursiers au chapitre 3. Nous concluons ce chapitre sur un léger écart du modèle de volatilité stochastique canonique utilisé ainsi que des simulations de Monte Carlo portant sur le modèle résultant. Enfin, nous nous efforçons dans les chapitres 4 et 5 à fournir les assises théoriques et pratiques de l'extension des méthodes Monte Carlo séquentielles notamment le filtrage et le lissage particulaire lorsque la structure markovienne est plus prononcée. En guise d'illustration, nous donnons l'exemple d'un modèle de volatilité stochastique dégénéré dont une approximation présente une telle propriété de dépendance.

In this thesis, we propose a new and simple approach of extending the single-factor Heston stochastic volatility model to a more flexible one in solving option pricing problems.  In this approach, the volatility process for the underlying asset dynamics depends on the time to maturity of the option. As this idea is inspired by the Heath-Jarrow-Morton framework which models the evolution of the full dynamics of forward rate curves for various maturities, we name this approach as the HJM-type stochastic volatility (HJM-SV)  model. We conduct an empirical analysis by calibrating this model to real-market option data for underlying assets including an equity  (ABB stock) and a market index (EURO STOXX 50), for two separated time spans from Jan 2017 to Dec 2017 (before the COVID-19 pandemic) and from Nov 2019 to Nov 2020 (after the start of COVID-19 pandemic). We investigate the optimal way of dividing the set of option maturities into three classes, namely, the short-maturity, middle-maturity, and long-maturity classes. We calibrate our HJM-SV model to the data in the following way, for each class a single-factor Heston stochastic volatility model is calibrated to the corresponding market data. We address the question that how well the new HJM-SV model captures the feature of implied volatility surface given by the market data.

Merton’s Optimal Investment-Consumption Problem is a classic optimization problem in finance. It aims to find the optimal controls for a portfolio with both risky and risk-less assets, inorder to maximize an investor’s utility function. One of the controls is the optimal allocationof wealth invested in a risky asset and the other control is the consumption rate. The problemis solved by using Dynamic Programming and the related Hamilton-Jacobi-Bellman equation.One of the disadvantages of the original problem is the consideration of constant volatility. Inthis thesis, we extend Merton’s problem considering the Grzelak-Oosterlee-Van Veeren modelthat describes the dynamics of a risky asset with stochastic volatility and stochastic interestrate. We derive the related Hamilton-Jacobi-Bellman for Merton’s problem considering theGrzelak-Oosterlee-Van Veeren model. We simulate the controls from Merton’s problem intwo different cases, one case where the volatility and interest rate are stochastic, following theGOVV-model. In the other case, the volatility and interest rate are assumed to be constant, asin Merton’s problem. The results obtained from simulations show that the case with stochasticvolatility and interest gave the same results as the case where the volatility and the interest ratewere assumed to be constant.

In the commodity and energy markets, there are two kinds of risk that traders and analysts are concerned a lot about: multiple underlying risk and average price risk. Spread options, swaps and swaptions are widely used to hedge multiple underlying risks and Asian (average price) options can deal with average price risk. But when those two risks are combined together, then we need to consider Asian spread options and Asian-European spread options for hedging purposes. For an Asian or Asian-European spread call option, its payoff depends on the difference of two underlyings' average price or of one average price and one final (at expiration) price. Asian and Asian-European spread option pricing is challenging work. Even under the basic assumption that each underlying price follows a log-normal distribution, the average price does not have a distribution with a simple form. In this dissertation, for the first time, a systematic analysis of Asian spread option and Asian-European spread option pricing is proposed, several original approaches for the Black-Scholes-Merton model and a special stochastic volatility model are developed and some numerical computation tests are conducted as well.

Understanding financial asset return correlation is a key facet in asset allocation and investor’s portfolio optimization strategy. For the last decades, several studies have investigated this relationship between stock and bond returns. But, fewer studies have dealt with multi-asset return dynamics. While initial literature attempted to understand the fundamental pattern of comovements, later studies model the economic state variables influencing such time-varying comovements of primarily stock and bond returns. Research widely acknowledges that return distributions of financial assets are non-normal. When the joint distributions of the asset returns follow a non-elliptical structure, linear correlation fails to provide sufficient information of their dependence structure. In particular two issues arise from this existing empirical evidence. The first is to propose a more reliable alternative density specification for a higher-dimensional case. The second is to formulate a measure of the variables’ dependence structure which is more instructive than linear correlation. In this work I use a time-varying conditional multivariate elliptical and non-elliptical copula to examine the return comovements of three different asset classes: financial assets, commodities and real estate in the US market. I establish the following stylized facts about asset return comovements. First, the static measures of asset return comovements overestimate the asset return comovements in the economic expansion phase, while underestimating it in the periods of economic contraction. Second, Student t-copulas outperform both elliptical and non-elliptical copula models, thus confirming the ii dominance of Student t-distribution. Third, findings show a significant increase in asset return comovements post August 2007 subprime crisis ... [cont.].

This thesis is concerned with the small-time asymptotics and expansions of call option prices, when the log-return processes of the underlying stock prices follow several Levy-based models. To be specific, we derive the time-to-maturity asymptotic behavior for both at-the-money (ATM), out-of-the-money (OTM) and in-the-money (ITM) call-option prices under several jump-diffusion models and stochastic volatility models with Levy jumps. In the OTM and ITM cases, we consider a general stochastic volatility model with independent Levy jumps, while in the ATM case, we consider the pure-jump CGMY model with or without an independent Brownian component. An accurate modeling of the option market and asset prices requires a mixture of a continuous diffusive component and a jump component. In this thesis, we first model the log-return process of a risk asset with a jump diffusion model by combining a stochastic volatility model with an independent pure-jump Levy process. By assuming smoothness conditions on the Levy density away from the origin and a small-time large deviation principle on the stochastic volatility model, we derive the small-time expansions, of arbitrary polynomial order, in time-t, for the tail distribution of the log-return process, and for the call-option price which is not at-the-money. Moreover, our approach allows for a unified treatment of more general payoff functions. As a consequence of our tail expansions, the polynomial expansion in t of the transition density is also obtained under mild conditions. The asymptotic behavior of the ATM call-option prices is more complicated to obtain, and, in general, is given by fractional powers of t, which depends on different choices of the underlying log-return models. Here, we focus on the CGMY model, one of the most popular tempered stable models used in financial modeling. A novel second-order approximation for ATM option prices under the pure-jump CGMY Levy model is derived, and then extended to a model with an additional independent Brownian component. The third-order asymptotic behavior of the ATM option prices as well as the asymptotic behavior of the corresponding Black-Scholes implied volatilities are also addressed.

Understanding financial asset return correlation is a key facet in asset allocation and investor’s portfolio optimization strategy. For the last decades, several studies have investigated this relationship between stock and bond returns. But, fewer studies have dealt with multi-asset return dynamics. While initial literature attempted to understand the fundamental pattern of comovements, later studies model the economic state variables influencing such time-varying comovements of primarily stock and bond returns. Research widely acknowledges that return distributions of financial assets are non-normal. When the joint distributions of the asset returns follow a non-elliptical structure, linear correlation fails to provide sufficient information of their dependence structure. In particular two issues arise from this existing empirical evidence. The first is to propose a more reliable alternative density specification for a higher-dimensional case. The second is to formulate a measure of the variables’ dependence structure which is more instructive than linear correlation. In this work I use a time-varying conditional multivariate elliptical and non-elliptical copula to examine the return comovements of three different asset classes: financial assets, commodities and real estate in the US market. I establish the following stylized facts about asset return comovements. First, the static measures of asset return comovements overestimate the asset return comovements in the economic expansion phase, while underestimating it in the periods of economic contraction. Second, Student t-copulas outperform both elliptical and non-elliptical copula models, thus confirming the ii dominance of Student t-distribution. Third, findings show a significant increase in asset return comovements post August 2007 subprime crisis ... [cont.].