Parameter Estimation for Bivariate Shock Models with Singular Distribution for Censored Data with Concomitant Order Statistics

When two‐component parallel systems are tested, the data consist of Type‐II censored data X_(i), i= 1, n, from one component, and their concomitants Y _[i] randomly censored at X_(r), the stopping time of the experiment. Marshall & Olkin's (1967) bivariate exponential distribution is used to illustrate statistical inference procedures developed for this data type. Although this data type is motivated practically, the likelihood is complicated, and maximum likelihood estimation is difficult, especially in the case where the parameter space is a non‐open set. An iterative algorithm is proposed for finding maximum likelihood estimates. This article derives several properties of the maximum likelihood estimator (MLE) including existence, uniqueness, strong consistency and asymptotic distribution. It also develops an alternative estimation method with closed‐form expressions based on marginal distributions, and derives its asymptotic properties. Compared with variances of the MLEs in the finite and large sample situations, the alternative estimator performs very well, especially when the correlation between X and Y is small.     

On some study of skew-t distributions

Abstract

In this note, through ratio of independent random variables, new families of univariate and bivariate skew-t distributions are introduced. Probability density function for each skew-t distribution will be given. We also derive explicit forms of moments of the univariate skew-t distribution and recurrence relations for its cumulative distribution function. Finally we illustrate the flexibility of this class of distributions with applications to a simulated data and the volcanos heights data.     

A new class of bivariate lifetime distributions

This paper introduces a new class of bivariate lifetime distributions. Let {X_i}_{i ? 1} and {Y_i}_{i ? 1} be two independent sequences of independent and identically distributed positive valued random variables. Define T₁ = min?(X₁, …, X_M) and T₂ = min?(Y₁, …, Y_N), where (M, N) has a discrete bivariate phase-type distribution, independent of {X_i}_{i ? 1} and {Y_i}_{i ? 1}. The joint survival function of (T₁, T₂) is studied.     

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.     

Subsampling-extrapolation bandwidth selection in bivariate kernel density estimation

This paper focuses on bivariate kernel density estimation that bridges the gap between univariate and multivariate applications. We propose a subsampling-extrapolation bandwidth matrix selector that improves the reliability of the conventional cross-validation method. The proposed procedure combines a U-statistic expression of the mean integrated squared error and asymptotic theory, and can be used in both cases of diagonal bandwidth matrix and unconstrained bandwidth matrix. In the subsampling stage, one takes advantage of the reduced variability of estimating the bandwidth matrix at a smaller subsample size m (m < n); in the extrapolation stage, a simple linear extrapolation is used to remove the incurred bias. Simulation studies reveal that the proposed method reduces the variability of the cross-validation method by about 50% and achieves an expected integrated squared error that is up to 30% smaller than that of the benchmark cross-validation. It shows comparable or improved performance compared to other competitors across six distributions in terms of the expected integrated squared error. We prove that the components of the selected bivariate bandwidth matrix have an asymptotic multivariate normal distribution, and also present the relative rate of convergence of the proposed bandwidth selector.     

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     

Weak and strong uniform consistency rates of kernel density estimates for randomly censored data

Let X₁,., X_n, be i.i.d. random variables with distribution function F, and let Y₁,.,.,Y_n be i.i.d. with distribution function G. For i = 1, 2,.,., n set δ_i, = 1 if X_i ≤ Y_i, and 0 otherwise, and X_i, = min{X_i, K_i}. A kernel-type density estimate of f, the density function of F w.r.t. Lebesgue measure on the Borel o-field, based on the censored data (δ_i, X_i), i = 1,.,.,n, is considered. Weak and strong uniform consistency properties over the whole real line are studied. Rates of convergence results are established under higher-order differentiability assumption on f. A procedure for relaxing such assumptions is also proposed.     

Distribution of the correlation coefficient for the class of bivariate elliptical models

We consider n pairs of random variables (X₁₁,X₂₁),(X₁₂,X₂₂),… (X_1n,X_2n) having a bivariate elliptically contoured density of the form where θ₁ θ₂ are location parameters and Δ = ((λ_ik)) is a 2 × 2 symmetric positive definite matrix of scale parameters. The exact distribution of the Pearson product-moment correlation coefficient between X₁ and X₂ is obtained. The usual case when a sample of size n is drawn from a bivariate normal population is a special case of the abovementioned model.     

An Affine-Invariant Generalization of the Wilcoxon Signed-Rank Test for the Bivariate Location Problem

This paper proposes an affine‐invariant test extending the univariate Wilcoxon signed‐rank test to the bivariate location problem. It gives two versions of the null distribution of the test statistic. The first version leads to a conditionally distribution‐free test which can be used with any sample size. The second version can be used for larger sample sizes and has a limiting χ₂² distribution under the null hypothesis. The paper investigates the relationship with a test proposed by Jan & Randles (1994). It shows that the Pitman efficiency of this test relative to the new test is equal to 1 for elliptical distributions but that the two tests are not necessarily equivalent for non‐elliptical distributions. These facts are also demonstrated empirically in a simulation study. The new test has the advantage of not requiring the assumption of elliptical symmetry which is needed to perform the asymptotic version of the Jan and Randles test.     

Distribution of order statistics from selected subsets of concomitants

