Search

In this paper we obtain some equivalent conditions and sufficient conditions for the local and nonlocal asymptotics of the φ-moments of the overshoot and undershoot of a random walk, where φ is a nonnegative, long-tailed function. By the strong Markov property, it can be shown that the moments of the overshoot and undershoot and the moments of the first ascending ladder height of a random walk satisfy some renewal equations. Therefore, in this paper we first investigate the local and nonlocal asymptotics for the moments of the first ascending ladder height of a random walk, and then give some equivalent conditions and sufficient conditions for the asymptotics of the solutions to some renewal equations. Using the above results, the main results of this paper are obtained.

In this paper we investigate the ‘local’ properties of a random mapping model, T n D̂, which maps the set {1, 2, …, n} into itself. The random mapping T n D̂ , which was introduced in a companion paper (Hansen and Jaworski (2008)), is constructed using a collection of exchangeable random variables D̂ 1, …, D̂ n which satisfy In the random digraph, G n D̂ , which represents the mapping T n D̂ , the in-degree sequence for the vertices is given by the variables D̂ 1, D̂ 2, …, D̂ n , and, in some sense, G n D̂ can be viewed as an analogue of the general independent degree models from random graph theory. By local properties we mean the distributions of random mapping characteristics related to a given vertex v of G n D̂ - for example, the numbers of predecessors and successors of v in G n D̂ . We show that the distribution of several variables associated with the local structure of G n D̂ can be expressed in terms of expectations of simple functions of D̂ 1, D̂ 2, …, D̂ n . We also consider two special examples of T n D̂ which correspond to random mappings with preferential and anti-preferential attachment, and determine, for these examples, exact and asymptotic distributions for the local structure variables considered in this paper. These distributions are also of independent interest.

Consider a random graph, having a pre-specified degree distribution F but other than that being uniformly distributed, describing the social structure (friendship) in a large community. Suppose one individual in the community is externally infected by an infectious disease and that the disease has its course by assuming that infected individuals infect their not yet infected friends independently with probability p. For this situation the paper determines R_0 and tau_0, the basic reproduction number and the asymptotic final size in case of a major outbreak. Further, the paper looks at some different local vaccination strategies where individuals are chosen randomly and vaccinated, or friends of the selected individuals are vaccinated, prior to the introduction of the disease. For the studied vaccination strategies the paper determines R_v: the reproduction number, and tau_v: the asymptotic final proportion infected in case of a major outbreak, after vaccinating a fraction v., 31 pages, 3 figures

Consider a sequence X = (Xn: n ≥ 1) of independent and identically distributed random variables, and an independent geometrically distributed random variable M with parameter p. The random variable SM = X1 + ∙ ∙ ∙ + XM is called a geometric sum. In this paper we obtain asymptotic expansions for the distribution of SM as p ↘ 0. If EX1 > 0, the asymptotic expansion is developed in powers of p and it provides higher-order correction terms to Renyi's theorem, which states that P(pSM > x) ≈ exp(-x/EX1). Conversely, if EX1 = 0 then the expansion is given in powers of √p. We apply the results to obtain corrected diffusion approximations for the M/G/1 queue. These expansions follow in a unified way as a consequence of new uniform renewal theory results that are also developed in this paper.

A large dam model is the object of study of this paper. The parametersLlowerandLupperdefine its lower and upper levels,L=Lupper-Lloweris large, and if the current level of water is between these bounds, the dam is assumed to be in a normal state. Passage across one or other of the levels leads to damage. LetJ1andJ2denote the damage costs of crossing the lower and, respectively, the upper levels. It is assumed that the input stream of water is described by a Poisson process, while the output stream is state dependent. LetLtdenote the dam level at timet, and letp1= limt→∞P{Lt=Llower} andp2= limt→∞P{Lt>Lupper} exist. The long-run average cost,J=p1J1+p2J2, is a performance measure. The aim of the paper is to choose the parameter controlling the output stream so as to minimizeJ.

This paper is concerned with the analysis of Markov decision processes in which a natural form of termination ensures that the expected future costs are bounded, at least under some policies. Whereas most previous analyses have restricted attention to the case where the set of states is finite, this paper analyses the case where the set of states is not necessarily finite or even countable. It is shown that all the existence, uniqueness, and convergence results of the finite-state case hold when the set of states is a general Borel space, provided we make the additional assumption that the optimal value function is bounded below. We give a sufficient condition for the optimal value function to be bounded below which holds, in particular, if the set of states is countable.

