Testing homogeneity of control and treatment populations — local optimality and related issues

Durairajan and Raman (1996 a, b) studied the robustness of Locally most powerful invariant (LMPI) tests for compound normal model in control and treatment populations. In the present paper, the Locally most powerful (LMP) tests are constructed for no contamination in normal mixture model through testing the parameter of mixture of distributions and the mixing proportion. The expected performance of LMP tests are compared using Efron's Statistical Curvature on the lines of Sen Gupta and Pal (1991). The Locally most powerful similar (LMPS) tests for the equality of control and treatment populations in the presence of nuisance parameters are also constructed. Further, the null and non-null distributions of the test statistics are derived and some power computations are made. Received: September 1, 1999; revised version: August 31, 2000     

Optimal Tests for No Contamination in Reliability Models

Inferences on mixtures of probability distributions, in general,and of life distributions, in particular, are receiving considerableimportance in recent years. The likelihood ratio procedure oftesting for the null hypothesis of no contamination is oftenvery cumbersome and lacks its usual asymptotic properties. Recently,SenGupta (1991) has introduced the notion of an `L-optimal' testfor such testing problems. The idea is to recast the originalseveral parametric hypotheses representation of the null hypothesisin terms of only a single hypothesis involving an appropriatelychosen parametric function. This approach is shown to be bothmathematically elegant and operationally simple for a quite generalclass of mixture distributions which contains, in particular,all mixtures of the one-parameter exponential family and alsoa very rich subclass of mixtures useful in life-testing and reliabilityanalysis. It is also illustrated through two examples—onebased on real-life data and the other on a simulated sample.     

SB-robustness of directional mean for circular distributions

In this paper we study the robustness of the directional mean (a.k.a. circular mean) for different families of circular distributions. We show that the directional mean is robust in the sense of finite standardized gross error sensitivity (SB-robust) for the following families: (1) mixture of two circular normal distributions, (2) mixture of wrapped normal and circular normal distributions and (3) mixture of two wrapped normal distributions. We also show that the directional mean is not SB-robust for the family of all circular normal distributions with varying concentration parameter. We define the circular trimmed mean and prove that it is SB-robust for this family. In general the property of SB-robustness of an estimator at a family of probability distributions is dependent on the choice of the dispersion measure. We introduce the concept of equivalent dispersion measures and prove that if an estimator is SB-robust for one dispersion measure then it is SB-robust for all equivalent dispersion measures. Three different dispersion measures for circular distributions are considered and their equivalence studied.     

Reliability inference for a general lower-truncated family of distributions under records

Based on record values, point and interval estimators are proposed in this paper for the parameters of a general lower-truncated family of distributions. Maximum likelihood and bias-corrected estimators are obtained for unknown model parameters. Based on a sufficient and complete statistic, the bias-corrected estimator is also shown to be uniformly minimum variance unbiased estimator. Different exact confidence intervals and exact confidence regions are constructed for the both model and truncated parameters, and other confidence interval estimates based on asymptotic distribution theory and bootstrap approaches are obtained as well. Finally, two real-life examples and a numerical study are presented to illustrate the performance of our methods.     

Locally optimal test for simultaneously testing the mean and variance of a normal distribution

A generalization of the locally most powerful unbiased (LMPU) test for the single parameter case to the k-parameter case was proposed by SenGupta and Vermeire (1986). In particular we defined a locally most mean power unbiased (LMMPU) test based on the mean curvature of the power hypersurface. Compared to the type C tests of Neyman and Pearson and the type D tests (Isaacson, 1951), LMMPU tests possess better theoretical properties and enjoy ease of construction of critical regions. In this paper we present an interesting example of a two-parameter univariate normal population for which Isaacson (1951, p. 233) was unsuccessful in finding a type D test. For the case of one observation, we prove that no Type D region exists but the LMMPU test is obtained - it is an example of a test with singular Hessian matrix for its power but is nevertheless a strictly locally unbiased (LU) test.     

Bayesian analysis of circular distributions based on non-negative trigonometric sums

ABSTRACT

Fernández-Durán [Circular distributions based on nonnegative trigonometric sums. Biometrics. 2004;60:499–503] developed a new family of circular distributions based on non-negative trigonometric sums that is suitable for modelling data sets that present skewness and/or multimodality. In this paper, a Bayesian approach to deriving estimates of the unknown parameters of this family of distributions is presented. Because the parameter space is the surface of a hypersphere and the dimension of the hypersphere is an unknown parameter of the distribution, the Bayesian inference must be based on transdimensional Markov Chain Monte Carlo (MCMC) algorithms to obtain samples from the high-dimensional posterior distribution. The MCMC algorithm explores the parameter space by moving along great circles on the surface of the hypersphere. The methodology is illustrated with real and simulated data sets.     

Testing for two states in a hidden Markov model

The authors consider hidden Markov models (HMMs) whose latent process has m ≥ 2 states and whose state‐dependent distributions arise from a general one‐parameter family. They propose a test of the hypothesis m = 2. Their procedure is an extension to HMMs of the modified likelihood ratio statistic proposed by Chen, Chen & Kalbfleisch (2004) for testing two states in a finite mixture. The authors determine the asymptotic distribution of their test under the hypothesis m = 2 and investigate its finite‐sample properties in a simulation study. Their test is based on inference for the marginal mixture distribution of the HMM. In order to illustrate the additional difficulties due to the dependence structure of the HMM, they show how to test general regular hypotheses on the marginal mixture of HMMs via a quasi‐modified likelihood ratio. They also discuss two applications.     

Robust Bayesian hypothesis testing in the presence of nuisance parameters

Robust Bayesian testing of point null hypotheses is considered for problems involving the presence of nuisance parameters. The robust Bayesian approach seeks answers that hold for a range of prior distributions. Three techniques for handling the nuisance parameter are studied and compared. They are (i) utilize a noninformative prior to integrate out the nuisance parameter; (ii) utilize a test statistic whose distribution does not depend on the nuisance parameter; and (iii) use a class of prior distributions for the nuisance parameter. These approaches are studied in two examples, the univariate normal model with unknown mean and variance, and a multivariate normal example.     

A general class of non parametric tests for comparing scale parameters

In this paper, a general class of non parametric tests is proposed for the two-sample scale problem. Testing of the scale parameter is very useful in real-life situations commonly faced in engineering, trade, cultivation, industries, medicine, etc. In all these fields, one will prefer the method that gives more consistent results. Thus, it is worthwhile to test the equality of scale parameters. The distribution of the proposed test is established. To assess the performance of the proposed test, the asymptotic efficacies are studied for some underlying distributions and the results are interpreted with useful information. To see the working of the proposed test, an illustrative example for the real-life data set is provided. The simulation study is also carried out to find the asymptotic power of the proposed test. An extension of the general class of tests to the multiple-sample problem is also discussed.     

Estimation of finite mixtures with symmetric components

It may sometimes be clear from background knowledge that a population under investigation proportionally consists of a known number of subpopulations, whose distributions belong to the same, yet unknown, family. While a parametric family is commonly used in practice, one can also consider some nonparametric families to avoid distributional misspecification. In this article, we propose a solution using a mixture-based nonparametric family for the component distribution in a finite mixture model as opposed to some recent research that utilizes a kernel-based approach. In particular, we present a semiparametric maximum likelihood estimation procedure for the model parameters and tackle the bandwidth parameter selection problem via some popular means for model selection. Empirical comparisons through simulation studies and three real data sets suggest that estimators based on our mixture-based approach are more efficient than those based on the kernel-based approach, in terms of both parameter estimation and overall density estimation.     

Estimation of a scale parameter in mixture models with unknown location

Estimation of the scale parameter in mixture models with unknown location is considered under Stein's loss. Under certain conditions, the inadmissibility of the “usual” estimator is established by exhibiting better estimators. In addition, robust improvements are found for a specified submodel of the original model. The results are applied to mixtures of normal distributions and mixtures of exponential distributions. Improved estimators of the variance of a normal distribution are shown to be robust under any scale mixture of normals having variance greater than the variance of that normal distribution. In particular, Stein's (Ann. Inst. Statist. Math. 16 (1964) 155) and Brewster's and Zidek's (Ann. Statist. 2 (1974) 21) estimators obtained under the normal model are robust under the t model, for arbitrary degrees of freedom, and under the double-exponential model. Improved estimators for the variance of a t distribution with unknown and arbitrary degrees of freedom are also given. In addition, improved estimators for the scale parameter of the multivariate Lomax distribution (which arises as a certain mixture of exponential distributions) are derived and the robustness of Zidek's (Ann. Statist. 1 (1973) 264) and Brewster's (Ann. Statist. 2 (1974) 553) estimators of the scale parameter of an exponential distribution is established under a class of modified Lomax distributions.     

Bayesian inference of the inverse Weibull mixture distribution using type-I censoring

A large number of models have been derived from the two-parameter Weibull distribution including the inverse Weibull (IW) model which is found suitable for modeling the complex failure data set. In this paper, we present the Bayesian inference for the mixture of two IW models. For this purpose, the Bayes estimates of the parameters of the mixture model along with their posterior risks using informative as well as the non-informative prior are obtained. These estimates have been attained considering two cases: (a) when the shape parameter is known and (b) when all parameters are unknown. For the former case, Bayes estimates are obtained under three loss functions while for the latter case only the squared error loss function is used. Simulation study is carried out in order to explore numerical aspects of the proposed Bayes estimators. A real-life data set is also presented for both cases, and parameters obtained under case when shape parameter is known are tested through testing of hypothesis procedure.     

Simultaneous tests for homogeneity of two zero-inflated (beta) populations

ABSTRACT

Motivated by an example in marine science, we use Fisher’s method to combine independent likelihood ratio tests (LRTs) and asymptotic independent score tests to assess the equivalence of two zero-inflated Beta populations (mixture distributions with three parameters). For each test, test statistics for the three individual parameters are combined into a single statistic to address the overall difference between the two populations. We also develop non parametric and semiparametric permutation-based tests for simultaneously comparing two or three features of unknown populations. Simulations show that the likelihood-based tests perform well for large sample sizes and that the statistics based on combining LRT statistics outperforms the ones based on combining score test statistics. The permutation-based tests have overall better performance in terms of both power and type I error rate. Our methods are easy to implement and computationally efficient, and can be expanded to more than two populations and to other multiple parameter families. The permutation tests are entirely generic and can be useful in various applications dealing with zero (or other) inflation.     

Testing the Homogeneity of Two Survival Functions Against a Mixture Alternative Based on Censored Data

When the survival distribution in a treatment group is a mixture of two distributions of the same family, traditional parametric methods that ignore the existence of mixture components or the nonparametric methods may not be very powerful. We develop a modified likelihood ratio test (MLRT) for testing homogeneity in a two sample problem with censored data and compare the actual type I error and power of the MLRT with that nonparametric log-rank test and parametric test through Monte-Carlo simulations. The proposed test is also applied to analyze data from a clinical trial on early breast cancer.     

A new generalized Lindley distribution

A new generalized Lindley distribution, based on weighted mixture of two gamma distributions, is proposed. This model includes the Lindley, gamma and exponential distributions as and other forms of Lindley distributions as special cases. Lindley distribution based on two gamma with two consecutive shape parameter is investigated in some details. Statistical and reliability properties of this model are derived. The size-biased, the length-biased and Lorenze curve are established. Estimation of the underlying parameters via the moment method and maximum likelihood has been investigated and their values are simulated. Finally, fitting this model to a set of real-life data is discussed.     

Score test for homogeneity of dispersion in generalized Poisson mixed models with excess zeros

In many applications, the clustered count data often contain excess zeros and the zero-inflated generalized Poisson mixed (ZIGPM) regression model may be suitable. However, dispersion in ZIGPM is often treated as fixed unknown parameter, and this assumption may be not appropriate in some situations. In this article, a score test for homogeneity of dispersion parameter in ZIGPM regression model is developed and corresponding test statistic is obtained. Sampling distribution and power of the score test statistic are investigated through Monte Carlo simulation. Finally, results from a biological example illustrate the usefulness of the diagnostic statistic.     

Construction of bivariate and multivariate weighted distributions via conditioning

Weighted distributions (univariate and bivariate) have received widespread attention over the last two decades because of their flexibility for analyzing skewed data. In this article, we propose an alternative method to construct a new family of bivariate and multivariate weighted distributions. For illustrative purposes, some examples of the proposed method are presented. Several structural properties of the bivariate weighted distributions including marginal distributions together with distributions of the minimum and maximum, evaluation of the reliability parameter, and verification of total positivity of order two are also presented. In addition, we provide some multivariate extensions of the proposed models. A real-life data set is used to show the applicability of these bivariate weighted distributions.     

Change-point problems for the von Mises distribution

A generalized likelihood ratio procedure and a Bayes procedure are considered for change-point problems for the mean direction of the von Mises distribution, both when the concentration parameter is known and when it is unknown. These tests are based on sample resultant lengths. Tables that list critical values of these test statistics are provided. These tests are shown to be valid even when the data come from other similar unimodal circular distributions. Some empirical studies of powers of these test procedures are also incorporated.     

On a Relation Between a Family of Distributions Attaining the Bhattacharyya Bound and That of Linear Combinations of the Distributions from an Exponential Family

Abstract

An unbiased estimation problem of a function g(θ) of a real parameter is considered. A relation between a family of distributions for which an unbiased estimator of a function g(θ) attains the general order Bhattacharyya lower bound and that of linear combinations of the distributions from an exponential family is discussed. An example on a family of distributions involving an exponential and a double exponential distributions with a scale parameter is given. An example on a normal distribution with a location parameter is also given.     

Discrimination among bivariate beta-generated distributions

Some properties of the general families of bivariate distributions generated by beta dependent random variables are derived and discussed here. Some classic measures of dependence and information are derived, and their behaviours and properties are discussed as well. Finally, a discrimination procedure within this general family of bivariate distributions is proposed based on Shannon entropy. A real-life example is presented to illustrate the model as well as the inferential results developed here.