Order Statistics from the Generalized Exponential Distribution and Applications

With linear dispersion effects, the standard factorial designs are not optimal estimation of a mean model. A sequential two-stage experimental design procedure has been proposed that first estimates the variance structure, and then uses the variance estimates and the variance optimality criterion to develop a second stage design that efficiency estimates the mean model. This procedure has been compared to an equal replicate design analyzed by ordinary least squares, and found to be a superior procedure in many situations.

However with small first stage sample sizes the variance estiamtes are not reliable, and hence an alternative procedure could be more beneficial. For this reason a Bayesian modification to the two-stage procedure is proposed which will combine the first stage variance estiamtes with some prior variance information that will produce a more efficient procedure. This Bayesian procedure will be compared to the non-Bayesian twostage procedure and to the two one-stage alternative procedures listed above. Finally, a recommendation will be made as to which procedure is preferred in certain situations.     

Dynamic sampling plans in on-line control charts

Three simple dynamic sampling plans for detecting the change point are investigated in the discrete-time case. The first is a two-rate sampling CUSUM procedure. The second is a two-rate sampling Shiryayev-Roberts procedure. The third is a periodic sequential testing procedure. Two problems are discussed. First, simple design methods are provided for practical use. Second, a comparison between the three plans is made in the continuous-time case, which shows that by properly choosing the design parameters, the three plans can be made equally efficient in certain senses.     

Refining binomial confidence intervals

A method for refining an equivariant binomial confidence procedure is presented which, when applied to an existing procedure, produces a new set of equivariant intervals that are uniformly superior. The family of procedures generated from this method constitute a complete class within the class of all equivariant procedures. In certain cases it is shown that this class is also minimal complete. Also, an optimally property, monotone minimaxity, is investigated, and monotone minimax procedures are constructed.     

Bootstrapping quantiles in a fixed design regression model with censored data

We consider the problem of estimating the quantiles of a distribution function in a fixed design regression model in which the observations are subject to random right censoring. The quantile estimator is defined via a conditional Kaplan-Meier type estimator for the distribution at a given design point. We establish an a.s. asymptotic representation for this quantile estimator, from which we obtain its asymptotic normality. Because a complicated estimation procedure is necessary for estimating the asymptotic bias and variance, we use a resampling procedure, which provides us, via an asymptotic representation for the bootstrapped estimator, with an alternative for the normal approximation.     

Two-sample nonparametric stochastic order inference with an application in plant physiology

In this paper, a new nonparametric methodology is developed for testing whether the changing pattern of a response variable over multiple ordered sub-populations from one treatment group differs with the one from another treatment group. The question is formalized into a nonparametric two-sample comparison problem for the stochastic order among subsamples, through U-statistics with accommodations for zero-inflated distributions. A novel bootstrap procedure is proposed to obtain the critical values with given type I error. Following the procedure, bootstrapped p-values are obtained through simulated samples. It is proven that the distribution of the test statistics is independent from the underlying distributions of the subsamples, when certain sufficient statistics provided. Furthermore, this study also develops a feasible framework for power studies to determine sample sizes, which is necessary in real-world applications. Simulation results suggest that the test is consistent. The methodology is illustrated using a biological experiment with a split-plot design, and significant differences in changing patterns of seed weight between treatments are found with relative small subsample sizes.     

Uncorrelated residuals and an exact test for two variance components in experimental design

The error contrasts from an experimental design can be constructed from uncorrelated residuals normally associated with the linear model. In this paper uncorrelated residuals are defined for the linear model that has a design matrix which is less than full rank, typical of many experimental design representations. It transpires in this setting, that for certain choices of uncorrelated residuals, corresponding to recursive type residuals, there is a natural partition of information when two variance components are known to be present. Under an assumtion of normality of errors this leads to construction of appropriate F-tests for testing heteroscedasticity. The test, which can be optimal, is applied to two well known data sets to illustrate its usefullness.     

Wilcoxon–Mann–Whitney test for stratified samples and Efron's paradox dice

Two-treatment multi-center clinical trials are the most common type of clinical trials in practice. The aim of this paper is to discuss a curious property of certain standard nonparametric procedures used in the analysis of such clinical trials. Different analyses of a simulated data example are presented, which lead to contrasting and surprising results. The source of the potentially misleading outcome is then explored while relating the simulated data with the concept of Efron's paradox dice and the notion of nontransitivity. With the root of the problem established, an alternate nonparametric method from the literature is shown to address the problem. Finally, pointing out an interpretational concern of using the alternate procedure, a modification to this procedure is also suggested and corresponding theoretical results are presented.     

Designs for discrimination between binary response models

We propose a test statistic for discrimination between alternative univariate binary response models which is asymptotically equivalent to the likelihood ratio statistic and Pearson's goodness of fit statistic. We propose an optimal design procedure. Under certain conditions we prove that the maximum value of the power can be obtained when the degrees of freedom of the test statistic is one. Several mathematical properties of the incomplete gamma function ratio and the non-central chi-squared distribution are required in the discussion and these are established.     

