Explosion time for some Laplace transforms of the Wishart process

In this article, we focus upon a family of matrix valued stochastic processes and study the problem of determining the smallest time such that their Laplace transforms become infinite. In particular, we concentrate upon the class of Wishart processes, which have proved to be very useful in different applications by their ability in describing non-trivial dependence. Thanks to this remarkable property we are able to explain the behavior of the explosion times for the Laplace transforms of the Wishart process and its time integral in terms of the relative importance of the involved factors and their correlations.     

Laplace transform approach to option pricing for time-changed Brownian models

In this article, we consider European option pricing for time-changed Brownian models using Laplace transform. We obtain a general formula for the option price as the integral of a real-valued function involving the Laplace transform of the random time change. Unlike the usual Fourier transform technique, our method does not suffer from difficulties specific to complex integration, such as the evaluation of multiple-valued functions, and allows for a model-independent analysis of the truncation error. In the numerical analysis part, we compare option prices in variance gamma (VG), normal inverse Gaussian (NIG), and generalized hyperbolic (GH) models obtained by Laplace transform with those obtained by the Fourier transform method introduced by Carr and Madan in 1999. The results show that our method converges faster than the Fourier approach when the Laplace transforms of the subordinators decay exponentially, for examples like NIG and GH models.     

A Note on Continuous Spatial-Temporal Dynamics of Stochastic Processes

We propose a general form to analyze the space-time interdependency of continuous space-time stochastic processes. We present a new space-time approach based on the intensity function of the underlying point process. These formulations can be, to some extent, analytically solved to obtain explicit formulae of interest. We define a general function that controls the space-time interaction and allows for closed forms depending on the particular choice of several mathematical tools playing a role in this interaction function. In particular, we make use of copulas and Laplace transforms to provide interesting examples of the dynamics of the random intensity function and, in turn, of the number of points contained in a given region.     

The discrete delta and nabla Mittag-Leffler distributions

In this paper, we extend Bernstein theorem by using basic tools of calculus on time scales, and, as a further application of it, the discrete nabla and delta Mittag-Leffler distributions are introduced here with respect to their Laplace transforms on the discrete time scale. For these discrete distributions, infinite divisibility and geometric infinite divisibility are proved along with some statistical properties. The delta and nabla Mittag-Leffler processes are defined.     

Ruin Probabilities for the Perturbed Compound Poisson Risk Process with Investment

In this article, we consider the perturbed compound Poisson risk process with investment incomes. The risk reserve process is perturbed by an independent Brownian motion and the surplus is invested at a constant force of interest. We investigate the asymptotic behavior of the ruin probability as the initial reserve goes to infinity. Bounds and time-dependent bounds are derived for the ultimate ruin probability and the probabilities of ruin within finite time, respectively. We also obtain an explicit expression for the Laplace transform of the ultimate ruin probability.     

A note on the passage time of finite-state Markov chains

Consider a Markov chain with finite state {0, 1, …, d}. We give the generation functions (or Laplace transforms) of absorbing (passage) time in the following two situations: (1) the absorbing time of state d when the chain starts from any state i and absorbing at state d; (2) the passage time of any state i when the chain starts from the stationary distribution supposed the chain is time reversible and ergodic. Example shows that it is more convenient compared with the existing methods, especially we can calculate the expectation of the absorbing time directly.     

The Pearson Diffusions: A Class of Statistically Tractable Diffusion Processes

Abstract.  The Pearson diffusions form a flexible class of diffusions defined by having linear drift and quadratic squared diffusion coefficient. It is demonstrated that for this class explicit statistical inference is feasible. A complete model classification is presented for the ergodic Pearson diffusions. The class of stationary distributions equals the full Pearson system of distributions. Well-known instances are the Ornstein–Uhlenbeck processes and the square root (CIR) processes. Also diffusions with heavy-tailed and skew marginals are included. Explicit formulae for the conditional moments and the polynomial eigenfunctions are derived. Explicit optimal martingale estimating functions are found. The discussion covers GMM, quasi-likelihood, non-linear weighted least squares estimation and likelihood inference too. The analytical tractability is inherited by transformed Pearson diffusions, integrated Pearson diffusions, sums of Pearson diffusions and Pearson stochastic volatility models. For the non-Markov models, explicit optimal prediction-based estimating functions are found. The estimators are shown to be consistent and asymptotically normal.     

A new stochastic order based upon Laplace transform with applications

