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.     

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

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.     

Characterization of some types of completeness resp. Total completeness and their conservation under direct products

Some notions ofL _p (μ)-completeness resp. totally L _p (μ)-completeness (1≦p≦∞) are characterized for families of probability distributions dominated by aσ-finite measureμ and their conservation with respect to direct products is proved. Furthermore, it is shown that totallyL _∞(μ)-completeness does not implyL ₁(μ)-completeness and that there are families of probability distributions in the i.i.d. case induced by the order statistic, which are L₁(μ)-complete but not totallyL _∞(μ)-complete.     

Sequential point estimation of the powers of a normal scale parameter

We consider the sequential point estimation problem of the powers of a normal scale parameter σ^r with r≠ 0 when the loss function is squared error plus linear cost. It is shown that the regret due to using our fully sequential procedure in ignorance of σ is asymptotically minimized for estimating σ⁻². We also propose a bias-corrected procedure to reduce the risk and show that the larger the distance between r and −2 is, the more effective our bias-corrected procedure is. Received August 2000     

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.     

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.     

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.     

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

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

Estimation of the parameters of log-gamma distribution using order statistics

In this work we propose a technique of estimating the location parameter μ and scale parameter σ of log-gamma distribution by U-statistics constructed by taking best linear functions of order statistics as kernels. The efficiency comparison of the proposed estimators with respect to maximum likelihood estimators is also made.     

Nonparametric regression function estimation using interaction least squares splines and comlexity regularization

Let (X, Y) be a pair of random variables withsupp(X)⊆[0,1] ^l andEY ²<∞. Letm ^* be the best approximation of the regression function of (X, Y) by sums of functions of at mostd variables (1≤d≤l). Estimation ofm ^* from i.i.d. data is considered. For the estimation interaction least squares splines, which are defined as sums of polynomial tensor product splines of at mostd variables, are used. The knot sequences of the tensor product splines are chosen equidistant. Complexity regularization is used to choose the number of the knots and the degree of the splines automatically using only the given data. Without any additional condition on the distribution of (X, Y) the weak and strongL ₂-consistency of the estimate is shown. Furthermore, for everyp≥1 and every distribution of (X, Y) withsupp(X)⊆[0,1] ^l ,y bounded andm ^* p-smooth, the integrated squared error of the estimate achieves up to a logarithmic factor the (optimal) rate      

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.     

Bounds on expectations of small order statistics from decreasing density populations

For i > (n + 1)/2, Danielak (Statistics 37:305–324, 2003) established an optimal positive upper mean-variance bound on the expectation of ith order statistic based on the i.i.d. sample of size n from the decreasing density population. We show that the best bounds on the expected deviation of the ith order statistics from the population mean, i ≤ (n + 1)/2, expressed in more general scale units generated by pth absolute central moments with p > 1 amount to zero. We also determine the respective strictly negative bounds in the mean absolute deviation units.     

On maximum likelihood prediction based on Type II doubly censored exponential data

