Characterization of the multivariate normal distribution by conditional normal distributions

Summary LetX andY be two random vectors with values in ℝ ^k and ℝ∝, respectively. IfZ=(X ^T,Y ^T) ^T is multivariate normal thenX givenY=y andY givenX=x are (multivariate) normal; the converse is wrong. In this paper simple additional conditions are stated such that the converse is true, too. Furthermore, the case is treated that the random vectorZ=(X ₁ ^T , …,X _t ^T ) ^T is splitted intot≥3 partsX ₁, …,X _t.     

A consistent test for multivariate normality based on the empirical characteristic function

LetX ₁,X ₂, …,X _n be independent identically distributed random vectors in IR ^d ,d ⩾ 1, with sample mean and sample covariance matrixS _n. We present a practicable and consistent test for the composite hypothesisH _d: the law ofX ₁ is a non-degenerate normal distribution, based on a weighted integral of the squared modulus of the difference between the empirical characteristic function of the residualsS _n ^−1/2 (X _j − ) and its pointwise limit exp (−1/2|t|²) underH _d. The limiting null distribution of the test statistic is obtained, and a table with critical values for various choices ofn andd based on extensive simulations is supplied.     

Locally minimax test of independence in elliptically symmetrical distributions with additional observations

Summary LetX=(X _ij )=(X ₁, ...,X _n )’,X’ _i =(X _i1, ...,X _ip )’,i=1,2, ...,n be a matrix having a multivariate elliptical distribution depending on a convex functionq with parameters, 0,σ. Let ϱ²=ϱ ₂ ^-2 be the squared multiple correlation coefficient between the first and the remainingp ₂+p ₃=p−1 components of eachX _i . We have considered here the problem of testingH ₀:ϱ²=0 against the alternativesH ₁:ϱ ₁ ^-2 =0, ϱ ₂ ^-2 >0 on the basis ofX andn ₁ additional observationsY ₁ (n ₁×1) on the first component,n ₂ observationsY ₂(n ₂×p ₂) on the followingp ₂ components andn ₃ additional observationsY ₃(n ₃×p ₃) on the lastp ₃ components and we have derived here the locally minimax test ofH ₀ againstH ₁ when ϱ ₂ ^-2 →0 for a givenq. This test, in general, depends on the choice ofq of the familyQ of elliptically symmetrical distributions and it is not optimality robust forQ.     

Testing whetherF is more IFR thanG

A distributionF is said to be “more IFR” than another distributionG ifG ⁻¹ F is convex. WhenF(0) =G(0) = 0, the problem of testingH ₀ :F(x) =G (θx) for someθ > 0 andx ⩾ 0, against the alternativeH _A:F is more IFR thanG, is considered in this paper. Both cases, whenG is completely specified (one-sample case) and when it is not specified but a random sample form it is available (two-sample case) are considered. The proposed tests are based onU-statistics. The asymptotic relative efficiency of the tests are compared with several other tests and the test statistics remain asymptotically normal under certain dependency assumptions. Research supported in part by a grant from the US Air Force Office of Scientific Research.     

One-Sample Bayesian Predictive Interval of Future Ordered Observations for the Pareto Distribution

Nigm et al. (2003, statistics 37: 527–536) proposed Bayesian method to obtain predictive interval of future ordered observation Y _(j) (r < j≤ n ) based on the right type II censored samples Y ₍₁₎ < Y ₍₂₎ < ... < Y _(r) from the Pareto distribution. If some of Y ₍₁₎ < ... < Y _(r-1) are missing or false due to artificial negligence of typist or recorder, then Nigm et al.’s method may not be an appropriate choice. Moreover, the conditional probability density function (p.d.f.) of the ordered observation Y _(j) (r < j ≤ n ) given Y ₍₁₎ <Y ₍₂₎ < ... < Y _(r) is equivalent to the conditional p.d.f. of Y _(j) (r < j ≤ n ) given Y _(r). Therefore, we propose another Bayesian method to obtain predictive interval of future ordered observations based on the only ordered observation Y _(r), then compares the length of the predictive intervals when using the method of Nigm et al. (2003, statistics 37: 527–536) and our proposed method. Numerical examples are provided to illustrate these results.     

