Search

The article presents discussion of a paper by statistician I.J. Good. The author supposes that most statisticians approaching the problem of testing goodness of fit would have been very much less ambitious than Good has been. Most of the people, if spend some time looking at the first planet, Mercury and the last two, Neptune and Pluto, do not agree with Bode's law, whereas all the intermediate planets from Venus to Uranus seem to fit the law remarkably well. Those statisticians who were sufficiently knowledgeable would presumably meditate on the matters that Good has meditated on; but, having done this, most people would make a merely informal judgment and completely shirk the challenge that Good has met, to evaluate the weight of evidence. Good remarks that if one is going to test goodness of fit of a theory to some data, using either a Neyman-Pearson approach or a Bayesian approach, one have to frame explicit hypotheses. Over a good many years the author has, from time to time tried to think about the general problem of testing goodness of fit, and remain skeptical about whether goodness of fit can be adequately discussed either in Neyrnan-Pearson or in Bayesian terms because of this necessity to have, not only an explicit null hypothesis, but also a fully defined explicit alternative hypothesis or set of alternative hypotheses.

The article presents summaries of papers delivered at the 125th annual meeting of the American Statistical Association. The meeting took place on September 1965 in Philadelphia, Pennsylvania. In the paper "On the Construction of Stochastic Working Life Tables," by Haskel Benishay, stochastic working life tables are constructed and evaluated. The essence of the construction of stochastic rather than the traditionally deterministic working life tables is in the utilization of the realistic assumptions that the number of births, the number of entries into the labor force, the duration of life, and the duration of working life, are random variables. The end result is a working life table with a double entry in each cell rather than the single traditional one. The implications of this approach for the evaluation of fluctuations in various subcomponents of the labor force are evaluated. Another paper "Measuring Elasticities of Air Travel," by Samuel L. Brown, summarizes research on the response of demand for air travel to differences of fares, incomes, and elapsed travel time. There are three kinds of studies. Multiple regression studies of time series of air traffic, fares and income in the national domestic market yield short-period fare elasticities and income-elasticities. A quarterly first-difference model stands well the usual statistical tests.

The article presents abstracts and titles of abstracts submitted for papers presented at the 1969 Annual Meetings of the American Statistical Association and Biometric Society in New York. Some of the abstracts are Financial and Non-Financial Factors Affecting Post-High School Plans, 1939-1965," by Walter Adams, "Experimental Designs and Mathematical Models," by Sidney Addelman, "Consistency of Linear Least Squares Estimate in a Regression Model with Lagged Variable," by M. Ashanullah, "The Use of Prior Information in Testing of New Drugs," by David M. Allen, "A Comparison of Probability Density Estimates," by Gary Anderson, "Stochastic Models of City Population Density," by Robert T. Amsden, "ldentification and Input Signal Synthesis Problems in Discrete-Time Stochastic Control Systems," by Masanao Akoi, "A Computer Language for the Analysis of Variance," by David J. Aemor, "A Generalized Mathematical Model for Biological Recovery with Applications in Radiation Biology," by George M. Angelton."

The article presents an alphabetical list of abstracts and titles of abstracts submitted for papers presented at meetings of the American Statistical Association and Biometric Society in Washington D.C., in December 1967. The abstract "Indians and Smudges on the Census Schedule," by Donald S. Akers and Elizabeth A. Larmon, is an explanation of how certain anomalies in the reported characteristics of American Indians in the 1960 Census are to be explained by systematic smudging of the census schedules. The abstract "Input-Output and the Economy," by Clopper Almon explains that input-output and regression analysis were once considered rival tools for analysis of the economy. In those days, the constancy of the coefficients seemed a crucial issue in deciding whether or not to use input-output. Today, thanks to the progress in computer hardware and model building technique, regression analysis and input-output complement one another and the constancy of coefficients need be assumed only where we know of no better assumption.

The article presents the author's comments on an article "Optimum Sample Size for a Problem in Choosing the Population With the Largest Mean," by statistician Paul N. Somerville, published in the June 1970 issue of "Journal of the American Statistical Association." In the article Somerville described a procedure for selecting the population with the largest mean from normal populations with unknown means and a common known variance. The author claims that his procedure gives a unique minimax solution and that the minimax is a reasonable solution to the problem. According to the author, Somerville has not made it clear that his procedure gives only the local minimax solution. This makes practical applications of his procedure severely limited for it is difficult to think of a situation in which an experimenter would like to determine a local minimax value.

