On testing hypotheses with divergence statistics

Using divergence measures based on entropy functions, a procedure to test statistical hypotheses is proposed. Replacing the parameters by suitable estimators in the expresion of the divergence measure, the test statistics are obtained. Asymptotic distributions for these statistics are given in several cases when maximum likelihood estimators are considered, so they can be used to construct confidence intervals and to test statistical hypotheses based on one or more samples. These results can also be applied to multinomial populations. Tests of goodness of fit and tests of homogeneity can be constructed.     

Robust confidence intervals for log-location-scale models with right censored data

In this paper, we combine empirical likelihood and estimating functions for censored data to obtain robust confidence regions for the parameters and more generally for functions of the parameters of distributions used in lifetime data analysis. The proposed method works with type I, type II or randomly censored data. It is illustrated by considering inference for log-location-scale models. In particular, we focus on the log-normal and the Weibull models and we tackle the problem of constructing robust confidence regions (or intervals) for the parameters of the model, as well as for quantiles and values of the survival function. The usefulness of the method is demonstrated through a Monte Carlo study and by examples on two lifetime data sets.     

Some aspects of the reparametrization of statistical models

Definitions are given for orthogonal parameters in the context of Bayesian inference and likelihood inference. The exact orthogonalizing transformations are derived for both cases, and the connection between the two settings is made precise. These parametrizations simplify the interpretation of likelihood functions and posterior distributions. Further, they make numerical maximization and integration procedures easier to apply. Several applications are studied.     

Estimation in two-parameter exponential distributions

Exponential distributions are used extensively in the field of life-testing. Estimation of parameters is revisited in two-parameter exponential distributions. A comparison study between the maximum likelihood method, the unbiased estimates which are linear functions of the maximum likelihood method, the method of product spacings, and the method of quantile estimates are presented. Finally, a simulation study is given to demonstrate the small sample properties     

Exact confidence regions for scale and location parameters based on sample quartiles

In this paper we present relatively simple (ruler, paper, and pencil) nonparametric procedures for constructing joint confidence regions for (i) the median and the inner quartile range for the symmetric one-sample problem and (ii) the shift and ratio of scale parameters for the two-sample case. Both procedures are functions of the sample quartiles and have exact confidence levels when the populations are continuous. The one-sample case requires symmetry of first and third quartiles about the median.

The confidence regions we propose are always convex, nested for decreasing confidence levels and are compact for reasonably large sample sizes. Both exact small sample and approximate large sample distributions are given.     

Discriminating Between the Log-Normal and Log-Logistic Distributions

Log-normal and log-logistic distributions are often used to analyze lifetime data. For certain ranges of the parameters, the shape of the probability density functions or the hazard functions can be very similar in nature. It might be very difficult to discriminate between the two distribution functions. In this article, we consider the discrimination procedure between the two distribution functions. We use the ratio of maximized likelihood for discrimination purposes. The asymptotic properties of the proposed criterion are investigated. It is observed that the asymptotic distributions are independent of the unknown parameters. The asymptotic distributions are used to determine the minimum sample size needed to discriminate between these two distribution functions for a user specified probability of correct selection. We perform some simulation experiments to see how the asymptotic results work for small sizes. For illustrative purpose, two data sets are analyzed.     

Subroutines for computing minimum variance unbiased estimators of functions of the parameters in the normal and gamma distributions

In this paper subroutines are given which calculate the uniformly minimum variance unbiased estimators (UMVUE’ s) of a broad class of functions of the parameters of the normal and gamma distributions. These subroutines employ the new expressions for the UMVUE’ s given recently by Gray, Watkins, and Schucany (1973), Woodward and Gray (1975), and Gray, Schucany, and Woodward (1976). In order to employ the subroutines here the user need only be able to provide a FORTRAN function subprogram to calculate derivatives of the function, either analytically or numerically.     

Lifetime distributions motivated by series and parallel structures

According to Ross, any system can be represented either as a series arrangement of parallel structures or as a parallel arrangement of series structures. Motivated by this, we propose new three-parameter lifetime distributions by compounding geometric, power series, and exponential distributions. The distributions can allow for decreasing, increasing, bathtub-shaped, and upside down bathtub-shaped hazard rates. A mathematical treatment of the new distributions is provided including expressions for their density functions, Shannon and Rényi entropies, mean residual life functions, hazard rate functions, quantiles, and moments. The method of maximum likelihood is used for estimating parameters. Five of the new distributions are studied in detail. Finally, two illustrative data examples and a sensitivity analysis are presented.     

