A note on the ratio of normal and Laplace random variables

The normal and Laplace are the two earliest known continuous distributions in statistics and the two most popular models for analyzing symmetric data. In this note, the exact distribution of the ratio | X / Y | is derived when X and Y are respectively normal and Laplace random variables distributed independently of each other. A MAPLE program is provided for computing the associated percentage points. An application of the derived distribution is provided to a discriminant problem.     

Determination of sample size for selecting the smallest of k poisson population means

A procedure for selecting a Poisson population with smallest mean is considered using an indifference zone approach. The objective is to determine the smallest sample size n required from k ≥ 2 populations in order to attain the desired probability of correct selection. Since the means procedure is not consistent with respect to the difference or ratio alone, two distance measures are used simultaneously to overcome the difficulty in obtaining the smallest probability of correct selection that is greater than some specified limit. The constants required to determine n are computed and tabulated. The asymptotic results are derived using a normal approximation. A comparison with the exact results indicates that the proposed approximation works well. Only in the extreme cases small increases in n are observed. An example of industrial accident data is used to illustrate this procedure.     

A modification of the partial sequential probability ratio test

For the problem of discriminating between two simple hypoth¬eses concerning a Koopman - Darmois parameter, a modification of the partial sequential probability ratio test is proposed where instead of drawing only one fixed sample, two fixed samples are drawn and then Wald's SPRT is started. The OC and the ASN func¬tions are derived. Numerical comparisons are made with Wald's and Read's procedures for testing the normal mean with known variance. For some parameter values, the test procedure has a lower ASN than that of Read's procedure.     

On the bootstrap quantile-treatment-effect test

Let {X ₁, …, X _n } and {Y ₁, …, Y _m } be two samples of independent and identically distributed observations with common continuous cumulative distribution functions F(x)=P(X≤x) and G(y)=P(Y≤y), respectively. In this article, we would like to test the no quantile treatment effect hypothesis H ₀: F=G. We develop a bootstrap quantile-treatment-effect test procedure for testing H ₀ under the location-scale shift model. Our test procedure avoids the calculation of the check function (which is non-differentiable at the origin and makes solving the quantile effects difficult in typical quantile regression analysis). The limiting null distribution of the test procedure is derived and the procedure is shown to be consistent against a broad family of alternatives. Simulation studies show that our proposed test procedure attains its type I error rate close to the pre-chosen significance level even for small sample sizes. Our test procedure is illustrated with two real data sets on the lifetimes of guinea pigs from a treatment-control experiment.     

Assessing circular error probable when the errors are elliptical normal

An important problem of continuing interest to engineers is the need to assess the circular error probable (CEP), a measure of the impact accuracy of a projectile or a measure of GPS point positioning accuracy. One of the challenges in addressing this problem is to construct some accurate confidence bounds or intervals for CEP in the small sample settings, where certain amount of systematic biases exist in testing experiments. Currently there is no general method available to deal with this challenge due to the intractability of the distributions of the existing CEP estimators. In this paper, in order to meet this challenge, several new approximate formulas are derived for calculating CEP, which are more accurate than the existing ones but still simple to use. Both exact and empirical expressions for the derivatives of CEP with respect to the population means and variances are also given. Using these formulas, three kinds of confidence bounds or intervals for CEP are proposed, which are based on the parametric bootstrap, the asymptotic distribution, and the Cornish–Fisher expansion, respectively. Moreover, a bias-corrected estimator of CEP is provided. The performances of these procedures are evaluated based on some Monte Carlo simulation studies. Both the theoretical and simulation results show that the Cornish–Fisher expansion-based procedure performs slightly better than the other two procedures when the downrange and cross-range variances are assumed the same. However, when these two variances are different, the simulation demonstrates that the bootstrap approach can be superior to the Cornish–Fisher for the small samples (say n=10), and vice versa for the moderate samples (say n=20).     

Generalized Supremum Tests for the Equality of Cause Specific Hazard Rates

In this paper we propose two new classes of asymptotically distribution-free Renyi-type tests for testing the equality of two risks in a competing risk model with possible censoring. This work extends the work of Aly, Kochar and McKeague [1994, Journal of American Statistical Association, 89, 994–999] and many of the existing tests for this problem belong to these newly proposed classes. The asymptotic properties of the proposed tests are investigated. Simulation studies are done to compare the performance with existing tests. A competing risks data set is analyzed to demonstrate the usefulness of the procedure.     