For a random sample of size n from an absolutely continuous bivariate population (X, Y), let X_i:n be the i th X-order statistic and Y_[i:n] be its concomitant. We study the joint distribution of (V_s:m, W_t:n−m), where V_s:m is the s th order statistic of the upper subset {Y_[i:n], i=n−m+1,…,n}, and W_t:n−m is the t th order statistic of the lower subset {Y_[j:n], j=1,…,n−m  } of concomitants. When m=⌈_np0⌉$m=\lceil {\mathit{np}}_0\rceil$, s=⌈_mp1⌉$s=\lceil {\mathit{mp}}_1\rceil$, and t=⌈(n−m)_p2⌉$t=\lceil (n-m)p_2\rceil$, 0<_pi<1,i=0,1,2$0<p_i<1,i=0,1,2$, and n→∞$n\to \infty$, we show that the joint distribution is asymptotically bivariate normal and establish the rate of convergence. We propose second order approximations to the joint and marginal distributions with significantly better performance for the bivariate normal and Farlie–Gumbel bivariate exponential parents, even for moderate sample sizes. We discuss implications of our findings to data-snooping and selection problems.     

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.     

Property of bivariate poisson distribution and its application to stochastic processes

Assuming that the random vectors X ₁ and X ₂ have independent bivariate Poisson distributions, the conditional distribution of X ₁ given X ₁?+?X ₂?=?n is obtained. The conditional distribution turns out to be a finite mixture of distributions involving univariate binomial distributions and the mixing proportions are based on a bivariate Poisson (BVP) distribution. The result is used to establish two properties of a bivariate Poisson stochastic process which are the bivariate extensions of the properties for a Poisson process given by Karlin, S. and Taylor, H. M. (1975). A First Course in Stochastic Processes, Academic Press, New York.     

A Test of Independence Based on the Likelihood of Cut-Points

We define a test statistic C _n based on the sum of the likelihood ratio statistics for testing independence in the 2 × 2 tables defined at n sample cut-points (X _i , Y _i ). The asymptotic distribution of C _n , given the cut-points, is sum of dependent χ² variables with one degree of freedom. We use the bootstrap to obtain the distribution of C _n . We compare the performance of several tests of bivariate independence, including Pearson, Spearman, and Kendall correlations, Blum-Kiefer-Rosenblatt statistic, and C _n under several copulas and given marginal distributions.     

The ruin probabilities of a discrete time risk model with one-sided linear claim sizes and dependent risks

This article investigates the ruin probabilities of a discrete time risk model with dependent claim sizes and dependent relation between insurance risks and financial risks. The risk-free and risky investments of an insurer lead to stochastic discount factors {θ_n}_{n ? 1}. The claim sizes are assumed to follow a one-sided linear process with independent and identically distributed (i.i.d.) innovations {?_n}_{n ? 1}. The i.i.d. random pairs {(?_n, θ_n)}_{n ? 1} follow a common bivariate Sarmanov-dependent distribution. When the common distribution of the innovations is heavy tailed, we establish some asymptotic estimates for the ruin probabilities of this discrete time risk model.     

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.     

Efficient computation of location depth contours by methods of computational geometry

The concept of location depth was introduced as a way to extend the univariate notion of ranking to a bivariate configuration of data points. It has been used successfully for robust estimation, hypothesis testing, and graphical display. The depth contours form a collection of nested polygons, and the center of the deepest contour is called the Tukey median. The only available implemented algorithms for the depth contours and the Tukey median are slow, which limits their usefulness. In this paper we describe an optimal algorithm which computes all bivariate depth contours in O(n ²) time and space, using topological sweep of the dual arrangement of lines. Once these contours are known, the location depth of any point can be computed in O(log² n) time with no additional preprocessing or in O(log n) time after O(n ²) preprocessing. We provide fast implementations of these algorithms to allow their use in everyday statistical practice.     

On some moving-average processes with exponentially distributed innovations

Summary In this paper we discusse the stationary sequence of random variables which are formed from an independent identically distributed sequence, according to the moving-average model of ordern. Some properties of the process are considered. The joint bivariate exponential distribution is given, as well as the distribution of the sum.     

Improved bootstrap confidence intervals in certain toxicological experiments

The bootstrap, the jackknife, and classical methods are compared through their confidence intervals for the proportion of affected fetuses in a common type of animal experiment. Specifically, suppose that for the ith of M pregnant animals, there are x _i affected fetuses out of n _i total in the litter. The conditional distribution of x _i given n _i is sometimes modeled as binomial (n _i p _i ), where p _i is a realization from some unknown continuous density. The p _i are not observable and it is of interest in some toxicological experiments to find confidence intervals for E(p). Theory suggests that the proposed parametric bootstrap should produce higher order agreement between the nominal and actual coverage than that exhibited by the usual nonparametric bootstrap. Some simulation results provide additional evidence of this superiority of the modified parametric bootstrap over the jack-knife and classical approaches. The proposed resampling is flexible enough to handle a more general model allowing correlation between p _i and n _i .     

Reliability analysis of systems with components having two dependent subcomponents

In this article, a system that consists of n independent components each having two dependent subcomponents (A_i, B_i), i = 1, …, n is considered. The system is assumed to compose of components that have two correlated subcomponents (A_i, B_i), and functions iff both systems of subcomponents A₁, A₂, …, A_n and B₁, B₂, …, B_n work under certain structural rules. The expressions for reliability and mean time to failure of such systems are obtained. A sufficient condition to compare two systems of bivariate components in terms of stochastic ordering is also presented.     

A bivariate density estimation method based on level crossings

This article introduces a method of nonparametric bivariate density estimation based on a bivariate sample level crossing function, which leads to the construction of a bivariate level crossing empirical distribution function (BLCEDF). An efficiency function for this BLCEDF relative to the classical empirical distribution function (EDF), is derived. The BLCEDF gives more efficient estimates than the EDF in the tails of any underlying continuous distribution, for both small and large sample sizes. On the basis of BLCEDF we define a bivariate level crossing kernel density estimator (BLCKDE) and study its properties. We apply the BLCEDF and BLCKDE for various distributions and provide results of simulations that confirm the theoretical properties. A real-world example is given.