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

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.     

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

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     

Bounds for the mean residual life function of a <Emphasis Type="Italic">k</Emphasis>-out-of-<Emphasis Type="Italic">n</Emphasis> system

In the reliability studies, k-out-of-n systems play an important role. In this paper, we consider sharp bounds for the mean residual life function of a k-out-of-n system consisting of n identical components with independent lifetimes having a common distribution function F, measured in location and scale units of the residual life random variable X _t  = (X−t|X > t). We characterize the probability distributions for which the bounds are attained. We also evaluate the so obtained bounds numerically for various choices of k and n.     

Asymptotic normality of M-estimators in a semiparametric model with longitudinal data

Suppose that the longitudinal observations (Y _ij , X _ij , t _ij ) for i = 1, . . . ,n; j = 1, . . . ,m _i are modeled by the semiparamtric model where β ₀ is a k × 1 vector of unknown parameters, g(·) is an unknown estimated function and e _ij are unobserved disturbances. This article consider M-type regressions which include mean, median and quantile regressions. The M-estimator of the slope parameter β ₀ is obtained through piecewise local polynomial approximation of the nonparametric component. The local M-estimator of g(·) is also obtained by replacing β ₀ in model with its M-estimator and using local linear approximation. The asymptotic distribution of the estimator of β ₀ is derived. The asymptotic distributions of the local M-estimators of g(·) at both interior and boundary points are also established. Various applications of our main results are given. The research is supported in part by National Natural Science Foundation of China (Grant No. 10671089).     

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 new interpretation of optimality forE-optimal designs in linear regression models

The optimal design problem for the estimation of several linear combinationsc′ _l ϑ (l=1, …,m) is considered in the usual linear regression modely=f′(x)ϑ (f(x) ∈ ℝ ^k ,ϑ ∈ ℝ ^k ). An optimal design minimizes a (weighted)p-norm of the variances of the least squares estimates for the different linear combinationsc′ _l ϑ. A generalized Elfving theorem is used to derive the relation of the new optimality criterion to theE-optimal design problem. It is shown that theE-optimal design for the parameterϑ minimizes such a (weighted)p-norm whenever the vectorc=(c′ ₁, …, c′_k)′ is an inball vector of a symmetric convex and compact “Elfving set” in.     

Tests und Konfidenzintervalle für exponentielle Verteilungen und deren Anwendung auf einige diskrete Verteilungen

Zusammenfassung Es werden Verteilungen betrachtet, die (bezüglich irgendeines Ma?es) eine Dichte der GestaltC(ϑ) exp [ϑ x] besitzen. Für solche Verteilungen werden (ein- und zweiseitige) Tests und Konfidenzintervalle mit gewissen Optimalit?tseigenschaften entwickelt, und zwar fürϑ, für die Differenzϑ ₁ −ϑ ₂, sowie für einige Versionen desk-Stichproben Problems. Sodann werden einige Hilfss?tze über den bedingten Erwartungswert und die bedingte Varianz von zweiparametrigen Verteilungen abgeleitet, die bezüglich des einen Parameters reproduktiv sind und eine bezüglich des zweiten Parameters ersch?pfende und vollst?ndige Funktion besitzen. Schlie?lich werden die allgemeinen Ergebnisse auf einige diskrete Verteilungen (Binomial, Poisson, negativ Binomial, Pascal) angewendet und der Zusammenhang mit verschiedenen bekannten Tests diskutiert.

---

Summary Probability distributions are considered which (with respect to any measure) possess a density function of the typeC(ϑ) exp [ϑ x]. For distributions of this type (one and twosided) tests and confidence intervals with some optimal properties are given, namely forϑ, for the differenceϑ ₁ −ϑ ₂, and for several versions of thek-sample problem. Furthermore, some lemmas concerning the conditional expectation and the conditional variance are proved for two-parameter families of distributions which are reproductive in one parameter and possess a complete statistic, sufficient for the second parameter. Finally the general results are applied to some discrete distributions (binomial, Poisson, negative binomial, Pascal) and the relationship to several fairly known tests is discussed.

     

Optimal process parameters under LINEX loss function with general input quality characteristic

In this paper we generalize the quality and cost trade-off problem of Chang and Hung (Qual Quant 41: 291–301, 2007) under the LINEX loss function. We consider the general input characteristic given by the random variable X with moment generating function m _X (t) and output characteristic given by the deterministic transformation Y  =  g(X). The two cases we consider are when g(X) is an affine function of X and X follows (1) the gamma distribution, and (2) the double exponential distribution.     

On characterizing discrete distributions via conditional expectations of weak record values

Let (W _n ,n ≥ 0) denote the sequence of weak records from a distribution with support S = { α₀,α₁,...,α _N }. In this paper, we consider regression functions of the form ψ _n (x) = E(h(W _n ) |W _n+1 = x), where h(·) is some strictly increasing function. We show that a single function ψ _n (·) determines F uniquely up to F(α₀). Then we derive an inversion formula which enables us to obtain F from knowledge of ψ _n (·), ψ _n-1(·), h(·) and F(α₀).     

Two-sample rank tests with truncated populations