The lack of memory property of the exponential distribution plays an important part in the branch of applied probability. This property assumes the information regarding the probability distribution. In the present paper we give a characteristic property of the exponential distribution based on the means of the conditional distributions. Considering a random variable T with finite mean and such that P(T > O) > 0, and denoting by y a positive number such that P(T > y) > 0, we show that T has an exponential distribution if and only if the mean of the conditional distribution, given T > y, exceeds the mean of the unconditional distribution by the quantity y for all y. Also given in the paper is a similar characterization for the geometric distribution. Since the information concerning the expected values is easily accessible, we expect these properties to be useful in dealing with the practical problems. [ABSTRACT FROM AUTHOR]

The model considered in this paper is simple linear regression (Ey[sub i] = beta[sub 1] + beta[sub 2] x[sub 1], i = 1, ..., n), and the problem is to make statistical inferences about an unknown value of x corresponding to one or more additional observed values of y. The maximum likelihood estimator x of x and the classical (l - alpha) 100% confidence set S for x have some undesirable properties. For example, x has infinite mean square error and P {S = (- infinity, + infinity)} > 0. The purpose of this paper is to demonstrate that insight and understanding, as well as a useful class of solutions, can be obtained by looking at the problem from a Bayesian point of view. A result which follows from a general Bayes solution is that the inverse estimator [4] is Bayes with respect to a particular informative prior. [ABSTRACT FROM AUTHOR]

The article discusses moments of the distribution of sample size, N in a sequential probability ratio test (SPRT). The present paper provides variance, the third and the fourth moments of N. The details are worked out in five common applications of the SPRIT. The relation of the variance of N to the truncation of a SPRT is discussed is also discussed in the paper. Scholar A. Wald indicated in passing how one can obtain the moments of N, but the only published literature where the author encountered a general expression for the variance of N. However, their expression is incorrect. Using scholar J. Wolfowitz's results, which they do, or differentiating Wald's, fundamental identity twice one gets provided. In many practical applications of the SPRT, μ and moments in an equation derived are differentiable functions of a real-valued parameter. The limiting expressions for the moments can then be determined by standard methods of mathematical analysis. However, for the third and fourth moments the actual technique may involve an excessive amount of arithmetic.

In an earlier paper [5] the authors compared the power of the BLUS test with the probability of a correct decision of the Durbin-Watson bounds test. A method to compute the distribution of the Von Neumann ratio under the null hypothesis and under the alternative hypothesis was given. In the present paper the latter method is used to tabulate the BLUS-test statistic and to compute the exact significance points of the Durbin-Watson test for several examples. Powers of both tests are computed and compared. It appears that, for the cases considered, the power of the exact Durbin-Watson test exceeds that of the BLUS procedure, while the latter is greater than the probability of a correct decision in the Durbin-Watson bounds test. [ABSTRACT FROM AUTHOR]

In a previous paper [4] an estimate of the 4-variate product moment E[{x[sub I] - E(x[sub I])} {x[sub j] - E(x[sub j])} {x[sub k] - E(x[sub k])} {x[sub h] - E(x[sub h])}] was obtained. This estimate had a slight bias, as we pointed out. In this paper an unbiased estimate of the 4-variate product moment is obtained. Asymptotic results for the 3-variate and 4-variate product moments and their estimates are also obtained. [ABSTRACT FROM AUTHOR]

The purpose of this paper is to present an easily applied test useful in determining whether or not a supplier's data have been truncated. The proposed test has the following desirable properties: (I) it is the uniformly most powerful rank test, (ii) it is asymptotically uniformly most powerful, (iii) power computations can easily be made for arbitrary sample sizes, formulas for such computations being given in the paper. Although formulated in the context of verifying a supplier's data, the test can be applied to other situations in which false representation of data in the form of truncation is important. Such is the case, for example, in reliability demonstrations or legal suits involving physical measurements. [ABSTRACT FROM AUTHOR]