Let X = (X1, X2, …, Xn) be an exchangeable random vector, and write X(1:i) = min{X1, X2, …, Xi}, 1 ≤ i ≤ n. In this paper we obtain conditions under which X(1:i) decreases in i in the hazard rate order. A result involving more general (that is, not necessarily exchangeable) random vectors is also derived. These results are applied to obtain the limiting behaviour of the hazard rate function of the lifetimes of various coherent systems in reliability theory. The notions of the Samaniego signatures and the minimal signatures of such systems are extensively used in the paper. An interesting relationship between these two signatures is obtained. The results are illustrated in a series of examples.

We study the uncertain dichotomous choice model. In this model, a group of expert decision makers is required to select one of two alternatives. The applications of this model are relevant to a wide variety of areas. A decision rule translates the individual opinions of the members into a group decision, and is optimal if it maximizes the probability of the group making a correct choice. In this paper, we assume the correctness probabilities of the experts to be independent random variables selected from some given distribution. Moreover, the ranking of the members in the group is (at least partly) known. Thus, one can follow rules based on this ranking. The extremes are the expert rule and the majority rule. The probabilities of the two extreme rules being optimal were compared in a series of early papers, for a variety of distributions. In most cases, the asymptotic behaviours of the probabilities of the two extreme rules followed the same patterns. Do these patterns hold in general? If not, what are the ranges of possible asymptotic behaviours of the probabilities of the two extreme rules being optimal? In this paper, we provide satisfactory answers to these questions.

This thesis consists of four papers.In paper 1, we prove central limit theorems for Markov chains under (local) contraction conditions. As a corollary we obtain a central limit theorem for Markov chains associated with iterated function systems with contractive maps and place-dependent Dini-continuous probabilities.In paper 2, properties of inverse subordinators are investigated, in particular similarities with renewal processes. The main tool is a theorem on processes that are both renewal and Cox processes.In paper 3, distributional properties of supercritical and especially immortal branching processes are derived. The marginal distributions of immortal branching processes are found to be compound geometric.In paper 4, a description of a dynamic population model is presented, such that samples from the population have genealogies as given by a Lambda-coalescent with mutations. Depending on whether the sample is grouped according to litters or families, the sampling distribution is either regenerative or non-regenerative.

In this paper, we describe a class of Wiener functionals that are ‘indeterminate by their moments’, that is, whose distributions are not uniquely determined by their moments. In particular, it is proved that the integral of a geometric Brownian motion is indeterminate by its moments and, moreover, shown that previous proofs of this result are incorrect. The main result of this paper is based on geometric inequalities in Gauss space and on a generalization of the Krein criterion due to H. L. Pedersen.

This paper studies the subexponential properties of the stationary workload, actual waiting time and sojourn time distributions in work-conserving single-server queues when the equilibrium residual service time distribution is subexponential. This kind of problem has been previously investigated in various queueing and insurance risk settings. For example, it has been shown that, when the queue has a Markovian arrival stream (MAS) input governed by a finite-state Markov chain, it has such subexponential properties. However, though MASs can approximate any stationary marked point process, it is known that the corresponding subexponential results fail in the general stationary framework. In this paper, we consider the model with a general stationary input and show the subexponential properties under some additional assumptions. Our assumptions are so general that the MAS governed by a finite-state Markov chain inherently possesses them. The approach used here is the Palm-martingale calculus, that is, the connection between the notion of Palm probability and that of stochastic intensity. The proof is essentially an extension of the M/GI/1 case to cover ‘Poisson-like’ arrival processes such as Markovian ones, where the stochastic intensity is admitted.

This paper studies a discrete-time single-server queue with two independent inputs and service interruptions. One of the inputs to the queue is an independent and identically distributed process. The other is a much more general process and it is not required to be Markov nor is it required to be stationary. The service interruption process is also general and it is not required to be Markov or to be stationary. This paper shows that a stochastic decomposition property for the virtual waiting-time process holds in the discrete-time single-server queue with service interruptions. To the best of the author's knowledge, no stochastic decomposition results for virtual waiting-time processes in non-work-conserving queues, such as queues with service interruptions, have been obtained before and only work-conserving queues have been studied in the literature.

The full-information best-choice problem, as posed by Gilbert and Mosteller in 1966, asks us to find a stopping rule which maximizes the probability of selecting the largest of a sequence of n i.i.d. standard uniform random variables. Porosiński, in 1987, replaced a fixed n by a random N, uniform on {1,2,…,n} and independent of the observations. A partial-information problem, imbedded in a 1980 paper of Petruccelli, keeps n fixed but allows us to observe only the sequence of ranges (max - min), as well as whether or not the current observation is largest so far. Recently, Porosiński compared the solutions to his and Petruccelli's problems and found that the two problems have identical optimal rules as well as risks that are asymptotically equal. His discovery prompts the question: why? This paper gives a good explanation of the equivalence of the optimal rules. But even under the lens of a planar Poisson process model, it leaves the equivalence of the asymptotic risks as somewhat of a mystery. Meanwhile, two other problems have been shown to have the same limiting risks: the full-information problem with the (suboptimal) Porosiński-Petruccelli stopping rule, and the full-information ‘duration of holding the best’ problem of Ferguson, Hardwick and Tamaki, which turns out to be nothing but the Porosiński problem in disguise.

