A robust multivariate measurement error model with skew-normal/independent distributions and Bayesian MCMC implementation

Skew-normal/independent distributions are a class of asymmetric thick-tailed distributions that include the skew-normal distribution as a special case. In this paper, we explore the use of Markov Chain Monte Carlo (MCMC) methods to develop a Bayesian analysis in multivariate measurement errors models. We propose the use of skew-normal/independent distributions to model the unobserved value of the covariates (latent variable) and symmetric normal/independent distributions for the random errors term, providing an appealing robust alternative to the usual symmetric process in multivariate measurement errors models. Among the distributions that belong to this class of distributions, we examine univariate and multivariate versions of the skew-normal, skew-t, skew-slash and skew-contaminated normal distributions. The results and methods are applied to a real data set.     

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.     

The shape of the hazard function: Does the generalized gamma have the last word?

This paper compares the five-parameter beta generalized gamma (BGG) distribution to the three-parameter generalized gamma (GG). Both distributions include the four standard hazard shapes that we believe is an important property for any parametric family. For several BGG distributions, we select matching GGs and compute the Kullback-Liebler distance, observing remarkable agreement. We explore the beta parameters' influence on the matched GG parameters, detecting a strong connection between the distributions. Lastly, we compare the distributions using two real-data examples. We conclude from these comparisons that the BGG is not likely to be more useful for analytical purposes than the simpler GG.     

Properties of Coherent Systems with Dependent Components

We introduce two new families of univariate distributions that we call hyperminimal and hypermaximal distributions. These families have interesting applications in the context of reliability theory in that they contain that of coherent system lifetime distributions. For these families, we obtain distributions, bounds, and moments. We also define the minimal and maximal signatures of a coherent system with exchangeable components which allow us to represent the system distribution as generalized mixtures (i.e., mixtures with possibly negative weights) of series and parallel systems. These results can also be applied to order statistics (k-out-of-n systems). Finally, we give some applications studying coherent systems with different multivariate exponential joint distributions.     

A strategy for constructing multivariate distributions

Given a random vector (X₁,…, X_n) for which the univariate and bivariate marginal distributions belong to some specified families of distributions, we present a procedure for constructing families of multivariate distributions with the specified univariate and bivariate margins. Some general properties of the resulting families of multivariate distributions are reviewed. This procedure is illustrated by generalizing the bivariate Plackett (1965) and Clayton (1978) distributions to three dimensions. In addition to providing rich families of models for data analysis, this method of construction provides a convenient way of simulating observations from multivariate distributions with specific types of univariate and bivariate marginal distributions. A general algorithm for simulating random observations from these families of multivariate distributions is presented     

The generalized exponential family of distributions and its characteristics

In this paper we generalize the exponential family (EF) of distributions into a wider family which includes important distributions such as the normal, log-normal, Student-t, Cauchy, logistic and Birnbaum–Saunders distributions. Furthermore, we derive several characteristics of the proposed family. The importance of such family is also discussed.     

On proportional odds models

Recently, Marshall and Olkin (Biometrika 84(3):641–652 1997) introduced a family of distributions by adding a new parameter to a survival function. In this paper, we give physical interpretation of the family using odds function. It is shown that the family of distributions satisfies the property of proportional odds function. We, then, develop a generalized family and study its properties. Further, we give various definitions of proportional odds model in the bivariate set up. Based on these, we introduce new families of bivariate distributions and study their properties.     

On the normality of transformed beta and unit- gamma random variables

In this paper we consider properties of the logarithmic and Tukey's lambda-type transformations of random variables that follow beta or unit-gamma distributions. Beta distributions often arise as models for random proportions, and unit-gamma distributions, although not well- known, may serve the same purpose. The latter possess many properties similar to those of beta distributions. Some transformations of random variables that follow a beta distribution are considered by Johnson (1949) and Johnson and Kotz (1970,1973). These are used to obtain a -new"random variable that potentially approximately follows a normal distribution, so that practical analyses become possible. We study normality -related properties of the above transformations. This is done for the first time for unit-gamma distributions. Under the logarithmic transformation the beta and unit-gamma distributions become, respectively, the logarithmic F and generalized logistic distributions. The distributions of the transformed beta and unit-gamma distributions after application of Tukey's lambda-type transformations cannot be derived easily; however, we obtain the first four moments and expressions for the skewness and kudos is of the transformed variables. Values of skewness and kurtosis for a variety of different parameter values are calculated, and in consequence, the near (or not near) normality of the transformed variables is evaluated. Comments on the use of the various transformations are provided..     

Random coefficient regressions: parametric goodness-of-fit tests

Random coefficient regression models have been applied in different fields during recent years and they are a unifying frame for many statistical models. Recently, Beran and Hall (Ann. Statist. 20 (1992) 1970) raised the question of the nonparametric study of the coefficients distribution. Nonparametric goodness-of-fit tests were considered in Delicado and Romo (Ann. Inst. Statist. Math. 51 (1999) 125). In this nonparametric framework, the study of parametric families for the coefficient distributions was started by Beran (Ann. Inst. Statist. Math. (1993) 639). Here we propose statistics for parametric goodness-of-fit tests and we obtain their asymptotic distributions. Moreover, we construct bootstrap approximations to these distributions, proving their validity. Finally, a simulation study illustrates our results.     

Multivariate log-skewed distributions with normal kernel and their applications

