Multiple comparisons with the average for exponential location parameters based on doubly censored sample under heteroscedasticity

In this paper, a one-stage multiple comparison procedures with the average for exponential location parameters based on the doubly censored sample under heteroscedasticity is proposed. These intervals can be used to identify a subset which includes all no-worse-than-the-average treatments in an experimental design and to identify better-than-the-average, worse-than-the-average and not-much-different-from-the-average products in agriculture, emerging market, pharmaceutical industries. The critical values are tabulated in a table for practical use. A simulation study on the confidence length and coverage probabilities is done. At last, an example of comparing four drugs in the treatment of leukemia is given to demonstrate the proposed procedures.     

Approximate one-sided tolerance limits in random effects model and in some mixed models and comparisons

An approximate closed-form one-sided tolerance limit (TL) in a general mixed model is proposed. One-sided TLs for the distribution of observable random variable and for the distribution of unobservable random variable in one-way random model are obtained as special cases from the one for the general mixed model. Applications to a two-way nested random model are also given. The merits of the TLs are evaluated using Monte Carlo simulation and compared with the existing ones. Our comparison studies indicate that the approximate TLs are quite satisfactory for all parameter and sample size configurations, and better than the existing ones in some cases. Approximate confidence intervals for exceedance probabilities in one-way random effects model are also proposed. The procedures are illustrated using three examples.     

A Goodness-of-Fit Test for Location-Scale Max-Stable Distributions

In this article, a technique based on the sample correlation coefficient to construct goodness-of-fit tests for max-stable distributions with unknown location and scale parameters and finite second moment is proposed. Specific details to test for the Gumbel distribution are given, including critical values for small sample sizes as well as approximate critical values for larger sample sizes by using normal quantiles. A comparison by Monte Carlo simulation shows that the proposed test for the Gumbel hypothesis is substantially more powerful than some other known tests against some alternative distributions with positive skewness coefficient.     

EMPIRICAL EVIDENCE FOR ADAPTIVE CONFIDENCE INTERVALS AND IDENTIFICATION OF OUTLDSRS USING METHODS OF TRIMMING

In an attempt to apply robust procedures, conventional t-tables are used to approximate critical values of a Studentized t-statistic which is formed from the ratio of a trimmed mean to the square root of a suitably normed Winsorized sum of squared deviations. It is shown here that the approximation is poor if the proportion of trimming is chosen to depend on the data. Instead a data dependent alternative is given which uses adaptive trimming proportions and confidence intervals based on trimmed likelihood statistics. Resulting statistics have high efficiency at the normal model, proper coverage for confidence intervals, yet retain breakdown point one half. Average lengths of confidence intervals are competitive with those of recent Studentized confidence intervals based on the biweight over a range of underlying distributions. In addition, the adaptive trimming is used to identify potential outliers. Evidence in the form of simulations and data analysis support the new adaptive trimming approach.     

Inference on unknown parameters of a Burr distribution under hybrid censoring

Based on hybrid censored data, the problem of making statistical inference on parameters of a two parameter Burr Type XII distribution is taken up. The maximum likelihood estimates are developed for the unknown parameters using the EM algorithm. Fisher information matrix is obtained by applying missing value principle and is further utilized for constructing the approximate confidence intervals. Some Bayes estimates and the corresponding highest posterior density intervals of the unknown parameters are also obtained. Lindley’s approximation method and a Markov Chain Monte Carlo (MCMC) technique have been applied to evaluate these Bayes estimates. Further, MCMC samples are utilized to construct the highest posterior density intervals as well. A numerical comparison is made between proposed estimates in terms of their mean square error values and comments are given. Finally, two data sets are analyzed using proposed methods.     

Multiple comparisons with a control for exponential location parameters under heteroscedasticity

