Discussion on the common threshold of several exponential distributions with different rates

The integration of results of independent studies in order to make inferences about a common threshold is an important problem with many practical applications. In this article, we apply the generalized variable method to make inferences on the common threshold of several exponential distributions when the scale (or rate) parameters are unknown and unequal. The merits of the proposed method are computed numerically and compared with other existing methods. Numerical results of both simulation studies and real data analyses show that the proposed method is applicable and its performance is better than other methods even when sample sizes are small.     

Testing the equality means of several log-normal distributions

Researches propose various methods for comparing the means of two log-normal distributions. Some of these methods have been recently extended to test the equality means of several log-normal populations. Investigations show that none of the established methods is satisfactory. In this article, we provide three methods based on the computational approach test, which is a parametric bootstrap approach, for testing the means of several log-normal distributions. Further, we compare our methods with the existing methods through Monte Carlo simulation. The numerical results show that the Type I errors of these procedures are satisfactory regardless of the sample size, number of populations, and the true parameters. Finally, we explain the considered methods by real examples.     

Simultaneous confidence intervals for comparing several inverse Gaussian means under heteroscedasticity

Recently, Zhang [Simultaneous confidence intervals for several inverse Gaussian populations. Stat Probab Lett. 2014;92:125–131] proposed simultaneous pairwise confidence intervals (SPCIs) based on the fiducial generalized pivotal quantity concept to make inferences about the inverse Gaussian means under heteroscedasticity. In this paper, we propose three new methods for constructing SPCIs to make inferences on the means of several inverse Gaussian distributions when scale parameters and sample sizes are unequal. One of the methods results in a set of classic SPCIs (in the sense that it is not simulation-based inference) and the two others are based on a parametric bootstrap approach. The advantages of our proposed methods over Zhang’s (2014) method are: (i) the simulation results show that the coverage probability of the proposed parametric bootstrap approaches is fairly close to the nominal confidence coefficient while the coverage probability of Zhang’s method is smaller than the nominal confidence coefficient when the number of groups and the variance of groups are large and (ii) the proposed set of classic SPCIs is conservative in contrast to Zhang’s method.     

Comparing the mean vectors of two independent multivariate log-normal distributions

The multivariate log-normal distribution is a good candidate to describe data that are not only positive and skewed, but also contain many characteristic values. In this study, we apply the generalized variable method to compare the mean vectors of two independent multivariate log-normal populations that display heteroscedasticity. Two generalized pivotal quantities are derived for constructing the generalized confidence region and for testing the difference between two mean vectors. Simulation results indicate that the proposed procedures exhibit satisfactory performance regardless of the sample sizes and heteroscedasticity. The type I error rates obtained are consistent with expectations and the coverage probabilities are close to the nominal level when compared with the other method which is currently available. These features make the proposed method a worthy alternative for inferential analysis of problems involving multivariate log-normal means. The results are illustrated using three examples.     

Exponentiated Weibull regression for time-to-event data

The Weibull, log-logistic and log-normal distributions are extensively used to model time-to-event data. The Weibull family accommodates only monotone hazard rates, whereas the log-logistic and log-normal are widely used to model unimodal hazard functions. The increasing availability of lifetime data with a wide range of characteristics motivate us to develop more flexible models that accommodate both monotone and nonmonotone hazard functions. One such model is the exponentiated Weibull distribution which not only accommodates monotone hazard functions but also allows for unimodal and bathtub shape hazard rates. This distribution has demonstrated considerable potential in univariate analysis of time-to-event data. However, the primary focus of many studies is rather on understanding the relationship between the time to the occurrence of an event and one or more covariates. This leads to a consideration of regression models that can be formulated in different ways in survival analysis. One such strategy involves formulating models for the accelerated failure time family of distributions. The most commonly used distributions serving this purpose are the Weibull, log-logistic and log-normal distributions. In this study, we show that the exponentiated Weibull distribution is closed under the accelerated failure time family. We then formulate a regression model based on the exponentiated Weibull distribution, and develop large sample theory for statistical inference. We also describe a Bayesian approach for inference. Two comparative studies based on real and simulated data sets reveal that the exponentiated Weibull regression can be valuable in adequately describing different types of time-to-event data.     

Rank-based empirical likelihood inference on medians of k populations

We propose a nonparametric method, called rank-based empirical likelihood (REL), for making inferences on medians and cumulative distribution functions (CDFs) of k populations. The standard distribution-free approach to testing the equality of k medians requires that the k population distributions have the same shape. Our REL-ratio (RELR) test for this problem requires fewer assumptions and can effectively use the symmetry information when the distributions are symmetric. Furthermore, our RELR statistic does not require estimation of variance, and achieves asymptotic pivotalness implicitly. When the k populations have equal medians we show that the REL method produces valid inferences for the common median and CDFs of k populations. Simulation results show that the REL approach works remarkably well in finite samples. A real data example is used to illustrate the proposed REL method.     

Testing goodness-of-fit for some lifetime distributions with conventional Type-I censoring