In this paper, the maximum likelihood predictor (MLP) of the kth ordered observation, t _k, in a sample of size n from a two-parameter exponential distribution as well as the predictive maximum likelihood estimators (PMLE's) of the location and scale parameters, θ and β, based on the observed values t _r, …, t _s (1≤r≤s<k≤n), are obtained in closed forms, contrary to the belief they cannot be so expressed. When θ is known, however, the PMLE of β and MLP of t _k do not admit explicit expressions. It is shown here that they exist and are unique; sharp lower and upper bounds are also provided. The derived predictors and estimators are reasonable and also have good asymptotic properties. As applications, the total duration time in a life test and the failure time of a k-out-of-n system may be predicted. Finally, an illustrative example is included. Received: August 1999     

Multivariate gamma distributions with one-factorial accompanying correlation matrices and applications to the distribution of the multivariate range

Summary A new representation for the characteristic function of the joint distribution of the Mahalanobis distances betweenk independentN(μ, Σ)-distributed points is given. Especially fork=3 the corresponding distribution function is obtained as a special case of multivariate gamma distributions whose accompanying normal distribution has a positive semidefinite correlation matrix with correlationsϱ _ij=−a _i a _j. These gamma distribution functions are given here by one-dimensional parameter integrals. With some further trivariate gamma distributions third order Bonferroni inequalities are derived for the upper tails of the distribution function of the multivariate range ofk independentN(μ, I)-distributed points. From these inequalities very accurate (conservative) approximations to upperα-level bounds can also be computed for studentized multivariate ranges.     

On the upper and lower semicontinuity of the Aumann integral

Let (T,τ,μ) be a finite measure space, X be a Banach space, P be a metric space and let L₁(μ,X) denote the space of equivalence classes of X-valued Bochner integrable functions on (T,τ,μ). We show that if φ:T×P→2^X is a set-valued function such that for each fixed pεP, φ(·,p) has a measurable graph and for each fixed tεT, φ(t,·) is either upper or lower semicontinuous then the Aumann integral of φ, i.e.,∫_Tφ(t,p)dμ(t)= {∫_Tx(t)dμ(t):xεS_φ(p)}, where S_φ(p)= {yεL₁(μ,X):y(t)εφ(t,p)μ−a.e.}, is either upper or lower semicontinuous in the variable p as well. Our results generalize those of Aumann (1965, 1976) who has considered the above problem for X=Rⁿ, and they have useful applications in general equilibrium and game theory.     

Die messende Prüfung bei Abweichung von der Normalverteilungsannahme

Summary For sampling inspection by variables in the one-sided case (item bad if variablex>a) under the usual assumption of normality with known variance ₂ the operating characteristic is given by , wherep denotes the fraction defective. If instead of a normal distribution ((·–a–)/) there is a distributionF((·–a–)/) whereF is sufficiently regular and normed like , one has the approximative operating characteristic . It is shown that for arbitrarily fixed parametersn andc the function takes the valueL _n,c ⁽⁾ (p) at the pointp _F (p)=1–F(–^–1(p)). Sufficient conditions for a simple behavior of the differencep _F (p)–p are given. In the cases of rectangular and symmetrically truncated normal distribution these conditions are shown to be fulfilled.     

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.     

Vergleich des verallgemeinerten nichtzentralent-Tests mit dem Binomialtest

Zusammenfassung Es sei {F _,(x); –<<, >0} mitF _,(x)=F((x–)/)–F(x) eine standardisierte Verteilungsfunktion — die Familie der zulässigen Verteilungsfunktionen. Der (früher eingeführte) verallgemeinerte nichtzentralet-Test für die Hypothese {PP ₀} mitP:=F _,(x ₀) gegen die Alternative {P>P ₀} zum Niveau wird mit dem entsprechenden nichtparametrischen Test (Test für die Hypothese {pP ₀} über den Parameterp einer Binomialverteilung gegen die Alternative {p>P ₀}) verglichen. Für dent-Test wird die relative asymptotische Effizienz bestimmt.Beide Tests lassen sich als Tests für das zur WahrscheinlichkeitP ₀ gehörende Quantil einer Verteilungsfunktion interpretieren. Der klassische zentrale Student-Test ergibt sich als Spezialfall (F(x)=(x),P ₀=0,5).

---

Summary Let {F _,(x);–<<, >0} withF _,(x ₀):=F((x–)/–F(x) a standarized distribution function — the family of admissible distribution functions. The (earlier introduced) generalized noncentralt-test for the hypothesis {PP ₀} withP:=F _,(x ₀) against the alternative {P>P ₀} at level of significance is compared with the corresponding nonparametric test (Binomial test). The relative asymptotic efficiency of thet-test is determined. Both kinds of tests can be interpreted as quantiltests. In caseF(x)=(x),P ₀=0,5 one gets the classical central Student-test.