Estimation of the Population Size in the Classical Occupancy Problem

In the classical occupancy problem where the random variable X is the number of N elements selected by K individuals when each element is equally likely to be chosen by any of the individuals, it is desired to estimate N. Three estimators given in the literature are compared with three estimators derived in this article, two of which are based on Bayesian methods, utilizing a simulation study. One of the Bayes estimators appears to perform the best along with one proposed in the literature in 1968. The estimators are also compared using data obtained from a cemetery in Ohio.     

Unbiased Estimation of the Distribution Function of an Exponential Population Using Order Statistics with Application in Ranked Set Sampling

In this paper we consider the problem of unbiased estimation of the distribution function of an exponential population using order statistics based on a random sample. We present a (unique) unbiased estimator based on a single, say ith, order statistic and study some properties of the estimator for i = 2. We also indicate how this estimator can be utilized to obtain unbiased estimators when a few selected order statistics are available as well as when the sample is selected following an alternative sampling procedure known as ranked set sampling. It is further proved that for a ranked set sample of size two, the proposed estimator is uniformly better than the conventional nonparametric unbiased estimator, further, for a general sample size, a modified ranked set sampling procedure provides an unbiased estimator uniformly better than the conventional nonparametric unbiased estimator based on the usual ranked set sampling procedure.     

On level and powee of anderson and hauck's peocedure for testing equivalence in compative bioavailability

The paper refers to a new procedure of Anderson and Hauck for the well-known problem of deciding whether two means of normal distributions with common but unknown variance differ within an admissible range. Analytical expressions for level and power of this method are derived and some numerical examples are given. A bound for the difference between real and no¬minal level is obtained which does not involve the variance.     

The optimal allocation combination for the two-sided sequential screening procedure based on the individual misclassification error

In this paper, the expected total costs (ETCs) of three kinds of quality cost functions for the two-sided sequential screening procedure (SQSP) based on the individual misclassification error are obtained, where the ETC is the sum of the expected cost of inspection, the expected cost of rejection and the expected cost of quality. The general formulas for all the desired probabilities and three ETCs when k screening variables are allocated into r-stages are derived. The optimal allocation combination for each ETC is determined based on the criterion of minimum ETC. Finally, we give two examples to illustrate the selection of the optimal allocation combination for the SQSP.     

One‐Way anova for Functional Data via Globalizing the Pointwise F‐test

In this paper, we propose and study a new global test, namely, GPF test, for the one‐way anova problem for functional data, obtained via globalizing the usual pointwise F‐test. The asymptotic random expressions of the test statistic are derived, and its asymptotic power is investigated. The GPF test is shown to be root‐n consistent. It is much less computationally intensive than a parametric bootstrap test proposed in the literature for the one‐way anova for functional data. Via some simulation studies, it is found that in terms of size‐controlling and power, the GPF test is comparable with two existing tests adopted for the one‐way anova problem for functional data. A real data example illustrates the GPF test.     

On testing correlation between growth curves

Fitting the growth curves by polynomials, this paper is intended to test whether or not there is any correlation between two characters. The likelihood ratio statistic is derived and is shown to be distributed under the null hypothesis as the product of three independent U statistics as defined in Anderson (1958). Box's procedure is then applied to approximate the critical region. A numerical example is given to illustrate the procedure.     

Exact calculation of power and sample size in bioequivalence studies using two one‐sided tests

The number of subjects in a pharmacokinetic two‐period two‐treatment crossover bioequivalence study is typically small, most often less than 60. The most common approach to testing for bioequivalence is the two one‐sided tests procedure. No explicit mathematical formula for the power function in the context of the two one‐sided tests procedure exists in the statistical literature, although the exact power based on Owen's special case of bivariate noncentral t‐distribution has been tabulated and graphed. Several approximations have previously been published for the probability of rejection in the two one‐sided tests procedure for crossover bioequivalence studies. These approximations and associated sample size formulas are reviewed in this article and compared for various parameter combinations with exact power formulas derived here, which are computed analytically as univariate integrals and which have been validated by Monte Carlo simulations. The exact formulas for power and sample size are shown to improve markedly in realistic parameter settings over the previous approximations. Copyright © 2014 John Wiley & Sons, Ltd.     

Optimum allocation in multivariate stratified sampling design in the presence of nonresponse with Gamma cost function

