Inflated beta distributions

This paper considers the issue of modeling fractional data observed on [0,1), (0,1] or [0,1]. Mixed continuous-discrete distributions are proposed. The beta distribution is used to describe the continuous component of the model since its density can have quite different shapes depending on the values of the two parameters that index the distribution. Properties of the proposed distributions are examined. Also, estimation based on maximum likelihood and conditional moments is discussed. Finally, practical applications that employ real data are presented.     

A moment-generating function with application to the birnbaum-saunders distribution

The Birnbaum-Saunders distribution is a fatigue life distribution that was derived from a model assuming that failure is due to the development and growth of a dominant crack. This distribution has been shown to be applicable not only for fatigue analysis but also in other areas of engineering science. Because of its increasing use, it would be desirable to obtain expressions for the expected value of different powers of this distribution.

In this article, the moment-generating function for the sinh-normal distribution is derived. It is shown that this moment-generating function can be used to obtain both integer and fractional moments for the Birnbaum-Saunders distribution. Thus it is now possible to obtain an expression for the expected value of the square root of a Birnbaum-Saunders random variable. A general expression for integer noncentral moments for the Birnbaum-Saunders distribution is derived using the moment-generating function of the sinh-normal distribution. Also included is an approximation of the moment-generating function that can be used fcx small values of the shape parameter.     

A new approximation to the distribution function of the studentized maximum root

The Studentized maximum root (SMR) distribution is useful for constructing simultaneous confidence intervals around product interaction contrasts in replicated two-way ANOVA. A three-moment approximation to the SMR distribution is proposed. The approximation requires the first three moments of the maximum root of a central Wishart matrix. These values are obtained by means of numerical integration. The accuracy of the approximation is compared to the accuracy of a two-moment approximation for selected two-way table sizes. Both approximations are reasonably accurate. The three-moment approximation is generally superior.     

Estimation for the Parameters of Generalized Half-normal Distribution Based on Progressive Type-I Interval Censoring

Based on progressively Type-I interval censored sample, the problem of estimating unknown parameters of a two parameter generalized half-normal(GHN) distribution is considered. Different methods of estimation are discussed. They include the maximum likelihood estimation, midpoint approximation method, approximate maximum likelihood estimation, method of moments, and estimation based on probability plot. Several Bayesian estimates with respect to different symmetric and asymmetric loss functions such as squared error, LINEX, and general entropy is calculated. The Lindley’s approximation method is applied to determine Bayesian estimates. Monte Carlo simulations are performed to compare the performances of the different methods. Finally, analysis is also carried out for a real dataset.     

Comparison of improved regression estimators with and without moments

For estimating the coefficients in a linear regression model, the double k–class estimators are considered and the small disturbance asymptotic approximation for their density function is obtained. Then employing the criterion of concentration probability around the true parameter values, a comparison is made between the estimators possessing finite moments and the estimators having no finite moments.     

On an application of the diffusion approximation of the generalized birth and death process for stochastic population projection

Assuming that both birth and death rates are density and time dependent, a diffusion approximation of the generalized birth and death process has been considered in this paper to obtain a suitable stochastic population model describing the population size and its moments. A simple method of estimating the parameters of the model Is discussed. The predictions of the expected size of the population, and the variance are made and compared with the corresponding census figures as well as with another deterministic projection series made for the corresponding period.     

DENSITY FUNCTION APPROXIMATION WITH PARTICULAR REFERENCE TO THE DISTRIBUTION OF OVERLAP OF CIRCLES

Given that two circles overlap, the area in common is a function of the distance between their centres. This paper adopts a suitable random distribution for the intercentre distance and then derives the distribution of the area of overlap. An approximation is sought for the density function using a criterion which enables bounds to be placed on the difference between the moments of the density function and those of the approximation. This is an approach of general applicability. The importance of matching the end-point behaviour of the density and the approximation is stressed. It is shown that the distribution of the area of overlap may be well approximated by a mixture of beta distributions in which the parameters change smoothly with the ratio of radii.     

A note on maximum entropy spectrum estimation

This article studies the maximum entropy spectrum estimation. After a bnei discussion on iiow co select ciic appropriate constraints and tiie objec¬tive functions, we decide to choose the constraints containing only the first four sample moments and, consequently, to employ the second order spectral entropy as the objective function. The resulting (maximum entropy) spectral estimate is the power spectral density of an ARMA sequence. Examples for comparing our proposal with the traditional maximum entropy spectral estimate follow at the end.     

A new compound class of log-logistic Weibull–Poisson distribution: model,properties and applications

