A note on the ratio of normal and Laplace random variables

The normal and Laplace are the two earliest known continuous distributions in statistics and the two most popular models for analyzing symmetric data. In this note, the exact distribution of the ratio | X / Y | is derived when X and Y are respectively normal and Laplace random variables distributed independently of each other. A MAPLE program is provided for computing the associated percentage points. An application of the derived distribution is provided to a discriminant problem.     

Study of incompatibility or near compatibility of bivariate discrete conditional probability distributions through divergence measures

Consider a two-dimensional discrete random variable (X, Y) with possible values 1, 2, …, I for X and 1, 2, …, J for Y. For specifying the distribution of (X, Y), suppose both conditional distributions, of X given Y and of Y given X, are provided. Under this setting, we present here different ways of measuring discrepancy between incompatible conditional distributions in the finite discrete case. In the process, we also suggest different ways of defining the most nearly compatible distributions in incompatible cases. Many new divergence measures are discussed along with those that are already known for determining the most nearly compatible joint distribution P. Finally, a comparative study is carried out between all these divergence measures as some examples.     

A compound beta distribution with applications in finance

If X and Y are gamma distributed independent random variables then it is well known that the ratio X / (X + Y) has the beta distribution. In this note, the distribution of W = X / (X + Y) is considered when X and Y have the compound gamma distribution. We refer to the distribution of W as compound beta and describe an application to consumer price indices to show that compound beta is a better model than one based on the standard beta distribution. We derive various properties of W, including its probability density function, cumulative distribution function, hazard rate function and moments.     

DISTRIBUTIONAL CHARACTERIZATIONS THROUGH SCALING RELATIONS

Investigated here are aspects of the relation between the laws of X and Y where X is represented as a randomly scaled version of Y. In the case that the scaling has a beta law, the law of Y is expressed in terms of the law of X. Common continuous distributions are characterized using this beta scaling law, and choosing the distribution function of Y as a weighted version of the distribution function of X, where the weight is a power function. It is shown, without any restriction on the law of the scaling, but using a one‐parameter family of weights which includes the power weights, that characterizations can be expressed in terms of known results for the power weights. Characterizations in the case where the distribution function of Y is a positive power of the distribution function of X are examined in two special cases. Finally, conditions are given for existence of inverses of the length‐bias and stationary‐excess operators.     

Estimation of δ=P(X<Y) for Burr XII distribution based on the progressively first failure-censored samples

Let X and Y have two-parameter Burr XII distributions. The maximum-likelihood estimator of δ=P(X<Y) is studied under the progressively first failure-censored samples. Three confidence intervals of δ are constructed by using an asymptotic distribution of the maximum-likelihood estimator of δ and two bootstrapping procedures, respectively. Some computational results from intensive simulations are presented. An illustrative example is provided to demonstrate the application of the proposed method.     

CHARACTERIZATIONS OF SHIFT-INVARIANT DISTRIBUTIONS BASED ON SUMMATION MODULO ONE

Let X₁Y_1,…, Y_n be independent random variables. We characterize the distributions of X and Y_j satisfying the equation {X+Y₁++Y_n}=_dX, where {Z} denotes the fractional part of a random variable Z. In the case of full generality, either X is uniformly distributed on [0,1), or Y_j has.a shifted lattice distribution and X is shift-invariant. We also give a characterization of shift-invariant distributions. Finally, we consider some special cases of this equation.     

New sensitivity analysis subordinated to a contrast

ABSTRACT

In a model of the form Y = h(X₁, …, X_d) where the goal is to estimate a parameter of the probability distribution of Y, we define new sensitivity indices which quantify the importance of each variable X_i with respect to this parameter of interest. The aim of this paper is to define goal oriented sensitivity indices and we will show that Sobol indices are sensitivity indices associated to a particular characteristic of the distribution Y. We name the framework we present as Goal Oriented Sensitivity Analysis (GOSA).     

Bayesian analysis of head and neck cancer data using generalized inverse Lindley stress–strength reliability model

In this article, we present the analysis of head and neck cancer data using generalized inverse Lindley stress–strength reliability model. We propose Bayes estimators for estimating P(X > Y), when X and Y represent survival times of two groups of cancer patients observed under different therapies. The X and Y are assumed to be independent generalized inverse Lindley random variables with common shape parameter. Bayes estimators are obtained under the considerations of symmetric and asymmetric loss functions assuming independent gamma priors. Since posterior becomes complex and does not possess closed form expressions for Bayes estimators, Lindley’s approximation and Markov Chain Monte Carlo techniques are utilized for Bayesian computation. An extensive simulation experiment is carried out to compare the performances of Bayes estimators with the maximum likelihood estimators on the basis of simulated risks. Asymptotic, bootstrap, and Bayesian credible intervals are also computed for the P(X > Y).     

