Multivariate multistage methodologies for simultaneous all pairwise comparisons

We consider the problem of constructing simultaneous fixed-width confidence intervals for all pairwise treatment differences μ₁−μ _J , in the presence ofk(≥2) independent populationsN _p (μ₁,Σ), 1≤i≠j≤k. Appropriate purely sequential, accelerated sequential and three-stage sampling strategies have been developed and variousfirst-order asymptotic properties are then derived when Σ _pxp is completely unknown, but positive definite (p.d.). In the two special cases when the largest component variance in Σ is a known multiple of one of the variances or Σ=σ² H where σ(>0) is unknown, butH _pxp is known and p.d., the original multistage sampling strategies are specialized. Under such special circumstances, associatedsecond-order characteristics are then developed. It is to be noted that our present formulation and the methodologies fill important voids in the context of multivariate multiple comparisons which is a challenging area that has not yet been fully explored. Moderate sample performances of the proposed techniques were very encouraging and detailed remarks on these were included in Mukhopadhyay and Aoshima (1997).     

A note on the maximum deviation of the scale-contaminated normal to the best normal distribution

In this paper we consider the case of the scale-contaminated normal (mixture of two normals with equal mean components but different component variances: (1−p)N(μ,σ²)+pN(μ,τ²) with σ and τ being non-negative and 0≤p≤1). Here is the scale error and p denotes the amount with which this error occurs. It's maximum deviation to the best normal distribution is studied and shown to be montone increasing with increasing scale error. A closed-form expression is derived for the proportion which maximizes the maximum deviation of the mixture of normals to the best normal distribution. Implications to power studies of tests for normality are pointed out. Received May 2001     

Distribution of correlation coefficient under mixture normal model

The present paper obtains the nonnull distribution of the product moment correlation coefficient r when sample is drawn from a mixture of two bivariate Gaussian distributions. The moments of 1−r ² have been used to derive the nonnull density of r. Received September 2000     

Multi-stage point and interval estimation of the largest mean ofK normal populations and the associated second-order properties

Summary We havek independent normal populations with unknown meansμ ₁, …,μ _k and a common unknown varianceσ ². Both point and interval estimation procedures for the largest mean are proposed by means of sequential and three-stage procedures. For the point estimation problem, we require that the maximal risk be at mostW, a preassigned positive number. For the other problem, we wish to construct a fixed-width confidence interval having the confidence coefficient at least 1-α, a preassigned number between zero and one. Asymptotic second order expansions are provided for various characteristics, such as average sample size, associated risks etc., for the suggested multi-stage estimation procedures.     

Estimating a restricted normal mean

Let X ₁, X ₂, ..., X _n be a random sample from a normal distribution with unknown mean μ and known variance σ ². In many practical situations, μ is known a priori to be restricted to a bounded interval, say [−m, m] for some m > 0. The sample mean , then, becomes an inadmissible estimator for μ. It is also not minimax with respect to the squared error loss function. Minimax and other estimators for this problem have been studied by Casella and Strawderman (Ann Stat 9:870–878, 1981), Bickel (Ann Stat 9:1301–1309, 1981) and Gatsonis et al. (Stat Prob Lett 6:21–30, 1987) etc. In this paper, we obtain some new estimators for μ. The case when the variance σ ² is unknown is also studied and various estimators for μ are proposed. Risk performance of all estimators is numerically compared for both the cases when σ ² may be known and unknown.     

Estimates of technical inefficiency in stochastic frontier models with panel data: generalized panel jackknife estimation

Estimates of technical inefficiency based on fixed effects estimation of the stochastic frontier model with panel data are biased upward. Previous work has attempted to correct this bias using the bootstrap, but in simulations the bootstrap corrects only part of the bias. The usual panel jackknife is based on the assumption that the bias is of order T ⁻¹ and is similar to the bootstrap. We show that when there is a tie or a near tie for the best firm, the bias is of order T ^−1/2, not T ⁻¹, and this calls for a different form of the jackknife. The generalized panel jackknife is quite successful in removing the bias. However, the resulting estimates have a large variance.     

Two families of kurtosis measures

