Comparisons of ten estimation methods for the parameters of Marshall–Olkin extended exponential distribution

The aim of this article is to compare via Monte Carlo simulations the finite sample properties of the parameter estimates of the Marshall–Olkin extended exponential distribution obtained by ten estimation methods: maximum likelihood, modified moments, L-moments, maximum product of spacings, ordinary least-squares, weighted least-squares, percentile, Crámer–von-Mises, Anderson–Darling, and Right-tail Anderson–Darling. The bias, root mean-squared error, absolute and maximum absolute difference between the true and estimated distribution functions are used as criterion of comparison. The simulation study reveals that the L-moments and maximum products of spacings methods are highly competitive with the maximum likelihood method in small as well as in large-sized samples.     

A goodness-of-fit test for overdispersed binomial (or multinomial) models

Overdispersion or extra variation is a common phenomenon that occurs when binomial (multinomial) data exhibit larger variances than that permitted by the binomial (multinomial) model. This arises when the data are clustered or when the assumption of independence is violated. Goodness-of-fit (GOF) tests available in the overdispersion literature have focused on testing for the presence of overdispersion in the data and hence they are not applicable for choosing between the several competing overdispersion models. In this paper, we consider a GOF test proposed by Neerchal and Morel [1998. Large cluster results for two parametric multinomial extra variation models. J. Amer. Statist. Assoc. 93(443), 1078–1087], and study its distributional properties and performance characteristics. This statistic is a direct analogue of the usual Pearson chi-squared statistic, but is also applicable when the clusters are not necessarily of the same size. As this test statistic is for testing model adequacy against the alternative that the model is not adequate, it is applicable in testing two competing overdispersion models.     

One and two sample confidence intervals for estimating the mean of skewed populations: an empirical comparative study

In this paper we consider and propose some confidence intervals for estimating the mean or difference of means of skewed populations. We extend the median t interval to the two sample problem. Further, we suggest using the bootstrap to find the critical points for use in the calculation of median t intervals. A simulation study has been made to compare the performance of the intervals and a real life example has been considered to illustrate the application of the methods.     

Small confidence sets for the mean of a spherically symmetric distribution

Summary.  Suppose that X has a k -variate spherically symmetric distribution with mean vector θ and identity covariance matrix. We present two spherical confidence sets for θ , both centred at a positive part Stein estimator     . In the first, we obtain the radius by approximating the upper α -point of the sampling distribution of     by the first two non-zero terms of its Taylor series about the origin. We can analyse some of the properties of this confidence set and see that it performs well in terms of coverage probability, volume and conditional behaviour. In the second method, we find the radius by using a parametric bootstrap procedure. Here, even greater improvement in terms of volume over the usual confidence set is possible, at the expense of having a less explicit radius function. A real data example is provided, and extensions to the unknown covariance matrix and elliptically symmetric cases are discussed.     

Comparing Several Regression Models with Unequal Variances

We consider the problem of comparing k regression models, when the variances are not assumed to be equal. For this problem, the classical F test can lead to misleading results, and there is no simple test which adequately controls the size when the sample sizes are small. For k = 2, the most widely used test is the “weighted F test,” also known as the “asymptotic Chow test.” But this test does not work well for small samples, and various modifications have been proposed in the literature. For k > 2, few tests are available and only the parametric-bootstrap (PB) test of Tian et al. (2009) Tian, L., Ma, C., Vexler, A. (2009). A parametric bootstrap test for comparing heteroscedastic regression models. Communications in Statistics—Simulation and Computation, 38, 1026–1036.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar] controls the size fairly adequately. In this article, we propose three fairly simple F tests which can easily be applied in the general case, k ? 2, and avoid the complications of the PB test. Our simulations indicate that these tests have satisfactory performance. Also, our simulations confirm that the power properties of our proposed tests are similar to the PB test. Therefore, our proposed tests provide simple alternatives to the PB test, which can easily be used by practitioners who may not be familiar with the PB.     

