63 research outputs found
Simulation of multivariate diffusion bridge
We propose simple methods for multivariate diffusion bridge simulation, which
plays a fundamental role in simulation-based likelihood and Bayesian inference
for stochastic differential equations. By a novel application of classical
coupling methods, the new approach generalizes a previously proposed simulation
method for one-dimensional bridges to the multi-variate setting. First a method
of simulating approximate, but often very accurate, diffusion bridges is
proposed. These approximate bridges are used as proposal for easily
implementable MCMC algorithms that produce exact diffusion bridges. The new
method is much more generally applicable than previous methods. Another
advantage is that the new method works well for diffusion bridges in long
intervals because the computational complexity of the method is linear in the
length of the interval. In a simulation study the new method performs well, and
its usefulness is illustrated by an application to Bayesian estimation for the
multivariate hyperbolic diffusion model.Comment: arXiv admin note: text overlap with arXiv:1403.176
Point processes with finite-dimensional conditional probabilities
AbstractWe study the structure of point processes N with the property that the P(θtNâ¡|Ft) vary in a finite-dimensional space where θt is the shift and Ft the Ď-field generated by the counting process up to time t. This class of point processes is strictly larger than Neutsâ class of Markovian arrival processes. On the one hand, it allows for more general features like interarrival distributions which are matrix-exponential rather than phase type, on the other the probabilistic interpretation is a priori less clear. Nevertheless, the properties are very similar. In particular, finite-dimensional distributions of interarrival times, moments, Laplace transforms, Palm distributions, etc., are shown to be given by two fundamental matrices C,D just as for the Markovian arrival process. We also give a probabilistic interpretation in terms of a piecewise deterministic Markov process on a compact convex subset of Rp, whose jump times are identical to the epochs of N
Mortality modeling and regression with matrix distributions
In this paper we investigate the flexibility of matrix distributions for the
modeling of mortality. Starting from a simple Gompertz law, we show how the
introduction of matrix-valued parameters via inhomogeneous phase-type
distributions can lead to reasonably accurate and relatively parsimonious
models for mortality curves across the entire lifespan. A particular feature of
the proposed model framework is that it allows for a more direct interpretation
of the implied underlying aging process than some previous approaches.
Subsequently, towards applications of the approach for multi-population
mortality modeling, we introduce regression via the concept of proportional
intensities, which are more flexible than proportional hazard models, and we
show that the two classes are asymptotically equivalent. We illustrate how the
model parameters can be estimated from data by providing an adapted EM
algorithm for which the likelihood increases at each iteration. The practical
feasibility and competitiveness of the proposed approach are illustrated for
several sets of mortality data
Multivariate matrix-exponential distributions
We review what is currently known about one-dimensional distributions on
the non-negative reals with rational Laplace transform, also known as
matrix-exponential distributions. In particular we discuss a flow
interpreation which enables one to mimic certain probabilisticly
inspired arguments which are known from the theory of phase-type distributions.
We then move on to present ongoing research for higher dimensions.
We discuss a characterization result, some closure properties, and
a number of examples. Finally we present open problems and future
perspectives
On the construction of bivariate exponential distributions with an arbitrary correlation coefficient
Aggregate Markov models in life insurance: properties and valuation
In multi-state life insurance, an adequate balance between analytic
tractability, computational efficiency, and statistical flexibility is of great
importance. This might explain the popularity of Markov chain modelling, where
matrix analytic methods allow for a comprehensive treatment. Unfortunately,
Markov chain modelling is unable to capture duration effects, so this paper
presents aggregate Markov models as an alternative. Aggregate Markov models
retain most of the analytical tractability of Markov chains, yet are
non-Markovian and thus more flexible. Based on an explicit characterization of
the fundamental martingales, matrix representations of the expected accumulated
cash flows and corresponding prospective reserves are derived for
duration-dependent payments with and without incidental policyholder behaviour.
Throughout, special attention is given to a semi-Markovian case. Finally, the
methods and results are illustrated in a numerical example
An Alternative Formula to Price American Options.
We give a new way to price American options, using Samuelson´s formula. We first obtain the option price corresponding to a European option at time t, weighting it by the probability that the underlying asset takes the value S at time t. This factor is given by the solution of the Fokker-Planck (Kolmogorov) equation for the transition probability density. The main advantage of this approach is that we can introduce systematically the effect of macroeconomic factors. If a macroeconomic framework is given by a dynamic system in the form of a set of ordinary differential equations we only have to solve a partial differential equation, for the transition probability density. In this context, we verify, for the sake of consistency, that this formula is consistent with the Black-Scholes model.American options, Fokker-Planck, Black-Scholes, Samuelson, density probability function.
- âŚ