In this paper, we study the on-line parameter estimation problem for a partially observable system subject to deterioration and random failure. The state of the system evolves according to a continuous-time homogeneous Markov process with a finite state space. The state of the system is hidden except for the failure state. When the system is operating, only the information obtained by condition monitoring, which is related to the working state of the system, is available. The condition monitoring observations are assumed to be in continuous range, so that no discretization is required. A recursive maximum likelihood (RML) algorithm is proposed for the on-line parameter estimation of the model. The new RML algorithm proposed in the paper is superior to other RML algorithms in the literature in that no projection is needed and no calculation of the gradient on the surface of the constraint manifolds is required. A numerical example is provided to illustrate the algorithm.

We study the asymptotic properties of a minimal spanning tree formed by n points uniformly distributed in the unit square, where the minimality is amongst all rooted spanning trees with a direction of growth. We show that the number of branches from the root of this tree, the total length of these branches, and the length of the longest branch each converges weakly. This model is related to the study of record values in the theory of extreme-value statistics and this relation is used to obtain our results. The results also hold when the tree is formed from a Poisson point process of intensity n in the unit square. In this paper we introduce a notion of a minimal directed spanning tree. To illustrate this, consider the following model of transmission of radio waves by transmitters, receivers and amplifiers. Suppose a source transmitter located at the origin transmits messages which can be received only by receivers placed in the positive quadrant. These receivers in turn amplify the message and transmit it to the receivers lying in the positive quadrant with respect to these amplifiers. In this way each receiver lying in the positive quadrant with respect to the transmitter receives the message. The graph which represents this model of transmission may be viewed as a directed spanning tree. For transmission of radio waves in this model, the required strength of the source transmitter is clearly related to the number of receivers which receive the message directly from the transmitter at the origin, as well as the sum of the distances of these receivers from the origin and the distance of the receiver farthest away from the origin which receives the message directly from the origin. In this paper we investigate these three factors when the receivers are located at random in the unit square (say) and we study the asymptotics as the number of transmitters goes to infinity. This model of transmission of radio waves is in contrast to that introduced by Gilbert (1961) where the radio waves can travel in any direction from the point of their origin. However, they may be received only by receivers located within a certain fixed distance from the source of transmission.

This thesis, consisting of six papers and a summary, studies the area of continuous time financial mathematics. A unifying theme for many of the problems studied is the implications of possible mis-specifications of models. Intimately connected with this question is, perhaps surprisingly, convexity properties of option prices. We also study qualitative behavior of different optimal stopping boundaries appearing in option pricing.In Paper I a new condition on the contract function of an American option is provided under which the option price increases monotonically in the volatility. It is also shown that American option prices are continuous in the volatility.In Paper II an explicit pricing formula for the perpetual American put option in the Constant Elasticity of Variance model is derived. Moreover, different properties of this price are studied.Paper III deals with the Russian option with a finite time horizon. It is shown that the value of the Russian option solves a certain free boundary problem. This information is used to analyze the optimal stopping boundary.A study of perpetual game options is performed in Paper IV. One of the main results provides a condition under which the value of the option is increasing in the volatility.In Paper V options written on several underlying assets are considered. It is shown that, within a large class of models, the only model for the stock prices that assigns convex option prices to all convex contract functions is geometric Brownian motion.Finally, in Paper VI it is shown that the optimal stopping boundary for the American put option is convex in the standard Black-Scholes model.

Suppose that {X s , 0 ≤ s ≤ T} is an m-dimensional geometric Brownian motion with drift, μ is a bounded positive Borel measure on [0,T], and ϕ : ℝ m → [0,∞) is a (weighted) l q (ℝ m )-norm, 1 ≤ q ≤ ∞. The purpose of this paper is to study the distribution and the moments of the random variable Y given by the L p (μ)-norm, 1 ≤ p ≤ ∞, of the function s ↦ ϕ(X s ), 0 ≤ s ≤ T. By using various geometric inequalities in Wiener space, this paper gives upper and lower bounds for the distribution function of Y and proves that the distribution function is log-concave and absolutely continuous on every open subset of the distribution's support. Moreover, the paper derives tail probabilities, presents sharp moment inequalities, and shows that Y is indetermined by its moments. The paper will also discuss the so-called moment-matching method for the pricing of Asian-styled basket options.