A new class of distributions called the log-logistic Weibull–Poisson distribution is introduced and its properties are explored. This new distribution represents a more flexible model for lifetime data. Some statistical properties of the proposed distribution including the expansion of the density function, quantile function, hazard and reverse hazard functions, moments, conditional moments, moment generating function, skewness and kurtosis are presented. Mean deviations, Bonferroni and Lorenz curves, Rényi entropy and distribution of the order statistics are derived. Maximum likelihood estimation technique is used to estimate the model parameters. A simulation study is conducted to examine the bias, mean square error of the maximum likelihood estimators and width of the confidence interval for each parameter and finally applications of the model to real data sets are presented to illustrate the usefulness of the proposed distribution.     

A Laguerre polynomial approximation for a goodness-of-fit test for exponential distribution based on progressively censored data

This paper proposes an approximation to the distribution of a goodness-of-fit statistic proposed recently by Balakrishnan et al. [Balakrishnan, N., Ng, H.K.T. and Kannan, N., 2002, A test of exponentiality based on spacings for progressively Type-II censored data. In: C. Huber-Carol et al. (Eds.), Goodness-of-Fit Tests and Model Validity (Boston: Birkhäuser), pp. 89–111.] for testing exponentiality based on progressively Type-II right censored data. The moments of this statistic can be easily calculated, but its distribution is not known in an explicit form. We first obtain the exact moments of the statistic using Basu's theorem and then the density approximants based on these exact moments of the statistic, expressed in terms of Laguerre polynomials, are proposed. A comparative study of the proposed approximation to the exact critical values, computed by Balakrishnan and Lin [Balakrishnan, N. and Lin, C.T., 2003, On the distribution of a test for exponentiality based on progressively Type-II right censored spacings. Journal of Statistical Computation and Simulation, 73 (4), 277–283.], is carried out. This reveals that the proposed approximation is very accurate.     

Simulation-assisted saddlepoint approximation

A general saddlepoint/Monte Carlo method to approximate (conditional) multivariate probabilities is presented. This method requires a tractable joint moment generating function (m.g.f.), but does not require a tractable distribution or density. The method is easy to program and has a third-order accuracy with respect to increasing sample size in contrast to standard asymptotic approximations which are typically only accurate to the first order.

The method is most easily described in the context of a continuous regular exponential family. Here, inferences can be formulated as probabilities with respect to the joint density of the sufficient statistics or the conditional density of some sufficient statistics given the others. Analytical expressions for these densities are not generally available, and it is often not possible to simulate exactly from the conditional distributions to obtain a direct Monte Carlo approximation of the required integral. A solution to the first of these problems is to replace the intractable density by a highly accurate saddlepoint approximation. The second problem can be addressed via importance sampling, that is, an indirect Monte Carlo approximation involving simulation from a crude approximation to the true density. Asymptotic normality of the sufficient statistics suggests an obvious candidate for an importance distribution.

The more general problem considers the computation of a joint probability for a subvector of random T, given its complementary subvector, when its distribution is intractable, but its joint m.g.f. is computable. For such settings, the distribution may be tilted, maintaining T as the sufficient statistic. Within this tilted family, the computation of such multivariate probabilities proceeds as described for the exponential family setting.     

Pareto Poisson–Lindley distribution with applications

A new lifetime distribution is introduced based on compounding Pareto and Poisson–Lindley distributions. Several statistical properties of the distribution are established, including behavior of the probability density function and the failure rate function, heavy- and long-right tailedness, moments, the Laplace transform, quantiles, order statistics, moments of residual lifetime, conditional moments, conditional moment generating function, stress–strength parameter, Rényi entropy and Song's measure. We get maximum-likelihood estimators of the distribution parameters and investigate the asymptotic distribution of the estimators via Fisher's information matrix. Applications of the distribution using three real data sets are presented and it is shown that the distribution fits better than other related distributions in practical uses.     

The McDonald arcsine distribution: a new model to proportional data

We propose a new three-parameter continuous model called the McDonald arcsine distribution, which is a very competitive model to the beta, beta type I and Kumaraswamy distributions for modelling rates and proportions. We provide a mathematical treatment of the new distribution including explicit expressions for the density function, moments, generating and quantile functions, mean deviations, two probability measures based on the Bonferroni and Lorenz curves, Shannon entropy, Rényi entropy and cumulative residual entropy. Maximum likelihood is used to estimate the model parameters and the expected information matrix is determined. An application of the proposed model to real data shows that it can give consistently a better fit than other important statistical models.     