We introduce two classes of multivariate log-skewed distributions with normal kernel: the log canonical fundamental skew-normal (log-CFUSN) and the log unified skew-normal. We also discuss some properties of the log-CFUSN family of distributions. These new classes of log-skewed distributions include the log-normal and multivariate log-skew normal families as particular cases. We discuss some issues related to Bayesian inference in the log-CFUSN family of distributions, mainly we focus on how to model the prior uncertainty about the skewing parameter. Based on the stochastic representation of the log-CFUSN family, we propose a data augmentation strategy for sampling from the posterior distributions. This proposed family is used to analyse the US national monthly precipitation data. We conclude that a high-dimensional skewing function lead to a better model fit.     

On the monotonic properties of discrete failure rates

As is well known, the monotonicity of failure rate of a life distribution plays an important role in modeling failure time data. In this paper, we develop techniques for the determination of increasing failure rate (IFR) and decreasing failure rate (DFR) property for a wide class of discrete distributions. Instead of using the failure rate, we make use of the ratio of two consecutive probabilities. The method developed is applied to various well known families of discrete distributions which include the binomial, negative binomial and Poisson distributions as special cases. Finally, a formula is presented to determine explicitly the failure rate of the families considered. This formula is used to determine the failure rate of various classes of discrete distributions. These formulas are explicit but complicated and cannot normally be used to determine the monotonicity of the failure rates.     

Shared gamma frailty models based on reversed hazard rate for modeling Australian twin data

In this paper, we consider shared gamma frailty model with the reversed hazard rate (RHR) with two different baseline distributions, namely the generalized inverse Rayleigh and the exponentiated Gumbel distributions. With these two baseline distributions we propose two different shared frailty models. We develop the Bayesian estimation procedure using Markov Chain Monte Carlo technique to estimate the parameters involved in these models. We present a simulation study to compare the true values of the parameters with the estimated values. A search of the literature suggests that currently no work has been done for these two baseline distributions with a shared gamma frailty with the RHR so far. We also apply these two models by using a real life bivariate survival data set of Australian twin data given by Duffy et a1. (1990) and a better model is suggested for the data.     

Simulating from polynomial-normal distributions

The family of polynomial-normal distributions includes a number of widely used distributions, such as the Gram–Charlier and Edgeworth distributions. In this note, we present three simple algorithms, (i) CDF Inversion, (ii) Acceptance–Rejection, (iii) and Ratio–of–Uniforms, for simulating variates from a polynomial-normal distribution. Details on the efficiency of the Acceptance–Rejection and the Ratio–of–Uniforms algorithms and a comparison across the various implementations are provided.     

Inverse Weibull power series distributions: properties and applications

Inverse Weibull (IW) distribution is one of the widely used probability distributions for nonnegative data modelling, specifically, for describing degradation phenomena of mechanical components. In this paper, by compounding IW and power series distributions we introduce a new lifetime distribution. The compounding procedure follows the same set-up carried out by Adamidis and Loukas [A lifetime distribution with decreasing failure rate. Stat Probab Lett. 1998;39:35–42]. We provide mathematical properties of this new distribution such as moments, estimation by maximum likelihood with censored data, inference for a large sample and the EM algorithm to determine the maximum likelihood estimates of the parameters. Furthermore, we characterize the proposed distributions using a simple relationship between two truncated moments and maximum entropy principle under suitable constraints. Finally, to show the flexibility of this type of distributions, we demonstrate applications of two real data sets.     

Phase Type Distributions with Finite Support

This research is motivated by the fact that many random variables of practical interest have a finite support. For fixed a < b, we consider the distribution of a random variable X = (a + Ymod(b ? a)), where Y is a phase type (PH) random variable. We demonstrate that as we traverse for Y the entire set of PH distributions (or even any subset thereof like Coxian that is dense in the class of distributions on [0, ∞)), we obtain a class of matrix exponential distributions dense in (a, b). We call these Finite Support Phase Type Distributions (FSPH) of the first kind. A simple example shows that though dense, this class by itself is not very efficient for modeling; therefore, we introduce (and derive the EM algorithms for) two other classes of finite support phase type distributions (FSPH). The properties of denseness, connection to Markov chains, the EM algorithm, and ability to exploit matrix-based computations should all make these classes of distributions attractive not only for applied probability but also for a much wider variety of fields using statistical methodologies.     

Prediction from a normal model using a generalized inverse Gaussian prior

In this paper, we derive prediction distribution of future response(s) from the normal distribution assuming a generalized inverse Gaussian (GIG) prior density for the variance. The GIG includes as special cases the inverse Gaussian, the inverted chi-squared and gamma distributions. The results lead to Bessel-type prediction distributions which is in contrast with the Student-t distributions usually obtained using the inverted chi-squared prior density for the variance. Further, the general structure of GIG provides us with new flexible prediction distributions which include as special cases most of the earlier results obtained under normal-inverted chi-squared or vague priors.     

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.     

Variational Bayes for Phase-Type Distribution

This article develops an algorithm for estimating parameters of general phase-type (PH) distribution based on Bayes estimation. The idea of Bayes estimation is to regard parameters as random variables, and the posterior distribution of parameters which is updated by the likelihood function provides estimators of parameters. One of the advantages of Bayes estimation is to evaluate uncertainty of estimators. In this article, we propose a fast algorithm for computing posterior distributions approximately, based on variational approximation. We formulate the optimal variational posterior distributions for PH distributions and develop the efficient computation algorithm for the optimal variational posterior distributions of discrete and continuous PH distributions.     

Identification of models using failure rate and mean residual life of doubly truncated random variables

In this paper, we study the relationship between the failure rate and the mean residual life of doubly truncated random variables. Accordingly, we develop characterizations for exponential, Pareto II and beta distributions. Further, we generalize the identities for the Pearson and the exponential family of distributions given respectively in Nair and Sankaran (1991) and Consul (1995). Applications of these measures in the context of lengthbiased models are also explored.