Limit theorems for multivariate discrete distributions

According to the usual law of small numbers a multivariate Poisson distribution is derived by defining an appropriate model for multivariate Binomial distributions and examining their behaviour for large numbers of trials and small probabilities of marginal and simultaneous successes. The weak limit law is a generalization of Poisson's distribution to larger finite dimensions with arbitrary dependence structure. Compounding this multivariate Poisson distribution by a Gamma distribution results in a multivariate Pascal distribution which is again asymptotically multivariate Poisson. These Pascal distributions contain a class of multivariate geometric distributions. Finally the bivariate Binomial distribution is shown to be the limit law of appropriate bivariate hypergeometric distributions. Proving the limit theorems mentioned here as well as understanding the corresponding limit distributions becomes feasible by using probability generating functions.     

Parametric inference from system lifetime data under a proportional hazard rate model

In this paper, we discuss the statistical inference of the lifetime distribution of components based on observing the system lifetimes when the system structure is known. A general proportional hazard rate model for the lifetime of the components is considered, which includes some commonly used lifetime distributions. Different estimation methods—method of moments, maximum likelihood method and least squares method—for the proportionality parameter are discussed. The conditions for existence and uniqueness of method of moments and maximum likelihood estimators are presented. Then, we focus on a special case when the lifetime distributions of the components are exponential. Computational formulas for point and interval estimations of the unknown mean lifetime of the components are provided. A Monte Carlo simulation study is used to compare the performance of these estimation methods and recommendations are made based on these results. Finally, an example is provided to illustrate the methods proposed in this paper.     

A unified approach to bivariate discrete distributions

Here we introduce a bivariate generalized hypergeometric factorial moment distribution (BGHFMD) through its probability generating function (p.g.f.) whose marginal distributions are the generalized hypergeometric factorial moment distributions introduced by Kemp and Kemp (Bull Int Stat Inst 43:336–338,1969). Well-known bivariate versions of distributions such as binomial, negative binomial and Poisson are special cases of this distribution. A genesis of the distribution and explicit closed form expressions for the probability mass function of the BGHFMD, its factorial moments and the p.g.f.’s of its conditional distributions are derived here. Certain recurrence relations for probabilities, moments and factorial moments of the bivariate distribution are also established.     

Progressive stress accelerated life tests under finite mixture models

In this paper, progressive stress accelerated life tests are considered when the lifetime of a product under use condition follows a finite mixture of distributions. The experiment is performed when each of the components in the mixture follows a general class of distributions which includes, among others, the Weibull, compound Weibull, power function, Gompertz and compound Gompertz distributions. It is assumed that the scale parameter of each component satisfies the inverse power low, the progressive stress is directly proportional to time and the cumulative exposure model for the effect of changing stress holds. Based on type-I censoring, the maximum likelihood estimates (MLEs) of the parameters under consideration are obtained. A special attention is paid to a mixture of two Rayleigh components. Simulation results are carried out to study the precision of the MLEs and to obtain confidence intervals for the parameters involved.     

On a class of discrete distributions arising from the birth-death-with-immigration process

This paper considers a class of distributions which may be regarded as the convolution of a negative binomial and a stopped-sum generalized hypergeometric factorial-moment random variables. Some properties are derived and it is shown that this class of distributions is a subset of distributions for the birth-and-death process with immigration (also reversible counter system). Formulations by mixing, limiting distributions and maximum likelihood equations are also discussed.     

Optimal rules and robust Bayes estimation of a Gamma scale parameter

For estimating an unknown scale parameter of Gamma distribution, we introduce the use of an asymmetric scale invariant loss function reflecting precision of estimation. This loss belongs to the class of precautionary loss functions. The problem of estimation of scale parameter of a Gamma distribution arises in several theoretical and applied problems. Explicit form of risk-unbiased, minimum risk scale-invariant, Bayes, generalized Bayes and minimax estimators are derived. We characterized the admissibility and inadmissibility of a class of linear estimators of the form $cX\,{+}\,d$ , when $X\sim \varGamma (\alpha ,\eta )$ . In the context of Bayesian statistical inference any statistical problem should be treated under a given loss function by specifying a prior distribution over the parameter space. Hence, arbitrariness of a unique prior distribution is a critical and permanent question. To overcome with this issue, we consider robust Bayesian analysis and deal with Gamma minimax, conditional Gamma minimax, the stable and characterize posterior regret Gamma minimax estimation of the unknown scale parameter under the asymmetric scale invariant loss function in detail.     

Compound gamma,beta and F distributions