In this paper, we study the classification of matrix GI/M/1-type Markov chains with a tree structure. We show that the Perron–Frobenius eigenvalue of a Jacobian matrix provides information for classifying these Markov chains. A fixed-point approach is utilized. A queueing application is presented to show the usefulness of the classification method developed in this paper.

This paper studies the law of any power of the integral of geometric Brownian motion over any finite time interval. As its main results, two integral representations for this law are derived. This is by enhancing the Laplace transform ansatz of Yor (1992) with complex analytic methods, which is the main methodological contribution of the paper. The one of our integrals has a similar structure to that obtained by Yor, while the other is in terms of Hermite functions as those of Dufresne (2001). Performing or not performing a certain Girsanov transformation is identified as the source of these two forms of the laws. While our results specialize for exponents equal to 1 to those obtained by Yor, they yield on specialization representations for the exponent equal to minus 1 laws which are markedly different from those obtained by Dufresne., 18 pages, 2 figures

In this paper, we consider denumerable-state continuous-time Markov decision processes with (possibly unbounded) transition and reward rates and general action space under the discounted criterion. We provide a set of conditions weaker than those previously known and then prove the existence of optimal stationary policies within the class of all possibly randomized Markov policies. Moreover, the results in this paper are illustrated by considering the birth-and-death processes with controlled immigration in which the conditions in this paper are satisfied, whereas the earlier conditions fail to hold.

This paper presents an algorithmic procedure to calculate the delay distribution of a type k customer in a first-come-first-served (FCFS) discrete-time queueing system with multiple types of customers, where each type has different service requirements (the MMAP[K]/PH[K]/1 queue). First, we develop a procedure, using matrix analytical methods, to handle arrival processes that do not allow batch arrivals to occur. Next, we show that this technique can be generalized to arrival processes that do allow batch arrivals to occur. We end the paper by presenting some numerical examples.

The purpose of this paper is to study ageing properties of first-passage times of increasing Markov chains. We extend the literature to some new ageing classes, such as the IFR(2), NBU(2), DRLLt and NBULt classes. We also give sufficient conditions in the finite case, that are more efficient computationally, just in terms of the transition matrix K, in the discrete case, or the generator matrix Q, in the continuous case. For the uniformizable, continuous-time Markov processes, we derive conditions in terms of the discrete uniformized Markov chain for the NBU(2) and the NBULt classes. In the last section, a review of the main results in this direction in the literature is given, and we compare some of the conditions stated in this paper with others given in the literature about some other ageing classes. Some examples where these results are applied are given.

In [20] Krajíček and Pudlák discovered connections between problems in computational complexity and the lengths of first-order proofs of finite consistency statements. Later Pudlák [25] studied more statements that connect provability with computational complexity and conjectured that they are true. All these conjectures are at least as strong as $\mathsf {P}\neq \mathsf {NP}$ [23–25].One of the problems concerning these conjectures is to find out how tightly they are connected with statements about computational complexity classes. Results of this kind had been proved in [20, 22].In this paper, we generalize and strengthen these results. Another question that we address concerns the dependence between these conjectures. We construct two oracles that enable us to answer questions about relativized separations asked in [19, 25] (i.e., for the pairs of conjectures mentioned in the questions, we construct oracles such that one conjecture from the pair is true in the relativized world and the other is false and vice versa). We also show several new connections between the studied conjectures. In particular, we show that the relation between the finite reflection principle and proof systems for existentially quantified Boolean formulas is similar to the one for finite consistency statements and proof systems for non-quantified propositional tautologies.

Siegmund and Worsley (1995) considered the problem of testing for signals with unknown location and scale in a Gaussian random field defined on ℝN. The test statistic was the maximum of a Gaussian random field in anN+1 dimensional ‘scale space’,Ndimensions for location and 1 dimension for the scale of a smoothing filter. Scale space is identical to a continuous wavelet transform with a kernel smoother as the wavelet, though the emphasis here is on signal detection rather than image compression or enhancement. Two methods were used to derive an approximate null distribution forN=2 andN=3: one based on the method of volumes of tubes, the other based on the expected Euler characteristic of the excursion set. The purpose of this paper is two-fold: to show how the latter method can be extended to higher dimensions, and to apply this more general result to χ2fields. The result of Siegmund and Worsley (1995) then follows as a special case. In this paper the results are applied to the problem of searching for activation in brain images obtained by functional magnetic resonance imaging (fMRI).