The linear combination,product and ratio of Laplace random variables

The distributions of linear combinations, products and ratios of random variables arise in many areas of engineering. In this paper, the exact distributions of the linear combination α X+β Y, the product |X Y| and the ratio |X/Y| are derived when X and Y are independent Laplace random variables. The Laplace distribution, being the oldest model for continuous data, has been one of the most popular models for measurement errors in engineering.     

On Censored Bivariate Random Variables: Copula,Characterization, and Estimation

Let (X, Y) be a bivariate random vector with joint distribution function F_{X, Y}(x, y) = C(F(x), G(y)), where C is a copula and F and G are marginal distributions of X and Y, respectively. Suppose that (X_i, Y_i), i = 1, 2, …, n is a random sample from (X, Y) but we are able to observe only the data consisting of those pairs (X_i, Y_i) for which X_i ? Y_i. We denote such pairs as (X*_i, Y_i*), i = 1, 2, …, ν, where ν is a random variable. The main problem of interest is to express the distribution function F_{X, Y}(x, y) and marginal distributions F and G with the distribution function of observed random variables X* and Y*. It is shown that if X and Y are exchangeable with marginal distribution function F, then F can be uniquely determined by the distributions of X* and Y*. It is also shown that if X and Y are independent and absolutely continuous, then F and G can be expressed through the distribution functions of X* and Y* and the stress–strength reliability P{X ? Y}. This allows also to estimate P{X ? Y} with the truncated observations (X*_i, Y_i*). The copula of bivariate random vector (X*, Y*) is also derived.     

Stress-Strength Reliability of a Two-Parameter Bathtub-shaped Lifetime Distribution Based on Progressively Censored Samples

Based on progressively Type-II censored samples, this article deals with inference for the stress-strength reliability R = P(Y < X) when X and Y are two independent two-parameter bathtub-shape lifetime distributions with different scale parameters, but having the same shape parameter. Different methods for estimating the reliability are applied. The maximum likelihood estimate of R is derived. Also, its asymptotic distribution is used to construct an asymptotic confidence interval for R. Assuming that the shape parameter is known, the maximum likelihood estimator of R is obtained. Based on the exact distribution of the maximum likelihood estimator of R an exact confidence interval of that has been obtained. The uniformly minimum variance unbiased estimator are calculated for R. Bayes estimate of R and the associated credible interval are also got under the assumption of independent gamma priors. Monte Carlo simulations are performed to compare the performances of the proposed estimators. One data analysis has been performed for illustrative purpose. Finally, we will generalize this distribution to the proportional hazard family with two parameters and derive various estimators in this family.     

The exact distribution function of the ratio of two dependent quadratic forms

A closed-form representation of the distribution function of the ratio of two linear combinations of Chi-squared variables is derived. The ratio is of the following form R = (X + aY)/(bY + Z), where X, Y, Z are independent Chi-square variables and a, b > 0. Two methods of obtaining the distribution function of this ratio are used. The exact density function of such a ratio is then obtained by differentiation. Two numerical examples are provided.     

Information Properties for Concomitants of Order Statistics in Farlie–Gumbel–Morgenstern (FGM) Family

Let (X _i , Y _i ), i = 1, 2,…, n be independent and identically distributed random variables from some continuous bivariate distribution. If X _(r) denotes the rth-order statistic, then the Y's associated with X _(r) denoted by Y _[r] is called the concomitant of the rth-order statistic. In this article, we derive an analytical expression of Shannon entropy for concomitants of order statistics in FGM family. Applying this expression for some well-known distributions of this family, we obtain the exact form of Shannon entropy, some of the information properties, and entropy bounds for concomitants of order statistics. Some comparisons are also made between the entropy of order statistics X _(r) and the entropy of its concomitants Y _[r]. In this family, we show that the mutual information between X _(r) and Y _[r], and Kullback–Leibler distance among the concomitants of order statistics are all distribution-free. Also, we compare the Pearson correlation coefficient between X _(r) and Y _[r] with the mutual information of (X _(r), Y _[r]) for the copula model of FGM family.     

Adaptive Deconvolution on the Non‐negative Real Line