Summary In this paper a compound gamma distribution has been derived by compounding a gamma distribution with another gamma distribution. The resulting compound gamma distribution has been reduced to the Beta distributions of the first kind and the second kind and to theF distribution by suitable transformations. This includes theLomax distribution as a special case which enjoys a useful property. Moment estimators for two of its parameters are explicitly obtained, which tend to a bivariate normal distribution. The paper contains expressions for a bivariate probability density function, its conditional expectation, conditional variance and the product moment correlation coefficient. Finally, all the parameters of the compound gamma distribution are explicitly expressed in terms of the functions of the moments of the functions of random variables in two different ways. This note is based on a technical report prepared by the author while he was with the Procter and Gamble Company.     

The UMVUE ofP(X<Y) based on type-II censored samples from Weinman multivariate exponential distributions

Based on two independent samples from Weinman multivariate exponential distributions with unknown scale parameters, uniformly minimum variance unbiased estimators ofP(X<Y) are obtained for both, unknown and known common location parameter. The samples are permitted to be Type-II censored with possibly different numbers of observations. Since sampling from two-parameter exponential distributions is contained in the model as a particular case, known results for complete and censored samples are generalized. In the case of an unknown common location parameter with a certain restriction of the model, the UMVUE is shown to have a Gauss hypergeometric distribution, which is further examined. Moreover, explicit expressions for the variances of the estimators are derived and used to calculate the relative efficiency.     

Adversarial and Amiable Inference in Medical Diagnosis,Reliability and Survival Analysis

In this paper, we develop a family of bivariate beta distributions that encapsulate both positive and negative correlations, and which can be of general interest for Bayesian inference. We then invoke a use of these bivariate distributions in two contexts. The first is diagnostic testing in medicine, threat detection and signal processing. The second is system survivability assessment, relevant to engineering reliability and to survival analysis in biomedicine. In diagnostic testing, one encounters two parameters that characterize the efficacy of the testing mechanism: test sensitivity and test specificity. These tend to be adversarial when their values are interpreted as utilities. In system survivability, the parameters of interest are the component reliabilities, whose values when interpreted as utilities tend to exhibit co‐operative (amiable) behavior. Besides probability modeling and Bayesian inference, this paper has a foundational import. Specifically, it advocates a conceptual change in how one may think about reliability and survival analysis. The philosophical writings of de Finetti, Kolmogorov, Popper and Savage, when brought to bear on these topics constitute the essence of this change. Its consequence is that we have at hand a defensible framework for invoking Bayesian inferential methods in diagnostics, reliability and survival analysis. Another consequence is a deeper appreciation of the judgment of independent lifetimes. Specifically, we make the important point that independent lifetimes entail at a minimum, a two‐stage hierarchical construction.     

Stochastic comparisons of coherent systems

The study of stochastic comparisons of coherent systems with different structures is a relevant topic in reliability theory. Several results have been obtained for specific distributions. The present paper is focused on distribution-free comparisons, that is, orderings which do not depend on the component distributions. Different assumptions for the component lifetimes are considered which lead us to different comparison techniques. Thus, if the components are independent and identically distributed (IID) or exchangeable, the orderings are obtained by using signatures. If they are just ID (homogeneous components), then ordering results for distorted distributions are used. In the general case or in the case of independent (heterogeneous) components, a similar technique based on generalized distorted distributions is applied. In these cases, the ordering results may depend on the copula used to model the dependence between the component lifetimes. Some illustrative examples are included in each case.     

Special functions for the study of economic dynamics: The case of the Lucas-Uzawa model

The special functions are intensively used in mathematical physics to solve differential systems. We argue that they should be most useful in economic dynamics, notably in the assessment of the transition dynamics of endogenous economic growth models. We illustrate our argument on the famous Lucas-Uzawa model, which we solve by the means of Gaussian hypergeometric functions. We show how the use of Gaussian hypergeometric functions allows for an explicit representation of the equilibrium dynamics of all variables in level. The parameters of the involved hypergeometric functions are identified using the Pontryagin conditions arising from the underlying optimization problems. In contrast to the pre-existing approaches, our method is global and does not rely on dimension reduction.     

Tests of independence in a bivariate exponential distribution

The problem of testing independence in a two component series system is considered. The joint distribution of component lifetimes is modeled by the Pickands bivariate exponential distribution, which includes the widely used Marshall and Olkins distribution and the Gumbels type II distribution. The case of identical components is first addressed. Uniformly most powerful unbiased test (UMPU) and likelihood ratio test are obtained. It is shown that inspite of a nuisance parameter, the UMPU test is unconditional and this test turns out to be the same as the likelihood ratio test. The case of nonidentical components is also addressed and both UMPU and likelihood ratio tests are obtained. A UMPU test is obtained to test the identical nature of the components and extensions to the type II censoring scheme and multiple component systems are also discussed. Some modifications to account for the difference in parameters under test and use conditions are also discussed.     