In this paper, the simple estimate formed by the ratio of the sample means, the estimate obtained by the statistician M.H. Quenouille's method of bias reduction, the unbiased estimate of the statisticians L.A. Goodman and H.O. Hartley and the estimate proposed by H.O. Hartley are compared. These estimates are respectively denoted by t1, t2, t3 and t4. Mean square errors of these estimates are compared under J. Durbin's model. In practice g has often been found to lie between 0 and 2. In this paper expressions for the Mean Square Errors of these estimates are obtained for general values of n, g and h, and they are compared for finite values of n when g = 0, 1 and 2. When g = 0. The method of obtaining exact expressions for the Mean Square Errors (MSE's) of these estimates is similar to that of J.N.K. Rao and J.T. Webster. After a considerable amount of simplification, the following expressions for the biases and MSE's are obtained. When t=0, Durbin compared t1 with t2, and Rao compared t1 with t4.

This paper is written with three objectives in mind. First, to point out that the problem of grouping, where a larger number of elements n are combined into m mutually exclusive groups (m ~n) should be recognized as a problem in Integer Programming, and that such recognition can help us in avoiding complete enumeration of stages in grouping from n to n - 1 ... to m, and of alternative possibilities in each stage. Second, to formulate mathematically some simple versions of the relevant Integer Programming Problem, so that the available computer codes can solve it. When the grouping attempts to minimize the within groups sums of squares, the so-called string property is proved to be necessary for the minimum (Lemma 1). It is shown that the string property can be exploited to write the non-linear within group sums of squares as a linear function (Lemma 2). An attempt is made to generalize the string property for the higher dimensional case. Third, to give some numerical examples to clarify the mathematical models and to illustrate the advantages of Integer Programming formulation by comparing with the numerical examples appearing in the literature [8], [13 ]. Zero-one programming where the integer variables can take values of zero or one only appears to be most efficient from the limited experience of the author. The paper is divided into three sections which closely correspond to the three objectives noted above. [ABSTRACT FROM AUTHOR]

The paper is concerned with testing an application of a theory of investment behavior based on the neo-classical theory of optimal capital accumulation. The demand for capital is determined to maximize net worth under the neo-classical conditions of production. With the assumption of a lag in the completion of investment projects, gross investment is a distributed lag function of period-to-period changes in desired capital stock. The investment process is divided into a series of separate stages, the first stage being a change in the demand for capital services, intermediate stages being for letting contracts, and so on, and the final stage the actual investment expenditure. To extend the theory to include anticipated investment expenditures, these are treated as a special kind of intermediate stage. The purposes of the paper are to test hypotheses about the internal structure of the investment process when viewed in this way, and to investigate whether forecasting actual investment expenditures by treating the investment process in this way can improve upon the method used by the OBE-SEC. [ABSTRACT FROM AUTHOR]

This paper concerns a new Normal approximation to the beta distribution and its relatives, in particular, the binomial, Pascal, negative binomial, F, t, Poisson, gamma, and chi square distributions. The approximate Normal deviates are expressible in terms of algebraic functions and logarithms, but for desk calculation it is preferable in most cases to use an equivalent expression in terms of a function specially tabulated here. Graphs of the error arc provided. They show that the approximation is good even in the extreme tails except for beta distributions which are J or U shaped or nearly so, and they permit correction to obtain still more accuracy. For use beyond the range of the graphs, some standard recursive relations and some classical continued fractions are listed, with some properties of the latter which seem to be partly new. Various Normal approximations are compared, with further graphs. The new approximation appears far more accurate than the others. Everything an ordinary user of the approximation might want to know is included in this paper. The theory behind the approximation and most proofs are postponed to a second paper immediately following this one. [ABSTRACT FROM AUTHOR]

Several writers have studied the problem of obtaining confidence bands for a straight line regression under the requirement that the bands be at a specified confidence level, (1 - α) say, for values of the independent variable restricted to a pre-specified closed interval. The bands studied by other writers have been either bands of equal width throughout the interval or trapezoidal bands. This paper studies confidence bands of the classical hyperbolic type under the restriction noted above. Relevant distribution results differ according to whether the pre-specified interval on the independent variable is symmetric or asymmetric about the mean of the independent variable values in the experiment. In the former case the problem under study, is shown to be equivalent to a problem studied by Halperin et al (JASA, Sept. 1967) for which distribution theory and tables are already available. The symmetrical case is generalized to obtain confidence bands in multiple linear regression at level (1 - α) for an ellipsoidal region on the independent variables centered at the point of means of the independent variable values used in the experiment. For straight line regression some numerical comparisons are made with bands of the other types mentioned above. These comparisons suggest that the bands proposed in this paper are uniformly superior to bands of equal width in the sense of having smaller average width; the proposed bands appear to be superior to trapezoidal bands in the same sense for intervals of practical interest, i.e. within or reasonably close to the range on the independent variables used in the experiment. [ABSTRACT FROM AUTHOR]

This paper presents spectral criteria to evaluate seasonal adjustment procedures, and applies them to the estimated spectral properties of the monthly unemployment series, as adjusted by past and present methods of the Bureau of Labor Statistics (BLS) and of the Census Bureau. It also re-examines some conclusions reported in the Gordon Committee Report (1962) and in Nerlove (1964) on past methods of adjustment. The paper concludes that (1) the loss in spectral power between the unadjusted and adjusted series spectra over most of the frequency range could arise with a satisfactory adjustment, (2) the excessive loss of spectral power at seasonal and trend-cycle frequencies present in past methods of adjustment has been reduced in present methods, (3) the effect of deviations from desired spectral properties on the uses of seasonally adjusted data must be examined in the time domain; for certain applications they have been found to be unimportant. [ABSTRACT FROM AUTHOR]

In this paper a method is developed for estimating the parameters in the multivariate normal distribution in which the missing observations are not restricted to follow certain patterns as in most previous papers. The large sample properties of the estimators are discussed. Equivalence with maximum likelihood estimators has been established for a subclass of problems. The results of some simulation studies are provided to support the theoretical development. [ABSTRACT FROM AUTHOR]

The 1967 Committee on Publications, chaired by David L. Wallace, found that many American Statistical Association members desired more review and survey papers. These have been hard to come by and so as the author's last act before leaving office, decided to provide a short survey paper on some related ideas in a field where nearly all of the statisticians sometimes work — that of contingency tables. These ideas are largely available in literature and yet they have not often been put together, though statisticians I. J. Good's monograph and Leo Goodman's many papers form good sources. But the author's paper is not intended as a review of the literature, only as a survey of one set of ideas about estimation in the analysis of contingency tables. The author fears that the first act of most social scientists upon seeing a contingency table is to compute chi-square for it. Sometimes this process is enlightening, sometimes wasteful, but sometimes it does not go quite far enough. The author collected 500 samples of writing published about 1961.

This paper describes a method of estimating the number of marriages during an inter-censal period among a group of never-married persons. Two alternative procedures are suggested: the first based on data on the never-married category and the second based on data on the ever-married category. The essential data needed are marital distributions by quinquennial age-groups from two successive censuses. The marital data from the 1951 and 1961 censuses of India are used to discuss the methodological problems involved in the estimation of marriage frequencies and in the construction of nuptiality tables based on them. Marriage and death are taken as the two attrition factors for decrement of single population. Gross and Net Attrition Tables for the single population of India are constructed for each sex. Section 5 describes various nuptiality measures derived from the Attrition Tables. It is suggested that the approach developed in this paper can be utilized for reliable marital and, subsequently, fertility projections. [ABSTRACT FROM AUTHOR]

Experimental design considerations have been virtually ignored by economists who have conducted computer simulation experiments with models of economic systems. The objective of this paper is to spell out in detail the relationship between existing experimental design techniques and techniques of data analysis and the design of simulation experiments with economic systems. We begin by defining the problem of experimental design as applied to computer simulation experiments. With the aid of an example model, we explore several techniques of data analysis and a number of specific experimental design problems. Although this paper is oriented towards the design of computer simulation experiments in economics, the techniques which are discussed are of a general nature and should be applicable to the design of simulation experiments in other disciplines. [ABSTRACT FROM AUTHOR]

A growing interest has been witnessed in recent, years in multivariate discrete probability models. The multivariate logarithmic series distribution (LSD) is a multivariate analogue of the univariate LSD. This paper investigates the marginal and the conditional distributions of the multivariate LSD. It. records its moment properties and also provides some regression and correlation analysis. An interesting modal property of the distribution is discovered. The paper also discusses the estimation of the parameters by the methods of maximum likelihood and unbiased minimum variance. At the end is discussed in detail an application of the multivariate LSD to the field of population and community ecology. [ABSTRACT FROM AUTHOR]

This paper considers the problem of obtaining efficient estimates for the parameters of a system of M regression equations. The disturbance terms of this system are assumed to be related by both serial and contemporaneous correlation. Under the further assumption that the serial correlation is a first order autoregressive process, the paper develops an estimator that is consistent and has the same asymptotic normal distribution as the Aitken estimator which assumes the covariance matrix to be known. The paper concludes with a discussion of some alternative covariance specifications and points out certain difficulties with the standard single equation procedures for handling autoregressive schemes. [ABSTRACT FROM AUTHOR]

In this paper confidence bands are given for the entire line for simple linear regressions. The conventional bands that have been given in the past are curvilinear. In this paper we consider confidence bands that are straight lines. It is shown that under certain conditions the "width" of these straight line bands is less than the width of the conventional curvilinear bands. [ABSTRACT FROM AUTHOR]

The Bradley-Terry model for a paired-comparison experiment with t treatments postulates a set of t 'true' treatment ratings π1, π2, ..., π1 such that π1 is ≥ 0, σpi;i = 1 and the probability for preferring treatment i to treatment j is π2 (π2 + πj)-1. Thus, according to this model, every comparison of two treatments results in a definite preference for one of the two. This is an unrealistic restriction since when there is no difference between the responses due to two treatments, any method of expressing preference for one over the other is somewhat arbitrary. This paper considers a modification of the Bradley-Terry model by introducing an additional parameter, called threshold parameter, into the model. This permits 'ties' in the model. The problem of estimation and tests of hypotheses for the parameters of the modified model is also dealt with in the paper. [ABSTRACT FROM AUTHOR]

An earlier paper [3] provides axioms for the seasonal adjustment of economic time series. Seasonally adjusted data should be minimum variance, unbiased, and linear in an appropriate sense. These axioms yield a unique method for seasonal adjustment. The seasonally adjusted data may be obtained by an application of ordinary least squares regression. The problem of seasonal adjustment of economic time series in [3] is not the same as the problem of seasonal adjustment of data for econometric analysis. The first problem is completely resolved by appeal to the axioms of minimum variance, unbiasedness, and linearity. The second problem requires formulation as a standard problem in econometrics: estimation of the parameters of a single equation in a system of simultaneous equations. Lovell [4] has proposed to solve the problem of seasonal adjustment of data for econometric analysis by applying ordinary least squares directly to a structural equation in a system of simultaneous equations. Ordinary least squares estimation of the parameters of a structural equation usually results in "least squares bias." However, under certain special assumptions such a procedure can be justified, as Wold [6] has pointed out. Under these assumptions Lovell's procedure is valid. In this paper the seasonal adjustment of data for econometric analysis is formulated as a general simultaneous equations problem. Conditions for identification of the parameters to be estimated and methods for constructing consistent and asymptotically unbiased and efficient estimates are derived. Extensions to multivariate problems of seasonal adjustment for econometric analysis are sketched. [ABSTRACT FROM AUTHOR]

In a recent paper in this Journal, Rothenberg, Fisher and Tilanus [1] discuss a class of estimators of tile location parameter of the Cauchy distribution, taking the form of the arithmetic average of a central subset, of the sample order statistics. They show that the average of roughly the middle quarter of the ordered sample has minimum asymptotic variance within this class, and that asymptotically it eliminates about. 36 per cent of the efficiency loss of the median (the most commonly used estimate,r) in comparison to the maximum likelihood estimator (m.l.e.). Of course both the m.l.e, and the best linear unbiased estimator based on the order statistics (BLUE) achieve full asymptotic efficiency in the Cramer-Rao sense and there can be no dispute about the relative merits of the three estimators asymptotically, or about the inferiority of the median (with asymptotic efficiency 8/pi[sup 2] * This character cannot be converted in ASCII text) 0.8 compared with about 0.88 for the estimator of Rothenberg et al.). In any practical situation however, we will be concerned with estimation from samples of finite size and asymptotic properties will not necessarily give any guidance here. We are essentially concerned with two points in assessing the relative merits of estimators in small samples, their case of application and "small-sample efficiency" which is conveniently measured as the ratio of the Cramer-Rao lower bound to the variance of the estimator. In this paper various estimators of the location of the Cauchy distribution arc compared in these two respects for samples of up to 20 observations. The small-sample properties of the m.l.e. have been extensively discussed elsewhere (Barnett [2]) and relevant results are summarized where necessary. The main purpose of the paper is to discuss general linear estimators based on the order statistics, and to assess their utility in the present context. Since this paper was prepared a further interesting 'quick estimator', b [ABSTRACT FROM AUTHOR]

