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.     

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.     

Logistic regression with a partially observed covariate

We present results of a Monte Carlo study comparing four methods of estimating the parameters of the logistic model logit (pr (Y = 1 | X, Z)) = α₀ + α ₁ X + α ₂ Z where X and Z are continuous covariates and X is always observed but Z is sometimes missing. The four methods examined are 1) logistic regression using complete cases, 2) logistic regression with filled-in values of Z obtained from the regression of Z on X and Y, 3) logistic regression with filled-in values of Z and random error added, and 4) maximum likelihood estimation assuming the distribution of Z given X and Y is normal. Effects of different percent missing for Z and different missing value mechanisms on the bias and mean absolute deviation of the estimators are examined for data sets of N = 200 and N = 400.     

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.     

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 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.     

Non-parametric Estimation of the Residual Distribution

Consider a heteroscedastic regression model Y = m ( X ) +σ( X )ε, where the functions m and σ are "smooth", and ε is independent of X . An estimator of the distribution of ε based on non-parametric regression residuals is proposed and its weak convergence is obtained. Applications to prediction intervals and goodness-of-fit tests are discussed.     

An Efficient Simulation Method for the Computation of a Class of Conditional Expectations

Suppose that the random vector X and the random variable Y are jointly continuous. Also suppose that an observation x of X can be easily simulated and that the probability density function of Y conditional on X = x is known. The paper presents an efficient simulation-based algorithm for estimating E{ g ( X , Y ) | h ( X , Y ) = r } where g and h are real-valued functions. This algorithm is applicable to time series problems in which X = ( X ₁, . . . , X _n−1) and Y = X_n where { x_t } is a discrete time stochastic process for which ( X₁ , . . . , X_n ) is a continuous random vector. A numerical example from time series analysis illustrates the algorithim, for prediction for an ARCH(1) process.     

A closed form solution for tre least squares regression problem with linear inequality constraints

In this paper we consider a linear model Y = Xβ+e with linear inequality constraints R'β≥r, where X and R are known and full column rank matrices. The closed form of the inequality constrained least squares (ICLS) estimator is given. We provide two examples which illustrate the use of this closed form in the computation of estimates.     

Effect on probabilities and quantiles of adding a quantity with small variance

This paper gives simple approximations for the distribution function and quantiles of the sum X + Y when X is a continuous variable and Y is an independent variable with variance small compared to that of X . The approximations are based around the distribution function or quantiles of X and require only the first two or three moments of Y to be known. Example evaluations with X having a normal, Student's t or chi-squared distribution suggest that the approximations are good in unbounded tail regions when the ratio of variances is less than 0.2.     

Measurement error in regression analysis

Consider the linear regression model Y = Xθ+ ε where Y denotes a vector of n observations on the dependent variable, X is a known matrix, θ is a vector of parameters to be estimated and e is a random vector of uncorrelated errors. If X'X is nearly singular, that is if the smallest characteristic root of X'X s small then a small perurbation in the elements of X, such as due to measurement errors, induces considerable variation in the least squares estimate of θ. In this paper we examine for the asymptotic case when n is large the effect of perturbation with regard to the bias and mean squared error of the estimate.     

Non-parametric tests for distributional treatment effect for randomly censored responses

Summary.  For a binary treatment ν =0, 1 and the corresponding 'potential response' Y ⁰ for the control group ( ν =0) and Y ¹ for the treatment group ( ν =1), one definition of no treatment effect is that Y ⁰ and Y ¹ follow the same distribution given a covariate vector X . Koul and Schick have provided a non-parametric test for no distributional effect when the realized response (1− ν ) Y ⁰+ ν Y ¹ is fully observed and the distribution of X is the same across the two groups. This test is thus not applicable to censored responses, nor to non-experimental (i.e. observational) studies that entail different distributions of X across the two groups. We propose ' X -matched' non-parametric tests generalizing the test of Koul and Schick following an idea of Gehan. Our tests are applicable to non-experimental data with randomly censored responses. In addition to these motivations, the tests have several advantages. First, they have the intuitive appeal of comparing all available pairs across the treatment and control groups, instead of selecting a number of matched controls (or treated) in the usual pair or multiple matching. Second, whereas most matching estimators or tests have a non-overlapping support (of X ) problem across the two groups, our tests have a built-in protection against the problem. Third, Gehan's idea allows the tests to make good use of censored observations. A simulation study is conducted, and an empirical illustration for a job training effect on the duration of unemployment is provided.     

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.     

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.     

Monotonicity of certain multivariate normal probabilities

Let (X, Y) have a (p+q)-dimensional normal distribution and let C, K be convex symmetric sets of dimensions p, q respectively. Under certain restrictions on the mean vector it is shown that P(X ε C, Y ε K) is a monotonically increasing function of the first canonical correlation coefficient between X and Y, provided the remaining coefficients are zero.     

Adjusted Kuk's model using two non sensitive characteristics unrelated to the sensitive characteristic

In this paper, we adjust the Kuk (1990 Kuk, A.Y.C. (1990). Asking sensitive questions indirectly. Biometrika 77(2):436–438.[Crossref], [Web of Science ®] , [Google Scholar]) model for both protection and efficiency by making use of proportions of two non sensitive characteristics which are unrelated to the main sensitive characteristic of interest. Various situations, where the proportions of the two non sensitive characteristics in the population of interest are known and that when these proportions are unknown, have been investigated. We compared the adjusted model and Kuk's model through a simulation study from both the protection and efficiency points of view.     

Linear Models in a general parametric form

The linear model Y - N(Xb, σ²∑) with a restriction R'b = M'u + c is considered, where X, R, M, ∑ and c are known. Explicit formulae are obtained for the best linear unbiased estimator of K'b, for the F-test of the hypothesis K'b = W'v + a, and for the simultaneous confidence intervals of the parameters K′_i b' s, where K = [K₁,K₂,…K_s], w, and a are known, none of the matrices X, ∑, R, M, K, and W is required to have full ranks, and the design X can be one - or multi-way,complete or incomplete, balanced or not balanced, connected or disconnected.     

Orthogonal Arrays for the Estimation of Global Sensitivity Indices Based on ANOVA High-Dimensional Model Representation

We present an efficient new approach to estimate global sensitivity indices based on ANOVA high-dimensional model representation. The method makes use of the properties of orthogonal arrays and extend the results of Wang et al. (2010 Wang , X. , Tang , Y. , Chen , X. , Zhang , Y. ( 2010 ). Design of experiment in global sensitivity analysis based on ANOVA high-dimensional model representation . Communications in Statistics—Simulation and Computation 39 : 1183 – 1195 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) to a more general situation. A theorem is given to illustrate that both the non influential and significant sensitivity indices are well estimated. This new method offers higher sampling efficiency than those of Sobol and Saltelli. We test its performance on two different models, which are widely used in engineering and statistical areas.     

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.     

Computation of the Generalized F Distribution

Exact expressions for the cumulative distribution function of a random variable of the form ( α ₁ X ₁+ α ₂ X ₂)/ Y are given where X ₁, X ₂ and Y are independent chi-squared random variables. The expressions are applied to the detection of joint outliers and Hotelling's mis-specified T ² distribution.