Siegmund and Worsley (1995) considered the problem of testing for signals with unknown location and scale in a Gaussian random field defined on ℝN. The test statistic was the maximum of a Gaussian random field in anN+1 dimensional ‘scale space’,Ndimensions for location and 1 dimension for the scale of a smoothing filter. Scale space is identical to a continuous wavelet transform with a kernel smoother as the wavelet, though the emphasis here is on signal detection rather than image compression or enhancement. Two methods were used to derive an approximate null distribution forN=2 andN=3: one based on the method of volumes of tubes, the other based on the expected Euler characteristic of the excursion set. The purpose of this paper is two-fold: to show how the latter method can be extended to higher dimensions, and to apply this more general result to χ2fields. The result of Siegmund and Worsley (1995) then follows as a special case. In this paper the results are applied to the problem of searching for activation in brain images obtained by functional magnetic resonance imaging (fMRI).

This paper raises the following question: let {Φn(A), A ⊂ ℝd} be a Poisson process with intensity nf(x), x ∈ ℝd and let c(Xi | Φn) be a Voronoi tile with nucleus Xi (a jump point of Φn). Let μ(.) denote Lebesgue measure in ℝd. Is it true that, for any bounded measurable subset B of ℝd, ∑Xi∈Bμ(c(Xi| Φn)) → μ(B) almost surely as n → ∞ only if f > 0 almost everywhere? This statement can be viewed as the strong law of large numbers for Voronoi tessellation. Though the positive answer may seem ‘obvious’, we could not find any such statement, especially for arbitrary measurable B and nonhomogeneous Poisson processes. For B with the boundary of Lebesgue measure 0 the proof is simple. We prove in this paper that the statement is true for ℝ1.

Consider a renewal process. The renewal events partition the process into i.i.d. renewal cycles. Assume that on each cycle, a rare event called 'success’ can occur. Such successes lend themselves naturally to approximation by Poisson point processes. If each success occurs after a random delay, however, Poisson convergence may be relatively slow, because each success corresponds to a time interval, not a point. In 1996, Altschul and Gish proposed a finite-size correction to a particular approximation by a Poisson point process. Their correction is now used routinely (about once a second) when computers compare biological sequences, although it lacks a mathematical foundation. This paper generalizes their correction. For a single renewal process or several renewal processes operating in parallel, this paper gives an asymptotic expansion that contains in successive terms a Poisson point approximation, a generalization of the Altschul-Gish correction, and a correction term beyond that.

This paper continues the work of Saikkonen (2001, Econometric Theory 17, 296–326) and develops an asymptotic theory of statistical inference in cointegrated vector autoregressive models with nonlinear time trends in cointegrating relations and general nonlinear parameter restrictions. Inference on parameters in cointegrating relations and short-run dynamics is studied separately. It is shown that Gaussian maximum likelihood estimators of parameters in cointegrating relations have mixed normal limiting distributions and that related Wald, Lagrange multiplier, and likelihood ratio tests for general nonlinear hypotheses have usual asymptotic chi-square distributions. These results are shown to hold even if parameters in the short-run dynamics are not identified. In that case suitable estimators of the information matrix have to be used to justify the application of Wald and Lagrange multiplier tests, whereas the likelihood ratio test is free of this difficulty. Similar results are also obtained when inference on parameters in the short-run dynamics is studied, although then Gaussian maximum likelihood estimators have usual normal limiting distributions. All results of the paper are proved without assuming existence of second partial derivatives of the likelihood function, and in some cases even differentiability with respect to nuisance parameters is not required.

In this paper, results on spectrally negative Lévy processes are used to study the ruin probability under some risk processes. These processes include the compound Poisson process and the gamma process, both perturbed by diffusion. In addition, the first time the risk process hits a given level is also studied. In the case of classical risk process, the joint distribution of the ruin time and the first recovery time is obtained. Some results in this paper have appeared before (e.g., Dufresne and Gerber (1991), Gerber (1990), dos Reis (1993)). We revisit them from the Lévy process theory's point of view and in a unified and simple way.

The paper considers one of the standard processes for modeling returns in finance, the stochastic volatility process with regularly varying innovations. The aim of the paper is to show how point process techniques can be used to derive the asymptotic behavior of the sample autocorrelation function of this process with heavy-tailed marginal distributions. Unlike other non-linear models used in finance, such as GARCH and bilinear models, sample autocorrelations of a stochastic volatility process have attractive asymptotic properties. Specifically, in the infinite variance case, the sample autocorrelation function converges to zero in probability at a rate that is faster the heavier the tails of the marginal distribution. This behavior is analogous to the asymptotic behavior of the sample autocorrelations of independent identically distributed random variables.