Given N points independently drawn from the uniform distribution on (0, 1), let p[sub n], be the size of the smallest interval that contains n out of the N points; let n[sub p], be the largest number of points to be found in any subinterval of (0, 1) of length p. This paper uses a result of Karlin, McGregor, Barton and Mallows to determine the distribution of n[sub p] for p = 1/k, k an integer. The paper gives simple determinations for the expectations and variances of p[sub n], for all fixed n > (N + 1)/2, and of n[sub 1/2]. The distribution and expectation of n[sub p] are estimated and tabulated for the cases p = 0.1(0.1)0.9, N =2(1)10. [ABSTRACT FROM AUTHOR]

This paper develops a stochastic model of accounts receivable, the debt arising from credit sales; and expresses the expectations and variances of (1) credit sales, (2) collections, (3) bad debts, (4) collections plus bad debts; and, (5) accounts receivable outstanding, in terms of parameters of the stochastic variables: (a) the number of credit sales made, (b) the size of sales; and, (c) the collection period or the effective term of loans. The paper discusses some interesting features of the results of the model and provides an empirical example to illustrate the model in context of actual credit sales made by a small wholesale food company during the years 1963 and 1964. [ABSTRACT FROM AUTHOR]

This paper discusses a procedure for finding robust; estimators of the location parameter of symmetric unimodal distributions. The estimators are based on robust rank tests and the methods used are applicable to other one parameter problems. To every density function there corresponds an asymptotically most powerful rank test (a.m.p.r.t.). For a set F of density functions the maximin rank test, R, maximizes the minimum limiting Pitman's efficiency of R relative to the a.m.p.r.t. for each member of F. This maximin test, R, can be used to construct estimators according to the proposal of Hodges and Lehman; it generates another estimator T in the following manner. If the test based on R is the a.m.p.r.t. for samples from a density function g(x - theta), then the estimator T will be the best linear unbiased estimate (b.l.u.e.) of the location parameter for samples from g(x). Unfortunately, the estimator T is not necessarily consistent for all members of F. A class of rank tests which generate linear combinations of a few order statistics is introduced and a simple estimator using the 33 1/3rd, 50th and 66 2/3rd percentiles is proposed. The relationship of the present paper to the work of Huber is discussed and it is shown that the b.l.u.e. corresponding to his least favorable distribution is the trimmed mean. [ABSTRACT FROM AUTHOR]