Largest excess of boundary crossings for martingales

Summary LetX ₁,X ₂, ... form a sequence of martingale differences and denote byZ(a, α) = sup n (S _n −an ^α)⁺ the largest excess forS _n =X ₁ + ... +X _n crossing the boundaryan ^α. We give a sufficient condition for the finiteness ofEZ(a, α)^β which is formulated in terms of bounds forE(X _i ^{+ p} andE(|X _i |γ|X ₁, ...,X _i-1), whereα, β, γ, p are suitably related. This general result is then applied to the case of independent random variables.     

On the power function of the exact test for ther×c contingency table

Summary LetN=[n _ij ] (i=1, …,r;j=1, …,c) be the matrix of observed frequencies in anr×c contingency table fromr possibly different multinomial populations with respective probabilitiesp _i =(p _i1, …,p _ic ).Freeman andHalton have proposed an exact conditional test for the hypothesisH ₀ :p _i =(p ₁, …p _c ) of the exact test is derived. Numerical values forβ(p) were previously computed for the special case:r=3,c=2 [Bennett andNakamura, 1964].     

Asymptotic results and tests for the choice of approximative models in nonlinear two-phases regression models, heteroscedatic case

This paper is devoted to the study of the least squares estimator of f for the classical, fixed design, nonlinear model X (t _i)=f(t _i)+ε(t _i), i=1,2,…,n, where the (ε(t _i))_i=1,…,n are independent second order r.v.. The estimation of f is based upon a given parametric form. In Brodeau (1993) this subject has been studied in the homoscedastic case. This time we assume that the ε(t _i) have non constant and unknown variances σ²(t _i). Our main goal is to develop two statistical tests, one for testing that f belongs to a given class of functions possibly discontinuous in their first derivative, and another for comparing two such classes. The fundamental tool is an approximation of the elements of these classes by more regular functions, which leads to asymptotic properties of estimators based on the least squares estimator of the unknown parameters. We point out that Neubauer and Zwanzig (1995) have obtained interesting results for connected subjects by using the same technique of approximation. Received: February 1996     

A characterization of the exponential distribution by higher order gap

Summary SupposeX is a non-negative random variable with an absolutely continuous (with respect to Lebesgue measure) distribution functionF (x) and the corresponding probability density functionf(x). LetX ₁,X ₂,...,X _n be a random sample of sizen fromF andX _i,n is thei-th smallest order statistics. We define thej-th order gapg _i,j(n) asg _i,j(n)=X _i+j,n–X_i,n 1i<n, 1nn–i. In this paper a characterization of the exponential distribution is given by considering a distribution property ofg _i,j(n).     

A bound on the expected overshoot for some concave boundaries

LetX ₁,X ₂,… be i.i.d. with finite meanμ>0,S _n =X ₁+…+X _n . Forf(n)=n ^β ,c>0 we consider the stopping timesT _c =inf{n:S _n >c+f(n)} with overshootR _c =S _{T c} −(c+f(T _c )). For 0<β<1 we give a bound for sup _c≥0 ER _c in the spirit of Lorden’s well-known inequality forf=0.     

Strong consistency of the MLE under random censoring

LetX ₁, ...,X _n be an i.i.d. sample from some parametric family {_θ :θ (Θ} of densities. In the random censorship model one observesZ _i =min (X _i ,Y _i ) andδ _i =1_{ x _i ≤Y _i}, whereY _i is a censoring variable being independent ofX _i . In this paper we investigate the strong consistency ofθ _n maximizing the modified likelihood function based on (Z _i ,δ _i , 1≤i≤n. The main result constitutes an extension of Wald’s theorem for complete data to censored data. Work partially supported by the “Deutsche Forschungsgemeinschaft”.     

Linearity of regression for non-adjacent order statistics

Let X ₁,X ₂,…,X _n be a random sample from a continuous distribution with the corresponding order statistics X _1:n≤X _2:n≤…≤X _n:n. All the distributions for which E(X _{k+r: n}|X _k:n)=a X _k:n+b are identified, which solves the problem stated in Ferguson (1967). Received February 1998     

Distributions determined by conditioning on a pair of order statistics

LetX ₁,X ₂, ...,X _n (n≥3) be a random sample on a random variableX with distribution functionF having a unique continuous inverseF ⁻¹ over (a,b), −∞≤a<b≤∞ the support ofF. LetX _1:n <X _2:n <...<X _n:n be the corresponding order statistics. Letg be a nonconstant continuous function over (a,b). Then for some functionG over (a, b) and for some positive integersr ands, 1<r+1<s≤n

Locally minimax tests for a multinormal data problem

LetX be ap-normal random vector with unknown mean and unknown covariance matrix and letX be partitioned asX=(X ₍₁₎ ,X ₍₂₎ , ...,X _(r) ) whereX _(j) is a subvector of dimensionp _j such that _j=1 ^r p _j =p. We show that the tests, obtained by Dahel (1988), are locally minimax. These tests have been derived to confront H_o: =0 versusH ₁: 0 on the basis of sample of sizeN, X ₁, ..., X_N, drawn fromX andr additional samples of sizeN _j, U _i ^(j) , i=1, ..., N_j, drawn fromX ₍₁₎, ...X _(r) respectively. We assume that the (r+1) samples are independent and thatN _j>p _j forj=0, 1, ..., r (N _oN andp _op). Whenr=2 andp=2, a Monte Carlo study is performed to compare these tests with the likelihood ratio test (LRT) given by Srivastava (1985). We also show that no locally most powerful invariant test exists for this problem.     

On testing for independence against right tail increasing in bivariate models

A random variableY is right tail increasing (RTI) inX if the failure rate of the conditional distribution ofX givenY>y is uniformly smaller than that of the marginal distribution ofX for everyy0. This concept of positive dependence is not symmetric inX andY and is stronger than the notion of positive quadrant dependence. In this paper we consider the problem of testing for independence against the alternative thatY is RTI inX. We propose two distribution-free tests and obtain their limiting null distributions. The proposed tests are compared to Kendall's and Spearman's tests in terms of Pitman asymptotic relative efficiency. We have also conducted a Monte Carlo study to compare the powers of these tests.Research supported by an NSERC Canada operating grant at the University of Alberta.Part of this research was done while visiting the University of Alberta supported by the NSERC Canada grant of the first author.     

Biases,variances and covariances of raking ratio estimators for marginal and cell totals and averages of observed characteristics

Summary On the basis of a simple random sample from a population, on which a cross-classification is defined with known marginal frequenciesN _i. andN _.j , one wishes to estimate the cell frequenciesN _ij , as well as cell totalsY _ij , marginal totalsY _i. andY _.j , and the grand totalY for characteristics measured on the units. Various authors have discussed so-called raking ratio estimators, which are built up from the estimated cell values by addition. They have given the bias and variance of this estimator ofY. This paper derives biases, variances and covariances for the corresponding estimators of the cell and marginal totals and of the corresponding marginal averages.     

Optimal design of nonlinear experiments with parameter constraints

The nonlinear regression model with N observations y _i=η(x _i,θ) +ε_i, and with the parameter θ subject to q nonlinear constraints C _j (θ)=0; j=1, …,q, is considered. As an example, the spline regression with unknown nodes is taken. Expressions for the variances (variance matrices) of the LSE are discussed. Because of the complexity of these expressions, and the singularity of the variance matrix of the LSE for θ, the optimality criteria and their properties, in particular the convexity and the equivalence theorem are considered from different aspects. Also the possibility of restriction to designs with limited values of measures of nonlinearity is mentioned. Research supported by the VEGA-grant of the Slovak grant agency No. 1/7295/20.     

Moment consistency of estimators in partially linear models under NA samples

Consider the heteroscedastic regression model Y ^(j)(x _in , t _in ) = t _in β + g(x _in ) + σ _in e ^(j)(x _in ), 1 ≤ j ≤ m, 1 ≤ i ≤ n, where s_in²=f(u_in){\sigma_{in}^{2}=f(u_{in})}, (x _in , t _in , u _in ) are fixed design points, β is an unknown parameter, g(·) and f(·) are unknown functions, and the errors {e ^(j)(x _in )} are mean zero NA random variables. The moment consistency for least-squares estimators and weighted least-squares estimators of β is studied. In addition, the moment consistency for estimators of g(·) and f(·) is investigated.