Some competitors of the mood test for the two-sample scale problem

Let X_l,…,X_n (Y_l,…,Y_m) be a random sample from an absolutely continuous distribution with distribution function F(G).A class of distribution-free tests based on U-statistics is proposed for testing the equality of F and G against the alternative that X's are more dispersed then Y's. Let 2 ? C ? n and 2 ? d ? m be two fixed integers. Let ?_c,d(X_il,…,X_ic ; Y_jl,…,X_jd)=1(-1)when max as well as min of {X_il,…,X_ic ; Y_jl,…,Y_jd } are some X_i's (Y_j's)and zero oterwise. Let S_c,d be the U-statistic corresponding to ?_c,d.In case of equal sample sizes, S22 is equivalent to Mood's Statistic.Large values of S_c,d are significant and these tests are quite efficient     

Local asymptotic normality of intermediate and central order statistics

Consider n independent random variables Z_i,…, Z_n on R with common distribution function F, whose upper tail belongs to a parametric family F(t) = F_θ(t),t ≥ x₀, where θ ∈ ? ? R ^d. A necessary and sufficient condition for the family F_θ, θ ∈ ?, is established such that the k-th largest order statistic Z_n?k+1:n alone constitutes the central sequence yielding local asymptotic normality ( LAN ) of the loglikelihood ratio of the vector (Z_n?i+1:n)¹ _i=kof the k largest order statistics. This is achieved for k = k(n)→_n→∞∞ with k/n→_n→∞ 0.

In the case of vectors of central order statistics ( Z_r:n, Z_r+1:n,…, Z_s:n ), with r/n and s/n both converging to q ∈ ( 0,1 ), it turns out that under fairly general conditions any order statistic Z_m:n with r ≤ m ≤s builds the central sequence in a pertaining LAN expansion.These results lead to asymptotically optimal tests and estimators of the underlying parameter, which depend on single order statistics only     

Development and comparative study of two near-exact approximations to the distribution of the product of an odd number of independent Beta random variables