Two families of kurtosis measures are defined as K ₁(b)=E[ab ^−|z|] and K ₂(b)=E[a(1−|z|^b)] where z denotes the standardized variable and a is a normalizing constant chosen such that the kurtosis is equal to 3 for normal distributions. K ₂(b) is an extension of Stavig's robust kurtosis. As with Pearson's measure of kurtosis β₂=E[z ⁴], both measures are expected values of continuous functions of z that are even, convex or linear and strictly monotonic in ℜ⁻ and in ℜ⁺. In contrast to β₂, our proposed kurtosis measures give more importance to the central part of the distribution instead of the tails. Tests of normality based on these new measures are more sensitive with respect to the peak of the distribution. K ₁(b) and K ₂(b) satisfy Van Zwet's ordering and correlate highly with other kurtosis measures such as L-kurtosis and quantile kurtosis. RID="*" ID="*"  The authors thank the referees for their insightful comments that significantly improved the clarity of the article.     

A minimal characterization of the covariance matrix

Summary LetX be ak-dimensional random vector with mean vectorμ and non-singular covariance matrix Σ. We show that among all pairs (a, Δ),a ∈ IR ^k , Δ ∈ IR ^k×k positive definite and symmetric andE(X−a)′ Δ⁻¹(X−a)=k, (μ, Σ) is the unique pair which minimizes det Δ. This motivates certain robust estimators of location and scale. Research supported by the Nuffield Foundation.     

sk−p Fractional factorial designs in sb blocks

Draper and Guttman (1997) shows that for basic 2^k−p designs, p≥0, k − p replicates of blocks designs of size two are needed to estimate all the usual (estimable) effects. In this work, we provide an algebraic formal proof for the two-level blocks designs results and present results applicable to the general case; that is, for the case of s ^k factorial (p=0) or s ^k−p fractional factorial (p >0) designs in s ^b blocks, where 0<b<k− p, at least replicates are needed to clear up all possible effects. Through the theoretical development presented in this work, it can provide a clearer view on why those results would hold. We will also discuss the estimation equations given in Draper and Guttman (1997).  Research supported in part by the National Science Council of Taiwan, R.O.C., Grant No. NSC 89-2118-M110-010. Acknowledgement. The authors would like to thank the referee for very helpful comments.     

Choice of optimal initial designs in sequential experiments

Combined-optimal designs (Li and Lin, 2003) are obviously the best choices for the initial designs if we partition the experiment into two parts with equal size to obtain some information about the process, especially for the case not considering the blocking factor. In this paper, the definition of combined-optimal design is extended to the case when blocking factor is significant, and this new class of designs is called blocked combined-optimal designs. Some general results are obtained which relate 2^k−p_III initial designs with their complementary designs when , where n=2^k−p. By applying these results, we are able to characterize 2^k−p_III combined-optimal designs or blocked combined-optimal designs in terms of their complementary designs. It is also proved that both 2^k−p_III combined-optimal and blocked combined-optimal designs are not minimum aberration designs when and n−1−k > 2. And some combined-optimal and blocked combined-optimal designs with 16 and 32 runs are constructed for illustration. 2000 Mathematics Subject Classifications: 62K15, 62K05     

Tests of significance for the average square deviation from target in a finite lot

Consider lots of discrete items 1, 2, …,N with quality characteristicsx ₁,x ₂, …,x _N . Leta be a target value for item quality. Lot quality is identified with the average square deviation from target per item in the lot (lot average square deviation from target). Under economic considerations this is an appropriate lot quality indicator if the loss respectively the profit incurred from an item is a quadratic function ofx _i −a. The present paper investigates tests of significance on the lot average square deviationz under the following assumptions: The lot is a subsequence of a process of production, storage, transport; the random quality characteristics of items resulting from this process are i.i.d. with normal distributionN(μ, σ ²); the target valuea coincides with the process meanμ.     

Minimax estimators andΓ-minimax estimators for a bounded normal mean under the lossl p (θ, d)=|θ-d|p

Summary Let the random variableX be normal distributed with known varianceσ ²>0. It is supposed that the unknown meanθ is an element of a bounded intervalΘ. The problem of estimatingθ under the loss functionl _p (θ, d)=|θ-d| ^p p≥2 is considered. In case the length of the intervalθ is sufficiently small the minimax estimator and theΓ(β, τ)-minimax estimator, whereΓ(β, τ) represents special vague prior information, are given.     

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.     

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