This paper studies the dependence structure and bounds of several basic prototypical parallel queueing systems with correlated arrival processes to different queues. The marked feature of our systems is that each queue viewed alone is a standard single-server queuing system extensively studied in the literature, but those queues are statistically dependent due to correlated arrival streams. The major difficulty in analysing those systems is that the presence of correlation makes the explicit computation of a joint performance measure either intractable or computationally intensive. In addition, it is not well understood how and in what sense arrival correlation will improve or deteriorate a system performance measure. The objective of this paper is to provide a better understanding of the dependence structure of correlated queueing systems and to derive computable bounds for the statistics of a joint performance measure. In this paper, we obtain conditions on arrival processes under which a performance measure in two systems can be compared, in the sense of orthant and supermodular orders, among different queues and over different arrival times. Such strong comparison results enable us to study both spatial dependence (dependence among different queues) and temporal dependence (dependence over different time instances) for a joint performance measure. Further, we derive a variety of upper and lower bounds for the statistics of a stationary joint performance measure. Finally, we apply our results to synchronized queueing systems, using the ideas combined from the theory of orthant and supermodular dependence orders and majorization with respect to weighted trees (Xu and Li (2000)). Our results reveal how a performance measure can be affected, favourably or adversely, by different types of dependencies.

This paper arose from interest in assessing the quality of random number generators. The problem of testing randomness of a string of binary bits produced by such a generator gained importance with the wide use of public key cryptography and the need for secure encryption algorithms. All such algorithms are based on a generator of (pseudo) random numbers; the testing of such generators for randomness became crucial for the communications industry where digital signatures and key management are vital for information processing. The concept of approximate entropy has been introduced in a series of papers by S. Pincus and co-authors. The corresponding statistic is designed to measure the degree of randomness of observed sequences. It is based on incremental contrasts of empirical entropies based on the frequencies of different patterns in the sequence. Sequences with large approximate entropy must have substantial fluctuation or irregularity. Alternatively, small values of this characteristic imply strong regularity, or lack of randomness, in a sequence. Pincus and Kalman (1997) evaluated approximate entropies for binary and decimal expansions of e, π, √2 and √3 with the surprising conclusion that the expansion of √3 demonstrated much less irregularity than that of π. Tractable small sample distributions are hardly available, and testing randomness is based, as a rule, on fairly long strings. Therefore, to have rigorous statistical tests of randomness based on this approximate entropy statistic, one needs the limiting distribution of this characteristic under the randomness assumption. Until now this distribution remained unknown and was thought to be difficult to obtain. To derive the limiting distribution of approximate entropy we modify its definition. It is shown that the approximate entropy as well as its modified version converges in distribution to a χ2-random variable. The P-values of approximate entropy test statistics for binary expansions of e, π and √3 are plotted. Although some of these values for √3 digits are small, they do not provide enough statistical significance against the randomness hypothesis.

This paper focuses on the analyticity of the limiting behavior of a class of dynamical systems defined by iteration of non-expansive random operators. The analyticity is understood with respect to the parameters which govern the law of the operators. The proofs are based on contraction with respect to certain projective semi-norms. Several examples are considered, including Lyapunov exponents associated with products of random matrices both in the conventional algebra, and in the (max, +) semi-field, and Lyapunov exponents associated with non-linear dynamical systems arising in stochastic control. For the class of reducible operators (defined in the paper), we also address the issue of analyticity of the expectation of functionals of the limiting behavior, and connect this with contraction properties with respect to the supremum norm. We give several applications to queueing theory.

In this paper, we give necessary and sufficient conditions to ensure the validity of confidence intervals, based on the central limit theorem, in simulations of highly reliable Markovian systems. We resort to simulations because of the frequently huge state space in practical systems. So far the literature has focused on the property of bounded relative error. In this paper we focus on ‘bounded normal approximation’ which asserts that the approximation of the normal law, suggested by the central limit theorem, does not deteriorate as the reliability of the system increases. Here we see that the set of systems with bounded normal approximation is (strictly) included in the set of systems with bounded relative error.

On the basis of the Karlin and McGregor result, which states that the transition probability functions of a birth and death process can be expressed via the introduction of an orthogonal polynomial system and a spectral measure, we investigate in this paper how the Laplace transforms and the distributions of different transient characteristics related to excursions of a birth and death process can be expressed by means of the basic orthogonal polynomial system and the spectral measure. This allows us in particular to give a probabilistic interpretation of the series introduced by Stieltjes to study the convergence of the fundamental continued fraction associated with the system. Throughout the paper, we pay special attention to the case when the birth and death process is ergodic. Under the assumption that the spectrum of the spectral measure is discrete, we show how the distributions of different random variables associated with excursions depend on the fundamental continued fraction, the orthogonal polynomial system and the spectral measure.