In this paper, we consider the problem of adaptive density or survival function estimation in an additive model defined by Z=X+Y with X independent of Y, when both random variables are non‐negative. This model is relevant, for instance, in reliability fields where we are interested in the failure time of a certain material that cannot be isolated from the system it belongs. Our goal is to recover the distribution of X (density or survival function) through n observations of Z, assuming that the distribution of Y is known. This issue can be seen as the classical statistical problem of deconvolution that has been tackled in many cases using Fourier‐type approaches. Nonetheless, in the present case, the random variables have the particularity to be supported. Knowing that, we propose a new angle of attack by building a projection estimator with an appropriate Laguerre basis. We present upper bounds on the mean squared integrated risk of our density and survival function estimators. We then describe a non‐parametric data‐driven strategy for selecting a relevant projection space. The procedures are illustrated with simulated data and compared with the performances of a more classical deconvolution setting using a Fourier approach. Our procedure achieves faster convergence rates than Fourier methods for estimating these functions.     

On the Linear Combination of Laplace and Logistic Random Variables

The distribution of linear combinations of random variables arises explicitly in many areas of engineering. This has increased the need to have available the widest possible range of statistical results on linear combinations of random variables. In this note, the exact distribution of the linear combination α X+β Y is derived when X and Y are Laplace and logistic random variables distributed independently of each other. Extensive tabulations of the associated percentage points obtained by inverting the derived distribution are also given.     

Improved estimation of scale parameters of morgenstern type bivariate uniform distribution using ranked set sampling

ABSTRACT

In this article we suggest some improved version of estimators of scale parameter of Morgenstern-type bivariate uniform distribution (MTBUD) based on the observations made on the units of the ranked set sampling regarding the study variable Y which is correlated with the auxiliary variable X, when (X, Y) follows a MTBUD. We also suggest some linear shrinkage estimators of scale parameter of Morgenstern type bivariate uniform distribution (MTBUD). Efficiency comparisons are also made in this work.     

Laplace random variables with application to price indices

Laplace distributions are becoming increasingly popular models in economics and finance. In this note, the exact distribution of the ratio Z=|X/Y| is derived when X and Y are independent Laplace random variables. This distribution arises when one is interested in comparing the performances of two economic or financial entities. We consider estimation issues of the distribution and illustrate an application for consumer price indices from the six major economics. Several computer programs are given for implementation of the methods used.     

Multivariate normal distribution approaches for dependently truncated data

Many statistical methods for truncated data rely on the independence assumption regarding the truncation variable. In many application studies, however, the dependence between a variable X of interest and its truncation variable L plays a fundamental role in modeling data structure. For truncated data, typical interest is in estimating the marginal distributions of (L, X) and often in examining the degree of the dependence between X and L. To relax the independence assumption, we present a method of fitting a parametric model on (L, X), which can easily incorporate the dependence structure on the truncation mechanisms. Focusing on a specific example for the bivariate normal distribution, the score equations and Fisher information matrix are provided. A robust procedure based on the bivariate t-distribution is also considered. Simulations are performed to examine finite-sample performances of the proposed method. Extension of the proposed method to doubly truncated data is briefly discussed.     

Bayesian inference on P(X > Y) in bivariate Rayleigh model

In the literature, assuming independence of random variables X and Y, statistical estimation of the stress–strength parameter R = P(X > Y) is intensively investigated. However, in some real applications, the strength variable X could be highly dependent on the stress variable Y. In this paper, unlike the common practice in the literature, we discuss on estimation of the parameter R where more realistically X and Y are dependent random variables distributed as bivariate Rayleigh model. We derive the Bayes estimates and highest posterior density credible intervals of the parameters using suitable priors on the parameters. Because there are not closed forms for the Bayes estimates, we will use an approximation based on Laplace method and a Markov Chain Monte Carlo technique to obtain the Bayes estimate of R and unknown parameters. Finally, simulation studies are conducted in order to evaluate the performances of the proposed estimators and analysis of two data sets are provided.     

Inference on P(Y < X) in Bivariate Rayleigh Distribution

This paper deals with the estimation of reliability R = P(Y < X) when X is a random strength of a component subjected to a random stress Y, and (X, Y) follows a bivariate Rayleigh distribution. The maximum likelihood estimator of R and its asymptotic distribution are obtained. An asymptotic confidence interval of R is constructed using the asymptotic distribution. Also, two confidence intervals are proposed based on Bootstrap method and a computational approach. Testing of the reliability based on asymptotic distribution of R is discussed. Simulation study to investigate performance of the confidence intervals and tests has been carried out. Also, a numerical example is given to illustrate the proposed approaches.