In this paper, a new design-oriented two-stage two-sided simultaneous confidence intervals, for comparing several exponential populations with control population in terms of location parameters under heteroscedasticity, are proposed. If there is a prior information that the location parameter of k exponential populations are not less than the location parameter of control population, one-sided simultaneous confidence intervals provide more inferential sensitivity than two-sided simultaneous confidence intervals. But the two-sided simultaneous confidence intervals have advantages over the one-sided simultaneous confidence intervals as they provide both lower and upper bounds for the parameters of interest. The proposed design-oriented two-stage two-sided simultaneous confidence intervals provide the benefits of both the two-stage one-sided and two-sided simultaneous confidence intervals. When the additional sample at the second stage may not be available due to the experimental budget shortage or other factors in an experiment, one-stage two-sided confidence intervals are proposed, which combine the advantages of one-stage one-sided and two-sided simultaneous confidence intervals. The critical constants are obtained using the techniques given in Lam [9,10]. These critical constant are compared with the critical constants obtained by Bonferroni inequality techniques and found that critical constant obtained by Lam [9,10] are less conservative than critical constants computed from the Bonferroni inequality technique. Implementation of the proposed simultaneous confidence intervals is demonstrated by a numerical example.     

Statistical aspects of linkage analysis in interlaboratory studies

This paper investigates statistical issues that arise in interlaboratory studies known as Key Comparisons when one has to link several comparisons to or through existing studies. An approach to the analysis of such a data is proposed using Gaussian distributions with heterogeneous variances. We develop conditions for the set of sufficient statistics to be complete and for the uniqueness of uniformly minimum variance unbiased estimators (UMVUE) of the contrast parametric functions. New procedures are derived for estimating these functions with estimates of their uncertainty. These estimates lead to associated confidence intervals for the laboratories (or studies) contrasts. Several examples demonstrate statistical inference for contrasts based on linkage through the pilot laboratories. Monte Carlo simulation results on performance of approximate confidence intervals are also reported.     

Some Properties of the Exact and Score Methods for Binomial Proportion and Sample Size Calculation

In this article, we point out some interesting relations between the exact test and the score test for a binomial proportion p. Based on the properties of the tests, we propose some approximate as well as exact methods of computing sample sizes required for the tests to attain a specified power. Sample sizes required for the tests are tabulated for various values of p to attain a power of 0.80 at level 0.05. We also propose approximate and exact methods of computing sample sizes needed to construct confidence intervals with a given precision. Using the proposed exact methods, sample sizes required to construct 95% confidence intervals with various precisions are tabulated for p = .05(.05).5. The approximate methods for computing sample sizes for score confidence intervals are very satisfactory and the results coincide with those of the exact methods for many cases.     

Comparing Several Variances with a Control in Direction-mixed Families of Hypotheses

This article deals with simultaneous comparisons of k(k ? 1) populations with a control population with respect to variance in direction-mixed families of hypotheses. Computation of critical constants required for the implementation of the proposed procedures is discussed and selected values of the critical constants are tabulated. Power comparison of the proposed procedures with their existing competitors is carried out using Monte Carlo simulation. An example is given to illustrate the implementation of proposed procedures.     

Small sample asymptotic inference for the coefficient of variation: normal and nonnormal models

In applied statistics, the coefficient of variation is widely calculated and interpreted even when the sample size of the data set is very small. However, confidence intervals for the coefficient of variation are rarely reported. One of the reasons is the exact confidence interval for the coefficient of variation, which is given in Lehmann (Testing Statistical Hypotheses, 2nd Edition, Wiley, New York, 1996), is very difficult to calculate. Various asymptotic methods have been proposed in literature. These methods, in general, require the sample size to be large. In this article, we will apply a recently developed small sample asymptotic method to obtain approximate confidence intervals for the coefficient of variation for both normal and nonnormal models. These small sample asymptotic methods are very accurate even for very small sample size. Numerical examples are given to illustrate the accuracy of the proposed method.     

A Simple Approximate Procedure for Constructing Binomial and Poisson Tolerance Intervals

