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.     

Characterization theorems for some discrete distributions based on conditional structure

The joint distribution of (X,Y) is determined if the conditional expectation E {g(X)|Y = y} is given and the conditional distribution of Y|(X = x) is a conditional power series distribution, where g(·) is a function satisfying some minor conditions.     

Estimation of reliability of a component subjected to bivariate exponential stress

In this paper, we estimate the reliability of a component subjected to two different stresses which are independent of the strength of a component. We assume that the distribution of stresses follow a bivariate exponential (BVE) distribution. If X is the strength of a component subjected to two stresses (Y ₁,Y ₂), then the reliability of a component is given by R=P[Y ₁+Y ₂<X]. We estimate R when (Y ₁,Y ₂) follow different BVE models proposed by Marshall-Olkin (1967), Block-Basu-(1974), Freund (1961) and Proschan-Sullo (1974). The distribution of X is assumed to be exponential. The asymptotic normal (AN) distributions of these estimates of R are obtained.     

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.     

Efficient estimation in a nonlinear counting-process regression model

Suppose we observe i.i.d. copies of X, C, Y, where X is a counting process, C is a censoring process talcing only values 0 and 1, and Y is a covariate process. Assume that the intensity process of X is of the form C(s)a(s, Y(s)) with a unknown, but that the distribution of X, C, Y is unspecified otherwise. McKeague and Utikal proposed an estimator for the doubly cumulative hazard f f a(s, y) ds dy and determined its asymptotic distribution. We show that the estimator is regular and efficient in the sense of a Hájek-Inagaki convolution theorem for partially specified models.     

Nonlinear regression in a trivariate elliptical distribution via a linear combination of order statistics

In this paper, by assuming that (X, Y ₁, Y ₂)^T has a trivariate elliptical distribution, we derive the exact joint distribution of X and a linear combination of order statistics from (Y ₁, Y ₂)^T and show that it is a mixture of unified bivariate skew-elliptical distributions. We then derive the corresponding marginal and conditional distributions for the special case of t kernel. We also present these results for an exchangeable case with t kernel and illustrate the established results with an air-pollution data.     

Extreme behaviour for bivariate elliptical distributions

The authors examine the asymptotic behaviour of conditional threshold exceedance probabilities for an elliptically distributed pair (X, Y) of random variables. More precisely, they investigate the limiting behaviour of the conditional distribution of Y given that X becomes extreme. They show that this behaviour differs between regularly and rapidly varying tails.     

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.     

Exact and approximate distributions for the product of inverted Dirichlet components

It is well known that X/(X+Y) has the beta distribution when X and Y follow the inverted Dirichlet distribution. In this paper, we derive the exact distribution of the product P=XY (involving the Gauss hypergeometric function) and the corresponding moment properties. We also propose two approximations and show evidence of their goodness of fit. The work is motivated by real-life examples.     

Exponential models,brownian motion,and independence

Some examples of steep, reproductive exponential models are considered. These models are shown to possess a τ-parallel foliation in the terminology of Barndorff-Nielsen and Blaesild. The independence of certain functions follows directly from the foliation. Suppose X(t) is a Wiener process with drift where X(t) = W(t) + ct, 0 < t < T. Furthermore let Y = max [X(s), 0 < s < T]. The joint density of Y and X = X(T), the end value, is studied within the framework of an exponential model, and it is shown that Y(Y – X) is independent of X. It is further shown that Y(Y – X) suitably scaled has an exponential distribution. Further examples are considered by randomizing on T.     

Note on estimation of reliability under bivariate pareto stress-strength model

In this paper, we discuss the problem of estimating reliability (R) of a component based on maximum likelihood estimators (MLEs). The reliability of a component is given byR=P[Y<X]. Here X is a random strength of a component subjected to a random stress(Y) and (X,Y) follow a bivariate pareto(BVP) distribution. We obtain an asymptotic normal(AN) distribution of MLE of the reliability(R).     

Correlation estimation in Downton's bivariate exponential distribution using incomplete samples

Suppose (X, Y) has a Downton's bivariate exponential distribution with correlation ρ. For a random sample of size n from (X, Y), let X _r:n be the rth X-order statistic and Y _[r:n] be its concomitant. We investigate estimators of ρ when all the parameters are unknown and the available data is an incomplete bivariate sample made up of (i) all the Y-values and the ranks of associated X-values, i.e. (i, Y _[i:n]), 1≤i≤n, and (ii) a Type II right-censored bivariate sample consisting of (X _i:n , Y _[i:n]), 1≤i≤r<n. In both setups, we use simulation to examine the bias and mean square errors of several estimators of ρ and obtain their estimated relative efficiencies. The preferred estimator under (i) is a function of the sample correlation of (Y _i:n , Y _[i:n]) values, and under (ii), a method of moments estimator involving the regression function is preferred.     

Estimation of R=P[Y<X] for three-parameter generalized Rayleigh distribution

