We tackle the problem of simulating seat- and vote-shares for a party system of a given size. We show how these shares can be generated using unordered and ordered Dirichlet distributions. We show that a distribution with a mean vector given by the rule described in Taagepera and Allik (2006, Electoral Studies 25, 696–713) fits real-world data almost as well as a saturated model where there is a parameter for each rank/system size combination.

Given a finite strongly connected directed graph $G=(V, E)$ , we study a Markov chain taking values on the space of probability measures on V. The chain, motivated by biological applications in the context of stochastic population dynamics, is characterized by transitions between states that respect the structure superimposed by E: mass (probability) can only be moved between neighbors in G. We provide conditions for the ergodicity of the chain. In a simple, symmetric case, we fully characterize the invariant probability.

Consider a Pólya urn with balls of several colours, where balls are drawn sequentially and each drawn ball is immediately replaced together with a fixed number of balls of the same colour. It is well known that the proportions of balls of the different colours converge in distribution to a Dirichlet distribution. We show that the rate of convergence is $\Theta(1/n)$ in the minimal $L_p$ metric for any $p\in[1,\infty]$, extending a result by Goldstein and Reinert; we further show the same rate for the Lévy distance, while the rate for the Kolmogorov distance depends on the parameters, i.e. on the initial composition of the urn. The method used here differs from the one used by Goldstein and Reinert, and uses direct calculations based on the known exact distributions.

We study the limit behaviour of a class of random walk models taking values in the standard d-dimensional ( $d\ge 1$ ) simplex. From an interior point z, the process chooses one of the $d+1$ vertices of the simplex, with probabilities depending on z, and then the particle randomly jumps to a new location z′ on the segment connecting z to the chosen vertex. In some special cases, using properties of the Beta distribution, we prove that the limiting distributions of the Markov chain are Dirichlet. We also consider a related history-dependent random walk model in [0, 1] based on an urn-type scheme. We show that this random walk converges in distribution to an arcsine random variable.

Consider a haploid population of fixed finite size with a finite number of allele types and having Cannings exchangeable genealogy with neutral mutation. The stationary distribution of the Markov chain of allele counts in each generation is an important quantity in population genetics but has no tractable description in general. We provide upper bounds on the distributional distance between the Dirichlet distribution and this finite-population stationary distribution for the Wright–Fisher genealogy with general mutation structure and the Cannings exchangeable genealogy with parent independent mutation structure. In the first case, the bound is small if the population is large and the mutations do not depend too much on parent type; 'too much' is naturally quantified by our bound. In the second case, the bound is small if the population is large and the chance of three-mergers in the Cannings genealogy is small relative to the chance of two-mergers; this is the same condition to ensure convergence of the genealogy to Kingman's coalescent. These results follow from a new development of Stein's method for the Dirichlet distribution based on Barbour's generator approach and a probabilistic description of the semigroup of the Wright–Fisher diffusion due to Griffiths and Li (1983) and Tavaré (1984).

We characterise the class of distributions of random stochastic matrices X with the property that the products X(n)X(n − 1) · · · X(1) of independent and identically distributed copies X(k) of X converge almost surely as n → ∞ and the limit is Dirichlet distributed. This extends a result by Chamayou and Letac (1994) and is illustrated by several examples that are of interest in applications.

In a real n-1 dimensional affine space E, consider a tetrahedron T0, i.e. the convex hull of n points α1, α2, …, αn of E. Choose n independent points β1, β2, …, βn randomly and uniformly in T0, thus obtaining a new tetrahedron T1 contained in T0. Repeat the operation with T1 instead of T0, obtaining T2, and so on. The sequence of the Tk shrinks to a point Y of T0 and this note computes the distribution of the barycentric coordinates of Y with respect to (α1, α2, …, αn) (Corollary 2.3). We also obtain the explicit distribution of Y in more general cases. The technique used is to reduce the problem to the study of a random walk on the semigroup of stochastic (n,n) matrices, and this note is a geometrical application of a former result of Chamayou and Letac (1994).

A Bayesian approach for analyzing layered defense systems is presented. This approach incorporates the dependence of penetration probabilities on the size of attackers going into any layer. A general formula is developed for computing the predictive distribution of the number of attackers surviving any layer as well as the posterior distribution of the penetration probabilities under the a priori assumptions that: (i) the probabilities are dependent and their joint distribution is Dirichlet, and (ii) the probabilities are independent. Positive dependence of the penetration probabilities as well as the number of attackers surviving the different layers is also established.

We consider a general credibility model for the prediction of IBNR-claims which allows for random fluctuations in the underlying delay distribution. Such fluctuations always bring about decreasing credibility. It is shown that even negative credibility is achieved for more substantial fluctuations in the delay distribution. Special attention is paid to the mixed Poisson case for claim numbers including the discussion of parameter estimation.

An expansion in orthogonal polynomials is found for the transition density in a neutral multi-allele diffusion model where the mutation rates of allele types Ai → Aj are assumed to be uj(≥ O). The density is found when the mutation rate is positive for all allele types, and when some or all have zero mutation. The asymptotic conditional density is found for a mixture of positive and zero mutation rates.

The infinite alleles limit with equal mutation is studied. Eigenfunctions of the process are derived and the frequency spectrum found. An important result is that the first eigenfunction depends only on the homozygosity.

A density for the time to fixation with zero mutation is found for the K allele, and infinite alleles model.

An expression is found for the stationary density of the allele frequencies, in the infinitely-many alleles model. It is assumed that all alleles are neutral, that there is a constant mutation rate, and that the population is sufficiently large for the diffusion model to apply. Bounds on some moments are calculated from the density, and some applications are made to the problem of which allele is oldest in a population. A postulate made by Ewens (1972), concerning the distribution of allele numbers in a finite random sample from the neutral diffusion population, is shown to be correct.