A generalized maximum likelihood characterization of the normal distribution

LetP be a probability measure on ℝ andI _x be the set of alln-dimensional rectangles containingx. If for allx ∈ ℝⁿ and θ ∈ ℝ the inequality holds,P is a normal distributioin with mean 0 or the unit mass at 0. The result generalizes Teicher’s (1961) maximum likelihood characterization of the normal density to a characterization ofN(0, σ²) amongall distributions (including those without density). The m.l. principle used is that of Scholz (1980).     

A geometric property of indifference curves in the (mean-standard deviation) plane

We prove that the slope of the indifference curve through point (m σ) of the (meanstandard deviation) plane must be smaller than the ratio (1+m)/σ. We show that the indifference curves corresponding to a quadratic utility function do actually satisfy that inequality; and prove, by the way, that these curves are not quarters of circles (as it is usually asserted), but smaller parts of them.

---

Riassunto Si prova che, nel piano (media-scarto quadratico medio), l'inclinazione della curva d'indifferenza per il punto (m, σ) deve essere minore del rapporto (1+m)/gs. Si mostra che le curve d'indifferenza corrispondenti ad una funzione d'utilià quadratica soddisfano di fatto tale limitazione; e si prova, a questo proposito, che tali curve non sono — come generalmente asserito — quarti di circonferenze, ma archi più piccoli delle medesime.

---

---

The second part of this article was stimulated by an observation of E. Castagnoli. The author gratefully acknowledges this to him.     

A method to obtain new copulas from a given one

Given a strictly increasing continuous function φ from [0, 1] to [0, 1] and its pseudo-inverse φ^[−1], conditions that φ must satisfy for C_φ(x₁, . . . ,x_n)=φ^[−1](C(φ(x₁), . . . ,φ(x_n))) to be a copula for any copula C are studied. Some basic properties of the copulas obtained in this way are analyzed and several examples of generator functions φ that can be used to construct copulas C_φ are presented. In this manner, a method to obtain from a given copula C a variety of new copulas is provided. This method generalizes that used to construct Archimedean copulas in which the original copula C is the product copula, and it is related with mixtures     

A note on generalized wald’s method

Let {v _n(θ)} be a sequence of statistics such that whenθ =θ ₀,v _n(θ ₀) N _p(0,Σ), whereΣ is of rankp andθ εR ^d. Suppose that underθ =θ ₀, {Σ _n} is a sequence of consistent estimators ofΣ. Wald (1943) shows thatv _n ^T (θ ₀)Σ _n ⁻¹ v _n(θ ₀) x ²(p). It often happens thatv _n(θ ₀) N _p(0,Σ) holds butΣ is singular. Moore (1977) states that under certain assumptionsv _n ^T (θ ₀)Σ _n ⁻ v _n(θ ₀) x ²(k), wherek = rank (Σ) andΣ _n ⁻ is a generalized inverse ofΣ _n. However, Moore’s result as stated is incorrect. It needs the additional assumption that rank (Σ _n) =k forn sufficiently large. In this article, we show that Moore’s result (as corrected) holds under somewhat different, but easier to verify, assumptions. Research partly supported by the U.S. Army Research Office through the Mathematical Sciences Institute at Cornell University.     

Comparison of test vs. control treatments using distance optimality criterion

We consider the problem of comparison of one test treatment (τ₀) with a set of v control treatments (τ₁, τ₂, …, τ_v) using distance optimality [DS-optimality] criterion introduced by Sinha (1970) in some treatment-connected design settings. It turns out that the nature of DS-optimal designs is quite similar to that for the usual A−, D− and E− optimality criteria. However, the optimality problem is quite complicated in most situations. First we deal with the CRD model and derive DS-optimal allocations for a given set of treatments. The results are almost identical to the A-optimal allocations for such problems. Then we consider a block design set-up and examine the nature of DS-optimal designs. In the process, we introduce the method of weighted coverage probability and maximize the resulting expression to obtain an optimal design. Received: December 1999