Recently Professors T. J. Rothenberg, F. M. Fisher, and C. B. Tilanus published a paper proposing the class of trimmed means as estimators of the location parameter of the Cauchy distribution [5]. They showed that the asymptotic sampling variance of the estimators in this class is essentially minimized by using the middle 24% of the sample order statistics. The corresponding estimate has an asymptotic relative efficiency to the best estimator for complete samples (A.R.E.) of .87796 as compared to an A.R.E. of .81057 for the sample median. In this paper a few "quick estimators" are considered as estimators for the location parameter of the Cauchy Distribution. A "quick estimate" is a linear combination (a weighted average) of one or more order statistics. Our goal is to find a simple estimator, i.e. an estimator based on only a few order statistics, which has an A.R.E. of at least 90%. We found an estimator based on five order statistics which is considerably better than the optimum trimmed mean (using the middle 24% of the sample order statistics) and much better than the sample median. The A.R.E. of the optimum censored estimate with censored fractiles .38 and .62 is also found, and a comparison between the trimmed, censored, and proposed estimators is made. [ABSTRACT FROM AUTHOR]

Two recent papers have pointed out that living costs are not only affected by the prices and qualities of products but also by the number of forms in which each product is offered.¹ Moreover, a tentative procedure has been outlined for measuring the economic significance of changes in product variety and for adjusting cost-of-living indexes accordingly. The present article is intended to explore the points covered in the earlier papers more fully and to answer several questions that have been raised concerning them.² [ABSTRACT FROM AUTHOR]