In multivariate stratified sample survey with L strata, let p-characteristics are defined on each unit of the population. To estimate the unknown p-population means of each characteristic, a random sample is taken out from the population. In multivariate stratified sample survey, the optimum allocation of any characteristic may not be optimum for others. Thus the problem arises to find out an allocation which may be optimum for all characteristics in some sense. Therefore a compromise criterion is needed to workout such allocation. In this paper, the procedure of estimation of p-population means is discussed in the presence of nonresponse when the use of linear cost function is not advisable. A solution procedure is suggested by using lexicographic goal programming problem. The numerical illustrations are given for its practical utility.     

Estimation of regression equation with cauchy disturbances

In this paper we present two methods of estimating a linear regression equation with Cauchy disturbances. The first method uses the maximum likelihood principle and therefore the estimators obtained are consistent. The asymptotic covariance is derived which provides with the necessary statistics for the purpose of making inference in large samples. The second method is the method of least lines which minimizes the sum of absolute errors (MSAE) from the fitted regression. Then these two methods are compared through a Monte Carlo study. The maximum likelihood method emerges superior over the MSAE method. However, the MSAE procedure which does not depend on the distribution of the error term appears to be a close competitor to the maximum likelihood estimator.     

Prediction intervals for generalized-order statistics with random sample size

The problem of predicting future generalized-order statistics, by assuming the future sample size is a random variable, is discussed. A general expression for the coverage probability of the prediction intervals is derived. Since k-records and progressively type-II censored-order statistics are contained in the model of generalized-order statistics, the corresponding results for them can be deduced as special cases. When the future sample size has degenerate, binomial, Poisson and geometric distributions, numerical computations are given. The procedure for finding an optimal prediction interval is presented for each case. Finally, we apply our results to a real data set in life testing given in Lee and Wang [Statistical methods for survival data analysis. Hoboken, NJ: John Wiley and Sons; 2003, p. 58, Table 3.4] for illustrative the proposed procedure in this paper.     

Estimation and Simulation for the M-Wright Function

In this article, a structural form of an M-Wright distributed random variable is derived. The mixture representation then led to a random number generation algorithm. A formal parameter estimation procedure is also proposed. This procedure is needed to make the M-Wright function usable in practice. The asymptotic normality of the estimator is established as well. The estimator and the random number generation algorithm are then tested using synthetic data.     

Inferences on a linear combination of K multivariate normal mean vectors

In this paper, the hypothesis testing and confidence region construction for a linear combination of mean vectors for K independent multivariate normal populations are considered. A new generalized pivotal quantity and a new generalized test variable are derived based on the concepts of generalized p-values and generalized confidence regions. When only two populations are considered, our results are equivalent to those proposed by Gamage et al. [Generalized p-values and confidence regions for the multivariate Behrens–Fisher problem and MANOVA, J. Multivariate Aanal. 88 (2004), pp. 117–189] in the bivariate case, which is also known as the bivariate Behrens–Fisher problem. However, in some higher dimension cases, these two results are quite different. The generalized confidence region is illustrated with two numerical examples and the merits of the proposed method are numerically compared with those of the existing methods with respect to their expected areas, coverage probabilities under different scenarios.     

A Bayesian Approach for Selecting the Best Exponential Population with Interval Censored Samples

In this article, we study the problem of selecting the best population from among several exponential populations based on interval censored samples using a Bayesian approach. A Bayes selection procedure and a curtailed Bayes selection procedure are derived. We show that these two Bayes selection procedures are equivalent. A numerical example is provided to illustrate the application of the two selection procedure. We also use Monte Carlo simulation to study performance of the two selection procedures. The numerical results of the simulation study demonstrate that the curtailed Bayes selection procedure has good performance because it can substantially reduce the duration time of life test experiment.     

Some New Estimators for Small-Area Means With Application to the Assessment of Farmland Values

Regression models that account for main state effects and nested county effects are considered for the assessment of farmland values. Empirical predictors obtained by replacing the unknown variances in the formulas of the optimal predictors by maximum likelihood estimates are presented. The computations are carried out by simple iterations between two SAS procedures. Estimators for the prediction variances are derived, and a modification to secure the robustness of the predictors is proposed. The procedure is applied to data on nonirrigated cropland in the Corn Belt states and is shown to yield predictors with considerably lower prediction mean squared errors than the survey estimators and other regression-type estimators.