The problems of constructing tolerance intervals for the binomial and Poisson distributions are considered. Closed-form approximate equal-tailed tolerance intervals (that control percentages in both tails) are proposed for both distributions. Exact coverage probabilities and expected widths are evaluated for the proposed equal-tailed tolerance intervals and the existing intervals. Furthermore, an adjustment to the nominal confidence level is suggested so that an equal-tailed tolerance interval can be used as a tolerance interval which includes a specified proportion of the population, but does not necessarily control percentages in both tails. Comparison of such coverage-adjusted tolerance intervals with respect to coverage probabilities and expected widths indicates that the closed-form approximate tolerance intervals are comparable with others, and less conservative, with minimum coverage probabilities close to the nominal level in most cases. The approximate tolerance intervals are simple and easy to compute using a calculator, and they can be recommended for practical applications. The methods are illustrated using two practical examples.     

Tolerance intervals for symmetric location-scale families based on uncensored or censored samples

Exact methods for constructing two-sided tolerance intervals (TIs) and tolerance intervals that control percentages in both tails for a location-scale family of distributions are proposed. The proposed methods are illustrated by constructing TIs for a normal, logistic, and Laplace (double exponential) distributions based on type II singly censored samples. Factors for constructing one-sided and two-sided TIs for a logistic distribution are tabulated for the case of uncensored samples. Factors for constructing TIs based on censored samples for all three distributions are also tabulated. The factors for all cases are estimated by Monte Carlo simulation. An adjustment to the tolerance factors based on type II censored samples is proposed so that they can be used to find approximate TIs based on type I censored samples. Coverage studies of the approximate TIs based on type I censored samples indicate that the approximation is satisfactory as long as the proportion of censored observations is no more than 0.70. The methods are illustrated using some practical examples.     

βExpectation and β-content tolerance intervals for dependent observations

Methods of constructing exact tolerance intervals (β-expectation and β-content) for independent observations are well known. For the case of dependent observations, obtaining exact results is not possible. In this article we provide an approximate method of constructing β-expectation tolerance intervals via a Taylor series expansion. Examples of independent observations are considered to compare the intervals constructed with those obtained by the exact method. For the case of non-stationary type processes we have proposed a method of constructing approximate β-content tolerance intervals. Once again an example is given to illustrate the results.     

Multi-sample inference for the simple-tree alternative based on one-sample confidence intervals

The problem of testing homogeneity of several group means is considered against some patterned alternatives for the one-way classified data. The patterns of interest include the simple-tree and the trend alternatives. The approach is to begin with some suitably defined one-sample confidence intervals for the groups in a graphical display. Depending on the pattern of interest, orientation features of the display are examined, more formally, using proposed overall tests or rules. In the classical setup under normality, the case of known common variance is treated in detail; extensions to the case of unknown variance are indicated. When normality is in doubt, a nonparametric procedure based on the sign test is proposed. The necessary critical values are percentiles of either a multivariate normal distribution or a multivariate t-distribution. Although some existing tables can be used for the critical values (or the P-values) in some special cases, in general, the use of simulations is recommended and the steps are detailed in the appendix. An illustrative numerical example is provided.     

Short confidence intervals for variance components

Four approximate methods are proposed to construct confidence intervals for the estimation of variance components in unbalanced mixed models. The first three methods are modifications of the Wald, arithmetic and harmonic mean procedures, see Harville and Fenech (1985), while the fourth is an adaptive approach, combining the arithmetic and harmonic mean procedures. The performances of the proposed methods were assessed by a Monte Carlo simulation study. It was found that the intervals based on Wald's method maintained the nominal confidence levels across all designs and values of the parameters under study. On the other hand, the arithmetic (harmonic) mean method performed well for small (large) values of the variance component, relative to the error variance component. The adaptive procedure performed rather well except for extremely unbalanced designs. Further, compared with equal tails intervals, the intervals which use special tables, e.g., Table 678 of Tate and Klett (1959), provided adequate coverage while having much shorter lengths and are thus recommended for use in practice.     

