scispace - formally typeset
Search or ask a question

Showing papers in "Methodology and Computing in Applied Probability in 2022"




Journal ArticleDOI
TL;DR: In this article , the authors proposed several algorithms to solve McKean-Vlasov Forward Backward Stochastic Differential Equations (FBSDEs), which rely on the approximation power of neural networks to estimate the solution or gradient through minimization problems.
Abstract: We propose several algorithms to solve McKean-Vlasov Forward Backward Stochastic Differential Equations (FBSDEs). Our schemes rely on the approximating power of neural networks to estimate the solution or its gradient through minimization problems. As a consequence, we obtain methods able to tackle both mean-field games and mean-field control problems in moderate dimension. We analyze the numerical behavior of our algorithms on several multidimensional examples including non linear quadratic models.

4 citations


Journal ArticleDOI
TL;DR: In this article , the problem of the estimation of a parameter in Banach spaces, maximizing some criterion function which depends on an unknown nuisance parameter h, possibly infinite-dimensional, is considered.
Abstract: In the present paper, we consider the problem of the estimation of a parameter $$\varvec{\theta }$$ , in Banach spaces, maximizing some criterion function which depends on an unknown nuisance parameter h, possibly infinite-dimensional. The classical estimation methods are mainly based on maximizing the corresponding empirical criterion by substituting the nuisance parameter by a nonparametric estimator. We show that the M-estimators converge weakly to maximizers of Gaussian processes under rather general conditions. The conventional bootstrap method fails in general to consistently estimate the limit law. We show that the m out of n bootstrap, in this extended setting, is weakly consistent under conditions similar to those required for weak convergence of the M-estimators. The aim of this paper is therefore to extend the existing theory on the bootstrap of the M-estimators. Examples of applications from the literature are given to illustrate the generality and the usefulness of our results. Finally, we investigate the performance of the methodology for small samples through a short simulation study.

4 citations




Journal ArticleDOI
TL;DR: In this paper , the authors studied the optimal truncation level for a semiparametric tempered stable Lévy model and obtained a closed-form 2nd-order approximation of the optimal threshold in a high-frequency setting.
Abstract: Truncated realized quadratic variations (TRQV) are among the most widely used high-frequency-based nonparametric methods to estimate the volatility of a process in the presence of jumps. Nevertheless, the truncation level is known to critically affect its performance, especially in the presence of infinite variation jumps. In this paper, we study the optimal truncation level, in the mean-square error sense, for a semiparametric tempered stable Lévy model. We obtain a novel closed-form 2nd-order approximation of the optimal threshold in a high-frequency setting. As an application, we propose a new estimation method, which combines iteratively an approximate semiparametric method of moment estimator and TRQVs with the newly found small-time approximation for the optimal threshold. The method is tested via simulations to estimate the volatility and the Blumenthal-Getoor index of a generalized CGMY model and, via a localization technique, to estimate the integrated volatility of a Heston type model with CGMY jumps. Our method is found to outperform other alternatives proposed in the literature when working with a Lévy process (i.e., the volatility is constant), or when the index of jump intensity Y is larger than 3/2 in the presence of stochastic volatility.

2 citations


Journal ArticleDOI
TL;DR: In this article , it is shown that the optimal exercise times are the first times at which the underlying risky asset price process reaches some lower or upper stochastic boundaries depending on the current values of its running maximum or minimum as well as the maximum drawdown or maximum drawup.
Abstract: Abstract We present closed-form solutions to the problems of pricing of the perpetual American double lookback put and call options on the maximum drawdown and the maximum drawup with floating strikes in the Black-Merton-Scholes model. It is shown that the optimal exercise times are the first times at which the underlying risky asset price process reaches some lower or upper stochastic boundaries depending on the current values of its running maximum or minimum as well as the maximum drawdown or maximum drawup. The proof is based on the reduction of the original double optimal stopping problems to the appropriate sequences of single optimal stopping problems for the three-dimensional continuous Markov processes. The latter problems are solved as the equivalent free-boundary problems by means of the smooth-fit and normal-reflection conditions for the value functions at the optimal stopping boundaries and the edges of the three-dimensional state spaces. We show that the optimal exercise boundaries are determined as either the unique solutions of the associated systems of arithmetic equations or the minimal and maximal solutions of the appropriate first-order nonlinear ordinary differential equations.

2 citations


