A flexible family of density functions

The aim of this paper is to define a new family of probability density functions (MR pdf) based on the multiresolution analysis theory. Each function of this family can be seen as a particular type of density mixture.The MR pdf has advantages with regards to estimation over conventional mixtures and it is suitable to model a large variety of square integrable probability density functions.     

Mixtures of regression models with incomplete and noisy data

The estimation of the mixtures of regression models is usually based on the normal assumption of components and maximum likelihood estimation of the normal components is sensitive to noise, outliers, or high-leverage points. Missing values are inevitable in many situations and parameter estimates could be biased if the missing values are not handled properly. In this article, we propose the mixtures of regression models for contaminated incomplete heterogeneous data. The proposed models provide robust estimates of regression coefficients varying across latent subgroups even under the presence of missing values. The methodology is illustrated through simulation studies and a real data analysis.     

Mixtures of Farlie-Gumbel-Morgenstern Copulas

Some modifications of bivariate Farlie-Gumbel-Morgenstern (FGM) copulas are going to be explained in this article. These modifications are generated by using mixtures of bivariate FGM copula functions. The main goal of this study is to determine both the ranges of association parameter and the rate of correlation, and also observe the changes in local dependence function. An application, which is related with simulated data, is conducted and results are illustrated.     

Mixtures of Linear Regression with Measurement Errors

Existing research on mixtures of regression models are limited to directly observed predictors. The estimation of mixtures of regression for measurement error data imposes challenges for statisticians. For linear regression models with measurement error data, the naive ordinary least squares method, which directly substitutes the observed surrogates for the unobserved error-prone variables, yields an inconsistent estimate for the regression coefficients. The same inconsistency also happens to the naive mixtures of regression estimate, which is based on the traditional maximum likelihood estimator and simply ignores the measurement error. To solve this inconsistency, we propose to use the deconvolution method to estimate the mixture likelihood of the observed surrogates. Then our proposed estimate is found by maximizing the estimated mixture likelihood. In addition, a generalized EM algorithm is also developed to find the estimate. The simulation results demonstrate that the proposed estimation procedures work well and perform much better than the naive estimates.     

Models for Dependent Extremes Using Stable Mixtures

Abstract.  This paper unifies and extends results on a class of multivariate extreme value (EV) models studied by Hougaard, Crowder and Tawn. In these models, both unconditional and conditional distributions are themselves EV distributions, and all lower-dimensional marginals and maxima belong to the class. One interpretation of the models is as size mixtures of EV distributions, where the mixing is by positive stable distributions. A second interpretation is as exponential-stable location mixtures (for Gumbel) or as power-stable scale mixtures (for non-Gumbel EV distributions). A third interpretation is through a peaks over thresholds model with a positive stable intensity. The mixing variables are used as a modelling tool and for better understanding and model checking. We study EV analogues of components of variance models, and new time series, spatial and continuous parameter models for extreme values. The results are applied to data from a pitting corrosion investigation.     

Failure Time of Non Homogeneous Gamma Process

We consider the non homogeneous gamma process as a degradation model. The hitting time of a deterministic or random level is studied here. We provide its distribution (both cdf and pdf) explicitly in the first case and in the second case when the threshold is exponentially or gamma distributed. The general case for a random threshold can be approximated by considering mixtures of Erlang distributions. Aging properties are also discussed in this article.     

Identifiability of Finite Mixtures of Elliptical Distributions

Abstract.  We present general results on the identifiability of finite mixtures of elliptical distributions under conditions on the characteristic generators or density generators. Examples include the multivariate t -distribution, symmetric stable laws, exponential power and Kotz distributions. In each case, the shape parameter is allowed to vary in the mixture, in addition to the location vector and the scatter matrix. Furthermore, we discuss the identifiability of finite mixtures of elliptical densities with generators that correspond to scale mixtures of normal distributions.     

Interpreting DNA Mixtures with Related Contributors in Subdivided Populations

Abstract.  In this paper, we study the statistical interpretation of forensic DNA mixtures with related contributors in subdivided populations. Compact general formulae for match probabilities are obtained for two situations: a relative of one tested person is an unknown contributor of a DNA mixture; and two related unknowns are contributors. The effect of kinship and population structure is illustrated using a real case example.     

Characterizing and Approximating Infinite Scale Mixtures of Normals

ABSTRACT

The purpose of this study is to approximate and identify infinite scale mixtures of normals, SMN. A new method for approximating any infinite SMN with a known mixing measure by a finite SMN is presented. In the new method, the modulus of continuity of the normal family as a function of the scale is used to discretize the mixing measure. This method will be used to approximate univariate and multivariate SMN with mean 0. In the multivariate case, two different methods are used to approximate the infinite SMN. Several results related to SMN are proved and other known ones are presented. For example, SMN are characterized by their corresponding Laplace transforms.     

The distribution of the product of two independent generalized trapezoidal random variables

ABSTRACT

In this article, we derive the probability density function (pdf) of the product of two independent generalized trapezoidal random variables having different supports, in closed form, by considering all possible cases. We also show that the results for the product of two triangular and uniform random variables follow as special cases of our main result. As an illustration, we obtain pdf of product for a suitably constrained set of parameters and plot some graphs using MATLAB, which express variation in pdf with change in different parameters of the generalized trapezoidal distribution.     

Uniform Design for Experiments with Mixtures

The goal of uniform mixture design is to scatter the design points in the experimental region uniformly. The commonly used criteria, such as mean square distance, are based on the Euclidean distance. Based on the Lee distance, a new criterion is proposed in this article. And an algorithm, called NTLBG, is also proposed to refine the randomly generated design for the experimental design with mixtures. Some examples show that the design generated by the NTLBG algorithm has a lower criteria value.     