This paper deals with nonnegative random variables having Laplace transforms as their reliability functions. We study a new stochastic order based upon Laplace transform. Some applications in actuarial science, frailty models and reliability are presented as well.     

SOJOURN TIMES IN NON-HOMOGENEOUS QBD PROCESSES WITH PROCESSOR SHARING

We study sojourn times of customers in a processor sharing queue with a service rate that varies over time, depending on the number of customers and on the state of a random environment. An explicit expression is derived for the Laplace–Stieltjes transform of the sojourn time conditional on the state upon arrival and the amount of work brought into the system. Particular attention is paid to the conditional mean sojourn time of a customer as a function of his required amount of work, and we establish the existence of an asymptote as the amount of work tends to infinity. The method of random time change is then extended to include the possibility of a varying service rate. By means of this method, we explain the well-established proportionality between the conditional mean sojourn time and required amount of work in processor sharing queues without random environment. Based on numerical experiments, we propose an approximation for the conditional mean sojourn time. Although first presented for exponentially distributed service requirements, the analysis is shown to extend to phase-type services. The service discipline of discriminatory processor sharing is also shown to fall within the framework.     

A method for computing the autocovariance of renewal processes

In this paper we derive formulae for the autocovariance functions of renewal and renewal reward processes. The derivation is based on a Poissonization technique of a renewal process. The formulae are expressed in the form of Laplace transforms. In some cases we may invert the Laplace transforms analytically, but in general we have to invert them numerically.     

Estimation of the scale parameter of the half-logistic distribution with multiply type II censored sample

In this paper, we consider the maximum likelihood and Bayes estimation of the scale parameter of the half-logistic distribution based on a multiply type II censored sample. However, the maximum likelihood estimator(MLE) and Bayes estimator do not exist in an explicit form for the scale parameter. We consider a simple method of deriving an explicit estimator by approximating the likelihood function and discuss the asymptotic variances of MLE and approximate MLE. Also, an approximation based on the Laplace approximation (Tierney & Kadane, 1986) is used to obtain the Bayes estimator. In order to compare the MLE, approximate MLE and Bayes estimates of the scale parameter, Monte Carlo simulation is used.     

Goodness‐of‐Fit based on Downsampling with Applications to Linear Drift Diffusions

Abstract. A goodness‐of‐fit test for continuous‐time models is developed that examines if the parameter estimates are consistent with another for different sampling frequencies. The test compares parameter estimates obtained from estimating functions for downsamples of the data. We prove asymptotic results for stationary and ergodic processes, and apply the downsampling test to linear drift diffusions. Simulations indicate that the test is quite powerful in detecting non‐Markovian deviations from the linear drift diffusions.     

On the probability of ruin in a continuous risk model with two types of delayed claims

ABSTRACT

This article studies a risk model involving one type of main claims and two types of by-claims, which is an extension of the general risk model with delayed claims. We suppose that every main claim may not induce any by-claims or may induce one by-claim belonging to one of the two types of by-claims with a certain probability. In addition, assume that the by-claim and its associated main claim may occur at the same time and that the occurrence of the by-claim may be delayed. An integro-differential equation system for survival probabilities is derived by using two auxiliary risk models. The expression of the survival probability is obtained by applying Laplace transforms and Rouché theorem. Furthermore, we provide a method for solving the survival probability when the two by-claim amounts satisfy different exponential distributions. As a special case, an explicit expression of survival probability is given when all the claim amounts obey the same exponential distribution. Finally, numerical results are provided to examine the proposed method.     

On the umvu estimators after using a normalizing transformation

The conditional distribution given complete sufficient statistics is used along with the Rao-Blackwell theorem to obtain uniformly minimum variance unbiased (UMVU) estimators after a transformation to normality has been applied to data. The estimators considered are for the mean, the variance and the cumulative distribution of the original non-normal data. Previous procedures to obtain UMVU estimators have used Laplace transforms, Taylor expansions and the jackknife. An integration method developed in this paper requires only integrability of the normalizing transformation function. This method is easy to employ and it is always possible to obtain a numerical result.     

Inversion of the Renewal Density with Dead-time

Stationary renewal point processes are defined by the probability distribution of the distances between successive points (lifetimes) that are independent and identically distributed random variables. For some applications it is also interesting to define the properties of a renewal process by using the renewal density. There are well-known expressions of this density in terms of the probability density of the lifetimes. It is more difficult to solve the inverse problem consisting in the determination of the density of the lifetimes in terms of the renewal density. Theoretical expressions between their Laplace transforms are available but the inversion of these transforms is often very difficult to obtain in closed form. We show that this is possible for renewal processes presenting a dead-time property characterized by the fact that the renewal density is zero in an interval including the origin. We present the principle of a recursive method allowing the solution of this problem and we apply this method to the case of some processes with input dead-time. Computer simulations on Poisson and Erlang (2) processes show quite good agreement between theoretical calculations and experimental measurements on simulated data.     

