We establish heavy-traffic limits for nearly deterministic queues, such as the G/D/n many-server queue. Since waiting times before starting service in the G/D/n queue are equivalent to waiting times in an associated G n /D/1 model, where the G n interarrival times are the sum of n consecutive interarrival times in the original model, we focus on the G n /D/1 model and the generalization to G n /G n /1, where ‘cyclic thinning’ is applied to both the arrival and service processes. We establish different limits in two cases: (i) when (1 − ρ n )√n → β as n → ∞ and (ii) when (1 − ρ n )n → β as n → ∞, where ρ n is the traffic intensity in model n. The nearly deterministic feature leads to interesting nonstandard scaling.

Nonspatial stochastic models of populations subject to catastrophic events result in the common conclusion that the survival probability of the population is nondecreasing with respect to the random number of individuals removed at each catastrophe. The purpose of this paper is to prove that such a monotonic relationship is not true for simple spatial models based on Harris' contact processes, whose dynamics are described by hypergraph structures rather than traditional graph structures. More precisely, it is proved that, for a wide range of parameters, the destruction of (infinite) hyperplanes does not affect the existence of a nontrivial invariant measure, whereas the destruction of large (finite) cubes drives the population to extinction, a result that we depict by using the biological picture: forest fires are more devastating than tornadoes. This indicates that the geometry of the subsets struck by catastrophes is somewhat more important than their area, thus the need to consider spatial rather than nonspatial models in this context.

Sequential decision problems can often be modeled as Markov decision processes. Classical solution approaches assume that the parameters of the model are known. However, model parameters are usually estimated and uncertain in practice. As a result, managers are often interested in how estimation errors affect the optimal solution. In this paper we illustrate how sensitivity analysis can be performed directly for a Markov decision process with uncertain reward parameters using the Bellman equations. In particular, we consider problems involving (i) a single stationary parameter, (ii) multiple stationary parameters, and (iii) multiple nonstationary parameters. We illustrate the applicability of this work through a capacitated stochastic lot-sizing problem.

In this work we provide sufficient conditions under which a general counting process stopped at a random time independent from the process belongs to the reliability decreasing reversed hazard rate (DRHR) or increasing failure rate (IFR) class. We also give some applications of these results in generalized renewal and trend renewal processes stopped at a random time.

We consider a risk model with threshold strategy, where the insurance company pays off a certain percentage of the income as dividend whenever the current surplus is larger than a given threshold. We investigate the ruin time, ruin probability, and the total dividend, using methods and results from queueing theory.

We consider a multiclass single-server queueing network as a model of a packet switching network. The rates packets are sent into this network are controlled by queues which act as congestion windows. By considering a sequence of congestion controls, we analyse a sequence of stationary queueing networks. In this asymptotic regime, the service capacity of the network remains constant and the sequence of congestion controllers act to exploit the network's capacity by increasing the number of packets within the network. We show that the stationary throughput of routes on this sequence of networks converges to an allocation that maximises aggregate utility subject to the network's capacity constraints. To perform this analysis, we require that our utility functions satisfy an exponential concavity condition. This family of utilities includes weighted α-fair utilities for α > 1.

We consider the first passage percolation problem on the random graph with vertex set N x {0, 1}, edges joining vertices at a Euclidean distance equal to unity, and independent exponential edge weights. We provide a central limit theorem for the first passage times ln between the vertices (0, 0) and (n, 0), thus extending earlier results about the almost-sure convergence of ln / n as n → ∞. We use generating function techniques to compute the n-step transition kernels of a closely related Markov chain which can be used to explicitly calculate the asymptotic variance in the central limit theorem.

In this paper we deal with an M/G/1 vacation system with the sojourn time (wait plus service) limit and two typical vacation rules, i.e. multiple and single vacation rules. Using the level crossing approach, explicit expressions for the steady-state distributions of the virtual waiting times are obtained in vacation systems with exponential and constant service times, a general vacation time, and two vacation rules.

Enhanced by the global financial crisis, the discussion about an accurate estimation of regulatory (risk) capital a financial institution needs to hold in order to safeguard against unexpected losses has become highly relevant again. The presence of heavy tails in combination with small sample sizes turns estimation at such extreme quantile levels into an inherently difficult statistical issue. We discuss some of the problems and pitfalls that may arise. In particular, based on the framework of second-order extended regular variation, we compare different high-quantile estimators and propose methods for the improvement of standard methods by focusing on the concept of penultimate approximations.

We consider an insurance model, where the underlying point process is a Cox process. Using a martingale approach applied to diffusion processes, finite-time Lundberg inequalities are obtained. By change-of-measure techniques, Cramér–Lundberg approximations are derived.