On approximating the distribution of indefinite quadratic forms

This paper provides a simple methodology for approximating the distribution of indefinite quadratic forms in normal random variables. It is shown that the density function of a positive definite quadratic form can be approximated in terms of the product of a gamma density function and a polynomial. An extension which makes use of a generalized gamma density function is also considered. Such representations are based on the moments of a quadratic form, which can be determined from its cumulants by means of a recursive formula. After expressing an indefinite quadratic form as the difference of two positive definite quadratic forms, one can obtain an approximation to its density function by means of the transformation of variable technique. An explicit representation of the resulting density approximant is given in terms of a degenerate hypergeometric function. An easily implementable algorithm is provided. The proposed approximants produce very accurate percentiles over the entire range of the distribution. Several numerical examples illustrate the results. In particular, the methodology is applied to the Durbin–Watson statistic which is expressible as the ratio of two quadratic forms in normal random variables. Quadratic forms being ubiquitous in statistics, the approximating technique introduced herewith has numerous potential applications. Some relevant computational considerations are also discussed.     

General properties for the beta extended half-normal model

We formulate and study a four-parameter lifetime model called the beta extended half-normal distribution. This model includes as sub-models the exponential, extended half-normal and half-normal distributions. We derive expansions for the new density function which do not depend on complicated functions. We obtain explicit expressions for the moments and incomplete moments, generating function, mean deviations, Bonferroni and Lorenz curves and Rényi entropy. In addition, the model parameters are estimated by maximum likelihood. We provide the observed information matrix. The new model is modified to cope with possible long-term survivors in the data. The usefulness of the new distribution is shown by means of two real data sets.     

A Proposal for a New Bound for Discrete Distributions

In this work we re-examine some classical bounds for non negative integer-valued random variables by means of information theoretic or maxentropic techniques using fractional moments as constraints. The proposed new bound, no more analytically expressible in terms of moments or moment generating function (mgf), is built by mixing classical bounds and the Maximum Entropy (ME) approximant of the underlying distribution; such a new bound is able to exploit optimally all the information content provided by the sequence of given moments or by the mgf. Particular care will be devoted to obtain fractional moments from the available information given in terms of integer moments and/or moment generating function. Numerical examples show clearly that the bound improvement involving the ME approximant based on fractional moments is not trivial.     

Moments of nonparametric probability density functions of entropy estimators applied to testing the inverse Gaussian distribution

The inverse Gaussian (IG) distribution is widely used to model data and then it is important to develop efficient goodness of fit tests for this distribution. In this article, we introduce some new test statistics for examining the IG goodness of fit based on correcting moments of nonparametric probability density functions of entropy estimators. These tests are consistent against all alternatives. Critical points and power of the tests are explored by simulation. We show that the proposed tests are more powerful than competitor tests. Finally, the proposed tests are illustrated by real data examples.     

Moment density estimation for positive random variables

An unknown moment-determinate cumulative distribution function or its density function can be recovered from corresponding moments and estimated from the empirical moments. This method of estimating an unknown density is natural in certain inverse estimation models like multiplicative censoring or biased sampling when the moments of unobserved distribution can be estimated via the transformed moments of the observed distribution. In this paper, we introduce a new nonparametric estimator of a probability density function defined on the positive real line, motivated by the above. Some fundamental properties of proposed estimator are studied. The comparison with traditional kernel density estimator is discussed.     

On approximating the central and noncentral multivariate gamma distributions

In this paper a finite series approximation involving Laguerre polynomials is derived for central and noncentral multivariate gamma distributions. It is shown that if one approximates the density of any k nonnegative continuous random variables by a finite series of Laguerre polynomials up to the (n₁, …, n_k)th degree, then all the mixed moments up to the order (n₁, …, n_k) of the approximated distribution equal to the mixed moments up to the same order of the random variables. Some numerical results are given for the bivariate central and noncentral multivariate gamma distributions to indicate the usefulness of the approximations.     

A Viable Alternative to Resorting to Statistical Tables

It is shown in this article that, given the moments of a distribution, any percentage point can be accurately determined from an approximation of the corresponding density function in terms of the product of an appropriate baseline density and a polynomial adjustment. This approach, which is based on a moment-matching technique, is not only conceptually simple but easy to implement. As illustrated by several applications, the percentiles so obtained are in excellent agreement with the tabulated values. Whereas statistical tables, if at all available or accessible, can hardly ever cover all the potentially useful combinations of the parameters associated with a random quantity of interest, the proposed methodology has no such limitation.