LetX ₁,…,X _m andY ₁,…,Y _n be two independent samples from continuous distributionsF andG respectively. Using a Hoeffding (1951) type theorem, we obtain the distributions of the vector S=(S ₍₁₎,…,S _(n)), whereS _(j)=# (X _i ’s≤Y _(j)) andY _(j) is thej-th order statistic ofY sample, under three truncation models: (a)G is a left truncation ofF orG is a right truncation ofF, (b)F is a right truncation ofH andG is a left truncation ofH, whereH is some continuous distribution function, (c)G is a two tail truncation ofF. Exploiting the relation between S and the vectorR of the ranks of the order statistics of theY-sample in the pooled sample, we can obtain exact distributions of many rank tests. We use these to compare powers of the Hajek test (Hajek 1967), the Sidak Vondracek test (1957) and the Mann-Whitney-Wilcoxon test. We derive some order relations between the values of the probagility-functions under each model. Hence find that the tests based onS ₍₁₎ andS _(n) are the UMP rank tests for the alternative (a). We also find LMP rank tests under the alternatives (b) and (c).     

Multivariate t and beta distributions associated with the multivariate F distribution

Relationships between F, skew t and beta distributions in the univariate case are in this paper extended in a natural way to the multivariate case. The result is two new distributions: a multivariate t/skew t distribution (on ℜ^m) and a multivariate beta distribution (on (0,1)^m). A special case of the former distribution is a new multivariate symmetric t distribution. The new distributions have a natural relationship to the standard multivariate F distribution (on (ℜ⁺)^m) and many of their properties run in parallel. We look at: joint distributions, mathematically and graphically; marginal and conditional distributions; moments; correlations; local dependence; and some limiting cases. Received: March 2001     

A note on reliability properties of k-record statistics

LetY _k,n denote the nth (upper) k-record value of an infinite sequence of independent, identically distributed random variables with common continuous distribution function F. We show that if the nth k-record valueY _k,n has an increasing failure rate (IFR), thenY _l,n (l<k) andY _k+1,n+1(n≤k+1) also have IFR distributions. On the other hand, ifY _k,n has a decreasing failure rate (DFR), thenY _l,n (1>k) has also a DFR distribution. We also present some results concerning log convexity and log concavity ofY _k,n .     

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.     

Regression analysis and dependence

In this paper we study the relationship between regression analysis and a multivariate dependency measure. If the general regression model Y=f() holds for some function f, where 1i₁< i₂<···i_m k, and X₁,...,X_k is a set of possible explanatory random variables for Y. Then there exists a dependency relation between the random variable Y and the random vector (). Using the dependency statistic defined below, we can detect such dependency even if the function f is not linear. We present several examples with real and simulated data to illustrate this assertion. We also present a way to select the appropriate subset among the random variables X₁,X₂,...,X_k, which better explain Y.     

Notes and Comments: Sup-convolutions of HARA utilities in the affine term structure

Abstract In the financial literature, the problem of maximizing the expected utility of the terminal wealth has been investigated extensively (for a survey, see, e.g., Karatzas and Shreve (1998), p. 153, and references therein) by using different approaches. In this paper, we extend the existing literature in two directions. First, we let the utility function U(.) of the financial agent (who is a price taker) be implicitly defined through I(.)=(U ^′ (.))^–1, which is assumed to be additively separable, i.e., I(.)=∑ _k=1 ^N I _k (.). Second, we solve the investment problem in the general affine term structure model proposed by Duffie and Kan (1996) in which the functions I _k (.), k=1,...,N are associated to HARA utility functions (with possibly different risk aversion parameters), and we show that the utility maximization problem leads to a Riccati ODE. Moreover, we extend to the multi-factor framework the stability result proved in Grasselli (2003), namely, the almost-sure convergence of the solution with respect to the parameters of the utility function. Mathematics Subject Classification (2000): 91B28 Journal of Economic Literature Classification: G11     

On the application of the Minkowski-Farkas theorem to sampling designs

In the paper, we consider the following problem: Let {π_k} be a sequence satisfying 0π_kΣ1 (k=1,…, N) and π=n.Tben, is there an unordered sampling design such that, for each k=1,…N, the inclusion probability of unit k is equal to π? It is shown that it can be solved by the straightforward application of the Minkowski-Farkas theorem.     

Slow convergence of the number of near-maxima for Burr XII distributions

A near-maximum is an observation which falls within a distance a of the maximum observation in an independent and identically distributed sample of size n. Subject to some conditions on the tail thickness of the population distribution, the number K _n (a) of near-maxima is known to converge in probability to one or infinity, or in distribution to a shifted geometric law. In this paper we show that for all Burr XII distributions K _n (a) converges almost surely to unity, but this convergence property may not become clear under certain cases even for very large n. We explore the reason of such slow convergence by studying a distributional continuity between Burr XII and Weibull distributions. We have also given a theoretical explanation of slow convergence of K _n (a) for the Burr XII distributions by showing that the rate of convergence in terms of P{K _n (a) > 1} tending to zero changes very little with the sample size n. Illustrations of the limiting behaviour K _n (a) for the Burr XII and the Weibull distributions are given by simulations and real data. The study also raises an important issue that although the Burr XII provides overall better fit to a given data set than the Weibull distribution, cautions should be taken for the extrapolation of the upper tail behaviour in the case of slow convergence.