Adaptive bootstrap tests and their competitors in the c -sample location problem

This paper deals with a power comparison of different types of tests, parametric, nonparametric, robustified and adaptive ones for the two-sided c -sample location problem. A robustness study on level f in the case of heteroscedasticity and non-normal distributions is included in our study, too. First of all, we consider an adaptive test based on Hogg's concept and two adaptive Bootstrap tests using Hogg's principle. It turns out that the adaptive Hogg-test is the best one in the case of homoscedasticity but for heteroscedasticity, an adaptive Bootstrap test using Hogg's principle is preferable.     

Weibull extension of bivariate exponential regression model with different frailty distributions

We propose bivariate Weibull regression model with frailty in which dependence is generated by a gamma or positive stable or power variance function distribution. We assume that the bivariate survival data follows bivariate Weibull of Hanagal (Econ Qual Control 19:83–90, 2004; Econ Qual Control 20:143–150, 2005a; Stat Pap 47:137–148, 2006a; Stat Methods, 2006b). There are some interesting situations like survival times in genetic epidemiology, dental implants of patients and twin births (both monozygotic and dizygotic) where genetic behavior (which is unknown and random) of patients follows known frailty distribution. These are the situations which motivate to study this particular model. David D. Hanagal is on leave from Department of Statistics, University of Pune, Pune 411007, India.     

Control charts for generalized exponential distribution percentiles

The generalized exponential (GE) distribution, which was introduced by Mudholkar and Srivastava in 1993 Mudholkar, G. S., Srivastava, D. K. (1993). Exponentiated Weibull family for analyzing bathtub failure data. IEEE Transactions on Reliability 42:299–302. [Google Scholar], has been studied for various applications of lifetime modelings. In this article, five control charts, that comprise the Shewhart-type chart and four parametric bootstrap charts based on maximum likelihood estimation method, the moment estimation method, probability plot method, and least-square error method for the GE percentiles, are investigated. An extensive Monte Carlo simulation study is conducted to compare the performance among all five control charts in terms of average run length. Finally, an example is given for illustration.     

Mean-Squared errors of small area estimators under a multivariate linear model for repeated measures data

In this paper, we discuss the derivation of the first and second moments for the proposed small area estimators under a multivariate linear model for repeated measures data. The aim is to use these moments to estimate the mean-squared errors (MSE) for the predicted small area means as a measure of precision. At the first stage, we derive the MSE when the covariance matrices are known. At the second stage, a method based on parametric bootstrap is proposed for bias correction and for prediction error that reflects the uncertainty when the unknown covariance is replaced by its suitable estimator.     

Properties of Model-Averaged BMDLs: A Study of Model Averaging in Dichotomous Response Risk Estimation

Model averaging (MA) has been proposed as a method of accounting for model uncertainty in benchmark dose (BMD) estimation. The technique has been used to average BMD dose estimates derived from dichotomous dose-response experiments, microbial dose-response experiments, as well as observational epidemiological studies. While MA is a promising tool for the risk assessor, a previous study suggested that the simple strategy of averaging individual models' BMD lower limits did not yield interval estimators that met nominal coverage levels in certain situations, and this performance was very sensitive to the underlying model space chosen. We present a different, more computationally intensive, approach in which the BMD is estimated using the average dose-response model and the corresponding benchmark dose lower bound (BMDL) is computed by bootstrapping. This method is illustrated with TiO(2) dose-response rat lung cancer data, and then systematically studied through an extensive Monte Carlo simulation. The results of this study suggest that the MA-BMD, estimated using this technique, performs better, in terms of bias and coverage, than the previous MA methodology. Further, the MA-BMDL achieves nominal coverage in most cases, and is superior to picking the "best fitting model" when estimating the benchmark dose. Although these results show utility of MA for benchmark dose risk estimation, they continue to highlight the importance of choosing an adequate model space as well as proper model fit diagnostics.