Given a population of N units, it is required to draw a sample of n distinct units in such a way that the probability for the ith unit to be in the sample is proportional to its 'size' x. From the alternative methods of achieving this we consider here only the so-called systematic method which, to the best of our knowledge, was first developed by W. G. Madow (1949): The units in the population are listed in a 'particular' order, their x, accumulated and a systematic selection of n elements from a 'random start' is then made on the accumulation. In a more recent paper (H. O. Hartley and J. N. K. Rao (1962) ) an asymptotic estimation theory (for large N) associated with this procedure was developed for the case when the order of the listed units is random. In this paper we draw attention to certain properties of Madow's estimator: We utilize the fact that with systematic sampling the total number of different samples is N (rather than ([This eq. cannot be change in char.]) as with completely random sampling). This simplification in the definition of the variance of the estimator in repeated sampling enables us to identify the exact variance of Madow's estimator with a 'between sample mean square' in a special analysis of variance (see section 4) and compare it with the variance of the pps estimator in sampling with replacement as well as in other sampling procedures. We also develop two approximate methods of variance estimation (see section 5). We pay particular attention to the case when the units are listed in the order of their size. With this particular arrangement our method can be described as 'systematic with random start' and the gain in precision that we accomplish has of course, analogues in systematic sampling with equal probabilities employing ratio estimators in which there is a relation between the ratio ri =yi/Xi and xi Compared with other methods the present procedure combines the advantage of ease of systematic sample selection with the availability of exact variance formulas for any n and N. Moreover, it usually leads to a more efficient estimate. Its shortcoming resides in the fact that the estimation of the variance is based on certain assumptions. [ABSTRACT FROM AUTHOR]