Journal ArticleDOI
TL;DR: In this article , the authors extended the results to the case when a final penalty P is taken into consideration as well besides a proportional cost $$k>1$$ for capital injections, which requires amending the scale and Gerber-Shiu functions already introduced in Gajek and Kucinsky (Insur Math Econ 73:1-19, 2017).
Abstract: The recent papers Gajek and Kucinsky (Insur Math Econ 73:1–19, 2017) and Avram et al. (Mathematics 9(9):931, 2021) cost induced dichotomy for optimal dividends in the cramr-lundberg model. Avram et al. (Mathematics 9(9):931, 2021) investigated the control problem of optimizing dividends when limiting capital injections stopped upon bankruptcy. The first paper works under the spectrally negative Lévy model; the second works under the Cramér-Lundberg model with exponential jumps, where the results are considerably more explicit. The current paper has three purposes. First, it illustrates the fact that quite reasonable approximations of the general problem may be obtained using the particular exponential case studied in Avram et al. cost induced dichotomy for optimal dividends in the Cramér-Lundberg model (Avram et al. in Mathematics 9(9):931, 2021). Secondly, it extends the results to the case when a final penalty P is taken into consideration as well besides a proportional cost $$k>1$$ for capital injections. This requires amending the “scale and Gerber-Shiu functions” already introduced in Gajek and Kucinsky (Insur Math Econ 73:1–19, 2017). Thirdly, in the exponential case, the results will be made even more explicit by employing the Lambert-W function. This tool has particular importance in computational aspects and can be employed in theoretical aspects such as asymptotics.

2 citations



Journal ArticleDOI
TL;DR: In this paper , it is shown that the optimal stopping times are the first times at which the underlying process reaches some lower or upper stochastic boundaries depending on the current values of its running maximum or minimum.
Abstract: Abstract We present closed-form solutions to some double optimal stopping problems with payoffs representing linear functions of the running maxima and minima of a geometric Brownian motion. It is shown that the optimal stopping times are th first times at which the underlying process reaches some lower or upper stochastic boundaries depending on the current values of its running maximum or minimum. The proof is based on the reduction of the original double optimal stopping problems to sequences of single optimal stopping problems for the resulting three-dimensional continuous Markov process. The latter problems are solved as the equivalent free-boundary problems by means of the smooth-fit and normal-reflection conditions for the value functions at the optimal stopping boundaries and the edges of the three-dimensional state space. We show that the optimal stopping boundaries are determined as the extremal solutions of the associated first-order nonlinear ordinary differential equations. The obtained results are related to the valuation of perpetual real double lookback options with floating sunk costs in the Black-Merton-Scholes model.




Journal ArticleDOI
TL;DR: In this article , the authors focus on estimating reach probability of a closed unsafe set by a stochastic process and prove that the variance of the multi-level MC estimated reach probability under fixed assignment splitting is smaller or equal than under random assignment splitting methods.
Abstract: Abstract This paper focuses on estimating reach probability of a closed unsafe set by a stochastic process. A well-developed approach is to make use of multi-level MC simulation, which consists of encapsulating the unsafe set by a sequence of increasing closed sets and conducting a sequence of MC simulations to estimate the reach probability of each inner set from the previous set. An essential step is to copy (split) particles that have reached the next level (inner set) prior to conducting a MC simulation to the next level. The aim of this paper is to prove that the variance of the multi-level MC estimated reach probability under fixed assignment splitting is smaller or equal than under random assignment splitting methods. The approaches are illustrated for a geometric Brownian motion example.

Journal ArticleDOI
TL;DR: In this paper , the pool-adjacent-violators algorithm (PAVA) is modified such that the computation times are reduced substantially by studying the dependence of weighted least squares fits on the response vector to be approximated.
Abstract: Abstract In the context of estimating stochastically ordered distribution functions, the pool-adjacent-violators algorithm (PAVA) can be modified such that the computation times are reduced substantially. This is achieved by studying the dependence of antitonic weighted least squares fits on the response vector to be approximated.



Journal ArticleDOI
TL;DR: In this article , the authors investigated the valuation of annuity guarantees under a regime-switching model when the dynamics of the underlying stock price follow a self-exciting switching jump-diffusion process.
Abstract: This article investigates the valuation of annuity guarantees under a regime-switching model when the dynamics of the underlying stock price follow a self-exciting switching jump-diffusion process. In this framework, we add a jump component to a regime-switching geometric Brownian for large shocks on the stock price. The intensity of shock arrivals is a Hawkes process modulated by a continuous time hidden Markov chain with a finite number of states. The interest rate used for discounting is stochastic and correlated to the stock market. In an incomplete market, we define an equivalent martingale measure to price a variable annuity contract that guarantees a minimum living or death benefit. Under this equivalent martingale measure, we propose closed-form approximation formulas using the inverse Fourier transform technique. A numerical implementation highlights the impact of self-exciting jumps and economic regimes on the valuation of guarantees.


Journal ArticleDOI
TL;DR: In this paper , the distribution of the virtual and the actual waiting time in a GI/M/c queue with impatient customers is derived under the assumption that the distribution is Coxian.
Abstract: Abstract There are three contributions to the literature in this note. Firstly, we point out that under some weak conditions the result in Swensen (1986) on the remaining loads in a GI/M/ c queue with impatient customers, derived under the assumption that the distribution of inter-arrival times is Coxian, is also valid for the more general phase type distribution. From this result the distributions of the virtual and the actual waiting time can easily be obtained. Secondly, the relation to an alternative expression for the distribution of the virtual waiting time derived by Kawanishi and Takine (2015) is also clarified. Finally, we explain how the results on the remaining loads can be used to find distributions describing particular fixed servers and provide a couple of numerical examples of how this can be done.





