Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
|
|
- Ambrose Holmes
- 5 years ago
- Views:
Transcription
1 A Renewal Problem with Bulk Ordering of Components Author(s): D. R. Cox Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 21, No. 1 (1959), pp Published by: Wiley for the Royal Statistical Society Stable URL: Accessed: :15 UTC JSTOR is a not-for-profit service that helps scholars, researchers, and students discover, use, and build upon a wide range of content in a trusted digital archive. We use information technology and tools to increase productivity and facilitate new forms of scholarship. For more information about JSTOR, please contact support@jstor.org. Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at Royal Statistical Society, Wiley are collaborating with JSTOR to digitize, preserve and extend access to Journal of the Royal Statistical Society. Series B (Methodological)
2 180 [No. 1, A RENEWAL PROBLEM WITH BULK ORDERING OF COMPONENTS By D. R. Cox Birkbeck College, University of London [Received July, 1958] SUMMARY SUPPOSE that k components are in use simultaneously, that batches of n items are ordered, and that all items in use at one time must be from the same batch. After n - k + 1 failures have occurred the remaining k - 1 items are wasted. Properties of the system are investigated and different strategies of use compared. 1. INTRODUCTION Consider a system in which k similar components are in use together. Suppose further that a batch of n components, n > k, is obtained and that when n - k + 1 failures have occurred, so that only k - 1 unbroken components remain, a new batch must be secured. The remaining k - 1 components of the first batch are rejected. That is, the k components in use at one time must be alike, taken from the same batch. This set-up was suggested to me by Mr. A. Glaskin, then of the Hosiery and Allied Trades Research Association, in connection with women's stockings. Here k = 2, and it is assumed that any number, n, of identical stockings can be bought together and any two of them worn as a pair, but that when only one of the batch remains fit for use, it cannot be worn with stockings from a new batch and is therefore scrapped. Another application is to situations where, when the system is first manufactured, spare components must be ordered for the whole life of the system, and where, if not enough spares are ordered, a completely new batch of components must be obtained. The two aspects of this situation studied in this paper are first the appropriate choice of n, and second the advantages of trying to equalize wear on all unbroken components. More precisely we consider the following two strategies: A. Each component in turn is used for a very short period, the choice of components for use being such that at any time each component that has not failed has been in service the same total time. B. From the batch of n components, k are chosen randomly for use. When one fails, it is replaced by an unused component, and so on until no more remain. Of course, in practice it might very well be impossible to adhere at all strictly to A, but it is reasonable to consider it for theoretical work. It will be assumed that the life-times of different components are mutually independent random variables Xi,..., X., all with probability density function f(x) and mean [, and that the life-times are unaffected by the strategy of use; this last assumption may well be seriously wrong in some applications. Let T denote the life of the batch, i.e. the time up to the (n - k + 1)th failure.
3 1959] Bulk Ordering of Components EXPONENTIAL DISTRIBUTION OF LIFE-TIME When the distribution of life-time is exponential, f(x) = -le-x/, a very simple solution can be given. At any time, all components that have not failed are equivalent; their future life is independent of the length of time for which they have been used in the past. Therefore strategy of use is irrelevant. In any time interval (t, t + h), there is a constant chance kh/,u + o(h) of a failure on one of the k components in use. Thus T is the time up to the (n - k + 1)t event in a Poisson process with parameter k/,u. Therefore 2kT/u is distributed as X2 and in particular E(T)=-{ [(n -k + 1)/k. (1) In some applications, such as the second mentioned in section 1, we may wish to choose n so that T exceeds a given value to with high probability. If [t is known, a suitable n is easily found from tables of the x2 distribution. In other applications, such as the one to stockings, long-run behaviour is of most interest and (1) implies that components are used at a rate n/e(t)=-nk/[(n- k + 1)?] compared with a rate k/,u if the restriction to the use only of components from the batch were removed. Quite generally, for any3distribution of life-time with mean may define an index of efficiency I ke(t) (2) n,u measuring the ratio of the rate of use if the batch restriction were removed to that with batch size n. Equivalently I measures the fraction of the total life of all components that is utilized. For the exponential distribution, the index becomes Ie n-k 1 (3) n Equation (3) can be obtained directly from the consideration that at the end of a batch, the k - 1 components rejected are, with an exponential distribution of life-time, "as good as new". The loss measured by Ie has to be set against the advantage of spreading expenditure, the mean interval between purchases being (1), and against the cost of storing spares, the number of which varies between 0 and n - k, all values being equally frequent in the long run. In some applications the cost will depend on the number of batches bought as well as on the total number of components. If total costs can be assessed quantitatively, the optimum n can be found. 3. THE STRATEGY OF EQUALIZED WEAR Suppose now that the p.d.f. of life-time, f (x), is not necessarily exponential and that strategy A of section 1 is adopted. Denote by X(1),..., X(n) the order statistics of Xi,..., X., i.e. the n life-times rearranged in decreasing order. Now the first failure will occur after component life-time X(.), i.e. after clock-time nx(fl)/k. The second
4 182 Cox-A Renewal Problem with [No. 1, failure will occur after a further componen a further clock-time of (n - 1)(X(n-1) - X(n))/k. Finally the interval in life-time and the corresponding clock-time between penultimate and final failure are both X(k) - X(k+l). Therefore, if X(n+l) is defined to be 0 n-k n kt =, (n -j)(x(n-j) - X(n+l) X() + kx(k). (4) 1=0 1=k+1 The distribution o is a normal density function, the mean and variance of T can be found from the means, variances and covariances of the X(W), which are tabulated up to n = 20 (Teichroew, 1956; Greenberg and Sarhan, 1956); values of the means to two decimal places are given up to n 50 by Fisher and Yates (1957). If in a random sample of size n from a normal population of mean,t and standard deviation a, the expected value of X(W) is denoted by Z + mna, (5) then the index of efficienc where I ke(t) 1 n-g./., (6) n,u k-1 ng mn - (k -1) mkn. (7) ngk To obtain (7) the result n has been used. The entries (a) in Table 1 give gk for n < 8. The table is discussed later in the section. In order to examine how seriously the factors depend on the assumption of normality, it is worthwhile examining first a E(X(>) distribution of life-time +n-r+ that is rectangular 1 (8) over (- A, U + A). It is well known that in random samples of n from this distribution n and some simple calculation from (4), together with the fact that for the rectangular distribution a = V3 A, leads to the index of efficiency Ir, where Ir 1 -/V3 k(k - 1) a (9) n(n + 1) I- as the analogue of (6). A selection of values of the factor multiplying ai> in (9) is g in Table 1 (b).
5 1959] Bulk Ordering of Components 183 TABLE 1 Factor Determining the Index of Efficiency I for Various Distributions of Life-Time n k 2 (a) (b) (c) (d) (a) * (b) * - (c) (d) (a) - 1* (b). 1* * (c) (d) *375 5 (a) (b) (c) (d) (a) (b) * 1* (c)..... (d).... 0* *625 7 (a) (b) * *. (c) (d) * *750 8 (a) *. ** (b) (d) *- *- * To calculate the index I, (a) Normal distribution; factor is glk (b) Rectangular distribution; factor is V3 k(k - 1)/[n(n + 1)]. (c) Distribution a2xe-az; factor is h'. (d) Exponential distribution; factor is (k - 1)/n. Over the range tabulated, the factors for the normal and rectangular distributions are in quite close agreement. Bearing in mind the sort of accuracy that is likely to be required in applications, it seems reasonable to conclude that for these n and k, the exact shap of the frequency distribution of life-time, if roughly symmetrical, is for the present purpo not very important. A family of skew distributions, for which the expected values of order statistics are tabulated or easy to calculate does not seem to be available. If life-time X is log-normal, so that Y = log X is normal, it should be adequate, unless r or n - r is small and V(Y) large, to use E(X(r)) E(ey(r)) exp {E(Y(r)) + 1 V( Y(r))}, (10) where E(Y(r)) and V( Y(r)) refer to a normal sample and are
6 184 Cox-A Renewal Problem with [No. 1, (10) would be exact if Y(r) was normally distributed and is a good approximation in any case if V( Y(r)) < 1. This method is probably the best for numerical work with a particular application, because it is likely that many empirical distributions of life-time can be reasonably well approximated by a log-normal curve. It will, however, not be followed up further here. Instead, to obtain results directly comparable with those obtained in section 4 for strategy B, we consider first the p.d.f. a2xe-a, proportional to that of X2 and with me 2/a. In samples of n from this distribution, the ratio of the expected value of an statistic to the population mean is 00 E(X(r))/(2a-1) n 2(r-1)!(n-r)! 1. [1 - (x + 1) e-0]n-r (x + 1)r-le-rxx2dx 0 n! n-rs+r-1 2 * - E E (_ )s X 2( ) =0 t=o (r + s - 1)! (t + 1)(t + 2) (l ) s!(n- r -s)!(r+s- 1 -t)!(r+s)t+3 ( The series can be expressed in terms of confluent hypergeometric functions, but this does not seem helpful for numerical work, which- is straightforward from (11). Since the expected values may possibly be of use in some other application they are given for n < 5 in Table 2. The values are checked by noting that each column adds up to n. TABLE 2 Ratio to Population Mean of Expected Value of rth Order Statistic in Samples of n from a2xe-x n r * * * * * '351 The corresponding values of E(T) and hence the indices of efficiency, Ig, are easily calculated from (3). The indices themselves are given in Table 6 in section 5, but for comparison with the results for the normal and rectangular distributions, the facto such that Ig = 1 - hkj, (12) where a = 1 /V2, is given in Table responding factor, (k - 1)/n, for the exponential distribution. Again, there is moderately good agreement over this range with the other factors in the table, bearing in mind the wide range of distributional form. One further numerical result has been obtained for the special distribution (2e-x + -Lje-xI5), with n = 3, k = 2, I = This will be used in section 5.
7 1959] Bulk Ordering of Components 185 A final remark is that if n is large, the asymptotic theory of order statistics (Cramer, 1946, p. 367) can be used to approximate to L 4. THE STRATEGY OF REPLACEMENT ON FAILURE Suppose now that strategy B of section 1 is followed, components being replaced only on failure. We then have in effect k independent simple renewal processes in operation simultaneously, one arising from each of the k components in use originally. A result for such systems given by Cox and Smith (1954) can be used to obtain a first approximation to the index of efficiency* I' valid for large values of n/k. We consider the system at the (n - k + 1) th failure. For each of the k - 1 processes on which this last failure does not occur, there is a residual life-time, defined as the time remaining before the next failure on that process. Provided that the number of failures that have occurred on each process is not too small, these residual life-times will be those of an equilibrium renewal process, i.e. an ordinary renewal process considered at a random point a long time from the start of the process. The expectation of such a residual lifetime is ([2 + g2) /(2M). Thus the expected wasted time is (k - 1)([2 + cr2) /(2p) therefore I 1 (k- 1)(IL2 + g2) (13) 2nK2 (3 This formula is exact for exponentially distributed life-times. It seems likely that (13) is a reasonable approximation even when k/n is no a/i is fairly small. If so, and if n - k is large, there will be an appreciable gain from the use of strategy A. For then I is near unity, whereas I' is appreciably less than unity. To see what is involved in the exact evaluation of I', we use methods standard in renewal theory. Let F(x) be the distribution function and f(x) the density function of life-time, and let Fr(x) be the r-fold convolution of F(x) with itself. Let Nt be the number of failures in (0, t) in a single renewal process and let Zr denote the sum of a fixed number r of independent life-times. Then pr (Nt < r) = pr (Zr > t) = 1 Fr(t). (14) Thus the probability generating function of Nt, G(4, t), say is given by G(z; t) E- z rfr(t) + 1(15) r=o If an asterisk is used to denote a Laplace transform with respect to t, so that 00 G*(4, s) G(4, t) e-st dt, (16) we have the well-known equation that G*(, s) 1 - f*(s) (17) * A prime will be used for indices of eff
8 186 Cox-A Renewal Problem with [No. 1, on using the result that 0 e-stfr(t) dt = [f()]. (18) S Now let A(7) be the number of failures in (0, t) in the combined set of k renewal processes. Then G(k)(4, t), its probability generating function, is given by G(k)(4, t) = [G(4, t)]k, (19) since N(k) is the sum of k independent rando Denote by Tm the time in the combined sequence measured up to the mth failure; we require m - n - k + 1. Then by a reversal of the argument of (14) pr (Tm > t) = pr (N(7) < m), (20) from which it follows easily that 00 E Cm pr (Tm > t) = G(k)(4, t). (21) m=1 1 Now we are particularly interested in E(Tm) whxich can be written Go E(Tm) = f pr (Tm > t) dt. (22) Therefore E mnle m=1 0 an equation for the moment generating function of Tm can be obtained in a similar way. Thus, to determine E(Tn-k+l), and hence I', we need to invert (17) to obtain G(4, t), to raise this to the kth power, to integrate with respect to t and finally to pick out the appropriate power of C in (23). These formulae can be evaluated explicitly whenf*(s) is a simple rational function of s. Thus whenf(x) = xe-,f*(s) = 1/(1 + S)2, and by (17) from which and finally G*(4, s) -s+12_4 (24) G(4, t) =2Vg {(A/4 + 1) evl e-kt k G(k)(t (4)k/2 t = (r(v + )r (ly - I)k-r e(2rvg-kv\ E VnmE(T'm)= (k w +1) V -1)- (25) m=1 (4)Ok/2 (1 - ) r=o r k(l + C)2r,\,/ In particular, when k = 2, the right-hand side of (25) reduces to 4(5-4) /[4(1 -
9 1959] Bulk Ordering of Components 187 from which E(Tm) l m + + (26) It follows that the index I', which is given in general by f- ke(tn-k+l) (27) nflp where M is the mean life-time (and is equal to 2 for the above distribution), is ' 1-3. (28) 4n The approximate formula (13) with k = 2, g2/2 = -, is exact in this case, a rather surprising result. When k = 3, the corresponding generating function (25) becomes 2~(13-5) (29) 3(1-0)2 (9 - The approximate formula (13) reduces in this case to I' 1-3. (30) 2n The exact values of I' for n = 3, 4, 5 are , and , whereas the corresponding approximate values from (30) are 0 5, and 0 7 indicating excellent agreement even for the smallest value of n. The special case k = 2, n = 3, f(x) - Ie-x + -Lj-e-0/5 has also been investigated, the value of I' being 5/9. The general approximate formula (13) can be derived as a limit as n -o 00 with k fixed by arguments similar to those leading to (23), but the details are complicated and very heavy algebra seems necessary to get a closer approximation. If results are required for a particular distribution of life-time, and it is thought that the approximation (13) may not be adequate, the best procedure is either to fit a distribution with a rational Laplace transform and to work through the exact theory, or to use an empirical sampling method. The theoretical method is difficult if the fitted Laplace transform is of high degree, and it is therefore worth discussing very briefly suitable sampling methods. The best methods to use depend on the particular values of n and k and on the properties of the distribution being sampled. Consider, as an example, a normal distribution of life-times and suppose that n - k is sm4ll. Let X1,..., Xn be the n life-times i order of use, and let Tr = Tr(Xi,..., Xn) be the time up to the rth failure in the combined process. Now it will happen frequently that Tr is equal to the rth largest observation among X1,..., Xk, a random variable, Ur, of known expectation. We therefore determine by sampling the expectation of Zr = Ur - Tr. To do this, first take random sampling numbers for X1,..., Xk and find Ur. In simple cases, we can determine exactly pr (Zr = 0) and then can sample for Xk+?,..., Xn conditionally on Zr =# 0 Additional precision can be obtained, (a) by "splitting", i.e. by taking several sets of
10 188 Cox-A Renewal Problem with [No. 1, Xk+l,..., X. for those configurations for which a large contribution to E(Zr) is likely, and (b) by stratified sampling for X1,..., Xk, or alternatively by weighted sampling based on the range of X1,..., Xk. These remarks are intended simply as a general indication of the methods to be followed. As an example, E(T2) has been determined for a normal distribution of mean 10 and standard deviation 3, taking k = 2, n = 3. The mean of 20 determinations was 11 62? 0 032, giving I'-0 775? The value from the approximation (13) is The agreement is good considering that (13) is based on the assumption that n/k is large. A rough estimate is that the devices described above lead to a several hundredfold increase in precision, as compared with simple random sampling. 5. THE COMPARISON OF STRATEGIES The following results have been reached relevant to the comparison of strategies. (i) For the exponential distribution strategy does not matter. (ii) For the distribution a2xe-ax, the results collected in Table 3 show that for small n and k, strategy A gives a higher index of efficiency, but that the difference between strategies is rather small. (iii) The general approximate formulae of sections 3 and 4 indicate that strategy A should, in general, give the higher index of efficiency when g/u < 1. The index of efficiency is however near 1 when n> k, so that the main situation in which an appreciable difference is likely is when k/n is not small. (iv) The calculation for the distribution le-x + -L-je-xl5 shows that it is possible for strategy B to give the higher index. TABLE 3 Indices of Efficiency when the Distribution of Life-Time is a2xe-x n k 2 (a) (a) (a) (a) (b) (b) (b) (b) (c) 0625 (c) (c) 0812 (c) (a) (a) 0640 (a) (b) 0482 (b) (b) (c) 0 5 (c) (c) (a) 0402 (a) 0*557 (b) (b) (c) (c) 0* (a) (b) (c) 0400 (a) Strategy A: exact. (b) Strategy B: exact. (c) Strategy B: approximate. Now the optimum strategy need be neither A nor B. When the life of the system ends at the (n - k + 1)th failure, there are left the (k - 1) residual life-times for the items that have not failed. The best strategy will be the one that minimizes the expected value of
11 1959] Bulk Ordering of Components 189 the sum of these residual life-times. Now if the distribution function of life-time X is F(x) and the survivor curve, giving the probability of surviving to at least age x is FC(x) = - F(x), we can define the expectancy of life at age x to be 00 ef (x) E[X - x X > x] J' ydf(y) /Fc(x) -x x 00 = Fc(y) dy/fc(x), (31) x [ ~L {log je FJ(y) dy (32) An intuitively very reasonable strategy is to use at any instant those comp greatest expectancy of life. Now ef(x) is a decreasing function of x if and only if d2 00 j2 log F0(y) d < 0, dxx i.e. if and only if the integral of F0(y) is logarithmically concave. This condition is satisfied by many common distributions such as the normal and the x2 distributions with more than two degrees of freedom, so that strategy A would be expected to be optimum in these cases. The mixed exponential distribution that has been investigated has ef(x) an increasing function of x, so that strategy B would be expected to be the better in this case, as was indeed found. Finally, it must be repeated that all these comparisons depend on the assumption that the life of individual components is unaffected by strategy of use, so that, for example, there are no fatigue effects. REFERENCES Cox, D. R. & SMITH, W. L. (1954), "On the superposition of renewal processes", Biometrika, 41, CRMEmR, H. (1946), Mathematical methods of statistics. Princeton: University Press. FISHER, R. A. & YATES, F. (1957), Statistical tables for biological agricultural and medical research. Edinburgh: Oliver & Boyd. 5th ed. GREENBERG, B. G. & SARHAN, A. E. (1956), "Estimation of location and scale parameters by ord from singly and doubly censored samples, I", Ann. Math. Statist., 27, TEICHROEW, D. (1956), "Tables of expected values of order statistics and products of order st samples of twenty and less from the normal distribution", Ann. Math. Statist., 27, VOL. XXI. NO. 1. G
Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Regression Analysis when there is Prior Information about Supplementary Variables Author(s): D. R. Cox Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 22, No. 1 (1960),
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Some Applications of Exponential Ordered Scores Author(s): D. R. Cox Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 26, No. 1 (1964), pp. 103-110 Published by: Wiley
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Biometrika Trust Some Simple Approximate Tests for Poisson Variates Author(s): D. R. Cox Source: Biometrika, Vol. 40, No. 3/4 (Dec., 1953), pp. 354-360 Published by: Oxford University Press on behalf of
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Biometrika Trust Robust Regression via Discriminant Analysis Author(s): A. C. Atkinson and D. R. Cox Source: Biometrika, Vol. 64, No. 1 (Apr., 1977), pp. 15-19 Published by: Oxford University Press on
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
On the Estimation of the Intensity Function of a Stationary Point Process Author(s): D. R. Cox Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 27, No. 2 (1965), pp. 332-337
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
A Note on the Efficiency of Least-Squares Estimates Author(s): D. R. Cox and D. V. Hinkley Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 30, No. 2 (1968), pp. 284-289
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Biometrika Trust Some Remarks on Overdispersion Author(s): D. R. Cox Source: Biometrika, Vol. 70, No. 1 (Apr., 1983), pp. 269-274 Published by: Oxford University Press on behalf of Biometrika Trust Stable
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
On the Bound for a Pair of Consecutive Quartic Residues of a Prime Author(s): R. G. Bierstedt and W. H. Mills Source: Proceedings of the American Mathematical Society, Vol. 14, No. 4 (Aug., 1963), pp.
More informationEcological Society of America is collaborating with JSTOR to digitize, preserve and extend access to Ecology.
Measures of the Amount of Ecologic Association Between Species Author(s): Lee R. Dice Reviewed work(s): Source: Ecology, Vol. 26, No. 3 (Jul., 1945), pp. 297-302 Published by: Ecological Society of America
More informationInternational Biometric Society is collaborating with JSTOR to digitize, preserve and extend access to Biometrics.
400: A Method for Combining Non-Independent, One-Sided Tests of Significance Author(s): Morton B. Brown Reviewed work(s): Source: Biometrics, Vol. 31, No. 4 (Dec., 1975), pp. 987-992 Published by: International
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
The Interpretation of Interaction in Contingency Tables Author(s): E. H. Simpson Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 13, No. 2 (1951), pp. 238-241 Published
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
American Society for Quality A Note on the Graphical Analysis of Multidimensional Contingency Tables Author(s): D. R. Cox and Elizabeth Lauh Source: Technometrics, Vol. 9, No. 3 (Aug., 1967), pp. 481-488
More informationE. DROR, W. G. DWYER AND D. M. KAN
Self Homotopy Equivalences of Postnikov Conjugates Author(s): E. Dror, W. G. Dwyer, D. M. Kan Reviewed work(s): Source: Proceedings of the American Mathematical Society, Vol. 74, No. 1 (Apr., 1979), pp.
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Queues with Time-Dependent Arrival Rates: II. The Maximum Queue and the Return to Equilibrium Author(s): G. F. Newell Source: Journal of Applied Probability, Vol. 5, No. 3 (Dec., 1968), pp. 579-590 Published
More informationThe Periodogram and its Optical Analogy.
The Periodogram and Its Optical Analogy Author(s): Arthur Schuster Reviewed work(s): Source: Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character,
More informationINFORMS is collaborating with JSTOR to digitize, preserve and extend access to Management Science.
On the Translocation of Masses Author(s): L. Kantorovitch Source: Management Science, Vol. 5, No. 1 (Oct., 1958), pp. 1-4 Published by: INFORMS Stable URL: http://www.jstor.org/stable/2626967. Accessed:
More informationBiometrika Trust. Biometrika Trust is collaborating with JSTOR to digitize, preserve and extend access to Biometrika.
Biometrika Trust An Improved Bonferroni Procedure for Multiple Tests of Significance Author(s): R. J. Simes Source: Biometrika, Vol. 73, No. 3 (Dec., 1986), pp. 751-754 Published by: Biometrika Trust Stable
More informationB.N.Bandodkar College of Science, Thane. Random-Number Generation. Mrs M.J.Gholba
B.N.Bandodkar College of Science, Thane Random-Number Generation Mrs M.J.Gholba Properties of Random Numbers A sequence of random numbers, R, R,., must have two important statistical properties, uniformity
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More informationStatistics for scientists and engineers
Statistics for scientists and engineers February 0, 006 Contents Introduction. Motivation - why study statistics?................................... Examples..................................................3
More informationProperties of Rational and Irrational Numbers
Properties of Rational and Irrational Numbers September 8, 2016 Definition: The natural numbers are the set of numbers N = {1, 2, 3,...}, and the integers are the set of numbers Z = {..., 2, 1, 0, 1, 2,...}.
More informationBiometrika Trust. Biometrika Trust is collaborating with JSTOR to digitize, preserve and extend access to Biometrika.
Biometrika Trust A Stagewise Rejective Multiple Test Procedure Based on a Modified Bonferroni Test Author(s): G. Hommel Source: Biometrika, Vol. 75, No. 2 (Jun., 1988), pp. 383-386 Published by: Biometrika
More informationEstimation for Mean and Standard Deviation of Normal Distribution under Type II Censoring
Communications for Statistical Applications and Methods 2014, Vol. 21, No. 6, 529 538 DOI: http://dx.doi.org/10.5351/csam.2014.21.6.529 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation for Mean
More informationBiometrika Trust. Biometrika Trust is collaborating with JSTOR to digitize, preserve and extend access to Biometrika.
Biometrika Trust Discrete Sequential Boundaries for Clinical Trials Author(s): K. K. Gordon Lan and David L. DeMets Reviewed work(s): Source: Biometrika, Vol. 70, No. 3 (Dec., 1983), pp. 659-663 Published
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationk-protected VERTICES IN BINARY SEARCH TREES
k-protected VERTICES IN BINARY SEARCH TREES MIKLÓS BÓNA Abstract. We show that for every k, the probability that a randomly selected vertex of a random binary search tree on n nodes is at distance k from
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationDistributions of Functions of Random Variables. 5.1 Functions of One Random Variable
Distributions of Functions of Random Variables 5.1 Functions of One Random Variable 5.2 Transformations of Two Random Variables 5.3 Several Random Variables 5.4 The Moment-Generating Function Technique
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
On the Probability of Covering the Circle by Rom Arcs Author(s): F. W. Huffer L. A. Shepp Source: Journal of Applied Probability, Vol. 24, No. 2 (Jun., 1987), pp. 422-429 Published by: Applied Probability
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Biometrika Trust The Use of a Concomitant Variable in Selecting an Experimental Design Author(s): D. R. Cox Source: Biometrika, Vol. 44, No. 1/2 (Jun., 1957), pp. 150-158 Published by: Oxford University
More informationIt can be shown that if X 1 ;X 2 ;:::;X n are independent r.v. s with
Example: Alternative calculation of mean and variance of binomial distribution A r.v. X has the Bernoulli distribution if it takes the values 1 ( success ) or 0 ( failure ) with probabilities p and (1
More informationThe College Mathematics Journal, Vol. 16, No. 2. (Mar., 1985), pp
On Rearrangements of the Alternating Harmonic Series Fon Brown; L. O. Cannon; Joe Elich; David G. Wright The College Mathematics Journal, Vol. 16, No. 2. (Mar., 1985), pp. 135-138. Stable URL: http://links.jstor.org/sici?sici=0746-8342%28198503%2916%3a2%3c135%3aorotah%3e2.0.co%3b2-q
More informationPhilosophy of Science Association
Philosophy of Science Association Why Bohm's Theory Solves the Measurement Problem Author(s): Tim Maudlin Source: Philosophy of Science, Vol. 62, No. 3 (Sep., 1995), pp. 479-483 Published by: The University
More informationThe Econometric Society is collaborating with JSTOR to digitize, preserve and extend access to Econometrica.
On the Optimal Character of the (s, S) Policy in Inventory Theory Author(s): A. Dvoretzky, J. Kiefer, J. Wolfowitz Reviewed work(s): Source: Econometrica, Vol. 21, No. 4 (Oct., 1953), pp. 586-596 Published
More informationMath Review. for the Quantitative Reasoning measure of the GRE General Test
Math Review for the Quantitative Reasoning measure of the GRE General Test www.ets.org Overview This Math Review will familiarize you with the mathematical skills and concepts that are important for solving
More informationTHE COVARIANCE MATRIX OF NORMAL ORDER STATISTICS. BY C. S. DAVIS and M. A. STEPHENS TECHNICAL REPORT NO. 14 FEBRUARY 21, 1978
THE COVARIANCE MATRIX OF NORMAL ORDER STATISTICS BY C. S. DAVIS and M. A. STEPHENS TECHNICAL REPORT NO. 14 FEBRUARY 21, 1978 PREPARED UNDER GRANT DAAG29-77-G-0031 FOR THE U.S. ARMY RESEARCH OFFICE Reproduction
More informationCHAPTER 6 SOME CONTINUOUS PROBABILITY DISTRIBUTIONS. 6.2 Normal Distribution. 6.1 Continuous Uniform Distribution
CHAPTER 6 SOME CONTINUOUS PROBABILITY DISTRIBUTIONS Recall that a continuous random variable X is a random variable that takes all values in an interval or a set of intervals. The distribution of a continuous
More informationEstimation of Quantiles
9 Estimation of Quantiles The notion of quantiles was introduced in Section 3.2: recall that a quantile x α for an r.v. X is a constant such that P(X x α )=1 α. (9.1) In this chapter we examine quantiles
More informationFourier and Stats / Astro Stats and Measurement : Stats Notes
Fourier and Stats / Astro Stats and Measurement : Stats Notes Andy Lawrence, University of Edinburgh Autumn 2013 1 Probabilities, distributions, and errors Laplace once said Probability theory is nothing
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
The Variance of the Product of Random Variables Author(s): Leo A. Goodman Source: Journal of the American Statistical Association, Vol. 57, No. 297 (Mar., 1962), pp. 54-60 Published by: American Statistical
More informationHW7 Solutions. f(x) = 0 otherwise. 0 otherwise. The density function looks like this: = 20 if x [10, 90) if x [90, 100]
HW7 Solutions. 5 pts.) James Bond James Bond, my favorite hero, has again jumped off a plane. The plane is traveling from from base A to base B, distance km apart. Now suppose the plane takes off from
More information2. The CDF Technique. 1. Introduction. f X ( ).
Week 5: Distributions of Function of Random Variables. Introduction Suppose X,X 2,..., X n are n random variables. In this chapter, we develop techniques that may be used to find the distribution of functions
More information[313 ] A USE OF COMPLEX PROBABILITIES IN THE THEORY OF STOCHASTIC PROCESSES
[313 ] A USE OF COMPLEX PROBABILITIES IN THE THEORY OF STOCHASTIC PROCESSES BY D. R. COX Received 17 September 1954 ABSTRACT. The exponential distribution is very important in the theory of stochastic
More informationProblem 1 (20) Log-normal. f(x) Cauchy
ORF 245. Rigollet Date: 11/21/2008 Problem 1 (20) f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.0 0.2 0.4 0.6 0.8 4 2 0 2 4 Normal (with mean -1) 4 2 0 2 4 Negative-exponential x x f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.5
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
Selecting the Better Bernoulli Treatment Using a Matched Samples Design Author(s): Ajit C. Tamhane Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 42, No. 1 (1980), pp.
More informationSampling Distributions
In statistics, a random sample is a collection of independent and identically distributed (iid) random variables, and a sampling distribution is the distribution of a function of random sample. For example,
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
Fermat's Little Theorem: Proofs That Fermat Might Have Used Author(s): Bob Burn Source: The Mathematical Gazette, Vol. 86, No. 507 (Nov., 2002), pp. 415-422 Published by: The Mathematical Association Stable
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
6625 Author(s): Nicholas Strauss, Jeffrey Shallit, Don Zagier Source: The American Mathematical Monthly, Vol. 99, No. 1 (Jan., 1992), pp. 66-69 Published by: Mathematical Association of America Stable
More informationTable of z values and probabilities for the standard normal distribution. z is the first column plus the top row. Each cell shows P(X z).
Table of z values and probabilities for the standard normal distribution. z is the first column plus the top row. Each cell shows P(X z). For example P(X.04) =.8508. For z < 0 subtract the value from,
More information3 Continuous Random Variables
Jinguo Lian Math437 Notes January 15, 016 3 Continuous Random Variables Remember that discrete random variables can take only a countable number of possible values. On the other hand, a continuous random
More informationf X (x) = λe λx, , x 0, k 0, λ > 0 Γ (k) f X (u)f X (z u)du
11 COLLECTED PROBLEMS Do the following problems for coursework 1. Problems 11.4 and 11.5 constitute one exercise leading you through the basic ruin arguments. 2. Problems 11.1 through to 11.13 but excluding
More informationDiscrete Distributions
Discrete Distributions STA 281 Fall 2011 1 Introduction Previously we defined a random variable to be an experiment with numerical outcomes. Often different random variables are related in that they have
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 6
MATH 56A: STOCHASTIC PROCESSES CHAPTER 6 6. Renewal Mathematically, renewal refers to a continuous time stochastic process with states,, 2,. N t {,, 2, 3, } so that you only have jumps from x to x + and
More informationOPTIMUM DESIGN ON STEP-STRESS LIFE TESTING
Libraries Conference on Applied Statistics in Agriculture 1998-10th Annual Conference Proceedings OPTIMUM DESIGN ON STEP-STRESS LIFE TESTING C. Xiong Follow this and additional works at: http://newprairiepress.org/agstatconference
More informationThe Econometric Society is collaborating with JSTOR to digitize, preserve and extend access to Econometrica.
A Set of Independent Necessary and Sufficient Conditions for Simple Majority Decision Author(s): Kenneth O. May Source: Econometrica, Vol. 20, No. 4 (Oct., 1952), pp. 680-684 Published by: The Econometric
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
A Look at Some Data on the Old Faithful Geyser Author(s): A. Azzalini and A. W. Bowman Source: Journal of the Royal Statistical Society. Series C (Applied Statistics), Vol. 39, No. 3 (1990), pp. 357-365
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
Uncountably Many Inequivalent Analytic Actions of a Compact Group on Rn Author(s): R. S. Palais and R. W. Richardson, Jr. Source: Proceedings of the American Mathematical Society, Vol. 14, No. 3 (Jun.,
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
A Simple Non-Desarguesian Plane Geometry Author(s): Forest Ray Moulton Source: Transactions of the American Mathematical Society, Vol. 3, No. 2 (Apr., 1902), pp. 192-195 Published by: American Mathematical
More informationSurvival Analysis for Case-Cohort Studies
Survival Analysis for ase-ohort Studies Petr Klášterecký Dept. of Probability and Mathematical Statistics, Faculty of Mathematics and Physics, harles University, Prague, zech Republic e-mail: petr.klasterecky@matfyz.cz
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
On Runs of Residues Author(s): D. H. Lehmer and Emma Lehmer Source: Proceedings of the American Mathematical Society, Vol. 13, No. 1 (Feb., 1962), pp. 102-106 Published by: American Mathematical Society
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
More informationSTAT331. Cox s Proportional Hazards Model
STAT331 Cox s Proportional Hazards Model In this unit we introduce Cox s proportional hazards (Cox s PH) model, give a heuristic development of the partial likelihood function, and discuss adaptations
More information1 Delayed Renewal Processes: Exploiting Laplace Transforms
IEOR 6711: Stochastic Models I Professor Whitt, Tuesday, October 22, 213 Renewal Theory: Proof of Blackwell s theorem 1 Delayed Renewal Processes: Exploiting Laplace Transforms The proof of Blackwell s
More informationMixtures and Random Sums
Mixtures and Random Sums C. CHATFIELD and C. M. THEOBALD, University of Bath A brief review is given of the terminology used to describe two types of probability distribution, which are often described
More informationStatistical Inference: Estimation and Confidence Intervals Hypothesis Testing
Statistical Inference: Estimation and Confidence Intervals Hypothesis Testing 1 In most statistics problems, we assume that the data have been generated from some unknown probability distribution. We desire
More informationSTAT 6385 Survey of Nonparametric Statistics. Order Statistics, EDF and Censoring
STAT 6385 Survey of Nonparametric Statistics Order Statistics, EDF and Censoring Quantile Function A quantile (or a percentile) of a distribution is that value of X such that a specific percentage of the
More informationDiscrete Distributions
A simplest example of random experiment is a coin-tossing, formally called Bernoulli trial. It happens to be the case that many useful distributions are built upon this simplest form of experiment, whose
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationContinuous Distributions
A normal distribution and other density functions involving exponential forms play the most important role in probability and statistics. They are related in a certain way, as summarized in a diagram later
More informationUNIT 5:Random number generation And Variation Generation
UNIT 5:Random number generation And Variation Generation RANDOM-NUMBER GENERATION Random numbers are a necessary basic ingredient in the simulation of almost all discrete systems. Most computer languages
More informationExecutive Assessment. Executive Assessment Math Review. Section 1.0, Arithmetic, includes the following topics:
Executive Assessment Math Review Although the following provides a review of some of the mathematical concepts of arithmetic and algebra, it is not intended to be a textbook. You should use this chapter
More informationSIMPLIFIED CALCULATION OF PRINCIPAL COMPONENTS HAROLD HOTELLING
PSYCHOMETRIKA--VOL. 1, NO. 1 SIMPLIFIED CALCULATION OF PRINCIPAL COMPONENTS HAROLD HOTELLING The resolution of a set of n tests or other variates into components 7~, each of which accounts for the greatest
More informationShor Factorization Algorithm
qitd52 Shor Factorization Algorithm Robert B. Griffiths Version of 7 March 202 References: Mermin = N. D. Mermin, Quantum Computer Science (Cambridge University Press, 2007), Ch. 3 QCQI = M. A. Nielsen
More information5.2 Continuous random variables
5.2 Continuous random variables It is often convenient to think of a random variable as having a whole (continuous) interval for its set of possible values. The devices used to describe continuous probability
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
A Quincuncial Projection of the Sphere Author(s): C. S. Peirce Source: American Journal of Mathematics, Vol. 2, No. 4 (Dec., 1879), pp. 394-396 Published by: The Johns Hopkins University Press Stable URL:
More informationConvergence Rates for Renewal Sequences
Convergence Rates for Renewal Sequences M. C. Spruill School of Mathematics Georgia Institute of Technology Atlanta, Ga. USA January 2002 ABSTRACT The precise rate of geometric convergence of nonhomogeneous
More informationWEIBULL RENEWAL PROCESSES
Ann. Inst. Statist. Math. Vol. 46, No. 4, 64-648 (994) WEIBULL RENEWAL PROCESSES NIKOS YANNAROS Department of Statistics, University of Stockholm, S-06 9 Stockholm, Sweden (Received April 9, 993; revised
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationMind Association. Oxford University Press and Mind Association are collaborating with JSTOR to digitize, preserve and extend access to Mind.
Mind Association Response to Colyvan Author(s): Joseph Melia Source: Mind, New Series, Vol. 111, No. 441 (Jan., 2002), pp. 75-79 Published by: Oxford University Press on behalf of the Mind Association
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
Merging of Opinions with Increasing Information Author(s): David Blackwell and Lester Dubins Source: The Annals of Mathematical Statistics, Vol. 33, No. 3 (Sep., 1962), pp. 882-886 Published by: Institute
More informationNormal Random Variables and Probability
Normal Random Variables and Probability An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2015 Discrete vs. Continuous Random Variables Think about the probability of selecting
More informationChapter 2. Theory of Errors and Basic Adjustment Principles
Chapter 2 Theory of Errors and Basic Adjustment Principles 2.1. Introduction Measurement is an observation carried out to determine the values of quantities (distances, angles, directions, temperature
More informationQuadratic reciprocity (after Weil) 1. Standard set-up and Poisson summation
(December 19, 010 Quadratic reciprocity (after Weil Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ I show that over global fields k (characteristic not the quadratic norm residue symbol
More information1 Basic concepts from probability theory
Basic concepts from probability theory This chapter is devoted to some basic concepts from probability theory.. Random variable Random variables are denoted by capitals, X, Y, etc. The expected value or
More informationSampling Distributions
Sampling Distributions In statistics, a random sample is a collection of independent and identically distributed (iid) random variables, and a sampling distribution is the distribution of a function of
More informationJoint Probability Distributions and Random Samples (Devore Chapter Five)
Joint Probability Distributions and Random Samples (Devore Chapter Five) 1016-345-01: Probability and Statistics for Engineers Spring 2013 Contents 1 Joint Probability Distributions 2 1.1 Two Discrete
More informationεx 2 + x 1 = 0. (2) Suppose we try a regular perturbation expansion on it. Setting ε = 0 gives x 1 = 0,
4 Rescaling In this section we ll look at one of the reasons that our ε = 0 system might not have enough solutions, and introduce a tool that is fundamental to all perturbation systems. We ll start with
More informationChapter 6: Functions of Random Variables
Chapter 6: Functions of Random Variables We are often interested in a function of one or several random variables, U(Y 1,..., Y n ). We will study three methods for determining the distribution of a function
More informationEco517 Fall 2004 C. Sims MIDTERM EXAM
Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering
More informationThe Degree of the Splitting Field of a Random Polynomial over a Finite Field
The Degree of the Splitting Field of a Random Polynomial over a Finite Field John D. Dixon and Daniel Panario School of Mathematics and Statistics Carleton University, Ottawa, Canada fjdixon,danielg@math.carleton.ca
More informationa 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.
Chapter 1 LINEAR EQUATIONS 11 Introduction to linear equations A linear equation in n unknowns x 1, x,, x n is an equation of the form a 1 x 1 + a x + + a n x n = b, where a 1, a,, a n, b are given real
More informationAP Calculus BC Summer Assignment 2018
AP Calculus BC Summer Assignment 018 Name: When you come back to school, I will epect you to have attempted every problem. These skills are all different tools that we will pull out of our toolbo at different
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
A Note on Weighted Randomization Author(s): D. R. Cox Source: The Annals of Mathematical Statistics, Vol. 27, No. 4 (Dec., 1956), pp. 1144-1151 Published by: Institute of Mathematical Statistics Stable
More information1. If X has density. cx 3 e x ), 0 x < 0, otherwise. Find the value of c that makes f a probability density. f(x) =
1. If X has density f(x) = { cx 3 e x ), 0 x < 0, otherwise. Find the value of c that makes f a probability density. 2. Let X have density f(x) = { xe x, 0 < x < 0, otherwise. (a) Find P (X > 2). (b) Find
More informationPROBABILITY DISTRIBUTIONS: DISCRETE AND CONTINUOUS
PROBABILITY DISTRIBUTIONS: DISCRETE AND CONTINUOUS Univariate Probability Distributions. Let S be a sample space with a probability measure P defined over it, and let x be a real scalar-valued set function
More information04. Random Variables: Concepts
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 215 4. Random Variables: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationNotes on Random Variables, Expectations, Probability Densities, and Martingales
Eco 315.2 Spring 2006 C.Sims Notes on Random Variables, Expectations, Probability Densities, and Martingales Includes Exercise Due Tuesday, April 4. For many or most of you, parts of these notes will be
More informationTime: 1 hour 30 minutes
Paper Reference(s) 6684/0 Edexcel GCE Statistics S Silver Level S Time: hour 30 minutes Materials required for examination papers Mathematical Formulae (Green) Items included with question Nil Candidates
More informationOverview of Dispersion. Standard. Deviation
15.30 STATISTICS UNIT II: DISPERSION After reading this chapter, students will be able to understand: LEARNING OBJECTIVES To understand different measures of Dispersion i.e Range, Quartile Deviation, Mean
More informationSEQUENTIAL TESTS FOR COMPOSITE HYPOTHESES
[ 290 ] SEQUENTIAL TESTS FOR COMPOSITE HYPOTHESES BYD. R. COX Communicated by F. J. ANSCOMBE Beceived 14 August 1951 ABSTRACT. A method is given for obtaining sequential tests in the presence of nuisance
More information