This paper is concerned with the standard bivariate death process as well as with some Markovian modifications and extensions of the process that are of interest especially in epidemic modeling. A new and powerful approach is developed that allows us to obtain the exact distribution of the population state at any point in time, and to highlight the actual nature of the solution. Firstly, using a martingale technique, a central system of relations with two indices for the temporal state distribution will be derived. A remarkable property is that for all the models under consideration, these relations exhibit a similar algebraic structure. Then, this structure will be exploited by having recourse to a theory of Abel-Gontcharoff pseudopolynomials with two indices. This theory generalizes the univariate case examined in a preceding paper and is briefly introduced in the Appendix.

The paper is concerned with the final state and severity of a number of SIR epidemic models in finite populations. Two different classes of models are considered, namely the classical SIR Markovian models and the collective models introduced recently by the authors. First, by applying a simple martingale argument, it is shown that in both cases, there exists a common algebraic structure underlying the exact law of the final state and severity. Then, a unified approach to these statistics is developed by exploiting the theory of Abel-Gontcharoff pseudopolynomials (presented in a preceding paper).

Recently in the paper by Møller and Zuyev (1996), the following Gamma-type result was established. Given n points of a homogeneous Poisson process defining a random figure, its volume is Γ(n,λ) distributed, where λ is the intensity of the process. In this paper we give an alternative description of the class of random sets for which the Gamma-type results hold. We show that it corresponds to the class of stopping sets with respect to the natural filtration of the point process with certain scaling properties. The proof uses the martingale technique for directed processes, in particular, an analogue of Doob's optional sampling theorem proved in Kurtz (1980). As well as being compact, this approach provides a new insight into the nature of geometrical objects constructed with respect to a Poisson point process. We show, in particular, that in this framework the probability that a point is covered by a stopping set does not depend on whether it is a point of the process or not.