Log-concavity and monotonicity of hazard and reversed hazard functions of univariate and multivariate skew-normal distributions

This paper establishes the log-concavity property of several forms of univariate and multivariate skew-normal distributions. This property is then used to prove the monotonicity of the hazard as well as reversed hazard functions. The log-concavity and monotonicity of the hazard and reversed hazard functions of series and parallel systems of components is then discussed. The corresponding results for the multivariate normal distribution follow readily as special cases.     

Mixed systems with minimal and maximal lifetime variances

We consider the mixed systems composed of a fixed number of components whose lifetimes are i.i.d. with a known distribution which has a positive and finite variance. We show that a certain of the k-out-of-n systems has the minimal lifetime variance, and the maximal one is attained by a mixture of series and parallel systems. The number of the k-out-of-n system, and the probability weights of the mixture depend on the first two moments of order statistics of the parent distribution of the component lifetimes. We also show methods of calculating extreme system lifetime variances under various restrictions on the system lifetime expectations, and vice versa.     

On the residual lifetime of coherent systems with heterogeneous components

The residual lifetime is of significant interest in reliability and survival analysis. In this article, we obtain a mixture representation for the reliability function of the residual lifetime of a coherent system with heterogeneous components in terms of the reliability functions of residual lifetimes of order statistics. Some stochastic comparisons are made on the residual lifetimes of the systems. Some examples are also given to illustrate the main results.     

A product Pareto distribution

The Pareto distributions are becoming increasing prominent in several applied areas. In this note, a new Pareto distribution is introduced. It takes the form of the product of two Pareto probability density functions. Various structural properties of this distribution are derived, including its cumulative distribution function, moments, mean deviation about the mean, mean deviation about the median, entropy, asymptotic distribution of the extreme order statistics, method of moments estimates, maximum likelihood estimates and the Fisher information matrix. The calculations involve the use of several special functions.     

Stochastic comparisons of multivariate reversed frailty model

Multivariate frailty approaches are most commonly used to define distributions of random vectors, which represent lifetimes of individuals or components and stochastically compare them in terms of various multivariate orders. In this paper, we study a multivariate shared reversed frailty model and a general multivariate reversed frailty mixture model, and derive sufficient conditions for some of the stochastic orderings to hold among the random vectors. We also consider a particular case of a general multivariate mixture model in which the baseline distribution function is represented in terms of a copula and study stochastic comparisons (stochastic and lower orthant order) among the two random vectors.     

On universal admissibility of scale parameter estimators

A characterization result of Kushary (1998) regarding universal admissibility of equivariant estimators in the one parameter gamma distribution is generalized to a scale family of distributions with monotone likelihood ratio. New examples are given, among them the F-distribution with a scale parameter. In particular, universal admissibility is characterized within the class of location-scale equivariant estimators of the ratio of the variances of two normal distributions with unknown means. In this context the maximum likelihood estimator is shown to be universally inadmissible by virtue of a general sufficient condition for universal inadmissibility of a scale equivariant estimator. Received: January 2000     

Bayesian estimation of inefficiency heterogeneity in stochastic frontier models

Estimation of the one sided error component in stochastic frontier models may erroneously attribute firm characteristics to inefficiency if heterogeneity is unaccounted for. However, unobserved inefficiency heterogeneity has been little explored. In this work, we propose to capture it through a random parameter which may affect the location, scale, or both parameters of a truncated normal inefficiency distribution using a Bayesian approach. Our findings using two real data sets, suggest that the inclusion of a random parameter in the inefficiency distribution is able to capture latent heterogeneity and can be used to validate the suitability of observed covariates to distinguish heterogeneity from inefficiency. Relevant effects are also found on separating and shrinking individual posterior efficiency distributions when heterogeneity affects the location and scale parameters of the one-sided error distribution, and consequently affecting the estimated mean efficiency scores and rankings. In particular, including heterogeneity simultaneously in both parameters of the inefficiency distribution in models that satisfy the scaling property leads to a decrease in the uncertainty around the mean scores and less overlapping of the posterior efficiency distributions, which provides both more reliable efficiency scores and rankings.     

Distribution of personal income: Development of a new model and its application to U.S. income data

This paper proposes a four-parameter statistical model of the personalized distribution of income using the ‘income share elasticity’ approach suggested by Esteban (1986). Our proposed model includes the Singh–Maddala (1976) and Dagum (1977) distributions as special cases. The generalized beta II distribution of McDonald (1984) is also a variant of this model. It appears to give an excellent fit to US income data and its empirical performance turns out to be superior to those of the Singh–Maddala (1976), Dagum (1977), the five-parameter Champernowne (1953) distributions and the generalized beta II distribution of McDonald (1984) on some data.