Bayesian methods for mark-resighting surveys

We consider the problem of estimating the size of a closed population based on the results of a certain type of mark-resighting sampling design. The design is similar to the commonly used multiple capture-recapture design, yet in some cases economically more feasible and easy to use. Sampling is done by first tagging a number of randomly selected animals with visible markers and later randomly sighting them (for instance, for large animals by visually sampling from a helicopter) and counting the number of tagged animals. In this paper, we look at Bayesian methods for point and interval estimation of population size for this design. An example involving estimation of mountain sheep, a couple of simulated examples and simulation studies are given to demonstrate the advantages of the proposed procedure over the other available approximate procedures.     

A HEURISTIC TECHNIQUE FOR TESTING FOR THE PRESENCE OF INTERACTION IN NONREPLICATED FACTORIAL EXPERIMENTS

Abstract  This paper suggests univariate and multivariate techniques for investigating interaction in nonreplicated factorial experiments. The tests can detect certain types of interaction, but they are not powerful against all possible alternative hypotheses. The two-way factorial experiment is discussed in some detail and an example is used to demonstrate the procedure. The procedure is compared to other tests for interaction. These comparisons show that the procedure can detect some types of interaction which other tests cannot. Likewise other tests can detect interaction this procedure cannot.     

Discrimination between two binary data models: sequentially designed experiments

In this paper we propose a sequential procedure to design optimum experiments for discriminating between two binary data models. For the problem to be fully specified, not only the mode1link functions should be provided but also their associated linear predictor structures. Further, we suppose that one of the models is true, albeit it is not known which of them. Under these assumptions the procedure consists of making sequential choices of single experimental units to discriminate between the rival models as efficiently as possible. Depending on whether the models are nested or not, alternative methods are proposed.

To illustrate the procedure, a simulation study for the classical case of pro bit versus logit model is presented. It enables us to estimate the total sample sizes required to gain a certain power of discrimination and compare them to sample sizes for methods that were previously suggested in the literature.     

A method for screening active effects in supersaturated designs

A supersaturated design (SSD) is a design whose run size is not enough for estimating all the main effects. The goal in conducting such a design is to identify, presumably only a few, relatively dominant active effects with a cost as low as possible. However, data analysis of such designs remains primitive: traditional approaches are not appropriate in such a situation and several methods which were proposed in the literature in recent years are effective when used to analyze two-level SSDs. In this paper, we introduce a variable selection procedure, called the PLSVS method, to screen active effects in mixed-level SSDs based on the variable importance in projection which is an important concept in the partial least-squares regression. Simulation studies show that this procedure is effective.     

Taylor linearization sampling errors and design effects for poverty measures and other complex statistics

A systematic procedure for the derivation of linearized variables for the estimation of sampling errors of complex nonlinear statistics involved in the analysis of poverty and income inequality is developed. The linearized variable extends the use of standard variance estimation formulae, developed for linear statistics such as sample aggregates, to nonlinear statistics. The context is that of cross-sectional samples of complex design and reasonably large size, as typically used in population-based surveys. Results of application of the procedure to a wide range of poverty and inequality measures are presented. A standardized software for the purpose has been developed and can be provided to interested users on request. Procedures are provided for the estimation of the design effect and its decomposition into the contribution of unequal sample weights and of other design complexities such as clustering and stratification. The consequence of treating a complex statistic as a simple ratio in estimating its sampling error is also quantified. The second theme of the paper is to compare the linearization approach with an alternative approach based on the concept of replication, namely the Jackknife repeated replication (JRR) method. The basis and application of the JRR method is described, the exposition paralleling that of the linearization method but in somewhat less detail. Based on data from an actual national survey, estimates of standard errors and design effects from the two methods are analysed and compared. The numerical results confirm that the two alternative approaches generally give very similar results, though notable differences can exist for certain statistics. Relative advantages and limitations of the approaches are identified.     

Attaining specified reliability in component design: a stochastic approximation approach

In the design, manufacture and maintenance of components, particular attention is paid to component reliability R, the probability that the strength X of a component will exceed a stress Y to which it will be subjected. The problem addressed here is the design (or redesign) of a compoFent to meet a specified reliability R*. While certain characteristics of the random variables X and Y are assumed (symmetry of X about a unique median for example) it is not assumed that the form of the distribution of (X,Y) is known, nor that X and Y are independent. A design is recomnended based on a variation of the stochastic approximation procedure due to Dupac and Kral (1972) which in general estimates recursively the root of a regression curve assuming both independent and dependent regression variables are subject to experimental error.     

Adaptive nonparametric estimation in the presence of dependence

