Two Universality Properties Associated with the Monkey Model of Zipf's Law

The distribution of word probabilities in the monkey model of Zipf's law is associated with two universality properties: (1) the power law exponent converges strongly to $-1$ as the alphabet size increases and the letter probabilities are specified as the spacings from a random division of the unit interval for any distribution with a bounded density function on $[0,1]$; and (2), on a logarithmic scale the version of the model with a finite word length cutoff and unequal letter probabilities is approximately normally distributed in the part of the distribution away from the tails. The first property is proved using a remarkably general limit theorem for the logarithm of sample spacings from Shao and Hahn, and the second property follows from Anscombe's central limit theorem for a random number of i.i.d. random variables. The finite word length model leads to a hybrid Zipf-lognormal mixture distribution closely related to work in other areas.


Introduction
In his popular expository article on universality examples in mathematics and physics Tao [1] mentions Zipf's empirical power law of word frequencies as lacking any "convincing explanation for how the law comes about and why it is universal." By universality he means the idea where systems follow certain macroscopic laws that are largely independent of their microscopic details. In this article we drop down many levels from natural language to discuss two universality properties associated with the toy monkey model of Zipf's law. Both of these properties were considered in Perline [2], but the presentation there was incompletely developed, and here we clarify and greatly expand upon these two ideas by showing: (1) how the model displays a nearly universal tendency towards a −1 exponent in its power law behavior; and (2) the significance of the central limit theorem (CLT) for a random number of i.i.d. variables in the model with a finite word length cutoff. We will sometimes refer to the case with a finite word length cutoff as Monkey Twitter, as explained in Section 3.
Our paper is laid out as follows. In Section 2 we prove the strong tendency towards an approximate −1 exponent under very broad conditions by means of a limit theorem for the logarithms of random spacing due to Shao and Hahn [3]. A somewhat longer approach to this proof is in Perline [4], where an elementary derivation of power law behavior of the monkey model is first given before applying the Shao-Hahn limit result. In Section 3, we explain the underlying lognormal structure of the central part of the distribution of word probabilities in the case of a finite word length cutoff and show how it leads to a hybrid Zipf-lognormal distribution (called a lognormal-Pareto distribution in [2]). These are universality properties in exactly Tao's sense because the distribution of the word probabilities (the macro behavior) is within very broad bounds independent of the details of how the letter probabilities for the monkey's typewriter are selected (the micro behavior). In Section 4, we discuss the multiple connections between these results and well-known research from other areas where Pareto-Zipf type distributions have been investigated. In the remainder of this section we sketch some historical background.
It has been known for many years that the monkey-at-the-typewriter scheme for generating random "words" conforms to an inverse power law of word frequencies and therefore mimicks the inverse power form of Zipf's [5] statistical regularity for many natural languages. However, Zipf 's empirical word frequency law not only follows a power distribution, but as he emphasized, it also very frequently exhibits an exponent in the neighborhood of −1. That is, letting f 1 ≥ f 2 ≥ . . . f r · · · ≥ ... represent the observed ranked word frequencies, Zipf found f r ≈ Cr −1 , with C > 0 a constant depending on the sample size. Many qualifications have been raised regarding this approximation (including the issue of the divergence of the harmonic series); yet as I.J. Good [6] commented: "The Zipf law is unreliable but it is often a good enough approximation to demand an explanation." Figure 1 illustrates the fascinating universal character of this word frequency law using the text from four different authors writing in four different European languages in four different centuries. The plots are shown in the log-log form that Zipf employed, and they resemble countless other examples that researchers have studied since his time. The raw text files for the four books used in the graphs were downloaded from the Public Gutenberg databank [7]. Figure 1. Log-log plots of relative word frequencies by rank for four authors writing in four different European languages in four different centuries. The approximate -1 slope in all the graphs is an iconic feature of Zipf's word frequency law.
The monkey model has a convoluted history. The model can be thought of as actually embedded within the combinatoric logic of Mandelbrot's early work providing an informationtheoretic explanation of Zipf's law. For example, in Mandelbrot [8] there is an appendix entitled "Combinatorial Derivation of the Law p r = P r −B " (his notation) that effectively specifies a monkey model, including a Markov version, but is not explicitly stated as such. However, his derivation there is informal, and the first completely rigorous analysis of the general monkey model with independently combined letters was given only surprisingly recently by Conrad and Michenmacher [9]. Their analysis utilizes analytic number theory and is, as a result, somewhat complicated -a point they note at the end of their article. A simpler analysis using only elementary methods based on the Pascal pyramid has now been given by Bochkarev and Lerner [10]. They have also analyzed the more general Markov problem [11] and hidden Markov models [12]. Edwards, Foxall and Perkins [13] have provided a directly relevant analysis in the context of scaling properties for paths on graphs explaining how the Markov variation can generate both a power law or a weaker scaling law, depending on the nature of the transition matrix. Perline [4] also gives a simple approach to the demonstration of power law behavior.
A clarifying and important milestone in the history of this topic came from the cognitive psychologist, G. A. Miller [14,15]. Miller used the simple model with equal letter probabilities and an independence assumption to show that this version generates a distribution of word probabilities that conforms to an approximate inverse power law such that the probability P r of the r th largest word probability is (in a sense that needs to be made precise) of the form P r ≈ Cr −β , for −β < −1. In addition, he made the very interesting observation that by using a keyboard of K = 26 letters and one space character, with the space character having a probability of . 18 (about what is seen in empirical studies of English text), the value of −β in this case is approximately -1.06, very close to the nearly -1 value in samples of natural language text as illustrated in Figure 1. In his simplified model, it turns out that −β = −1 + ln(1 − s)/ ln K, where 0 < s < 1 is the space probability, so that −β approaches -1 from below as K increases. (Natural logarithms are denoted ln(x) and logarithms with any other radix will be explicitly indicated, as in log K (x). Note that Miller's β involves a ratio of logarithms, so that it is invariant with respect to the radix used in the numerator and denominator.) Consequently, Miller's model not only mimicks the form of Zipf's law for real languages, but with a large enough alphabet size it even mimicks the key parameter value. In this article we give a broad generalization of Miller's observation to the case of unequal letter probabilities. In light of our present results given in Section 2, the application of sample spacings from a uniform distribution in [2] to study the −1 behavior using an asymptotic regression line should now be viewed as just a step in a direction that ultimately led us to the Shao and Hahn limit used here.

2.
Proof that −β Tends Towards −1 from the Asymptotics of Log-Spacings For our analysis we specify a keyboard with an alphabet of K ≥ 2 letters {L 1 . . . , L K } and a space character S. The letter characters have non-zero probabilities q 1 , q 2 , . . . , q K . The space character has probability s, so that K i=1 q i + s = 1. A word is defined as any sequence of non-space letters terminating with a space. A word W of exactly n letters is a string such as W = L i 1 L i 2 . . . L in S and has a probability of the form P (W ) = P = q i 1 q i 2 . . . q in s because letters are struck independently. The space character with no preceding letter character will be considered a word of length zero. The rank ordered sequence of descending word probabilities in the ensemble of all possible words is written P 1 = s > P 2 ≥ . . . P r ≥ . . . (P 1 = s is always the first and largest word probability.) We break ties for words with equal probabilities by alphabetical ordering, so that each word probability has a unique rank r.
Conrad and Mitzenmacher [9] give a carefully constructed definition of power law behavior in the monkey model as the situation where there exist two positive constants C 1 and C 2 such that the inequality holds for sufficiently large r. Following the argument in Bochkarev and Lerner [10], β in the monkey model turns out to be simply the solution to with β > 1. 1 In the Miller model with equal letter probabilities, q 1 = · · · = q K = (1 − s)/K, so from equation (2), β in this case is found to be 1 − ln(1 − s)/ ln K. In the Fibonacci example given by Conrad and Mitzenmacher [9] and in Mitzenmacher [16] they use K = 2 letters with probabilities q 1 , q 2 = q 2 1 and q 1 < (−1 + To understand conditions leading to −β ≈ −1, we define spacings through a random division of the unit interval and then state the Shao-Hahn limit law. Let X 1 , X 2 , . . . , X K−1 be a sample of K − 1 i.i.d. random variables drawn from a distribution on [0, 1] with a bounded density function h(x) ≤ M < ∞. 2 Write the order statistics of the sample as X 1:K−1 ≥ X 2:K−1 · · · ≥ X K−1:K−1 . The K spacings D i are defined as the differences between the successive order statistics: We'll refer to this as a generalized broken stick process. By Shao and Hahn [3] Corollary 3.6, we have as K → ∞ and where a.s. signifies almost sure convergence, λ = .577 . . . is the Euler so dividing through by ln K gives as K → ∞. The right side of the limit in (6) goes to 0 because − 1 0 h(x) ln h(x)dx/ ln K goes to 0 by the boundedness of the density h(x). Expressing logarithms with a radix = K then leads to the limit Our universality property for β will now follow almost immediately from this. We use sample spacings to populate the K letter probabilities for the monkey keyboard. Since s is the probability of the space character, define Then from the limit (7), we have To see that this is the case, note: where the inequality in (10) follows from the geometric-arithmetic mean inequality. Therefore, from we see that m K ≤ −β and so −β a.s.
In the special case of Miller's model using all equal letter probabilities, m K = −β.
The broad generality of the Shao-Hahn limit leads to the near-universal characterization of the approximation −β ≈ −1 as K increases. Figure 2 presents graphical results using simulations of sample spacings from several distributions to illustrate this phenomenon. In Figure 2(a) we plot log 10 P r by log 10 r for the Miller model with exactly the parameter values he used, i.e., K = 26 letters, a space probability s = .18 and equal letter probabilities q 1 = · · · = q 26 = (1 − .18)/26. The graph is based on the ranks 1 ≤ r ≤ 475255 = 4 i=0 26 j , corresponding to all words of length ≤ 4 non-space letters. Figures 2(b) through 2(d) are based on generating K = 26 letter probabilities from three different continuous distributions with bounded densities h(x) on [0, 1] using a generalized broken stick method to obtain the sample spacings. Again s = .18 was used for the probability of the space character and the letter probabilities were populated with values from the spacings with q i = (1 − s)D i , 1 ≤ i ≤ 26, so that in each case, 26 i=1 q i = .82, as in Miller's example. For each graph in Figure 2(b) -2(d), we generated the largest 475255 word probability values in order to match the Miller example of Fig 2(a). The three continuous distributions, all defined on [0, 1], are: -a unform distribution with density h(x) = 1; -a beta B(3, 2) distribution with density -a triangular distribution with density The graphs in Figures 2(b)-2(d) illustrate our theoretical derivation, but we also note that the very linear plots indicate what appears to be an almost "immediate" convergence to power law behavior exceeding what might be expected from the asymptotics of Conrad-Mitzenmacher and Bochkarov-Lerner.

Monkey Twitter: Anscombe's CLT for the Model with a Finite Word Length Cutoff
The discussion of the finite word length version of the monkey model in [2] was incomplete because it did not explicitly provide the normalizing constants for the application of Anscombe's CLT for a random number of i.i.d. variables, as we do here. We also want to explain more clearly the nature of the hybrid Zipf-lognormal distribution that results from a finite length cutoff when letter probabilities are not identical.
The focus in Section 2 on the infinite ensemble of word probabilities P ∞ ranked to give P 1 > P 2 ≥ P 3 . . . has actually obscured the significant hierarchical structure of the monkey model -what Mandelbrot [17] (p. 345), called a lexicographic tree -which only becomes evident when word length is considered. To see this, write P len=n for the multiset of all word probabilities for monkey words of exactly length n. There are K n probabilities in P len=n having a total sum of (q 1 + q 2 + · · · + q K ) n s = (1 − s) n s and an average value of ((1 − s)/K) n s, declining geometrically with n. Now taking word length into account, Mandelbrot's lexicographic tree structure becomes clear with the simple case of K = 2 letters: the root node of the tree has a probability of s (for the "space" word of length 0); it has two branches to the next level consisting of nodes for words of length 1 with probabilities q 1 s and q 2 s; these each branch out to the next level with four nodes having probability values q 2 1 s, q 1 q 2 s, q 2 q 1 s, q 2 2 s and so on. The relevance of the CLT is seen first by representing any probability P ∈ P len=n as a product of i.i.d. random variables times the constant s: P = X 1 X 2 . . . X n s, where each X i takes on one of the letter probability values q 1 , q 2 , . . . , q K with probability 1/K (i.e., we use the natural counting measure to construct a probability space on P len=n ). Let µ 1 = K i=1 ln q i /K and σ 2 1 = K i=1 (ln q i − µ 1 ) 2 /K. Assume that σ 2 1 > 0, i.e., the letter probabilities are not all equal. Then µ 1 + ln s is the mean and σ 2 1 is the variance of ln P for P ∈ P len=1 , and it follows that nµ 1 +ln s and nσ 2 1 are the respective mean and variance of ln P = S n = n i=1 ln X i +ln s for P ∈ P len=n . Therefore, (S n − nµ 1 )/ nσ 2 1 is asymptotically normally distributed N (0, 1) (the term ln s is rendered negligible in the asymptics) so that for sufficiently large n, P itself will be approximately lognormal LN (nµ 1 , nσ 2 1 ). This approximate normality for the log-probabilities of words of fixed length is quite obvious and was noted in passing by Mandelbrot [18] (p. 210). However, he missed a stronger observation that the N n = n i=0 K i word probabilities for P ∈ P len≤n = ∪ n i=0 P len=i have a distribution that behaves in its central part away from the tails very much like P ∈ P len=n . That is, a version of the CLT can be applied to words of length ≤ n, not just to words of length n, and the two distributions, in one sense, are very close to each other.
To explain why this is so, for P ∈ P len≤n represent ln P = S Rn = Rn i=1 ln X i + ln s as a sum of a random number of i.i.d. random variables (plus the constant ln s), where R n is itself a random variable with a finite geometric distribution, Prob{R n = i} = K i /N n , (1 ≤ i ≤ n). We will not repeat the calculations of [2] here, but it is easily shown that R n /N n → 1 in probability. Because this limiting constant is 1, Anscombe's generalization of the CLT [19] (Theorem 3.1, p. 15) can be applied with the identical normalizing constants nµ 1 and nσ 2 1 as used above with ln P for P ∈ P len=n . Consequently, it is also true that for P ∈ P len≤n , the normalized sum (S Rn − nµ 1 )/ nσ 2 1 has an asymptotic N (0, 1) distribution. In other words the two random sums S n and S Rn behave so similarly in their centers that the word probabilities in P len=n and P len≤n are both approximately LN (nµ 1 , nσ 2 1 ). However, it is essential to remark that the two distributions have very different upper tail behavior. In this regard, Le Cam's [20] comment that French mathematicians use the term "central" referring to the CLT "because it refers to the center of the distribution as opposed to its tails" is particularly relevant.
The behavior of P ∈ P len≤n is illustrated in the graphs of Figure 3(a) and 3(b). The plot in Figure 3(a) was generated just as in Figure 2(b) except that a finite length cutoff of n ≤ 4 letters has been applied. To make this clear, in Figure 2(b) the plot shows the top 475255 word probabilities in P ∞ generated using K = 26 letters derived from uniform spacings. In contrast, using the same letter probabilities, the plot in Figure 3(a) shows the 475255 word probabilities generated with word length ≤ 4 letters, i.e., all the values of P ∈ P len≤4 . Word length in the case of unequal letter probabilities is certainly correlated, but not perfectly, with word probability: words of shorter length will, on average, have a higher probability than those of longer length, but except in Miller's degenerate case, there will always be reversals where a longer word will have a higher probability than a shorter word. 3 However, writing P 1:Nn > P 2:Nn · · · ≥ P Nn:Nn for the ranked values of P len≤n , it should be evident that for any given rank r, P r:Nn ≤ P r and that P r:Nn = P r when n is sufficiently large. In short, P len≤n inherits its upper tail power law behavior from P ∞ , which is illustrated by comparing the top part of the curve in Figure 3(a) with the corresponding part of Figure 2(b). Figure 3(b) uses the same data points from Figure 3(a) graphed as a normal quantile plot, and its roughly linear appearance for the logarithm of word probabilities for P ∈ P len≤n conforms to an approximate Guassian fit, although certainly the bending in the upper half of the distribution departs a bit from the linear trend of the lower half. The fact that a distribution can have a power law tail and lognormal central part and still look lognormal over essentially its entire range may seem surprising. The discussion in [21] about power law mimicry in the upper tail of lognormal distributions helps to explain why this can happen. We will call the distribution for P len≤n a Zipf-lognormal distribution.
The monkey model with a fixed word length cutoff can prove useful as a motivating idea. In the next section, we will discuss how models with the same branching tree structure have been proposed many times in the past, typically in settings where something like a finite word length is natural to consider. For the moment, think of the social networking service, Twitter, which allows members to exchange messages limited to at most 140 characters. Define Monkey Twitter with a finite limit of n + 1 characters. For convenience, in any implementation of a Monkey Twitter random experiment we will assume that monkeys would always fill up their allotted message space of n + 1 characters. Monkey words still require a terminating space character, so it is possible for a monkey to type (at most) one "non-word," which can vary in length from 1 to n + 1 characters, and will always be the last part of a message string. Non-words are discarded, and the probabilities for the legitimate monkey words, varying in length from 0 to n non-space characters plus a terminating space character, will correspond to the values in P len≤n . Figure 3. Figure 3(a) is a log-log plot of monkey word probabilities by rank for all words of length ≤ 4 non-space characters using letter probabilities from uniform spacings. The linear upper tail coincides closely with the previous Figure 2(b), but the power law clearly breaks down. Figure 3(b) shows the same word probabilities in a standard normal quantile plot. Its rough linearity confirms an approximate Gaussian fit over the whole distribution even though the upper tail is a power law as seen on the left in Figure 3(a). We refer to this distribution as a Zipf-lognormal hybrid, and it has many connections to other distributions discussed in the statistical literature.

Connections to Other Work
The subject of power law distributions is a vast topic of research reaching across diverse scientific disciplines [16,21,22,23]. Here we provide a brief sketch of how our results connect to some other work and how they can be considered in a much more general light.
The monkey model viewed in terms of its hierarchical tree structure is the starting point for understanding its general nature. Though Miller did not explicitly present his model as a branching tree, several researchers at almost the same time in the mid-to-late 1950s were highlighting this form using the equivalent of equal letter probabilities to motivate the occurrence of empirical power law distribution in other areas. For example, Beckman [24] introduced this idea in connection with the power law of city populations within a country ("Auerbach's law" [25]). Using essentially the same logic as Miller, but in a completely different setting, Beckman assumed that a given community will have K satellite communities, each with a constant decreasing fraction of the population at a higher level. That is, if A is the population of the largest city and 0 < p < 1, there will be K nearby satellite communities with population pA, K 2 smaller communities nearby to these with population p 2 A, and so on. Mandelbrot [26] (p. 226) had an apt expression for this pattern: he described it as "compensation between two exponentials" because the number of observations at each level increases geometrically while the mean value decreases geometrically. Beckman then went on to note that if instead of using a constant decreasing fraction p, one used a random variable X on (0, 1), the population at the n th level down would be a random variable of the form X 1 X 2 · · · X n A, leading to an approximate lognormal distribution of populations at this level. This corresponds to our discussion of monkey word probabilities in P len=n , although Beckman was not aware of the still stronger statement of lognormality for the probabilities of words of length ≤ n in P len≤n that we demonstrated using Anscombe's CLT.
Many more examples of a branching tree structure essentially equivalent to Miller's monkey model with equiprobable letters have been proposed over the years to motivate the occurrence of a huge variety of empirical power law distributions, including such size distributions as lake areas, island areas ("Korcak's law"), river lengths, etc. Indeed, Mandelbrot's [17] classic book on fractals is a rich source of these. This equiprobability case is so simple to analyze that it has been invoked over and over again to illustrate how power law behavior can arise. However, the more complicated case using unequal letter probabilities (or proportions) and a finite word length cutoff is what really underscores the close analogy between the monkey Zipf-lognormal distribution and other mixture models exhibiting hybrid power law tails and approximately lognormal central ranges.
Empirical distributions of this form have appeared with great frequency. In fact there is a long history of researchers, including Pareto, first discovering what appears to be an empirical power law over an entire distribution because they start out by looking at only the largest observations (cities, corporations, islands, etc.). However, when they extend their measurements to the part of the distribution below the upper tail, which is always more difficult to observe, the power law behavior typically breaks down. Several examples of this are chronicled in [21]. The power law for city sizes noted above is a perfect illustration. Auerbach in 1913 [25] looked at the 94 largest German cities from a 1910 census and showed a good power law fit, but because he did not include the thousands of smaller communities, he may not have been aware that this relationship breaks down. Modern studies of the full distribution of communities, such as Eeckhout's [27] analysis of the populations of 25,359 places from U.S. Census 2000, prove with high confidence that the bulk of the community populations fit an approximate lognormal distribution and that the power law behavior is confined to the upper tail.
Montroll and Shlesinger [28,29] explained how to generate a Pareto-lognormal hybrid for income distributions by using a hierarchical mixture model of lognormal distributions. This is motivated from several angles, including the notion that higher classes amplify their income by organizing in such a way as to benefit from the efforts of others. Reed and Hughes [30] have provided a far-ranging framework for understanding these hybrid distributions across the entire spectrum of disciplines where they have been discovered: "physics, biology, geography, economics, insurance, lexicography, internet ecology, etc." In [30] he and Hughes give a condensed summary showing that "if stochastic processes with exponential growth in expectation are killed (or observed) randomly, the distribution of the killed or observed state exhibits power-law behavior in one or both tails." This work encompasses: (1) geometric Brownian motion (GBM); (2) discrete multiplicative processes; (3) homogeneous birth and death processes; and (4) Galton-Watson branching processes. In one of many variations on this theme, in his GBM model [31] Reed specifies a stochastic process that uses lognormal distributions varying continuously in time with an exponential mixing distribution and shows that this generates a "Double Pareto-Lognormal Distribution." This has an asymptotic Pareto upper tail and a central part approximately lognormal; in addition, it exhibits interesting asymptotic behavior in the it lower tail where it is characterized by a direct (rather than an inverse) power law. Reed has gone to great effort to demonstrate the high quality of the fit of this distribution for all ranges of values (low, middle, high) across a wide variety of size distributions such as particle and oil field sizes [31], incomes [32], internet file sizes and the sizes of biological genera [30].
In today's nomenclature, the term Zipf 's law has come to mean any power law distribution with an exponent close to the value −1, not just the word frequency law. (To be clear here, when we refer to power law distributions, we also mean the hybrids with power law tails that we have been considering.) The subset of power laws with this restricted exponent value is surprisingly large and includes the distribution of firm sizes [33], city sizes [34], the famous Gutenberg-Richter law for earthquake magnitudes [35] and many others.
Gabaix [34] pointed out that while it has long been known that stochastic growth processes could generate power laws, "these studies stopped short of explaining why the Zipf exponent should be 1." To address this question, he has given a theoretical explanation of the genesis of a −1 exponent for the Zipf law for city populations, but in fact, his approach can be regarded more generally. His key idea is a variation of Gibrat's [36] Law of Proportion Effect: using a fixed number of normalized city population quantities that sum to 1 and assuming growth processes (expressed as percentages) randomly distributed as i.i.d. variables with a common mean and variance, he proves that a steady state will satisfy Zipf's law. This proof requires an additional strong assumption in order to reach a steady state, namely, a lower bound for the (normalized) population size. Gabaix goes on to discuss relaxed versions of his model and how it fits into the larger context of similar work done by others.
In their monograph on Zipf's law Saichev et al [37] modify and extend the core Gabaix idea in numerous ways that render it more realistic. For concreteness, they present their work using the terminology of financial markets and corporate asset values, but they are very clear on how their models are relevant to a broad range of "physical, biological, sociological and other applications." As with Reed, GBM plays an important role in their investigations, which incorporate birth and death processes, acquisitions and mergers, the subtleties of finite-size effects and other features. Notably, they focus on the sets of conditions that lead to an approximate −1 exponent.
In both Harrëmoes and Topsøe [38] and Corominas-Murtra and Solé [39] −1 emerges as a consequence of certain information theoretic ideas pertaining to the entropy function under growth assumptions. In [38] the focus is on language and the expansion of vocabulary size while simultaneously maintaining finite entropy. The perspective in [39] is broader, but still based on the idea of systems with an expanding number of possible states "characterized by special features on the behavior of the entropy." There is another topic area of statistical research that impacts on our work. Natural language word frequency distributions have been characterized as Large Number of Rare Events (LNRE) distributions because when collections of text are examined, no matter how big, there are always a large number of words that occur very infrequently in the sample [40]. LNRE behavior indicates a tip-of-the-iceberg phenomenon resulting from a bias that captures the most common words, but necessarily misses a vast quantity of rarely used words. This is a classic and much studied problem encountered in species abundance surveys in ecological research, where a key question becomes estimating the size of the zero abundance class -the typically great number of species not observed in a sample because of their rarity [41].
To get an idea of the significance of this issue, consider that Figure 3(a) graphs the population (or parent or theoretical) word probabilities for the Zipf-lognormal distribution, and not the sample frequencies as would be obtained from actually carrying out the Monkey Twitter experiment. Simple visual inspection of Figure 3(a) indicates that the linear part of the graph (i.e., the power law tail) holds for about 5 logarithmic decades or about the first 100000 word probabilities out of the total of 475255 probabilities plotted there. However, these 100000 probabilities comprise 97.4% of the total probability mass of all 475255 values. Intuitively, it should be clear that sampling from this Zipf-lognormal population distribution will be dominated by the Zipf part, not the lognormal part. Unless the sample size was astronomically large, so that the large number of low probability words showed up, the underlying structure of the parent distribution would not reveal itself. To carry this idea still further, imagine the situation if the word length cutoff was on the order of n = 140 characters, such as with real Twitter. No experiment could be run within any realistic time frame to ever hope to obtain a sample sufficiently large to uncover the true hybrid character of the population distribution -the sample would always appear as a Zipf law, not a Zipf-lognormal law. This kind of visibility bias has been a constant and recurring theme in all areas where Pareto-Zipf type distributions have been studied [2,21]. We believe its significance has been poorly appreciated in relation to this topic.
Finally, we will take a more bird's-eye view of matters and remark that our application of random spacings is very much in the spirit of the enormously fruitful study of random matrix theory and a class of stochastic systems referred to as the KPZ universality class. Along these lines Borodin and Gorin [42] have discussed a variety of probabilistic systems "that can be analyzed by essentially algebraic methods," yet are applicable to a broad array of topics. Miller's demonstration of power law behavior and an approximate −1 exponent with increasing K for the monkey model is in a similar vein. We regard this result as analogous to the Borodin and Gorin example of the De Moivre-Laplace proof of the CLT for a sequence of i.i.d. Bernoulli trials, which depends on having an explicit pre-limit distribution and then taking "the limit of the resulting expression." (We want to add, however, that Miller's model is ultimately too simple to reveal the approximate lognormal behavior of the monkey word probabilities -for that, we needed to assume non-identical letter probabilities.) Our point is to make clear that the monkey model fits into a much larger conceptual scheme than appears at first glance.

Conclusions
Sample spacings provide a natural way to populate the letter probabilities in the monkey model through a random division of the unit interval. The Shao and Hahn asymptotic limit law for the logarithms of spacings then leads to the result that the exponent in the power law of word frequencies will tend towards a −1 value under broad conditions as the number of letters in the alphabet increases. The monkey model can also be viewed as a branching tree structure. In that light, Anscombe's CLT reveals an underlying lognormal central part of the word frequency distribution when a finite word length cutoff is imposed. The resulting hybrid Zipf-lognormal distribution has many connections to other work. The visibility bias inherent in sampling from this distribution is similar to what has been noted as a characteristic of many different empirical power laws.