A goodness-of-fit test procedure is proposed for some lifetime distributions when the available data are subject to Type-I censoring. The proposed method extends the test procedure of Pakyari and Balakrishnan to other lifetime distributions. The extension to Weibull and log-normal models is studied in details. The new test recovers the nominal level of significance and exhibits more power in comparison to the existing tests for several alternative distributions by means of Monte Carlo simulations. Finally, a real dataset is considered for illustrative purposes.     

On asymptotic non-normality of null distributions of mrpp statistics

Severe departures from normality occur frequently for null distributions of statistics associated with applications of mulLi-response permutation procedures (MRPP) for either small or large finite populations. This paper describes the commonly encountered situation associated with asymptotic non-normality for null distributions of MRPP statistics which does not depend on the underlying multivariate distribution. In addition, this paper establishes the existence of a non-degenerate underlying distribution for which the null distributions of MRPP statistics are asymptotically non-normal for essentially all size structure configurations. It is known that MRPP statistics are symmetric versions of a broader class of statistics, most of which are asymmetric. Because of the non-normality associated with null distributions of MRPP statistics, this paper includes necessary results for inferences based on the exact first three moments of anv statistic in this broader class (analogous to existing results for MRPP statistics).     

Discriminating between the Weibull and log-normal distributions for Type-II censored data

Log-normal and Weibull distributions are the two most popular distributions for analysing lifetime data. In this paper, we consider the problem of discriminating between the two distribution functions. It is assumed that the data are coming either from log-normal or Weibull distributions and that they are Type-II censored. We use the difference of the maximized log-likelihood functions, in discriminating between the two distribution functions. We obtain the asymptotic distribution of the discrimination statistic. It is used to determine the probability of correct selection in this discrimination process. We perform some simulation studies to observe how the asymptotic results work for different sample sizes and for different censoring proportions. It is observed that the asymptotic results work quite well even for small sizes if the censoring proportions are not very low. We further suggest a modified discrimination procedure. Two real data sets are analysed for illustrative purposes.     

Modified Tukey's Control Chart

Phase I of control analysis requires large amount of data to fit a distribution and estimate the corresponding parameters of the process under study. However, when only individual observations are available, and no a priori knowledge exists, the presence of outliers can bias the analysis. A relatively recent and successful approach to address this situation is Tukey's Control Chart (TCC), a charting method that applies the Box Plot technique to estimate the control limits. This procedure has proven to be effective for symmetric distributions. However, when skewness is present the average run length performance diminishes significantly. This article proposes a modified version of TCC to consider skewness with minimum assumptions on the underlying distribution of observations. Using theoretical results and Monte Carlo simulation, the modified TCC is tested over several distributions proving a better representation of skewed populations, even in cases when only a limited number of observations are available.     

Inferences on the common mean of several heterogeneous log-normal distributions

We consider the problem of making inferences on the common mean of several heterogeneous log-normal populations. We apply the parametric bootstrap (PB) approach and the method of variance estimate recovery (MOVER) to construct confidence intervals for the log-normal common mean. We then compare the performances of the proposed confidence intervals with the existing confidence intervals via an extensive simulation study. Simulation results show that our proposed MOVER and PB confidence intervals can be recommended generally for different sample sizes and number of populations.     

Generalized gamma parameter estimation and moment evaluation

The generalized gamma distribution includes the exponential distribution, the gamma distribution, and the Weibull distribution as special cases. It also includes the log-normal distribution in the limit as one of its parameters goes to infinity. Prentice (1974) developed an estimation method that is effective even when the underlying distribution is nearly log-normal. He reparameterized the density function so that it achieved the limiting case in a smooth fashion relative to the new parameters. He also gave formulas for the second partial derivatives of the log-density function to be used in the nearly log-normal case. His formulas included infinite summations, and he did not estimate the error in approximating these summations.

We derive approximations for the log-density function and moments of the generalized gamma distribution that are smooth in the nearly log-normal case and involve only finite summations. Absolute error bounds for these approximations are included. The approximation for the first moment is applied to the problem of estimating the parameters of a generalized gamma distribution under the constraint that the distribution have mean one. This enables the development of a correspondence between the parameters in a mean one generalized gamma distribution and certain parameters in acoustic scattering theory.     

Inferences on correlation coefficients of bivariate log-normal distributions

This article considers inference on correlation coefficients of bivariate log-normal distributions. We developed generalized confidence intervals and hypothesis tests for the correlation coefficients, and extended the results to compare two independent correlations. Simulation studies show that the suggested methods work well. Two practical examples are used to illustrate the application of the proposed methods.     

A Unified Method for Checking Compatibility and Uniqueness for Finite Discrete Conditional Distributions