Using the concept of near-exact approximation to a distribution we developed two different near-exact approximations to the distribution of the product of an odd number of particular independent Beta random variables (r.v.'s). One of them is a particular generalized near-integer Gamma (GNIG) distribution and the other is a mixture of two GNIG distributions. These near-exact distributions are mostly adequate to be used as a basis for approximations of distributions of several statistics used in multivariate analysis. By factoring the characteristic function (c.f.) of the logarithm of the product of the Beta r.v.'s, and then replacing a suitably chosen factor of that c.f. by an adequate asymptotic result it is possible to obtain what we call a near-exact c.f., which gives rise to the near-exact approximation to the exact distribution. Depending on the asymptotic result used to replace the chosen parts of the c.f., one may obtain different near-exact approximations. Moments from the two near-exact approximations developed are compared with the exact ones. The two approximations are also compared with each other, namely in terms of moments and quantiles.     

Asymptotic properties of bivariate random extremes

The class of limit distribution functions (d.f.s) of bivariate extremes order statistics with random sample size, which is independent of all basic random variables (r.v.s), is fully characterized. Necessary and sufficient conditions, as well as, the domains of attraction of the limit d.f.s are obtained. Furthermore, when the interrelation of the random size and the basic r.v.s is not restricted, sufficient conditions of the convergence and the forms of the limit d.f.s are deduced.     

Joint distribution of simultaneous exposures to several carcinogens in a case–control study: sample size determination

Consider a population of individuals who are free of a disease under study, and who are exposed simultaneously at random exposure levels, say X,Y,Z,… to several risk factors which are suspected to cause the disease in the populationm. At any specified levels X=x, Y=y, Z=z, …, the incidence rate of the disease in the population ot risk is given by the exposure–response relationship r(x,y,z,…) = P(disease|x,y,z,…). The present paper examines the relationship between the joint distribution of the exposure variables X,Y,Z, … in the population at risk and the joint distribution of the exposure variables U,V,W,… among cases under the linear and the exponential risk models. It is proven that under the exponential risk model, these two joint distributions belong to the same family of multivariate probability distributions, possibly with different parameters values. For example, if the exposure variables in the population at risk have jointly a multivariate normal distribution, so do the exposure variables among cases; if the former variables have jointly a multinomial distribution, so do the latter. More generally, it is demonstrated that if the joint distribution of the exposure variables in the population at risk belongs to the exponential family of multivariate probability distributions, so does the joint distribution of exposure variables among cases. If the epidemiologist can specify the differnce among the mean exposure levels in the case and control groups which are considered to be clinically or etiologically important in the study, the results of the present paper may be used to make sample size determinations for the case–control study, corresponding to specified protection levels, i.e., size α and 1–β of a statistical test. The multivariate normal, the multinomial, the negative multinomial and Fisher's multivariate logarithmic series exposure distributions are used to illustrate our results.     

On the Influence of Record Terms in the Addition of Independent Random Variables

We discuss some problems connected with the role of record values and maximal values generated by sequences of random variables X₁, X₂,…, X _n in the process of the growth of sums X₁ +···+ X_n, n = 1, 2,….     

Identifiability of distributions under competing risks and complementary risks model

Let X₁,X₂,…,X_p be p random variables with cdf's F₁(x),F₂(x),…,F_p(x)respectively. Let U = min(X₁,X₂,…,X_p) and V = max(X₁,X₂,…,X_p).In this paper we study the problem of uniquely determining and estimating the marginal distributions F₁,F₂,…,F_p given the distribution of U or of V.

First the problem of competing and complementary risks are introduced with examples and the corresponding identification problems are considered when the X₁'s are independently distributed and U(V) is identified, as well as the case when U(V) is not identified. The case when the X₁'s are dependent is considered next. Finally the problem of estimation is considered.     

Incorporating support constraints into nonparametric estimators of densities

Let f_n ^? (x) be the usual Parzen-Rosenblatt kernel estimator of the pdf f of a random variable X based on a sample X₁,…,X_n from X.In many practical applications,it is knownt hat X>c and/or X<d for given constants c and d.Additionally, one might know the values of(c)and/or f(d).“mirrorimage”and“tieddown”modifications of f_n ^?incorporate this additional information into an estimator f_n which has support [c,d].This estimatoris interpreted in a manner which allows one to use most of the known convergence properties of kernel estimates in studying the behavior of f_n.     

On characterizations based on record values and order statistics

Consider an infinite sequence of independent random variables having common continuous c.d.f. F. For 1 ⩽ i ⩽ n, let X_i:n denote the ith order statistic of the first n random variables, and let {X(n), n ⩾ 1} be the sequence of upper record values. We examine the similarities and differences between the dependence structures of the X_i:n's and the X(n)'s, with an emphasis on the latter. We present an interesting situation involving a characterization of F using the moment sequence of records. We obtain characterizations based on the properties of certain regression functions associated with order statistics, record values, and the original observations. We discuss the resemblance between some known and some new characterizations based on order statistics, record values and those based on the properties of truncated F.     

Bounds on the expected maximum

Let (X₁,X₂, …,X_n) be jointly distributed random variables. Define X_n:n = max(X₁,X₂, …,X_n).Bounds on E(X_n:n), obtained by putting constraints on the distributions and/or dependence structure of the X_i's, are surveyed.     

Minimum-Distance Estimator for Stable Exponent

Assume that X ₁, X ₂,…, X _n is a sequence of i.i.d. random variables with α-stable distribution (α ∈ (0,2], the stable exponent, is the unknown parameter). We construct minimum distance estimators for α by minimizing the Kolmogorov distance or the Cramér–von-Mises distance between the empirical distribution function G _n , and a class of distributions defined based on the sum-preserving property of stable random variables. The minimum distance estimators can also be obtained by minimizing a U-statistic estimate of an empirical distribution function involving the stable exponent. They share the same invariance property with the maximum likelihood estimates. In this article, we prove the strong consistency of the minimum distance estimators. We prove the asymptotic normality of our estimators. Simulation study shows that the new estimators are competitive to the existing ones and perform very closely even to the maximum likelihood estimator.     

Mean intergrated squared error properties and optimal kernels when estimating a diatribution function

Let X₁,X₂,… X_n be a sample of independent identically distributed (i.i.d)random variables having an unknown absolutely continuous distribution function f with density f the twofold aim of his paper consists in, firstly deriving asymptotic expressions of the mean intergrated squared error (MISE) of a kernel estimator of F when f is either assumed to be continuous everywhere or problem of finding optimal kernels in these two cases is studied in detail.     

Probability inequalities for continuous unimodal random variables with finite support

A paramecer-free Bernstein-type upper bound is derived for the probability that the sum S of n i.i.d, unimodal random variables with finite support, X₁ ,X₂,…,X_n, exceeds its mean E(S) by the positive value nt. The bound for P{S - nμ ≥ nt} depends on the range of the summands, the sample size n, the positive number t, and the type of unimodality assumed for X_i. A two-sided Gauss-type probability inequality for sums of strongly unimodal random variables is also given. The new bounds are contrasted to Hoeffding's inequality for bounded random variables and to the Bienayme-Chebyshev inequality. Finally, the new inequalities are applied to a classic probability inequality example first published by Savage (1961).     

NON-EXISTENCE OF AN EXACT β-CONFIDENCE INTERVAL FOR A SCALE PARAMETER*

Abstract There are given k (≥22) independent distributions with c.d.f.'s F(x;θ_j) indexed by a scale parameter θ_j, j = 1,…, k. Let θ_[i] (i = 1,…, k) denote the ith smallest one of θ₁,…, θ_k. In this paper we wish to show that, under some regularity conditions, there does not exist an exact β-level (0≤β1) confidence interval for the ith smallest scale parameter θ_i based on k independent samples. Since the log transformation method may not yield the desired results for the scale parameter problem, we will treat the scale parameter case directly without transformation. Application is considered for normal variances. Two conservative one-sided confidence intervals for the ith smallest normal variance and the percentage points needed to actually apply the intervals are provided.     

Spectral analysis with replicated time series

A doubly stochastic process {x(b,t);b?B,t?Z} is considered, with (B,β,P_β) being a probability space so that for each b, {X(b,t);t ? Z} is a stationary process with an absolutely continuous spectral distribution. The population spectrum is defined as f(ω) = E_B[Q(b,ω)] with Q(b,ω) being the spectral density function of X(b,t). The aim of this paper is to estimate f(ω) by means of a random sample b₁,…,b_r from (B,β,P_β). For each b₁? B, the processes X(b₁,t) are observed at the same times t=1,…,N. Thus, r time series (x(b₁,t)} are available in order to estimate f(ω). A model for each individual periodogram, which involves f(ω), is formulated. It has been proven that a certain family of linear stationary processes follows the above model In this context, a kernel estimator is proposed in order to estimate f(ω). The bias, variance and asymptotic distribution of this estimator are investigated under certain conditions.     

On some joint distributions in fluctuation theory

For non-negative integral valued interchangeable random variables v₁, v₂,…,v_n, Takács (1967, 70) has derived the distributions of the statistics ?_n' ?¹_n' ?^(c)_n and ?^(-c)_n concerning the partial sums N_r = v₁ + v₂ + ··· + v_rr = 1,…,n. This paper deals with the joint distributions of some other statistics viz., (α^(c)_n, δ^(c)_n, Z_n), (β^(c)_n, Z_n) and (β^(-c)_n, Z_n) concerning the partial sums N_r = ε₁ + ··· + ε_rr = 1,2,…,n, of geometric random variables ε₁, ε₂,…,ε_n.     

On the admissibility of c[Xbar] + d with respect to the linex loss function

Let X₁, …,X_n be a random sample from a normal distribution with mean θ and variance σ² The problem is to estimate θ with loss function L(θ,e) = v(e-θ) where v(x) = b(exp(ax)-ax-l) and where a, b are constants with b>0, a¦0. Zellner (1986), showed that [Xbar] ? σ²a/2n dominates [Xbar] and hence [Xbar] is inadmissible. The question of what values of c and d render c[Xbar]+ d admissible is studied here.     

Moments of suprema of random variables

The supremum of random variables representing a sequence of rewards is of interest in establishing the existence of optimal stopping rules. Necessary and sufficient conditions are given for existence of moments of sup_n(X_n ? c_n) and sup_n(S_n ? c_n) where X₁, X₂, … are i.i.d. random variables, S_n = X₁ + … + X_n, and c_n = (nL(n))^1/r, 0 < r < 2, L = 1, L = log, and L = log log. Following Cohn (1974), “rates of convergence” results are used in the proof.     

Record values of exponentially distributed random variables

A sequence {X_n, n≥1} of independent and identically distributed random variables with absolutely continuous (with respect to Lebesque measure) cumulative distribution function F(x) is considered. X_j is a record value of this sequence if X_j>max(X₁,…,X_j?1), j>1. Let {X_L(n), n≥0} with L(o)=1 be the sequence of such record values and Z_n,n?1=X_L(n)–X_L(n?1). Some properties of Z_n,n?1 are studied and characterizations of the exponential distribution are discussed in terms of the expectation and the hazard rate of z_n,n?1.