Surles and Padgett [Inference for reliability and stress–strength for a scaled Burr type X distribution. Lifetime Data Anal. 2001;7:187–200] introduced a two-parameter Burr-type X distribution, which can be described as a generalized Rayleigh distribution. In this paper, we consider the estimation of the stress–strength parameter R=P[Y<X], when X and Y are both three-parameter generalized Rayleigh distributions with the same scale and locations parameters but different shape parameters. It is assumed that they are independently distributed. It is observed that the maximum-likelihood estimators (MLEs) do not exist, and we propose a modified MLE of R. We obtain the asymptotic distribution of the modified MLE of R, and it can be used to construct the asymptotic confidence interval of R. We also propose the Bayes estimate of R and the construction of the associated credible interval based on importance sampling technique. Analysis of two real data sets, (i) simulated and (ii) real, have been performed for illustrative purposes.     

Scale Checks in Censored Regression

Abstract. Suppose the random vector (X,Y) satisfies the regression model Y = m(X) + σ (X) ? , where m (?) and σ (?) are unknown location and scale functions and ? is independent of X. The response Y is subject to random right censoring, and the covariate X is completely observed. A new test for a specific parametric form of any scale function σ (?) (including the standard deviation function) is proposed. Its statistic is based on the distribution of the residuals obtained from the assumed regression model. Weak convergence of the corresponding process is obtained, and its finite sample behaviour is studied via simulations. Finally, characteristics of the test are illustrated in the analysis of a fatigue data set.     

Synthetic data method to incorporate external information into a current study

We consider the situation where there is a known regression model that can be used to predict an outcome, Y, from a set of predictor variables X . A new variable B is expected to enhance the prediction of Y. A dataset of size n containing Y, X and B is available, and the challenge is to build an improved model for Y| X ,B that uses both the available individual level data and some summary information obtained from the known model for Y| X . We propose a synthetic data approach, which consists of creating m additional synthetic data observations, and then analyzing the combined dataset of size n + m to estimate the parameters of the Y| X ,B model. This combined dataset of size n + m now has missing values of B for m of the observations, and is analyzed using methods that can handle missing data (e.g., multiple imputation). We present simulation studies and illustrate the method using data from the Prostate Cancer Prevention Trial. Though the synthetic data method is applicable to a general regression context, to provide some justification, we show in two special cases that the asymptotic variances of the parameter estimates in the Y| X ,B model are identical to those from an alternative constrained maximum likelihood estimation approach. This correspondence in special cases and the method's broad applicability makes it appealing for use across diverse scenarios. The Canadian Journal of Statistics 47: 580–603; 2019 © 2019 Statistical Society of Canada     

Non-central bivariate beta distribution

Let U, V and W be independent random variables, U and V having a gamma distribution with respective shape parameters a and b, and W having a non-central gamma distribution with shape and non-centrality parameters c and δ, respectively. Define X = U/(U + W) and Y = V/(V + W). Clearly, X and Y are correlated each having a non-central beta type 1 distribution, X ~ NCB1 (a,c;d){X \sim {\rm NCB1} (a,c;\delta)} and Y ~ NCB1 (b,c;d){Y \sim {\rm NCB1} (b,c;\delta)} . In this article we derive the joint probability density function of X and Y and study its properties.     

Bivariate Contours of Copula

Quantile functions associated with bivariate copulas are considered. Some of their structural properties are studied. Quantile functions allow one to express the cdf of the random variable C(X, Y), where (X, Y) is distributed as C(x, y) and where C is a copula. Quantile functions provide also a simple algorithm for simulating random observations.     

Propriétés statistiques des copules de valeurs extrêmes bidimensionnelles

Let (X, Y) be a bivariate random vector whose distribution function H(x, y) belongs to the class of bivariate extreme-value distributions. If F₁ and F₂ are the marginals of X and Y, then H(x, y) = C{F₁(x),F₂(y)}, where C is a bivariate extreme-value dependence function. This paper gives the joint distribution of the random variables Z = {log F₁(X)}/{log F₁(X)F₂(Y)} and W = C{F₁{(X),F₂(Y)}. Using this distribution, an algorithm to generate random variables having bivariate extreme-value distribution is présentés. Furthermore, it is shown that for any bivariate extreme-value dependence function C, the distribution of the random variable W = C{F₁(X),F₂(Y)} belongs to a monoparametric family of distributions. This property is used to derive goodness-of-fit statistics to determine whether a copula belongs to an extreme-value family.     

A LINEAR MODEL WITH ERRORS LACKING A VARIANCE1

The problem discussed is that of estimating β= (β¹, …, β_k) in the model Y=βX +ε when X has a specified multivariate distribution and the error ε does not necessarily have a finite second moment, for example, ε symmetric stable. We construct a moment estimator based on the empirical characteristic function and establish asymptotic unbiassedness and normality. Most of the paper is concerned with the case when X is normal. Forms of the suggested estimator are given in (2.5), (4.6) and (5.5).     

Sum,product and ratio of Pareto and gamma variables

The exact distributions of X+Y, X Y and X/(X+Y) are studied when X and Y are independent Pareto and gamma random variables. Applications are discussed, to real problems in clinical trials, computer networks and economics.