Checking compatibility for two given conditional distributions and identifying the corresponding unique compatible marginal distributions are important problems in mathematical statistics, especially in Bayesian inferences. In this article, we develop a unified method to check the compatibility and uniqueness for two finite discrete conditional distributions. By formulating the compatibility problem into a system of linear equations subject to constraints, it can be reduced to a quadratic optimization problem with box constraints. We also extend the proposed method from two-dimensional cases to higher-dimensional cases. Finally, we show that our method can be easily applied to checking compatibility and uniqueness for a regression function and a conditional distribution. Several numerical examples are used to illustrate the proposed method. Some comparisons with existing methods are also presented.     

Longitudinal data analysis in the presence of informative sampling: weighted distribution or joint modelling

ABSTRACT

Weighted distributions, as an example of informative sampling, work appropriately under the missing at random mechanism since they neglect missing values and only completely observed subjects are used in the study plan. However, length-biased distributions, as a special case of weighted distributions, remove the subjects with short length deliberately, which surely meet the missing not at random mechanism. Accordingly, applying length-biased distributions jeopardizes the results by producing biased estimates. Hence, an alternate method has to be used such that the results are improved by means of valid inferences. We propose methods that are based on weighted distributions and joint modelling procedure and compare them in analysing longitudinal data. After introducing three methods in use, a set of simulation studies and analysis of two real longitudinal datasets affirm our claim.     

Robust statistics for testing mean vectors of multivariate distributions

We develop a ‘robust’ statistic T² _R, based on Tiku's (1967, 1980) MML (modified maximum likelihood) estimators of location and scale parameters, for testing an assumed meam vector of a symmetric multivariate distribution. We show that T² _R is one the whole considerably more powerful than the prominenet Hotelling T² statistics. We also develop a robust statistic T² _D for testing that two multivariate distributions (skew or symmetric) are identical; T² _D seems to be usually more powerful than nonparametric statistics. The only assumption we make is that the marginal distributions are of the type (1/σ_k)f((x-μ_k)/σ_k) and the means and variances of these marginal distributions exist.     

A Bayesian Hierarchical Model for Multiple Comparisons in Mixed Models

We propose a Bayesian hierarchical model for multiple comparisons in mixed models where the repeated measures on subjects are described with the subject random effects. The model facilitates inferences in parameterizing the successive differences of the population means, and for them, we choose independent prior distributions that are mixtures of a normal distribution and a discrete distribution with its entire mass at zero. For the other parameters, we choose conjugate or vague priors. The performance of the proposed hierarchical model is investigated in the simulated and two real data sets, and the results illustrate that the proposed hierarchical model can effectively conduct a global test and pairwise comparisons using the posterior probability that any two means are equal. A simulation study is performed to analyze the type I error rate, the familywise error rate, and the test power. The Gibbs sampler procedure is used to estimate the parameters and to calculate the posterior probabilities.     

Comparisons of various types of normality tests

Normality tests can be classified into tests based on chi-squared, moments, empirical distribution, spacings, regression and correlation and other special tests. This paper studies and compares the power of eight selected normality tests: the Shapiro–Wilk test, the Kolmogorov–Smirnov test, the Lilliefors test, the Cramer–von Mises test, the Anderson–Darling test, the D'Agostino–Pearson test, the Jarque–Bera test and chi-squared test. Power comparisons of these eight tests were obtained via the Monte Carlo simulation of sample data generated from alternative distributions that follow symmetric short-tailed, symmetric long-tailed and asymmetric distributions. Our simulation results show that for symmetric short-tailed distributions, D'Agostino and Shapiro–Wilk tests have better power. For symmetric long-tailed distributions, the power of Jarque–Bera and D'Agostino tests is quite comparable with the Shapiro–Wilk test. As for asymmetric distributions, the Shapiro–Wilk test is the most powerful test followed by the Anderson–Darling test.     

Bayesian Conditional Mean Estimation in Log‐Normal Linear Regression Models with Finite Quadratic Expected Loss

Log‐normal linear regression models are popular in many fields of research. Bayesian estimation of the conditional mean of the dependent variable is problematic as many choices of the prior for the variance (on the log‐scale) lead to posterior distributions with no finite moments. We propose a generalized inverse Gaussian prior for this variance and derive the conditions on the prior parameters that yield posterior distributions of the conditional mean of the dependent variable with finite moments up to a pre‐specified order. The conditions depend on one of the three parameters of the suggested prior; the other two have an influence on inferences for small and medium sample sizes. A second goal of this paper is to discuss how to choose these parameters according to different criteria including the optimization of frequentist properties of posterior means.     

Multivariate CUSUM and EWMA Control Charts for Skewed Populations Using Weighted Standard Deviations

This article proposes a heuristic method of constructing multivariate cumulative sum and exponentially weighted moving average control charts for skewed populations based on the weighted standard deviation method which adjusts the variance–covariance matrix of quality characteristics and approximates the probability density function using several multivariate normal distributions. These control charts, however, reduce to the conventional control charts when the underlying distribution is symmetric. In-control and out-of-control average run lengths of the proposed control charts are compared with those of the conventional control charts for multivariate lognormal and Weibull distributions. Simulation results show that considerable improvements over the standard method can be achieved when the underlying distribution is skewed.