Aspects of Statistical Computing: What Packages for the 1980's Ought to Do

A wish list of desirable statistical computing capabilities is presented. This may help one question which of these capabilities can be satisfied by existing packages, which might be met through reasonable extensions to these packages, which might require substantial new development, and which ought to be supplied by the computing environment rather than the packages. These questions are explored, taking into account the nature of the statistical work and the choices presented by technology. Attention is given to the barriers to be overcome if future statistical packages are to take full advantage of new technology.     

A monte carlo evaluation of response surface analysis based on paired-comparison data

The use of the logit transformation on paired-comparison data in the weighted least squares analysis of response surfaces for aesthetic qualities of products is discussed. Monte Carlo simulations are employed to investigate the small sample properties of the estimators and test statistics. A secondary objective of the Monte Carlo simulations is the comparison of two transformation procedures. The simulations are of standard-item paired-compar-ison experiments in which ties are not allowed.     

On concurrence of several linear regressions and applications to problems of enzyme kinetics

In comparing several regressions E(y_ij) =α_i + β_ix_ij i = 1, 2, ..., k, j = 1,2, ..., n_i, researchers are generally interested in the following five problems: whether they have (1) equal slope, (2) equal intercept, (3) coincidence, (4) common intersection on X-axis, and (5) common intersection on (X,Y) - plane. Problems (1) - (3) can be put into the framework of the general linear hypothesis and the F-test can be used. However, problems (4) and (5) cannot be put into the general linear hypothesis because they are ratios of parameters. Hence, in this paper we consider the generalized likelihood ratio test for hypothesis testing. An application to an enzyme kinetics problem in Aniline Metabolism is demonstrated     

Time series analysis of stratospheric ozone

Time series analyses of monthly average total ozone measured at 37 stations throughout the world were used to estimate the extent to which the average ozone trend correlates with the depletion curve hypothesized as due to chlorofluorocarbons (CFCs). Statistical characteristics of stations in the ensemble were used to help define appropriate model and station selection criteria. The maximum likelihood procedure developed herein estimates the weighted average trend, its. variance, and the intra- and inter-station variance components of the trend. Correlations among trends at different stations are also taken into account. The models were subjected to much checking and criticism. Variations in statistical methodology are used to show that the results are insensitive to details of the model selection criteria. The method does not discriminate well between the hypothesized CFC trend and a linear trend. The trend estimates represent the sum of all long-term global effects. The variance includes all effects that differ from station-to-station. The estimated trend and 2α limits for 14 stations with 20-year records (1958-79) is an ozone increase through 1979 of (1.5+1.0) percent. At the 23 stations with shorter records, the trend is (1.0=1.7) percent. It is concluded that no significant depletion in stratospheric ozone has occurred from any cause through the end of 1979.     

On the asymptotic distribution of cramer-von mises one-sample test statistics under an alternative

The asymptotic normality of the Cramer-von Mises one-sample test statistic and one of its variants under an alternative cdf is demonstrated. The derivation herein is unique in that it does not require knowledge of the theory of weak convergence of probability measures defined on metrized function spaces, and thus is accessible to a broader class of students and practitioners.     

Practical bayesian approaches to the analystcs of a 2×2 contingency table with incompletely categorized data

To assess the efficacy of a treatment, patients are administered a pre-test, the treatment, and a post-test (identical to the pre-test). These patients are then categorized according to their outcomes observed on both tests,e.g., (S,S), (S,F), etc. Also, we observe "incomplete" information on the pre-tests' outcomes for some patients and the results of only the post-test being known for thers, A Bayesian framework is fit to the problem and Bayes factors, posterior odds ratios, and utility functions are given to evaluate th e treatment, A method of assessing the prior distribution is specified and a numerical example is worked.     

Estimation of the parameters of a regression model with a multivariate t error variable

This paper deals with a regression model for several vari¬ables under the assumption that the errors have a multivariate t-distribution. The parameters of the model, the regression parameters, as well as the scale parameters and the degress of freedom of the error variable are estimated and the estimation procedure is illustrated by a numerical example, Also, the prop¬erties of the estimators and tests for the regression parameters are discussed.     

On f-tests and linear hypothesis in a general linear model

A simple method of setting linear hypotheses testable by F-tests in a general linear model when the covariance matrix has a general form and is completely unknown, is provided. With some additional conditions imposed on the covariance matrix, there exist the UMP invariant tests of certain linear hypotheses. We derive them to compare the powers with those of F-tests obtained under no restrictions on the covariance matrix. The results are illustrated in a multiple regression model with some examples.     

On the asymptotic distribution of an estimate of the change point in a failure rate

The failure rate r(t) is assumed to have the shape of the"first"part of the"bathtub"model, i.e.r(t) is non-increasing for t<r and is constant for t> r. Asymptotic distribution of one of the estimates proposed earlier has been investigated in this paper. This leads to a test for the hypothesis HQ r<r ₀ vs H :r>r (where TQ > 0). Asymptotic expression for the power of this test under Pitman alternatives is derived. Some simulations are reported.