Page:EB1911 - Volume 22.djvu/406

Rh aggregate is given by a corresponding term in the expansion of (q + p), and by a well-known theorem this term is approximately equal to $1⁄√2npq$e− 2/2npq ; where is the number of integers by which the term is distant from np (or an integer close to np); provided that is of (or <) the order √n. Graphically, let the sortition made for each element be represented by the taking or not taking with respective frequency p and q a step of length i. If a body starting from zero takes successively n such steps, the point at which it will most probably come to a stop is at npi (measured from zero); the probability of its stopping at any neighbouring point within a range of ± √ni is given by the above-written law of frequency, i being the distance of the stopping-point from npi. Put and  then the probability may be written $$(1/\sqrt{\pi}c)\exp{-x^2/c^2}$$.

104. It is a short step, but a difficult one, from this case, in which the element is binomial—heads or tails—to the general case, in which the element has several values, according to the law of frequency—consists, for instance, of the member of points presented by a randomly-thrown die. According to the general theorem, if Q is the sum of numerous elements, each of which assumes different magnitudes according to a law of frequency, the function f being in general different for different elements, the number of times that Q assumes magnitudes between x and x + ∆x in the course of N trials is Nz∆x, if $$z=(1/\sqrt{2{\pi}k})\exp{-(x-a)^2/2k}$$; where a is the sum of the arithmetic means of all the elements, any one of which the square brackets denoting that the integrations extend between the extreme limits of the element's range, if the frequency-locus for each element is continuous, it being understood that [∫fr(x)dx] = 1; and k is the sum of the mean squares of error for each element, = ∑[∫2fr(ar + )d], if the frequency-locus for each element is continuous, where ar is the arithmetic mean of one of the elements, and the deviation of any value assumed by that element from ar, ∑ denoting summation over all the elements. When the frequency-locus for the element is not continuous, the integrations which give the arithmetic mean and mean square of error for the element must be replaced by summations. For example, in the case of the dice above instanced, the law of frequency for each element is that it assumes equally often each of the values 1, 2, 3, 4, 5, 6. Thus the arithmetic mean for each element is 3.5, and the mean square of error {(3.5 − 1)2 + (3.5 − 2)2 + &c.}/6 = 2.916. Accordingly, the sum of the points obtained by tossing a large number, n, of dice at random will assume a particular value x with a frequency which is approximately assigned by the equation

The rule equally applies to the case in which the elements are not similar; one might be the number of points on a die, another the number of points on a domino, and so on. Graphically, each element is no longer represented by a step which is either null or i, but by a step which may be, with an assigned probability, one or other of several degrees between those limits, the law of frequency and the range of i being different for the different elements.

105. Variant Proofs.—The evidence of these statements can only be indicated here. All the proofs which have been offered involve some postulate as to the deviation of the elements from their respective centres of gravity, their “errors.” If these errors extended to infinity, it might well happen that the law of error would not be fulfilled by a sum of such elements. The necessary and sufficient postulate appears to be that the mean powers of deviation for the elements, the second (above written) and the similarly formed third, fourth, &c., powers (up to some assigned power), should be finite.

106. (1) The proof which seems to flow most directly from this postulate proceeds thus. It is deduced that the mean powers of deviation for the proposed representative curve, the law of error (up to a certain power), differ from the corresponding powers of the actual locus by quantities which are negligible when the number of the elements is large. But loci which have their mean powers of deviation (up to some certain power) approximately equal may be considered as approximately coincident.

107. (2) The earliest and best-known proof is that which was

originated by Laplace and generalized by Poisson. Some idea of this celebrated theory may be obtained from the following free version, applied to a simple case. The case is that in which all the elements have one and the same locus of frequency, and that locus is symmetrical about the centre of gravity. Let the locus be represented by the equation =, where the centre of gravity is the origin, and the construction signifying that the probability of the element having a value (between say − ½∆ and + ½∆ is ∆. Square brackets denoting summation between extreme limits, put (a) for [Se√−1a∆] where  is an integer multiple of ∆ (or ∆x) = ∆x, say. Form the mth power of (a). The coefficient of e√−1ar∆x in ((a))m is the probability that the sum of the values of the m elements should be equal to r∆x; a probability which is equal to ∆xyr, where y is the ordinate of the locus representing the frequency of the compound quantity (formed by the sum of the elements). Owing to the symmetry of the function the value of yr, will not be altered if we substitute for e√−1ar∆x, e−√−1ar∆x, nor if we substitute ½(e+√−1ar∆x + e−√−1ar∆x), that is cos ar∆x. Thus ((a))m becomes a sum of terms of the form ∆xyr cos ar∆x, where y−r = y+r. Now multiply ((a))m thus expressed by cos t∆xa, where, t being an integer, t∆x =x, the abscissa of the “error” the probability of whose occurrence is to be determined. The product will consist of a sum of terms of the form ∆xyr ½(cos a(r + t)∆x + cos a(r − t)∆x). As every value of r − t (except zero) is matched by a value equal in absolute magnitude, −r + t, and likewise every value of r + t is matched by value −r − t, the series takes the form ∆xyr∑ cos qa∆x + ∆xyt, where q has all possible integer values from 1 to the largest value of |r| increased by |t|; and the term free from circular functions is the equivalent of ∆xyr cos a(r + t)∆x, when together with ∆xyr cos a(r − t)∆x, when Now substitute for a∆x a new symbol ; and integrate with respect to , the thus transformed ((a))m cos t∆xa between the limits = 0 and =. The integrals of all the terms which are of the form ∆xyr cos q will vanish, and there will be left surviving only ∆xyt. We thus obtain, as equal to ∆xyt, $$\int_0^\pi\{\chi(\beta/{\Delta}x)\}^{m}\cos{t\beta}\,d\beta$$. Now change the independent variable to a; then as

Replacing t∆x by x, and dividing both sides by ∆x, we have

Now expanding the cos ax which enters into the expression for (a), we obtain

Performing the summations indicated, we express (a) in terms of the mean powers of deviation for an element. Whence (a)m is expressible in terms of the mean powers of the compound locus. First and chief is the mean second power of deviation for the compound, which is the sum of the mean second powers of deviation for the elements, say k. It is found that the sought probability may be equated to $\int_0^{\pi/\Delta{x}}{dx}e^{-\frac{1}{2}a^2k}\cos{ax}+\frac{1}{4!}k_2\int_0^{\pi/\Delta{x}}{dx}a^4e^{-\frac{1}{2}a^2k}\cos{ax}$ -. . .,|undefined where k2 is the coefficient defined below. Here /∆x may be replaced by ∞, since the finite difference ∆x is small with respect to unity when the number of the elements is large; and thus the integrals involved become equateable to known definite integrals. If it were allowable to neglect all the terms of the series but the first the expression would reduce to $1⁄2!$e−u 2/k, the normal law of error. But it is allowable to neglect the terms after the first, in a first approximation, for values of x not exceeding a certain range, the number of the elements being large, and if the postulate above enunciated is satisfied. With these reservations it is proved that the sum of a number of similar and symmetrical elements conforms to the normal law of error. The proof is by parity extended to the case in which the elements have different but still symmetrical frequency functions; and, by a bolder use of imaginary quantities, to the case of unsymmetrical functions.