Test of Association Between Two Ordinal Variables While Adjusting for Covariates

We propose a new set of test statistics to examine the association between two ordinal categorical variables X and Y after adjusting for continuous and/or categorical covariates Z. Our approach first fits multinomial (e.g., proportional odds) models of X and Y, separately, on Z. For each subject, we then compute the conditional distributions of X and Y given Z. If there is no relationship between X and Y after adjusting for Z, then these conditional distributions will be independent, and the observed value of (X, Y) for a subject is expected to follow the product distribution of these conditional distributions. We consider two simple ways of testing the null of conditional independence, both of which treat X and Y equally, in the sense that they do not require specifying an outcome and a predictor variable. The first approach adds these product distributions across all subjects to obtain the expected distribution of (X, Y) under the null and then contrasts it with the observed unconditional distribution of (X, Y). Our second approach computes "residuals" from the two multinomial models and then tests for correlation between these residuals; we define a new individual-level residual for models with ordinal outcomes. We present methods for computing p-values using either the empirical or asymptotic distributions of our test statistics. Through simulations, we demonstrate that our test statistics perform well in terms of power and Type I error rate when compared to proportional odds models which treat X as either a continuous or categorical predictor. We apply our methods to data from a study of visual impairment in children and to a study of cervical abnormalities in human immunodeficiency virus (HIV)-infected women. Supplemental materials for the article are available online.     

Probability density function of quotient of order statistics from the pareto,power and weibull distributions

This paper deals with the probability density functions of quotient of order statistics. We use the Mellin transform technique, to find the distribution of the quotient Z= X/Xwhere X.,X(i < j) are the ith and jth order statistics from the Pareto, Power and Weibull distributions     

Confidence intervals on p(y < x) for small sample sizes

When F = G^a, confidence intervals are derived and presented in graphs for p = P(Y < X), when X and Y are independent and the sample sizes are at most 25. Also, it is demonstrated via a monte carlo simulation that this is a robust procedure, when the distributions of X and Y differ by a location parameter.     

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.     

Confidence intervals for discrete distributions based on a single observation

Confidence intervals are developed for the mode of a discrete unimodal distribution in the case where only a single observation is available. These intervals are centered on either the observation, X, or a weighted average of X with a constant, b, chosen by the investigator. Intervals are derived for nonrestricted unimodal distributions, for unimodal distributions with a symmetry property, and for a family of two-sided truncated geometric distributions.     

Some bivariate uniform distributions

Bivariate uniform distributions with dependent components are readily derived by distribution function transformations of the components of non-uniform dependent continuous bivariate random variables (X,Y). Contour plots of joint density functions show the various, and varying, forms of dependence which can arise from different distributional forms for (X,Y) and aids the choice of bivariate uniform distributions as empirical models.     

Bayesian Estimation of Survival Functions under Stochastic Precedence

When estimating the distributions of two random variables, X and Y, investigators often have prior information that Y tends to be bigger than X. To formalize this prior belief, one could potentially assume stochastic ordering between X and Y, which implies Pr(X < or = z) > or = Pr(Y < or = z) for all z in the domain of X and Y. Stochastic ordering is quite restrictive, though, and this article focuses instead on Bayesian estimation of the distribution functions of X and Y under the weaker stochastic precedence constraint, Pr(X < or = Y) > or = 0.5. We consider the case where both X and Y are categorical variables with common support and develop a Gibbs sampling algorithm for posterior computation. The method is then generalized to the case where X and Y are survival times. The proposed approach is illustrated using data on survival after tumor removal for patients with malignant melanoma.     

A property of the yule distribution and its applications

The relationship Y = RX between two random variables X and Y, where R is distributed independently of X in (0, l), is known to have important consequences in different fields such as income distribution analysis, Inventory decision models, etc.

In this paper it is shown that when X and Y are discrete random variables, relationships of similar nature lead to Yule-type distributions. The implications of the results are studied in connection with problems of income underreporting and inventory decision making.     

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.     

Moments for a ratio of correlated gamma variates

This paper considers the finite integral moments for the ratio, R = X/Y, where X and Y re correlated gamma distributed variables. An analytical and numerical comparison is given for two classes of underlying bivariate gamma distributions. It is shown that the two bivariate gamma structures provide indentical experessions for the mth unadjussted moment, E(R^m), if and only if either of the following conditions hold : 1) X and Y are uncorrelated of 2) m=1. A numerical evaluation is performed to determine the extent that the two methods differ whenever the variables are correlated     

Approximate Identifiability, Moments and Censored Data

It is well known that the joint distribution of a pair of random variables ( X,Y ) is not identifiable on the basis of the joint distribution of the function (min ( X,Y ), 1[ X < Y ]). This paper introduces the concept of approximate identifiability and studies its relevance to the function (min ( X,Y ), Y ). It shows that the distribution of ( X,Y ) is approximately identifiable on the basis of the distribution of (min ( X,Y ), Y ). The identification is explicitly executed by a method of moments. The method is applied to the analysis of censored distributions arising in the theory of clinical trials and is compared to the standard method of Kaplan and Meier.     