In educational and psychological applications as well as in other applications, it may be necessary to make certain comparisons of two regression lines when the variances are unequal. Such problems arise, for example, in studies comparing two alternative curriculums or two different teaching methods. By generalizing an idea which Scheffe used to obtain a test for the Behrens-Fisher problem, this paper develops some tests for comparing two regression lines when the two sets of error terms are normally distributed but with two different variances. Scheffe's test itself is a randomized test, but in this paper we present both randomized and non-randomized tests. Both simple and multiple regression are considered, but the simple regression tests are computationally easier than the multiple regression tests. The basic test statistic which is used is the ordinary t-statistic. Essentially two types of problems are dealt with: (A) determining whether the two regression lines are identical when they are known to be parallel; and (B) determining whether the two regression lines are parallel. Confidence bounds as well as tests of hypotheses are available. For Problem A, a minimax estimator of the distance between the two lines is obtained. In addition to the Scheffe-type tests, we also consider some tests based on an approach of Welch and Hajek. A numerical example is presented. [ABSTRACT FROM AUTHOR]

The article presents information about Samuel S. Wilks, a great contributor to statistics. Wilks' behavior toward applications was peculiarly split, he encouraged his students to work on applications, not always an easy thing to do in a strongly theoretical mathematics department; he often told students about applications that he regarded as "neat" or "cute" or "clever"; the statistical colloquium he guided often had speakers on practical applications of mathematical statistics; he himself wrote some practical papers in statistics, but in the classroom he rarely discussed applications. Repeatedly, however, practical problems explicitly influenced both his own publications and those of his students and Wilks often used applications as motivation for the discussion of distribution theory, usually going well beyond the needs of the original problem. The papers to be discussed exhibit incompletely and fragmentarily a major influence on the work of the man and his students. Wilks watched the development of the statistical theory of order statistics closely. Indeed he wrote a masterful summary of the literature. The general area involves the study of the statistical properties of ordered measurements.

The paper is concerned with randomized experiments with one treatment. Two randomization schemes are considered: randomized pairs and unrestricted randomization. If effective at all, the treatment is supposed to affect the conditional distribution of the "experimental" variable Y given another variable X, called "predictor". The distribution of X is not affected by the treatment. Using the general theory published elsewhere, the paper deduces the locally asymptotically optimal test of the hypothesis that the treatment has no effect. Apart from the usual difficulties connected with asymptotic tests (how large must N be?), the theory is easily applicable in many "live" cases even though the conditional distribution of Y given X may contain nuisance parameters and be of unusual form. [ABSTRACT FROM AUTHOR]

Various methods based on Student t variates have been suggested and used for obtaining simultaneous confidence intervals for several means, or for several contrasts among means. Determination of an overall confidence level for such intervals involves evaluating the probability mass of a multivariate t distribution over a hypercube centered at the origin, with sides paralleling the coordinate planes, or obtaining bounds for this probability mass. Since such distributions involve many nuisance parameters, an impossible number of tables would be necessary in order to make exact confidence intervals. In the virtual absence of tables, approximations and bounds become important. In this paper, an attempt has been made to investigate the adequacy of certain suggested approximations [2], [5], [8] by computing the exact distributions for some particular cases. These exact distributions have been compared with approximations. This paper is concerned with two-sided confidence intervals, rather than one-sided intervals. [ABSTRACT FROM AUTHOR]

This paper deals with the asymptotic relative efficiency (A.R.E.) and the asymptotic power of certain tests. These tests are distribution-free goodness-of-fit tests based on the empirical distribution function. They are the one and two sided Kolmogorov-Smirnov tests and a new one based on the absolute values of the order statistics. The method of obtaining the A.R.E. is due to Quade [5] and is outlined in Section 1. In this paper we are concerned with the case of scalar alternatives. The asymptotic power for the Kolmogorov-Smirnov tests is obtained in Section 2. Section 3 deals with a new test which allows a more facile use of a one-sided statistic while utilizing the same information normally included in a two-sided statistic. In Section 4, three examples are considered to illustrate the superiority of this new test. [ABSTRACT FROM AUTHOR]