Improved closed-form prediction intervals for binomial and Poisson distributions

The problems of constructing prediction intervals for the binomial and Poisson distributions are considered. Available approximate, exact and conditional methods for both distributions are reviewed and compared. Simple approximate prediction intervals based on the joint distribution of the past samples and the future sample are proposed. Exact coverage studies and expected widths of prediction intervals show that the new prediction intervals are comparable to or better than the available ones in most cases. The methods are illustrated using two practical examples.     

A Monte Carlo simulation study of two-sided tolerance intervals in balanced one-way random effects model for non-normal errors

Recently, Ong and Mukerjee [Probability matching priors for two-sided tolerance intervals in balanced one-way and two-way nested random effects models. Statistics. 2011;45:403–411] developed two-sided Bayesian tolerance intervals, with approximate frequentist validity, for a future observation in balanced one-way and two-way nested random effects models. These were obtained using probability matching priors (PMP). On the other hand, Krishnamoorthy and Lian [Closed-form approximate tolerance intervals for some general linear models and comparison studies. J Stat Comput Simul. 2012;82:547–563] studied closed-form approximate tolerance intervals by the modified large-sample (MLS) approach. We compare the performances of these two approaches for normal as well as non-normal error distributions. Monte Carlo simulation methods are used to evaluate the resulting tolerance intervals with regard to achieved confidence levels and expected widths. It turns out that PMP tolerance intervals are less conservative for data with large number of classes and small number of observations per class and the MLS procedure is preferable for smaller sample sizes.     

Function-based hypothesis testing in censored two-sample location-scale models

Function-based hypothesis testing in two-sample location-scale models has been addressed for uncensored data using the empirical characteristic function. A test of adequacy in censored two-sample location-scale models is lacking, however. A plug-in empirical likelihood approach is used to introduce a test statistic, which, asymptotically, is not distribution free. Hence for practical situations bootstrap is necessary for performing the test. A multiplier bootstrap and a model appropriate resampling procedure are given to approximate critical values from the null asymptotic distribution. Although minimum distance estimators of the location and scale are deployed for the plug-in, any consistent estimators can be used. Numerical studies are carried out that validate the proposed testing method, and real example illustrations are given.

     

Accurate computation of single and product moments of order statistics under progressive censoring

An accurate procedure is proposed to calculate approximate moments of progressive order statistics in the context of statistical inference for lifetime models. The study analyses the performance of power series expansion to approximate the moments for location and scale distributions with high precision and smaller deviations with respect to the exact values. A comparative analysis between exact and approximate methods is shown using some tables and figures. The different approximations are applied in two situations. First, we consider the problem of computing the large sample variance–covariance matrix of maximum likelihood estimators. We also use the approximations to obtain progressively censored sampling plans for log-normal distributed data. These problems illustrate that the presented procedure is highly useful to compute the moments with precision for numerous censoring patterns and, in many cases, is the only valid method because the exact calculation may not be applicable.     

The application of the durbin-watson test to the dynamic regression model under normal and non-normal errors

Until recently, a difficulty with applying the Durbin-Watson (DW) test to the dynamic linear regression model has been the lack of appropriate critical values. Inder (1986) used a modified small-disturbance distribution (SDD) to find approximate critical values. King and Wu (1991) showed that the exact SDD of the DW statistic is equivalent to the distribution of the DW statistic from the regression with the lagged dependent variables replaced by their means. Unfortunately, these means are unknown although they could be estimated by the actual variable values. This provides a justification for using the exact critical values of the DW statistic from the regression with the lagged dependent variables treated as non-stochastic regressors. Extensive Monte Carlo experiments are reported in this paper. They show that this approach leads to reasonably accurate critical values, particularly when two lags of the dependent variable are present. Robustness to non-normality is also investigated.