Journal ArticleDOI
TL;DR: In this paper , a framework for statistical causality detection in general classes of multivariate nonlinear time series models has been proposed, where causality may be present in either: trend, volatility or both structural components of the general multivariate Markov processes under study.
Abstract: Abstract The ability to test for statistical causality in linear and nonlinear contexts, in stationary or non-stationary settings, and to identify whether statistical causality influences trend of volatility forms a particularly important class of problems to explore in multi-modal and multivariate processes. In this paper, we develop novel testing frameworks for statistical causality in general classes of multivariate nonlinear time series models. Our framework accommodates flexible features where causality may be present in either: trend, volatility or both structural components of the general multivariate Markov processes under study. In addition, we accommodate the added possibilities of flexible structural features such as long memory and persistence in the multivariate processes when applying our semi-parametric approach to causality detection. We design a calibration procedure and formal testing procedure to detect these relationships through classes of Gaussian process models. We provide a generic framework which can be applied to a wide range of problems, including partially observed generalised diffusions or general multivariate linear or nonlinear time series models. We demonstrate several illustrative examples of features that are easily testable under our framework to study the properties of the inference procedure developed including the power of the test, sensitivity and robustness. We then illustrate our method on an interesting real data example from commodity modelling.

Journal ArticleDOI
TL;DR: In this paper , the authors studied the limiting behavior of the maximum of a bivariate (finite or infinite) moving average model, based on discrete random variables, assuming that the bivariate distribution of the iid innovations belong to the Anderson's class.
Abstract: We study the limiting behaviour of the maximum of a bivariate (finite or infinite) moving average model, based on discrete random variables. We assume that the bivariate distribution of the iid innovations belong to the Anderson's class (Anderson, 1970). The innovations have an impact on the random variables of the INMA model by binomial thinning. We show that the limiting distribution of the bivariate maximum is also of Anderson's class, and that the components of the bivariate maximum are asymptotically independent.

Journal ArticleDOI
TL;DR: In this paper , a fractional degree reference dependent stochastic dominance rule is developed which is a generalization of the integer degree reference-dependent Stochastic Dominance Rule.
Abstract: For addressing the Allis-type anomalies, a fractional degree reference dependent stochastic dominance rule is developed which is a generalization of the integer degree reference dependent stochastic dominance rules. This new rule can effectively explain why the risk comparison does not satisfy translational invariance and scaling invariance in some cases. The rule also has a good property that it is compatible with the endowment effect of risk. This rule can help risk-averse but not absolute risk-averse decision makers to compare risks relative to reference points. We present some tractable equivalent integral conditions for the fractional degree reference dependent stochastic dominance rule, as well as some practical applications for the rule in economics and finance.

Journal ArticleDOI
TL;DR: In this paper , a fractional degree reference dependent stochastic dominance rule is developed which is a generalization of the integer degree reference-dependent Stochastic Dominance Rule.
Abstract: For addressing the Allis-type anomalies, a fractional degree reference dependent stochastic dominance rule is developed which is a generalization of the integer degree reference dependent stochastic dominance rules. This new rule can effectively explain why the risk comparison does not satisfy translational invariance and scaling invariance in some cases. The rule also has a good property that it is compatible with the endowment effect of risk. This rule can help risk-averse but not absolute risk-averse decision makers to compare risks relative to reference points. We present some tractable equivalent integral conditions for the fractional degree reference dependent stochastic dominance rule, as well as some practical applications for the rule in economics and finance.

Journal ArticleDOI
TL;DR: In this paper , the authors define a complete framework of Hawkes processes with a Gamma density excitation function (i.e. estimation, simulation, goodness-of-fit) instead of an exponential-decaying function and demonstrate some mathematical properties about the transient regime of the process.
Abstract: Hawkes processes are temporal self-exciting point processes. They are well established in earthquake modelling or finance and their application is spreading to diverse areas. Most models from the literature have two major drawbacks regarding their potential application to insurance. First, they use an exponentially-decaying form of excitation, which does not allow a delay between the occurrence of an event and its excitation effect on the process and does not fit well on insurance data consequently. Second, theoretical results developed from these models are valid only when time of observation tends to infinity, whereas the time horizon for an insurance use case is of several months or years. In this paper, we define a complete framework of Hawkes processes with a Gamma density excitation function (i.e. estimation, simulation, goodness-of-fit) instead of an exponential-decaying function and we demonstrate some mathematical properties (i.e. expectation, variance) about the transient regime of the process. We illustrate our results with real insurance data about natural disasters in Luxembourg.