In a previous paper [12], an attempt was made to show how spectral techniques could be used to compare the effects of two seasonal adjustment procedures on the series to which they were applied. The two procedures compared were: (a) the technique currently used by the Bureau of Labor Statistics for seasonally adjusting employment, unemployment, and labor force monthly statistics, and (b) the so-called "residual" method, proposed by Brittain [2], Samuelson [16], and others. Spectra of the original and the seasonally adjusted series and the cross spectrum of the two were used to aid in the assessment of whether either procedure removed more than could be considered seasonal, introduced spurious regularities, and/or distorted temporal relationships. It was concluded that both techniques removed more than seasonal effects from, and produced some temporal distortion in, the series to which they were applied. Neither method appeared to be superior to the other. It is the purpose of this paper to carry the previous analysis one step further and to compare the BLS procedure with a modified version of the regression method of seasonal adjustment suggested by Cowden [3] and Mendershausen [11], and recently revived by Harman [9, 10] in an exceptionally sophisticated form. In addition to Hannan's work along the lines suggested, Nettheim [13] and Rosenblatt [15] have made studies. Rosenblatt [15] has carried out analyses similar to those reported here. [ABSTRACT FROM AUTHOR]

In this article, the author will talk about researcher R.A. Fisher's contributions to statistics. It is well-known that his work in genetics was of comparable status. It is largely represented by his book "The Genetical Theory of Natural Selection," though in his subsequent work his further association with ecological and experimental studies in evolutionary genetics, and his share in the development of studies in the human blood groups, might especially be recalled. Fisher's contributions to statistics began with a paper in 1912 advocating the method of maximum likelihood for fitting frequency curves, although the first paper of substance was his 1915 paper in the journal "Biometrika," on the sampling distribution of the correlation coefficient. Fisher tried to pose problems of analysis as the reduction and simplification of statistical data. He put forward his well-known concept of amount of information in estimation theory, such that information might be lost, but never gained, by analysis. His concept has been of great practical value, especially in large sample theory.

The principal part of this paper is devoted to the study of the distribution and density functions of the ratio of two normal random variables. It gives several representations of the distribution function in terms of the vivariate normal distribution and Nicholson's V function, both of which have been extensively studied, and for which tables and computational procedures are readily available. One of these representations leads to an easy derivation of the density function in terms of the Cauchy density and the normal density and integral. A number of graphs of the possible shapes of the density are given, together with an indication of when the density is unimodal or bimodal. The last part of the paper discusses the distribution of the ratio (u[sub 1] + ... + u[sub n])/(v[sub I] + ... + v[sub m]) where the u's and v's are independent, uniform variables. The exact distribution for all n and m is given, and some approximations discussed. [ABSTRACT FROM AUTHOR]

The use of simulation, as a technique for attacking difficult problems, has increased greatly with the availability of the digital computer. This is illustrated by the large number of references in Shubik's (1960) bibliography[sup 2] and in the large number of studies published since then. Simulation is essentially an extension of a technique known as empirical sampling, or distribution sampling, which has been used in the field of statistics for many years. The limitations of the technique, which are well known to statisticians, are apparently not as well known, or at least not as well recognized, by those using simulation today. The first part of this paper contains an historical survey of distribution sampling as used by statisticians. The material was originally prepared in 1953 and is reproduced here in slightly revised form to bring this history to the attention of present day simulators in order that the lessons that can be learned from this part can more readily be incorporated in the development of methodology today. The second part of this paper discusses the relevance of empirical sampling to the present day state of the art of simulation. The technique of generating random members, developed for empirical sampling can be applied directly to simulation. However in other aspects simulation is more difficult than empirical sampling and here the theory of distribution sampling does not have much to offer. The difficulties are due to lack of independence among time series, non-stationarity of the time series, and the large number of parameters. [ABSTRACT FROM AUTHOR]

In 1965 Stanley Warner [8] illustrated a technique whereby one could estimate from a sample the proportion of persons in a population possessing some characteristic X, without pointedly asking the question, "Do you possess characteristic X?". The present article uses a variation of these ideas suggested by Warner in a later paper [7] to estimate the distribution function of a continuous-type random variable. The technique is illustrated by estimating an income distribution from the responses of a sample of 500 individuals. The potential use of devices of this type in maintaining confidentiality of existing data files is apparent. [ABSTRACT FROM AUTHOR]

Given N points distributed at random on [0,1), let np, be the size of the largest number of points clustered within an interval of length p. Previous work finds Pr (np ≥ n), for n> N/2, and for n < N/2, p=1/L, L an integer. The formula for the case p=1/L is in terms of the sum of L × L determinants and is not computationally feasible for large L. The present paper derives such a computational formula. [ABSTRACT FROM AUTHOR]