Bayesian Goodness of Fit Testing with Mixtures of Triangular Distributions

Abstract.  We consider the consistency of the Bayes factor in goodness of fit testing for a parametric family of densities against a non-parametric alternative. Sufficient conditions for consistency of the Bayes factor are determined and demonstrated with priors using certain mixtures of triangular densities.     

Approximating Symmetric Distributions via Sampling and Coefficient of Variation

Abstract

The Coefficient of Variation is one of the most commonly used statistical tool across various scientific fields. This paper proposes a use of the Coefficient of Variation, obtained by Sampling, to define the polynomial probability density function (pdf) of a continuous and symmetric random variable on the interval [a, b]. The basic idea behind the first proposed algorithm is the transformation of the interval from [a, b] to [0, b-a]. The chi-square goodness-of-fit test is used to compare the proposed (observed) sample distribution with the expected probability distribution. The experimental results show that the collected data are approximated by the proposed pdf. The second algorithm proposes a new method to get a fast estimate for the degree of the polynomial pdf when the random variable is normally distributed. Using the known percentages of values that lie within one, two and three standard deviations of the mean, respectively, the so-called three-sigma rule of thumb, we conclude that the degree of the polynomial pdf takes values between 1.8127 and 1.8642. In the case of a Laplace (μ, b) distribution, we conclude that the degree of the polynomial pdf takes values greater than 1. All calculations and graphs needed are done using statistical software R.     

Role of concomitants of order statistics in determining parent bivariate distributions

In this paper, we establish the role of concomitants of order statistics in the unique identification of the parent bivariate distribution. From the results developed, we have illustrated by examples the process of determination of the parent bivariate distribution using a marginal pdf and the pdf of either of the concomitant of largest or smallest order statistic on the other variable. An application of the results derived in modeling of a bivariate distribution for data sets drawn from a population as well is discussed.     

Mixtures,embedding, and ancillarity

Ancillary statistics, proposed by Fisher (1925), can be constructed by forming a mixture model (Birnbaum 1962) or can be extracted or derived from a transformation-parameter model (Peisakoff 1951, Fraser 1961) or from the corresponding error-based structural model (Fraser 1968, 1979); these latter models involve an implicit mixture structure. Compound models with ancillaries can also be formed by a cross embedding, discussed from a technical viewpoint in this paper. Of the 25 examples in Buehler (1982), 22 are mixtures or implicit mixtures and 3 correspond to cross embedding. The cross embedding examples exemplify the nonuniqueness difficulties with ancillaries. This paper discusses a simple and two generalized versions of cross embedding but makes no general valuations of these for statistical inference; their role within inference is discussed in Evans, Fraser, and Monette (1984, 1985).     

Ordering Comparison of Logarithmic Series Random Variables with their Mixtures

The problem of comparing some known distributions in various types of stochastic orderings has been of interest to many authors. In particular, several authors have been recently concerned with the comparison of Poisson, binomial, and negative binomial distributions with their respective mixtures. Incidentally, these distributions are among the four well-known distributions of the family of generalized power series distributions (GPSD's). The remaining distribution is the logarithmic series distribution. In this paper, we shall be concerned with comparing this remaining distribution of the class GPSD with its mixture in terms of various types of stochastic orderings such as the simple stochastic, likelihood ratio, uniformly more variable, convex, hazard rate and expectation orderings. Derivation of the results in this case prove to be computationally trickier than the other three. The special case when the means of the two distributions are the same is also discussed. Finally, an illustrative explicit example is provided.     

Nonparametric Estimation of Distribution Functions of Nonstandard Mixtures

ABSTRACT

Nonstandard mixtures are those that result from a mixture of a discrete and a continuous random variable. They arise in practice, for example, in medical studies of exposure. Here, a random variable that models exposure might have a discrete mass point at no exposure, but otherwise may be continuous. In this article we explore estimating the distribution function associated with such a random variable from a nonparametric viewpoint. We assume that the locations of the discrete mass points are known so that we will be able to apply a classical nonparametric smoothing approach to the problem. The proposed estimator is a mixture of an empirical distribution function and a kernel estimate of a distribution function. A simple theoretical argument reveals that existing bandwidth selection algorithms can be applied to the smooth component of this estimator as well. The proposed approach is applied to two example sets of data.     

Relationships among moments of order statistics in samples from two related outlier models and some applications

Balakrishnan (1987a) has recently shown that the moments of order statistics in samples drawn from a continuous population with pdf f(x) symmetric about zero comprising a single outlier with pdf g(x) also symmetric about zero can be expressed in terms of the moments of order statistics in samples drawn from the population obtained by folding the pdf f(x) at zero and the moments of order statistics in samples drawn from the population obtained by folding the pdf f(x) at zero comprising a single outlier with pdf obtained by folding g(x) at zero. The cumulative round off error involved in the numerical evaluation of the moments of order statistics from the symmetric outlier model, using a table of the moments of order statistics from the folded population and the moments of order statistics from the folded outlier model, has also been studied by Balakrishnan (1987a) and shown to be not serious. Making use of these results we study here the robustness of some estimators of th location and scale parameters of a double exponential distribution.     

Application of Concomitants of Order Statistics of Independent Non Identically Distributed Random Variables in Estimation

In this article, we obtain expressions for the pdf of a single concomitant of order statistic and the joint pdf of a pair of concomitants of order statistics of independent non identically distributed random variables. Using these expressions, we find the means, variances and covariances of order statistics arising from independent non identically distributed bivariate Pareto distributions. A method of estimation of a common parameter involved in several bivariate Pareto distributions using concomitants of order statistics is also discussed.