A TEST OF INDEPENDENCE FOR BIVARIATE SYMMETRIC STABLE DISTRIBUTIONS

A test of independence in symmetric bivariate stable distributions is constructed using the empirical characteristic function as a test statistic. A particular class of distributions considered in detail is X = U + V, Y=V+W where U, V, and W are mutually independent symmetric stable distributions with the same index.     

Some robust tests of independence in symmetrical multivariate distributions

Let X =(x)^ij=(11¹, …, X,)_T, i = l, …n, be an n X random matrix having multivariate symmetrical distributions with parameters μ, Σ. The p-variate normal with mean μ and covariance matrix is a member of this family. Let be the squared multiple correlation coefficient between the first and the succeeding p₁ components, and let p² = + be the squared multiple correlation coefficient between the first and the remaining p₁ + p₂ =p – 1 components of the p-variate normal vector. We shall consider here three testing problems for multivariate symmetrical distributions. They are (A) to test p² =0 against; (B) to test against =0, 0; (C) to test against p² =0, We have shown here that for problem (A) the uniformly most powerful invariant (UMPI) and locally minimax test for the multivariate normal is UMPI and is locally minimax as p² 0 for multivariate symmetrical distributions. For problem (B) the UMPI and locally minimax test is UMPI and locally minimax as for multivariate symmetrical distributions. For problem (C) the locally best invariant (LBI) and locally minimax test for the multivariate normal is also LBI and is locally minimax as for multivariate symmetrical distributions.     

SOME CHARACTERIZATIONS OF THE CAUCHY DISTRIBUTION

If X and Y are independent standard Cauchy random variables then (i) Y and (X+Y)/(1-Xu) are independent, (ii) X and (X + Y)/(1 -XU) are identically distributed, and (iii) X and 2X/(1-X²) are identically distributed. Each of these three properties is shown to characterize the Cauchy distribution among absolutely continuous distributions. Some related uniform characterizations are discussed.     

Attribute control chart for some popular distributions

The binomial distribution is often used to display attribute control data. In this paper, a statistical model is settled for attribute control chart under truncated life test. By Burr X & XII, inverse Gaussian (IG), and exponential lifetime-truncated distributions, a Shewhart-type attribute control chart is built to display the data. The performance of attributed control chart constructed on truncated life test is evaluated by average run length, which compares the performance of all distributions. Our study arranges that IG is better distribution among all.     

Models for residual time to AIDS

The distributions of the time from Human Immunodeficiency Virus (HIV) infection to the onset of Acquired Immune Deficiency Syndrome (AIDS) and of the residual time to AIDS diagnosis are important for modeling the growth of the AIDS epidemic and for predicting onset of the disease in an individual. Markers such as CD4 counts carry valuable information about disease progression and therefore about the two survival distributions. Building on the framework set out by Jewell and Kalbfleisch (1992), we study these two survival distributions based on stochastic models for the marker process (X(t)) and a marker-dependent hazard (h()). We examine various plausible CD4 marker processes and marker-dependent hazard functions for AIDS proposed in recent literature. For a random effects plus Brownian motion marker process X(t)=(a+bt+BM(t))⁴, where a has a normal distribution, b<0 is an unknown parameter and BM(t) is Brownian motion, and marker-dependent hazard h(X(t)), we prove that, given CD4 cell count X(t), the residual time to AIDS distribution does not depend on the time since infection t. Using simulation and numerical integration, we find the marginal incubation period distribution, the marginal hazard and the residual time distribution for several combinations of marker processes and marker-dependent hazards. An example using data from the Multicenter AIDS Cohort Study is given. A simple regression model relating the cube root of residual time to AIDS to CD4 count is suggested.     

Proportions,sums and ratios

Data that are proportions arise frequently in all areas of the sciences and engineering. In this paper, the exact distributions of R=X+Y and W=X/(X+Y) and the corresponding moment properties are derived when X and Y are proportions and arise from the most flexible bivariate beta distribution known to date. The associated estimation procedures are developed. Finally, an application is illustrated to compositional data of lavas from Skye.     

Estimation and testing procedures for the reliability functions of a general class of distributions

We consider here the general class of distributions proposed by Sankaran and Gupta (2005) by zeroing in on two measures of reliability, R(t) = P(X > t) and P = P(X > Y). Thereafter, we develop point estimation for R(t) and ‘P’ and develop uniformly minimum variance unbiased estimators (UMVUES). Then we derive testing procedures for the hypotheses related to different parametric functions. Finally, we compare the results using the Monte Carlo simulation method. Using real data set, we illustrate the procedure clearly.     

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.     

Correlation analysis of ordered symmetrically dependent observations and their concomitants of order statistics

Given two jointly observed random vectors Y and Z of the same dimension, let Y be a reordered version of Y and Z the resulting vector of concomitants of order statistics. When X is a covariate of interest, also jointly observed with Y, the authors obtain the joint covariance structure of (X, y, Z) and the related correlation parameters explicitly, under the assumption that the vector (X, Y, Z) is normal and that its joint covariance structure is permutation symmetric. They also discuss extensions to elliptically contoured distributions.