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.     

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 a new transformation to normality

A Gaussian approximation to the distribution of the nonnegative random variable Y is developed using the Wilson and Hilferty (1931) approach. This approximation uses the symmetrizing transformation ((Y + b)/k₁)^h where k₁ is the first moment of Y and h and b are determined from the first three cumulants of Y. The approximation is illustrated in the case which Y is a non-central chi-square, where numerical evaluations indicate that the new transformation is an improvement over existing ones, especially for small values of k₁.     

An alternative proof for a version of cochran's theorem under elliptical settings

For an elliptically contoured n × p random matrix Y with mean μ and covariance proportional to ∑_Y, the necessary and sufficient conditions, under which (Y?μ)′W(Y?μ) with nonnegative definite W is generalized Wishart distributed, are obtained by using the higher moments of Y. This version of Cochran's theorem is general as the assumptions on ∑_Y=A?∑ with nonnegative definite A and ∑, P(Y?μ)=0, and P(Y≠μ)<1 have been relaxed. An example on two way balanced mixed models is given for illustration ot our main results.     

A NOTE ON THE MULTIVARIATE LINEAR MODEL WITH CONSTRAINTS ON THE DEPENDENT VECTOR

It is shown that the least squares estimators of B and Σ in the multivariate linear model {E Y _i= X ₁ B , D ( Y _i) =Σ, 1 ≤ i ≤ n , Y ₁ Y _n uncorrelated} subject to the constraints Y _i M = X _i N are just the usual least squares estimators = ( X'X )-¹ X'Y and ΣC = 1/n( Y-X )( Y-X ) in the unconstrained model where Σ has full rank. Tests of hypotheses concerning B are discussed for situations in which each Y _i has a multivariate normal distribution, and examples of the applicability of the model reviewed.     

Some inference results on pr(x < y) in the bivariate exponential model

This paper provides three different estimators for Pr(X < Y) when X and Y have a bivariate exponential distribution. The asymptotic variances of the three estimators are also derived. A test for the equality of the means of X and Y and confidence limits for the difference of the two means are presented. Our results are directly applicable in a reliability context with underlying bivariate exponential distribution.     

Vec and vech operators for matrices,with some uses in jacobians and multivariate statistics

The vec of a matrix X stacks columns of X one under another in a single column; the vech of a square matrix X does the same thing but starting each column at its diagonal element. The Jacobian of a one-to-one transformation X → Y is then ∣∣?(vecX)/?(vecY) ∣∣ when X and Y each have functionally independent elements; it is ∣∣ ?(vechX)/?(vechY) ∣∣ when X and Y are symmetric; and there is a general form for when X and Y are other patterned matrices. Kronecker product properties of vec(ABC) permit easy evaluation of this determinant in many cases. The vec and vech operators are also very convenient in developing results in multivariate statistics.     

On some properties of bivariate weighted distributions

Let X^w and Y^w be weighted random variables arising from the distribution of (X,Y). We explore implications of independence of X and Y on the dependence structure of (X^w, Y^w). We also show that when X and Y are independent and the weight function is symmetric, identical distribution of X^w and Y^w implies that of X and Y. We discuss application of these results to the study of a renewal process.     

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.     

Nonparametric estimation of survival functions for doubly censored observations when the lifetime hazard rate is proportionally related to the censoring hazard rates

Let T, X and Y be non-negative random variables, where T is the time of occurrence of an event of interest, X and Y being the lefl and right censoring variables respectively.

In this paper we propose a nonparametric estimator of the survival function, ST, when T, X and Y are supposed to be independent and their corresponding hazard rates are proportionally related. In this way, our results extend Ebrahimi's work (1985) to the doubly censored data case.     

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     

Gamma variates of fractional shape as functioials of a homogeneous multidimensional poisson process

If events are scattered in Rⁿ in accordance with a homogeneous Poisson process and if X is the location of the event with minimal [d]l^P norm, then in the case p = n the n^th absolute powers of the coordinates of X form a sample of size n from a gamma distribution with shape parameter 1/n. In an age of parallel computing, this fact may lead to some attractive simulation methods. One possibility is to generate R = [d]X[d] and U = Y/[d]X[d] independently, perhaps by setting U = Y/[d]Y[d] where Y has any p.d.f. which is a function only of ¦Y¦. We consider for example Y having the uniform distribution in an l^P ball.     

Bayesian tests for matched proportions

A sample of n subjects is observed in each of two states, S₁-and S₂. In each state, a subject is in one of two conditions, X or Y. Thus, a subject may be recorded as showing a change if its condition in the two states is ‘Y,X’ or ‘X,Y’ and, otherwise, the condition is unchanged. We consider a Bayesian test of the null hypothesis that the probability of an ‘X,Y’ change exceeds that of a ‘Y,X’ change by amount k_O. That is, we develop the posterior distribution of k_O, the difference between the two probabilities and reject the null hypothesis if k lies outside the appropriate posterior probability interval. The performance of the method is assessed by Monte Carlo and other numerical studies and brief tables of exact critical values are presented     

Approximations of distributions for some standardized partial sums in sequential analysis

In sequential analysis it is often necessary to determine the distributions of √t Y _t and/or √a Y _t , where t is a stopping time of the form t = inf{ n ≥ 1 : n+S_n+ξ_n> a }, Y _n is the sample mean of n independent and identically distributed random variables (iidrvs) Y_i with mean zero and variance one, S_n is the partial sum of iidrvs X_i with mean zero and a positive finite variance, and { ξ_n } is a sequence of random variables that converges in distribution to a random variable ξ as n →∞ and ξ_n is independent of ( X_n+1, Y_n+1), (X_n+2, Y_n+2), . . . for all n ≥ 1. Anscombe's (1952) central limit theorem asserts that both √t Y _t and √a Y _t are asymptotically normal for large a , but a normal approximation is not accurate enough for many applications. Refined approximations are available only for a few special cases of the general setting above and are often very complex. This paper provides some simple Edgeworth approximations that are numerically satisfactory for the problems it considers.     

CORRELATIONS AND CHARACTERIZATIONS OF THE UNIFORM DISTRIBUTION

Two characterizations of the uniform distribution on a suitable compact space are proved. These characterizations are applied to a number of particular examples of which the most interesting is the following: if X , Y and Z are independent n-vectors whose components are independent and identically distributed within a vector, then the pairwise independence of the product moment correlation coefficients between X , Y and Z implies that these vectors are normally distributed.     

Some remarks on classical and bayesian reliability estimation of binomial and poisson distributions

The problems of estimating the reliability function and Pr{X₁+...+X_k ≤ Y} are considered. The random variables X’s and Y are assumed to follow binomial and Poisson distributions. Classical estimators available in the literature are discussed and Bayes estimators are derived. In order to obtain the estimators of these parametric functions, the basic role is played by the estimators of factorial moments of the two distributions.     

Non-parametric Regression with Dependent Censored Data

Abstract.  Let ( X _i , Y _i ) ( i = 1 ,…, n ) be n replications of a random vector ( X , Y  ), where Y is supposed to be subject to random right censoring. The data ( X _i , Y _i ) are assumed to come from a stationary α -mixing process. We consider the problem of estimating the function m ( x ) = E ( φ ( Y ) |  X = x ), for some known transformation φ . This problem is approached in the following way: first, we introduce a transformed variable     , that is not subject to censoring and satisfies the relation     , and then we estimate m ( x ) by applying local linear regression techniques. As a by-product, we obtain a general result on the uniform rate of convergence of kernel type estimators of functionals of an unknown distribution function, under strong mixing assumptions.     

Tolerance regions for random vectors and best linear predictors

In the case that vectors X and Y have a joint multivariate normal distribution, tolerance regions are found for the best linear predictor of Y using X if samples are used to estimate the regression coeffierante. Tolerance regions are also found for Y. In addition, simultaneous tolerance intervals for all linear functions of Y or of the best linear predictor of Y using X are found.     

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.     

Score tests in the two - way layout of counts

For the balanced two-way layout of a count response variable Y classified by fixed or random factors A and B, we address the problems of (i) testing for individual and interactive effects on Y of two fixed factors, and (ii) testing for the effect of a fixed factor in the presence of a random factor and conversely. In case (i), we assume independent Poisson responses with µ_ij= E(Y| A=i,B=j) = α_iβ_jγ_ij corresponding respectively to the multiplicative

interactive and non-interactive cases. For case (ii) with factor A random, we derive a multivariate gamma-Poisson model by mixing on the random variable associated with each level of A. In each case Neyman C(α) score tests are derived. We present simulation results,and apply the interaction test to a data set, to evaluate and compare the size and power of the score test for interaction between two fixed factors, the competing Poisson-based likelihood ratio test, and the F-tests based on the assumptions that √Y+1 or log(Y+1) are approximately normal. Our results provide strong evidence that the normal-theory based F-tests typically are very far from nominal size, and that the likelihood ratio test is somewhat more liberal than the score test.