Optimal control of stochastic bandwidth-sharing networks is typically difficult. In order to facilitate the analysis, deterministic analogues of stochastic bandwidth-sharing networks, the so-called fluid models, are often taken for analysis, as their optimal control can be found more easily. The tracking policy translates the fluid optimal control policy back to a control policy for the stochastic model, so that the fluid optimality can be achieved asymptotically when the stochastic model is scaled properly. In this work we study the efficiency of the tracking policy, that is, how fast the fluid optimality can be achieved in the stochastic model with respect to the scaling parameter. In particular, our result shows that, under certain conditions, the tracking policy can be as efficient as feedback policies.

For a random walk with both downward and upward jumps (increments), the joint distribution of the exit time across a given level and the undershoot or overshoot at crossing is determined through its generating function, when assuming that the distribution of the jump in the direction making the exit possible has a Laplace transform which is a rational function. The expected exit time is also determined and the paper concludes with exact distribution results concerning exits from bounded intervals. The proofs use simple martingale techniques together with some classical expansions of polynomials and Rouché's theorem from complex function theory.

It is shown that the sample paths of Poisson shot noise with heavy-tailed semiexponential distributions satisfy a large deviation principle with a rate function that is insensitive to the shot shape. This demonstrates that, on the scale of large deviations, paths to rare events do not depend on the shot shape.

In this paper we analyze the so-called Parisian ruin probability, which arises when the surplus process stays below 0 longer than a fixed amount of time ζ > 0. We focus on a general spectrally negative Lévy insurance risk process. For this class of processes, we derive an expression for the ruin probability in terms of quantities that can be calculated explicitly in many models. We find its Cramér-type and convolution-equivalent asymptotics when reserves tend to ∞. Finally, we analyze some explicit examples.

Improved bounds on the copula of a bivariate random vector are computed when partial information is available, such as the values of the copula on a given subset of [0, 1]2, or the value of a functional of the copula, monotone with respect to the concordance order. These results are then used to compute model-free bounds on the prices of two-asset options which make use of extra information about the dependence structure, such as the price of another two-asset option.

In this paper we consider a single-server queue with Lévy input, and, in particular, its workload process (Qt)t≥0, focusing on its correlation structure. With the correlation function defined as r(t):= cov(Q0, Qt) / varQ0 (assuming that the workload process is in stationarity at time 0), we first study its transform ∫0∞r(t)e-ϑtdt, both for when the Lévy process has positive jumps and when it has negative jumps. These expressions allow us to prove that r(·) is positive, decreasing, and convex, relying on the machinery of completely monotone functions. For the light-tailed case, we estimate the behavior of r(t) for large t. We then focus on techniques to estimate r(t) by simulation. Naive simulation techniques require roughly (r(t))-2 runs to obtain an estimate of a given precision, but we develop a coupling technique that leads to substantial variance reduction (the required number of runs being roughly (r(t))-1). If this is augmented with importance sampling, it even leads to a logarithmically efficient algorithm.

In this paper we consider a structural form credit risk model with jumps. We investigate the credit spread, the price, and the fair premium of the zero-coupon bond for the proposed model. The price and the fair premium of the bond are associated with the Laplace transform of default time and the firm's expected present market value at default. We give sufficient conditions under which the Laplace transform and the expected present market value of a firm at default are twice continuously differentiable. We derive closed-form expressions for them when the jumps have a hyperexponential distribution. Using the closed-form expressions, we obtain numerical solutions for the default probability, the credit spread, and the fair premium of the bond.

We propose a method for estimating first passage time densities of one-dimensional diffusions via Monte Carlo simulation. Our approach involves a representation of the first passage time density as the expectation of a functional of the three-dimensional Brownian bridge. As the latter process can be simulated exactly, our method leads to almost unbiased estimators. Furthermore, since the density is estimated directly, a convergence of order 1 / √N, where N is the sample size, is achieved, which is in sharp contrast to the slower nonparametric rates achieved by kernel smoothing of cumulative distribution functions.

Suppose that N is a ℤ+-valued random variable and that X,X1,X2,… is a sequence of independent and identically distributed ℤ+ random variables independent of N. In this paper we are interested in properties of the conditional variable In particular, we want to know the asymptotic behavior of the conditional mean ENk or the conditional variance varNk as k→∞. We consider the cases when X is Poisson and when X is mixed Poisson. The problem is motivated by modeling loss reserves in nonlife insurance.

We study the leading term in the small-time asymptotics of at-the-money call option prices when the stock price process S follows a general martingale. This is equivalent to studying the first centered absolute moment of S. We show that if S has a continuous part, the leading term is of order √T in time T and depends only on the initial value of the volatility. Furthermore, the term is linear in T if and only if S is of finite variation. The leading terms for pure-jump processes with infinite variation are between these two cases; we obtain their exact form for stable-like small jumps. To derive these results, we use a natural approximation of S so that calculations are necessary only for the class of Lévy processes.