Lévy Processes,Phase-Type Distributions,and Martingales

Lévy processes are defined as processes with stationary independent increments and have become increasingly popular as models in queueing, finance, etc.; apart from Brownian motion and compound Poisson processes, some popular examples are stable processes, variance gamma processes, CGMY Lévy processes (tempered stable processes), NIG (normal inverse Gaussian) Lévy processes, and hyperbolic Lévy processes. We consider here a dense class of Lévy processes, compound Poisson processes with phase-type jumps in both directions and an added Brownian component. Within this class, we survey how to explicitly compute a number of quantities that are traditionally studied in the area of Lévy processes, in particular two-sided exit probabilities and associated Laplace transforms, the closely related scale function, one-sided exit probabilities and associated Laplace transforms coming up in queueing problems, and similar quantities for a Lévy process with reflection in 0. The solutions are in terms of roots to polynomials, and the basic equations are derived by purely probabilistic arguments using martingale optional stopping; a particularly useful martingale is the so-called Kella-Whitt martingale. Also, the relation to fluid models with a Brownian component is discussed.     

An RG-Factorization Approach for a BMAP/M/1 Generalized Processor-Sharing Queue

ABSTRACT

In this paper, we study a BMAP/M/1 generalized processor-sharing queue. We propose an RG-factorization approach, which can be applied to a wider class of Markovian block-structured processor-sharing queues. We obtain the expressions for both the distribution of the stationary queue length and the Laplace transform of the sojourn time distribution. From these two expressions, we develop an algorithm to compute the mean and variance of the sojourn time approximately.     

PIECEWISE POLYNOMIAL APPROXIMATIONS FOR HEAVY-TAILED DISTRIBUTIONS IN QUEUEING ANALYSIS

ABSTRACT

A basic difficulty in dealing with heavy-tailed distributions is that they may not have explicit Laplace transforms. This makes numerical methods that use the Laplace transform more challenging. This paper generalizes an existing method for approximating heavy-tailed distributions, for use in queueing analysis. The generalization involves fitting Chebyshev polynomials to a probability density function g(t) at specified points t ₁, t ₂, …, t _N . By choosing points t _i , which rapidly get far out in the tail, it is possible to capture the tail behavior with relatively few points, and to control the relative error in the approximation. We give numerical examples to evaluate the performance of the method in simple queueing problems.     

Sojourn time in an M/M/1 processor sharing queue with permanent customers

We explicitly compute the sojourn time distribution of an arbitrary customer in an M/M/1 processor sharing (PS) queue with permanent customers. We notably exhibit the orthogonal structure associated with this queuing system and we show how sieved Pollaczek polynomials and their associated orthogonality measure can be used to obtain an explicit representation for the complementary cumulative distribution function of the sojourn time of a customer. This explicit formula subsequently allows us to compute the two first moments of this random variable and to study the asymptotic behavior of its distribution. The most salient result is that the decay rate depends on the load of the system and the number K of permanent customers. When the load is above a certain threshold depending on K, the decay rate is identical to that of a regular M/M/1 PS queue.     

Laplace approximations to posterior moments and marginal distributions on circles,spheres, and cylinders

We extend recent work on Laplace approximations (Tierney and Kadane 1986; Tierney, Kass, and Kadane 1989) from parameter spaces that are subspaces of R^k to those that are on circles, spheres, and cylinders. While such distributions can be mapped onto the real line (for example, a distribution on the circle can be thought of as a function of an angle θ, 0 ? 0 ? 2π), that the end points coincide is not a feature of the real line, and requires special treatment. Laplace approximations on the real line make essential use of the normal integral in both the numerator and the denominator. Here that role is played by the von Mises integral on the circle, by the Bingham integrals on the spheres and hyperspheres, and by the normal-von Mises and normal-Bingham integrals on the cylinders and hypercylinders, respectively. We begin with a brief introduction to Laplace approximations and to previous Bayesian work on circles, spheres, and cylinders. We then develop the theory for parameter spaces that are hypercylinders, since all other shapes considered here are special cases. We compute some examples, which show reasonable accuracy even for small samples.