M-tests in multivariate models

Summary For testing homogeneity in multivariatek sample model, robust tests based onM-estimators are proposed and their asymptoticx ²-distributions are investigated. FurthermoreM-tests in multivariate regression models are discussed.     

Asymptotic properties of minimum contrast estimators for parameters of boolean models

Summary This paper presents a method for the estimation of parameters of random closed sets (racs’s) in ℝ ^d based on a single realization within a (large) convex sampling window. The essential idea first applied by Diggle (1981) in a special case consists in defining the estimation by minimizing a suitably defined distance (called contrast function) between the true and the empirical contact distribution function of the racs under consideration, where the most relevant case of Boolean models is discussed in details. The resulting estimates are shown to be strongly consistent (if the racs is ergodic) and asymptotically normal (if the racs is Boolean) when the sampling window expands unboundedly.     

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.     

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     

Prediction of <Emphasis Type="Italic">k</Emphasis>-records from a general class of distributions under balanced type loss functions

We study the problem of predicting future k-records based on k-record data for a large class of distributions, which includes several well-known distributions such as: Exponential, Weibull (one parameter), Pareto, Burr type XII, among others. With both Bayesian and non-Bayesian approaches being investigated here, we pay more attention to Bayesian predictors under balanced type loss functions as introduced by Jafari Jozani et al. (Stat Probab Lett 76:773–780, 2006a). The results are presented under the balanced versions of some well-known loss functions, namely squared error loss, Varian’s linear-exponential loss and absolute error loss or L ₁ loss functions. Some of the previous results in the literatures such as Ahmadi et al. (Commun Stat Theory Methods 34:795–805, 2005), and Raqab et al. (Statistics 41:105–108, 2007) can be achieved as special cases of our results. Partial support from Ordered and Spatial Data Center of Excellence of Ferdowsi University of Mashhad is acknowledged by J. Ahmadi. M. J. Jozani’s research supported partially by a grant of Statistical Research and Training Center. é. Marchand’s research supported by NSERC of Canada. A. Parsian’s research supported by a grant of the Research Council of the University of Tehran.     

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     

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     

An improved combinatorial approach for pricing Parisian options

Parisian options are path-dependent options whose payoff depends on whether the underlying asset’s price remains continuously at or above a given barrier over a given time interval. Costabile’s (Decis Econ Finance 25(2):111–125, 2002b) algorithm for pricing Parisian options based on a combinatorial approach in binomial tree has a time complexity of O( n³){O\left( {n^{3}}\right)}. We improve that algorithm to yield one with a time complexity of only O(n²){O\left({n^{2}}\right)}.     

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.     

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.     

Simultaneous estimation after selection and ranking and other procedures: The negative exponential case

We considerk (≥2) independent negative exponential populations with unknown location parameters and unknown but equal scale parameter. We incorporate the existing purely sequential and three-stage sampling procedures for selecting the “best” population and study the asymptotic second-order characteristics of the proposed fixed-size simultaneous confidence regions for the location parameters constructed after selection and ranking. Some direct estimation procedures have also been discussed.     

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.     

Generalized gauss-chebyshev inequalities for unimodal distributions

Letg be an even function on ℝ which is nondecreasing in |x|. Letk be a positive constant. Sharp inequalities relatingP(|X|≥k) toEg(X) are obtained for random variablesX which are unimodal with mode 0, and for random variablesX which are unimodal with unspecified mode. The bounds in the mode 0 case generalize an inequality due to Gauss (1823), whereg(x)=x ². The bounds in the second case generalize inequalities of Vysochanskiĭ and Petunin (1980, 1983) and Dharmadhikari and Joag-dev (1985).     

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.     

A fixed sample size selection procedure for negative binomial populations

Summary A fixed sample size procedure for selecting the ‘best’ ofk negative binomial populations is developed. Selection is made in such a way that the probability of correct selection is at leastP* whenever the distance between the probabilities of success is at leastδ*. The exponentr is assumed to be known and the same for all populations. Extensive computer calculations* were employed to obtain the exact least favorable configuration. The smallest sample sizes needed to meet specifications (P*,δ*) are tabulated forr=1 (1)5;δ*=0.05 (0.05) 0.55 andP*=0.75, 0.80, 0.90, 0.95, 0.98, 0.99 involvingk=3 (1) 6, 8, 10 populations. All the computations were carried out on the Alabama Supercomputer. Part of this work was completed when the authors were at the Department of Statistics, Oklahoma State University, Stillwater, OK 74078.     

Generalized entropies in coding theory

Equivalence of the generalized entropyH ^β (P, Φ _t ) defined in this paper andKapur’s entropy of orderα and typeβ, ie.H _α ^β (P), is established. The results given recently byCampbell follow as special cases. International Conference on System Sciences, Honolulu, January 1968.     

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     

Combining information for prediction in linear regression

This article deals with the prediction problem in linear regression where the measurements are obtained using k different devices or collected from k different independent sources. For the case of k=2, a Graybill-Deal type combined estimtor for the regression parameters is shown to dominate the individual least squares estimators under the covariance criterion. Two predictors ŷ _c and ŷ _p are proposed. ŷ _c is based on a combined estimator of the regression coefficient vector, and ŷ _p is obtained by combining the individual predictors from different models. Prediction mean square errors of both predictors are derived. It is shown that the predictor ŷ _p is better than the individual predictors for k≥2 and the predictor ŷ _c is better than the individual predictors for k=2. Numerical comparison between ŷ _c and ŷ _p shows that the former is superior to the latter for the case k=2.     

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.     

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.