In this paper, in work strongly related with that of Coffman et al. [5], Bruss and Robertson [2], and Rhee and Talagrand [15], we focus our interest on an asymptotic distributional comparison between numbers of ‘smallest’ i.i.d. random variables selected by either on-line or off-line policies. Let X 1,X 2,… be a sequence of i.i.d. random variables with distribution function F(x), and let X 1,n ,…,X n,n be the sequence of order statistics of X 1,…,X n . For a sequence (c n ) n≥1 of positive constants, the smallest fit off-line counting random variable is defined by N e (c n ) := max {j ≤ n : X 1,n + … + X j,n ≤ c n }. The asymptotic joint distributional comparison is given between the off-line count N e (c n ) and on-line counts N n τ for ‘good’ sequential (on-line) policies τ satisfying the sum constraint ∑ j≥1 X τ j I (τ j ≤n) ≤ c n . Specifically, for such policies τ, under appropriate conditions on the distribution function F(x) and the constants (c n ) n≥1, we find sequences of positive constants (B n ) n≥1, (Δ n ) n≥1 and (Δ' n ) n≥1 such that for some non-degenerate random variables W and W'. The major tools used in the paper are convergence of point processes to Poisson random measure and continuous mapping theorems, strong approximation results of the normalized empirical process by Brownian bridges, and some renewal theory.

Much of the literature in reliability and survival analysis considers failure models indexed by a single scale. There are situations which require that failure be described by several scales. An example from reliability is items under warranty whose failure is recorded by time and amount of use. An example from survival analysis is the death of a mine worker which is noted by age and the duration of exposure to dust. This paper proposes an approach for developing probabilistic models indexed by two scales: time, and usage, a quantity that is related to time. The relationship between the scales is described by an additive hazards model. The evolution of usage is described by stochastic processes like the Poisson, the gamma and the Markov additive. The paper concludes with an application involving the setting of warranties. Two features differentiate this work from related efforts: a use of specific processes for describing usage, and a use of Monte Carlo techniques for generating the models.

Target benefit (TB) plans that incorporate intergenerational risk sharing have been demonstrated to be welfare improving over the long term. However, there has been little discussion of the short-term benefits for members in a defined benefit (DB) plan that is transitioning to TB. In this paper, we adopt a two-step approach that is designed to ensure the long-term sustainability of the new plan, without unduly sacrificing the benefit security of current retirees. We propose a cohort-based transition plan for reducing intergenerational inequity. Our study is based on simulations using an economic scenario generator with some theoretical results under simplified settings.

Analyzing variation in treatment effects across subsets of the population is an important way for social scientists to evaluate theoretical arguments. A common strategy in assessing such treatment effect heterogeneity is to include a multiplicative interaction term between the treatment and a hypothesized effect modifier in a regression model. Unfortunately, this approach can result in biased inferences due to unmodeled interactions between the effect modifier and other covariates, and including these interactions can lead to unstable estimates due to overfitting. In this paper, we explore the usefulness of machine learning algorithms for stabilizing these estimates and show how many off-the-shelf adaptive methods lead to two forms of bias: direct and indirect regularization bias. To overcome these issues, we use a post-double selection approach that utilizes several lasso estimators to select the interactions to include in the final model. We extend this approach to estimate uncertainty for both interaction and marginal effects. Simulation evidence shows that this approach has better performance than competing methods, even when the number of covariates is large. We show in two empirical examples that the choice of method leads to dramatically different conclusions about effect heterogeneity.

Most of the real-life populations are heterogeneous and homogeneity is often just a simplifying assumption for the relevant statistical analysis. Mixtures of lifetime distributions that correspond to homogeneous subpopulations were intensively studied in the literature. Various distributional and stochastic properties of finite and continuous mixtures were discussed. In this paper, following recent publications, we develop further a mixture concept in the form of the generalized α-mixtures that include all mixture models that are widely explored in the literature. We study some main stochastic properties of the suggested mixture model, that is, aging and appropriate stochastic comparisons. Some relevant examples and counterexamples are given to illustrate our findings.

This paper establishes structural properties for the throughput of a large class of queueing networks with i.i.d. new-better-than-used service times. The main result obtained in this paper is applied to a wide range of networks, including tandems, cycles and fork-join networks with general blocking and starvation (as well as certain networks with splitting and merging of traffic streams), to deduce the concavity of their throughput as a function of system parameters, such as buffer and initial job configurations, and blocking and starvation parameters. These results have important implications for the optimal design and control of such queueing networks by providing exact solutions, reducing the search space over which optimization need be performed, or establishing the convergence of optimization algorithms. In order to obtain results for such disparate networks in a unified manner, we introduce the framework of constrained discrete event systems (CDES), which enables us to characterize any permutable and non-interruptive queueing network through its constraint set. The main result of this paper establishes comparison properties of the event occurrence processes of CDES as a function of the constraint sets, which are then translated into the above-mentioned concavity of the throughput as a function of system parameters in the context of queueing networks.

This paper deals with existence of bivariate Fréchet optimal lower bounds for two sets of events, and provides a practical approach to find this kind of bound. The main device used is linear programming ideas, coupled with construction of probability spaces. The highlight of this paper is that perturbation terms in the optimization process, even when a tie occurs, are not necessary in this practical implementation.

This paper consists of two parts. The first part provides a more elementary proof of the asymptotic theorem of the refusals stream for an M/GI/1/n queueing system discussed in Abramov (1991a). The central property of the refusals stream discussed in the second part of this paper is that, if the expectations of interarrival and service time of an M/GI/1/n queueing system are equal to each other, then the expectation of the number of refusals during a busy period is equal to 1. This property is extended for a wide family of single-server queueing systems with refusals including, for example, queueing systems with bounded waiting time.

In this paper we analyze a single server two-queue model with Bernoulli schedules. This discipline is very flexible and contains the exhaustive and 1-limited disciplines as special cases. We formulate the queueing system as a Riemann boundary value problem with shift. The boundary value problem is solved by exploring a Fredholm integral equation around the unit circle. Some numerical examples are presented at the end of the paper.

A number of recent papers have exhibited classes of queueing networks, with batches of customers served and routed through the network, which have generalised product-form equilibrium distributions. In this paper we look at these from a new viewpoint. In particular we show that, under standard assumptions, for a network to possess an equilibrium distribution that factorises into a product form over the nodes of the network for all possible transition rates, it is necessary and sufficient that it be equivalent to a suitably-defined single-movement network. We consider also the form of the state space for such networks.

This paper is concerned with the study of death processes with time-homogeneous non-linear death rates. An explicit formula is obtained for the joint distribution of the state XT and the variable ∫ T 0 g(Xt ), where g is any given real function and T corresponds to some appropriate stopping time. This is achieved by constructing a family of martingales and then by using a particular family of Abel–Gontcharoff pseudopolynomials (the theory of which has been introduced in a companion paper) and related Abelian-type expansions. Moreover, the distribution of the first crossing level of such a death process through a general upper boundary is also evaluated in terms of pseudopolynomials of that kind. The flexibility of the methods developed makes easy the extension to multidimensional death processes.

The paper is concerned with the distribution of the level N of the first crossing of a counting process trajectory with a lower boundary. Compound and simple Poisson or binomial processes, gamma renewal processes, and finally birth processes are considered. In the simple Poisson case, expressing the exact distribution of N requires the use of a classical family of Abel–Gontcharoff polynomials. For other cases convenient extensions of these polynomials into pseudopolynomials with a similar structure are necessary. Such extensions being applicable to other fields of applied probability, the central part of the present paper has been devoted to the building of these pseudopolynomials in a rather general framework.