Minimum variance unbiased estimation for modified power series distribution

In this paper we study the problem of finding the minimum variance unbiased (MVU) estimators of the functions of the para-meters of the modified power series distributions (MPSD). A theorem giving the necessary and sufficient conditions for the existence of the MVU estimators has been proved. Also, the estimators for a number of estimable functions of a parameter are obtained. Two other theorems dealing with the MVU estimation of the left truncated MPSD with unknown truncation point are also given. The particular case of the Lagrangian Poisson, the Lagrangian binomial and the Borel-Tanner distributions are considered and tables are also provided for the MVU estimators for some functions of the parameters. The variances of the estimators are also given for some cases.     

Generalized additive models for location, scale and shape

Summary.  A general class of statistical models for a univariate response variable is presented which we call the generalized additive model for location, scale and shape (GAMLSS). The model assumes independent observations of the response variable y given the parameters, the explanatory variables and the values of the random effects. The distribution for the response variable in the GAMLSS can be selected from a very general family of distributions including highly skew or kurtotic continuous and discrete distributions. The systematic part of the model is expanded to allow modelling not only of the mean (or location) but also of the other parameters of the distribution of y , as parametric and/or additive nonparametric (smooth) functions of explanatory variables and/or random-effects terms. Maximum (penalized) likelihood estimation is used to fit the (non)parametric models. A Newton–Raphson or Fisher scoring algorithm is used to maximize the (penalized) likelihood. The additive terms in the model are fitted by using a backfitting algorithm. Censored data are easily incorporated into the framework. Five data sets from different fields of application are analysed to emphasize the generality of the GAMLSS class of models.     

Evaluation of Tweedie exponential dispersion model densities by Fourier inversion

The Tweedie family of distributions is a family of exponential dispersion models with power variance functions V(μ)=μ ^p for . These distributions do not generally have density functions that can be written in closed form. However, they have simple moment generating functions, so the densities can be evaluated numerically by Fourier inversion of the characteristic functions. This paper develops numerical methods to make this inversion fast and accurate. Acceleration techniques are used to handle oscillating integrands. A range of analytic results are used to ensure convergent computations and to reduce the complexity of the parameter space. The Fourier inversion method is compared to a series evaluation method and the two methods are found to be complementary in that they perform well in different regions of the parameter space.     

Preliminary test estimation of the parameters of exponential and Pareto distributions for censored samples

This paper deals with the preliminary test estimation (PTE) of the parameters of the exponential and Pareto distributions in censored samples. The biases, risk functions, efficiency tables and the graphs for the relative efficiency for the proposed estimators for the parameters of the exponential and Pareto distributions are given. We find that the proposed estimators dominate the corresponding unrestricted (usual) estimators in the neighborhood of null hypothesis. The range of the parameters for which the proposed estimators dominate the corresponding usual estimators for different sample sizes and level of significance are given. The findings of the paper will be useful for the practitioners who are dealing with the censored samples in life testing experiments.     

A new class of probability distributions via cosine and sine functions with applications

ABSTRACT

In this paper, we introduce a new class of (probability) distributions, based on a cosine-sine transformation, obtained by compounding a baseline distribution with cosine and sine functions. Some of its properties are explored. A special focus is given to a particular cosine-sine transformation using the exponential distribution as baseline. Estimations of parameters of a particular cosine-sine exponential distribution are performed via the maximum likelihood estimation method. A simulation study investigates the performances of these estimates. Applications are given for four real data sets, showing a better fit in comparison to some existing distributions based on some goodness-of-fit tests.     

On an asymptotic relative efficiency concept based on expected volumes of confidence regions

ABSTRACT

The paper deals with an asymptotic relative efficiency concept for confidence regions of multidimensional parameters that is based on the expected volumes of the confidence regions. Under standard conditions the asymptotic relative efficiencies of confidence regions are seen to be certain powers of the ratio of the limits of the expected volumes. These limits are explicitly derived for confidence regions associated with certain plugin estimators, likelihood ratio tests and Wald tests. Under regularity conditions, the asymptotic relative efficiency of each of these procedures with respect to each one of its competitors is equal to 1. The results are applied to multivariate normal distributions and multinomial distributions in a fairly general setting.     

