We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Computationally, when we solve for the stationary probabilities for a countable-state Markov chain, the transition probability matrix of the Markov chain has to be truncated, in some way, into a finite matrix. Different augmentation methods might be valid such that the stationary probability distribution for the truncated Markov chain approaches that for the countable Markov chain as the truncation size gets large. In this paper, we prove that the censored (watched) Markov chain provides the best approximation in the sense that, for a given truncation size, the sum of errors is the minimum and show, by examples, that the method of augmenting the last column only is not always the best.
In a famous paper, Dwass (1967) proposed a method to deal with rank order statistics, which constitutes a unifying framework to derive various distributional results. In the present paper an alternative method is presented, which allows us to extend Dwass's results in several ways, namely arbitrary endpoints, horizontal steps and arbitrary probabilities for the three step types. Regarding these extensions the pertaining rank order statistics are extended as well to simple random walk statistics. This method has proved appropriate to generalize all results given by Dwass. Moreover, these discrete time results can be taken as a starting point to derive the corresponding results for randomized random walks by means of a limiting process.
In this paper various statistics for randomized random walks and their distributions are presented. The distributional results are derived by means of a limiting procedure applied to the pertaining discrete time process, which has been considered in part I of this work (Katzenbeisser and Panny 1996). This basic approach, originally due to Meisling (1958), seems to offer certain technical advantages, since it avoids the use of Laplace transforms and is even simpler than Feller's randomization technique.
Recursive trees have been used to model such things as the spread of epidemics, family trees of ancient manuscripts, and pyramid schemes. A tree Tn with n labeled nodes is a recursive tree if n = 1, or n > 1 and Tn can be constructed by joining node n to a node of some recursive tree Tn–1. For arbitrary nodes i < n in a random recursive tree we give the exact distribution of Xi,n, the distance between nodes i and n. We characterize this distribution as the convolution of the law of Xi,j+1 and n – i – 1 Bernoulli distributions. We further characterize the law of Xi,j+1 as a mixture of sums of Bernoullis. For i = in growing as a function of n, we show that is asymptotically normal in several settings.
We consider the application of importance sampling in steady-state simulations of finite Markov chains. We show that, for a large class of performance measures, there is a choice of the alternative transition matrix for which the ratio of the variance of the importance sampling estimator to the variance of the naive simulation estimator converges to zero as the sample path length goes to infinity. Obtaining this ‘optimal’ transition matrix involves computing the performance measure of interest, so the optimal matrix cannot be computed in precisely those situations where simulation is required to estimate steady-state performance. However, our results show that alternative transition matrices of the form Q = P + E/T, where P is the original transition matrix and T is the sample path length, can be expected to provide good results. Moreover, we provide an iterative algorithm for obtaining alternative transition matrices of this form that converge to the optimal matrix as the number of iterations increases, and present an example that shows that spending some computer time iterating this algorithm and then conducting the simulation with the resulting alternative transition matrix may provide considerable variance reduction when compared to naive simulation.
A Markov chain is used as a model for a sequence of random experiments. The waiting time for sequence patterns is considered. Recursive-type relations for the distribution of waiting times are obtained.
Consider a stationary Markov chain with state space consisting of the ξ -letter alphabet set Λ= {a1, a2, ···, aξ }. We study the variables M=M(n, k) and N=N(n, k), defined, respectively, as the number of overlapping and non-overlapping occurrences of a fixed periodic k-letter word, and use the Stein–Chen method to obtain compound Poisson approximations for their distribution.
The cumulative distribution of the finite sum of the binary sequence of order k is studied and some of its applications discussed. Certain properties of this sequence are investigated and uniformly superior bounds for the cumulative distribution under minimal information on the ‘success' probabilities are derived.
The paper deals with processes in (or ), one of whose components is skip-free. We obtain identities for distributions of hitting times for the components of the process generalizing the well-known one for the one-dimensional case. These relations reflect the fact that in this case spatial and time coordinates play, in some sense, symmetric roles. They turn out to be useful for solving several problems. For example, they allow us to find the distribution of the number of jumps of the process, which fall in a fixed set before the skip-free component of the process hits a fixed level. Examples are given showing how our results can be applied to models in branching processes, queueing, and risk theory.
Guided by analogy with Euler's spherical excess formula, we define a finite-additive functional on bounded convex polygons in ℝ2 (the Euler functional). Under certain smoothness assumptions, we find some sufficient conditions when this functional can be extended to a planar signed measure. A dual reformulation of these conditions leads to signed measures in the space of lines in ℝ2. In this way we obtain two sets of conditions which ensure that a segment function corresponds to a signed measure in the space of lines. The latter conditions are also necessary.
The random k-dimensional partial order Pk(n) on n points is defined by taking n points uniformly at random from [0,1]k. Previous work has concentrated on the case where k is constant: we consider the model where k increases with n.
We pay particular attention to the height Hk(n) of Pk(n). We show that k = (t/log t!) log n is a sharp threshold function for the existence of a t-chain in Pk(n): if k – (t/log t!) log n tends to + ∞ then the probability that Pk(n) contains a t-chain tends to 0; whereas if the quantity tends to − ∞ then the probability tends to 1. We describe the behaviour of Hk(n) for the entire range of k(n).
We also consider the maximum degree of Pk(n). We show that, for each fixed d ≧ 2, is a threshold function for the appearance of an element of degree d. Thus the maximum degree undergoes very rapid growth near this value of k.
We make some remarks on the existence of threshold functions in general, and give some bounds on the dimension of Pk(n) for large k(n).
At which (random) sample size will every population element have been drawn at least m times? This special coupon collector's problem is often referred to as the Dixie cup problem. Some asymptotic properties of the Dixie cup problem with unequal sampling probabilities are described.
We discuss the probability that among a number of random discrete values there is a unique maximal value. Some general inequalities are derived. As an application, we study the limit behavior of the probability of uniqueness for a subexponential number of binomially distributed variables.
Consider a number of events in a probability space. Let X be a random variable that is the number of events that occur. Given some of the moments of the distribution of X, it is possible to obtain bounds on the probability that at least one event occurs. The best possible bounds are given here. If there are many equiprobable events that are d- wise independent, and d is even, then the probability that at least one event happens is at least 1 — O(µ–d/2), where μ = E(X).
Some exact and asymptotic joint distributions are given for certain random variables defined on the excursions of a simple symmetric random walk. We derive appropriate recursion formulas and apply them to get certain expressions for the joint generating or characteristic functions of the random variables.
Let A1, A2, …, An and B1, B2,. ., BN be two sequences of events on the same probability space. Let mn(A) and mN(B), respectively, be the number of those Aj and Bj which occur. Let Si,j denote the joint ith binomial moment of mn(A) and jth binomial moment of mN(B), 0 ≤ i ≤ n, 0 ≤ j ≤ N. For fixed non-negative integers a and b, we establish both lower and upper bounds on the distribution P(mn(A) = r, mN(B) = u) by linear combinations of Si,j, 0 ≤ i ≤ a, 0 ≤ j ≤ b. When both a and b are even, all mentioned S¡,j are utilized in both the upper and the lower bound. In a set of remarks the results are analyzed and their relation to the existing literature, including the univariate case, is discussed.
The accuracy of the Poisson approximation to the distribution of the numbers of large and small m-spacings, when n points are placed at random on the circle, was analysed using the Stein–Chen method in Barbour et al. (1992b). The Poisson approximation for m≧2 was found not to be as good as for 1-spacings. In this paper, rates of approximation of these distributions to suitable compound Poisson distributions are worked out, using the CP–Stein–Chen method and an appropriate coupling argument. The rates are better than for Poisson approximation for m≧2, and are of order O((log n)2/n) for large m-spacings and of order O(1/n) for small m-spacings, for any fixed m≧2, if the expected number of spacings is held constant as n → ∞.
Using the electric network approach, we give very simple derivations for the expected first passage from the origin to the opposite vertex in the d-cube (i.e. the Ehrenfest urn model) and the Platonic graphs.
Let X1, X2, · ··, Xn be a sequence of n random variables taking values in the ξ -letter alphabet . We consider the number N = N(n, k) of non-overlapping occurrences of a fixed k-letter word under (a) i.i.d. and (b) stationary Markovian hypotheses on the sequence , and use the Stein–Chen method to obtain Poisson approximations for the same. In each case, results and couplings from Barbour et al. (1992) are used to show that the total variation distance between the distribution of N and that of an appropriate Poisson random variable is of order (roughly) O(kS(k)), where S(k) denotes the stationary probability of the word in question. These results vastly improve on the approximations obtained in Godbole (1991). In the Markov case, we also make use of recently obtained eigenvalue bounds on convergence to stationarity due to Diaconis and Stroock (1991) and Fill (1991).