We consider nonparametric estimation problems in the presence of dependent data, notably nonparametric regression with random design and nonparametric density estimation. The proposed estimation procedure is based on a dimension reduction. The minimax optimal rate of convergence of the estimator is derived assuming a sufficiently weak dependence characterised by fast decreasing mixing coefficients. We illustrate these results by considering classical smoothness assumptions. However, the proposed estimator requires an optimal choice of a dimension parameter depending on certain characteristics of the function of interest, which are not known in practice. The main issue addressed in our work is an adaptive choice of this dimension parameter combining model selection and Lepski's method. It is inspired by the recent work of Goldenshluger and Lepski [(2011), ‘Bandwidth Selection in Kernel Density Estimation: Oracle Inequalities and Adaptive Minimax Optimality’, The Annals of Statistics, 39, 1608–1632]. We show that this data-driven estimator can attain the lower risk bound up to a constant provided a fast decay of the mixing coefficients.     

Statistical inference for the accelerated failure time model under two-stage generalized case–cohort design

Abstract

In this article, we propose a two-stage generalized case–cohort design and develop an efficient inference procedure for the data collected with this design. In the first-stage, we observe the failure time, censoring indicator and covariates which are easy or cheap to measure, and in the second-stage, select a subcohort by simple random sampling and a subset of failures in remaining subjects from the first-stage subjects to observe their exposures which are different or expensive to measure. We derive estimators for regression parameters in the accelerated failure time model under the two-stage generalized case–cohort design through the estimated augmented estimating equation and the kernel function method. The resulting estimators are shown to be consistent and asymptotically normal. The finite sample performance of the proposed method is evaluated through the simulation studies. The proposed method is applied to a real data set from the National Wilm’s Tumor Study Group.     

Optimal design of the adaptive exponentially weighted moving average control chart over a range of mean shifts

The adaptive exponentially weighted moving average (AEWMA) control chart is a smooth combination of the Shewhart and exponentially weighted moving average (EWMA) control charts. This chart was proposed by Cappizzi and Masarotto (2003) to achieve a reasonable performance for both small and large shifts. Cappizzi and Masarotto (2003) used a pair of shifts in designing their control chart. In this study, however, the process mean shift is considered as a random variable with a certain probability distribution and the AEWMA control chart is optimized for a wide range of mean shifts according to that probability distribution and not just for a pair of shifts. Using the Markov chain technique, the results show that the new optimization design can improve the performance of the AEWMA control chart from an overall point of view relative to the various designs presented by Cappizzi and Masarotto (2003). Optimal design parameters that achieve the desired in-control average run length (ARL) are computed in several cases and formulas used to find approximately their values are given. Using these formulas, the practitioner can compute the optimal design parameters corresponding to any desired in-control ARL without the need to apply the optimization procedure. The results obtained by these formulas are very promising and would particularly facilitate the design of the AEWMA control chart for any in-control ARL value.     

Experimental design for response surface gradient estimation

In this paper, a new experimental design and gradient estimation procedure is presented for Phase I response surface optimization. The design is motivated by basic principles of differential calculus, which imply that if a point in Rⁿ has been reached by exactly minimizing a function along a given direction, then the gradient of the function at that point must be orthogonal to the search direction followed. While exact line search is not required for the new design to be effective, this principle implies that the dimension of the gradient estimation procedure may often be reduced from n to n-1 variables, and the experimenter is able to concentrate experimental effort within the most productive region around the center of the design. The new design and gradient estimation procedures are presented, and bias and variance properties are derived. The effectiveness of the new design is shown to depend on the experimenter's ability to terminate line search within a near-stationary region of the line search function A simple heuristic is presented which indicates whether the new design should be used at a given experimental region.     

Calibration estimation in dual-frame surveys

Survey statisticians make use of auxiliary information to improve estimates. One important example is calibration estimation, which constructs new weights that match benchmark constraints on auxiliary variables while remaining “close” to the design weights. Multiple-frame surveys are increasingly used by statistical agencies and private organizations to reduce sampling costs and/or avoid frame undercoverage errors. Several ways of combining estimates derived from such frames have been proposed elsewhere; in this paper, we extend the calibration paradigm, previously used for single-frame surveys, to calculate the total value of a variable of interest in a dual-frame survey. Calibration is a general tool that allows to include auxiliary information from two frames. It also incorporates, as a special case, certain dual-frame estimators that have been proposed previously. The theoretical properties of our class of estimators are derived and discussed, and simulation studies conducted to compare the efficiency of the procedure, using different sets of auxiliary variables. Finally, the proposed methodology is applied to real data obtained from the Barometer of Culture of Andalusia survey.     

Designing selection experiments with bernoulli populations

This paper is concerned with a fixed size subset selection problem for Bernoulli populations in the framework of the indifference zone approach. The goal is to select s populationswhich contain at least c of those with the t largest success probabilities. In order to control the probability of correct selection over the preference zone extensive tables of exact minimum sample sizes have been prepared to implement the single-stage procedure generalized from the well-known Sobel-Huyett procedure. It is shown how the tables can also be employed to design certain closedsequential procedures. These procedures curtail the sampling process of the single-stage procedureand may differ in their sampling rules. Two procedures working with play-the-winner rules are described in detail