A new lifetime distribution for a series-parallel system: properties,applications and estimations under progressive type-II censoring

A compound class of zero truncated Poisson and lifetime distributions is introduced. A specialization is paved to a new three-parameter distribution, called doubly Poisson-exponential distribution, which may represent the lifetime of units connected in a series-parallel system. The new distribution can be obtained by compounding two zero truncated Poisson distributions with an exponential distribution. Among its motivations is that its hazard rate function can take different shapes such as decreasing, increasing and upside-down bathtub depending on the values of its parameters. Several properties of the new distribution are discussed. Based on progressive type-II censoring, six estimation methods [maximum likelihood, moments, least squares, weighted least squares and Bayes (under linear-exponential and general entropy loss functions) estimations] are used to estimate the involved parameters. The performance of these methods is investigated through a simulation study. The Bayes estimates are obtained using Markov chain Monte Carlo algorithm. In addition, confidence intervals, symmetric credible intervals and highest posterior density credible intervals of the parameters are obtained. Finally, an application to a real data set is used to compare the new distribution with other five distributions.     

Estimation of Hazard, Density and Survivor Functions for Randomly Censored Data

Maximum likelihood estimation and goodness-of-fit techniques are used within a competing risks framework to obtain maximum likelihood estimates of hazard, density, and survivor functions for randomly right-censored variables. Goodness-of- fit techniques are used to fit distributions to the crude lifetimes, which are used to obtain an estimate of the hazard function, which, in turn, is used to construct the survivor and density functions of the net lifetime of the variable of interest. If only one of the crude lifetimes can be adequately characterized by a parametric model, then semi-parametric estimates may be obtained using a maximum likelihood estimate of one crude lifetime and the empirical distribution function of the other. Simulation studies show that the survivor function estimates from crude lifetimes compare favourably with those given by the product-limit estimator when crude lifetimes are chosen correctly. Other advantages are discussed.     

Discriminating among Weibull,log-normal,and log-logistic distributions

In this article, we consider the problem of the model selection/discrimination among three different positively skewed lifetime distributions. All these three distributions, namely; the Weibull, log-normal, and log-logistic, have been used quite effectively to analyze positively skewed lifetime data. In this article, we have used three different methods to discriminate among these three distributions. We have used the maximized likelihood method to choose the correct model and computed the asymptotic probability of correct selection. We have further obtained the Fisher information matrices of these three different distributions and compare them for complete and censored observations. These measures can be used to discriminate among these three distributions. We have also proposed to use the Kolmogorov–Smirnov distance to choose the correct model. Extensive simulations have been performed to compare the performances of the three different methods. It is observed that each method performs better than the other two for some distributions and for certain range of parameters. Further, the loss of information due to censoring are compared for these three distributions. The analysis of a real dataset has been performed for illustrative purposes.     

Construction of non-exchangeable bivariate distribution functions

A method is given for constructing bivariate distributions functions by means of the copula functions, and, hence, it is used for obtaining distribution functions that can describe the behaviour of non–exchangeable random vectors.     

A Marshall–Olkin family of heavy-tailed distributions which includes the lognormal one

Abstract

A new class of heavy-tailed distribution functions,, containing the lognormal distribution as a particular case is introduced. The class thus obtained depends on a set of three parameters, incorporating an additional distribution to the classical lognormal one. This new class of heavy-tailed distribution is presented as an alternative to other useful heavy-tailed distributions, such as the lognormal, Weibull, and Pareto distributions. The density and distribution functions of this new class are given by a closed expression which allows us to easily compute probabilities, quantiles, moments, and related measurements. Finally, some applications are shown as examples.     

Multivariate distributions defined in terms of contours

The paper deals with the problem of using contours as the basis for defining probability distributions. First, the most general probability densities with given contours are obtained and the particular cases of circular and elliptical contours are dealt with. It is shown that the so-called elliptically contoured distributions do not include all possible cases. Next, the case of contours defined by polar coordinates is analyzed including its simulation and parameter estimation. Finally, the case of cumulative distribution functions with given contours is discussed. Several examples are used for illustrative purposes.