A reinforcement-depletion urn model: A contiguity case

Summary An urn contains balls ofs different colors. The problem of the reinforcement of a specified color and random depletion of balls has been considered by Bernard (1977,Bull. Math. Biol.,39, 463–470) and Shenton (1981,Bull. Math. Biol.,43, 327–340), (1983,Bull. Math. Biol.,45, 1–9). Here we consider a special relation between a reinforcement and depletion, leading to a hypergeometric distribution. Research sponsored in part by the Applied Mathematical Sciences Research Program, Office of Energy Research, U.S. Department of Energy under contract DE-AC05-840R21400 with the Martin Marietta Energy Systems, Inc.     

On the estimation of entropy

Motivated by recent work of Joe (1989,Ann. Inst. Statist. Math.,41, 683–697), we introduce estimators of entropy and describe their properties. We study the effects of tail behaviour, distribution smoothness and dimensionality on convergence properties. In particular, we argue that root-n consistency of entropy estimation requires appropriate assumptions about each of these three features. Our estimators are different from Joe's, and may be computed without numerical integration, but it can be shown that the same interaction of tail behaviour, smoothness and dimensionality also determines the convergence rate of Joe's estimator. We study both histogram and kernel estimators of entropy, and in each case suggest empirical methods for choosing the smoothing parameter.     

A characterization of discrete unimodality with applications to variance upper bounds

Bertin and Theodorescu (1984,Statist. Probab. Lett.,2, 23–30) developed a characterization of discrete unimodality based on convexity properties of a discretization of distribution functions. We offer a new characterization of discrete unimodality based on convexity properties of a piecewise linear extension of distribution functions. This reliance on functional convexity, as in Khintchine's classic definition, leads to variance dilations and upper bounds on variance for a large class of discrete unimodal distributions. These bounds are compared to existing inequalities due to Muilwijk (1966,Sankhy, Ser. B,28, p. 183), Moors and Muilwijk (1971,Sankhy, Ser. B,33, 385–388), and Rayner (1975,Sankhy, Ser. B,37, 135–138), and are found to be generally tighter, thus illustrating the power of unimodality assumptions.     

Interval estimates for posterior probabilities in a multivariate normal classification model

This paper is devoted to the asymptotic distribution of estimators for the posterior probability that a p-dimensional observation vector originates from one of k normal distributions with identical covariance matrices. The estimators are based on training samples for the k distributions involved. Observation vector and prior probabilities are regarded as given constants. The validity of various estimators and approximate confidence intervals is investigated by simulation experiments.     

Distributions of orientations on Stiefel manifolds

The Riemann space whose elements are m × k (m k) matrices X, i.e., orientations, such that X′X = I_k is called the Stiefel manifold V_k,m. The matrix Langevin (or von Mises-Fisher) and matrix Bingham distributions have been suggested as distributions on V_k,m. In this paper, we present some distributional results on V_k,m. Two kinds of decomposition are given of the differential form for the invariant measure on V_k,m, and they are utilized to derive distributions on the component Stiefel manifolds and subspaces of V_k,m for the above-mentioned two distributions. The singular value decomposition of the sum of a random sample from the matrix Langevin distribution gives the maximum likelihood estimators of the population orientations and modal orientation. We derive sampling distributions of matrix statistics including these sample estimators. Furthermore, representations in terms of the Hankel transform and multi-sample distribution theory are briefly discussed.     

Estimation of the ratio of the scale parameters of two exponential distributions with unknown location parameters

We consider the estimation of the ratio of the scale parameters of two independent two-parameter exponential distributions with unknown location parameters. It is shown that the best affine equivariant estimator (BAEE) is inadmissible under any loss function from a large class of bowl-shaped loss functions. Two new classes of improved estimators are obtained. Some values of the risk functions of the BAEE and two improved estimators are evaluated for two particular loss functions. Our results are parallel to those of Zidek (1973, Ann. Statist., 1, 264–278), who derived a class of estimators that dominate the BAEE of the scale parameter of a two-parameter exponential distribution.     

Fisher information in window censored renewal process data and its applications

Suppose we have a renewal process observed over a fixed length of time starting from a random time point and only the times of renewals that occur within the observation window are recorded. Assuming a parametric model for the renewal time distribution with parameter θ, we obtain the likelihood of the observed data and describe the exact and asymptotic behavior of the Fisher information (FI) on θ contained in this window censored renewal process. We illustrate our results with exponential, gamma, and Weibull models for the renewal distribution. We use the FI matrix to determine optimal window length for designing experiments with recurring events when the total time of observation is fixed. Our results are useful in estimating the standard errors of the maximum likelihood estimators and in determining the sample size and duration of clinical trials that involve recurring events associated with diseases such as lupus.     

A Note on a Composition of Two Random Integral Mappings β and Some Examples

Abstract

The method of random integral representation, that is, the method of representing a given probability measure as the probability distribution of some random integral, was quite successful in the past few decades. In this note, we show that a composition of two random integral mappings ^β is again a random integral mapping. We illustrate our results with some examples.     

Likelihood ratio tests for comparingk populations —The two-parameter nonregular models

Summary The null and nonnull distributions of the likelihood ratio statistics for testing the homogeneity ofk given populations, each associated with a nonregular density depending on two truncation parameters, are investigated. This generalizes to the two-parameter case the work of Hogg (1956,Ann. Math. Statist.,27, 529–532), Barr (1966,J. Amer. Statist. Assoc.,61, 856–864) and Khatri and Jaiswal (1969,Aust. J. Statist.,11, 79–84; 1969, 1971,Ann. Inst. Statist. Math.,21, 127–136;23, 199–210).     

Shrinkage estimators of the location parameter for certain spherically symmetric distributions

We consider estimation of a location vector for particular subclasses of spherically symmetric distributions in the presence of a known or unknown scale parameter. Specifically, for these spherically symmetric distributions we obtain slightly more general conditions and larger classes of estimators than Brandwein and Strawderman (1991,Ann. Statist.,19, 1639–1650) under which estimators of the formX +ag(X) dominateX for quadratic loss, concave functions of quadratic loss and general quadratic loss.Research supported by NSF grant DMS-88-22622     

Asymptotic theorems for estimating the distribution function under random truncation

Representation theorem and local asymptotic minimax theorem are derived for nonparametric estimators of the distribution function on the basis of randomly truncated data. The convolution-type representation theorem asserts that the limiting process of any regular estimator of the distribution function is at least as dispersed as the limiting process of the product-limit estimator. The theorems are similar to those results for the complete data case due to Beran (1977, Ann. Statist., 5, 400–404) and for the censored data case due to Wellner (1982, Ann. Statist., 10, 595–602). Both likelihood and functional approaches are considered and the proofs rely on the method of Begun et al. (1983, Ann. Statist., 11, 432–452) with slight modifications.Division of Biostatistics, School of Public Health, Columbia Univ.     

Conditional limit results for type I polar distributions

Let (S ₁,S ₂) = (R cos(Θ), R sin(Θ)) be a bivariate random vector with associated random radius R which has distribution function F being further independent of the random angle Θ. In this paper we investigate the asymptotic behaviour of the conditional survivor probability when u approaches the upper endpoint of F. On the density function of Θ we impose a certain local asymptotic behaviour at 0, whereas for F we require that it belongs to the Gumbel max-domain of attraction. The main result of this contribution is an asymptotic expansion of , which is then utilised to construct two estimators for the conditional distribution function . Furthermore, we allow Θ to depend on u.      

Empirical Bayes estimation in a multiple linear regression model

Summary Estimation of the vector β of the regression coefficients in a multiple linear regressionY=Xβ+ε is considered when β has a completely unknown and unspecified distribution and the error-vector ε has a multivariate standard normal distribution. The optimal estimator for β, which minimizes the overall mean squared error, cannot be constructed for use in practice. UsingX, Y and the information contained in the observation-vectors obtained fromn independent past experiences of the problem, (empirical Bayes) estimators for β are exhibited. These estimators are compared with the optimal estimator and are shown to be asymptotically optimal. Estimators asymptotically optimal with rates nearO(n ⁻¹) are constructed. Supported in part by a Natural Sciences and Engineering Research Council of Canada grant.     

Estimation in Linear Models with Random Effects and Errors-in-Variables

The independent variables of linear mixed models are subject to measurement errors in practice. In this paper, we present a unified method for the estimation in linear mixed models with errors-in-variables, based upon the corrected score function of Nakamura (1990, Biometrika, 77, 127–137). Asymptotic normality properties of the estimators are obtained. The estimators are shown to be consistent and convergent at the order of n ^–1/2. The performance of the proposed method is studied via simulation and the analysis of a data set on hedonic housing prices.     

The expected values of invariant polynomials with matrix argument of elliptical distributions

ThisresearchissupportedbytheChineseAcademyofSciences.1.IntroductionTherearemailypublicatiollsonmatrixvariatedistributions(ormatrixdistributionforsimplicity),inparticular,abollttheirexpectedvaluesofzonalpolynomialsoftheirquadraticforms(of.[12],[151,[171and…     

Kernel Density and Hazard Rate Estimation for Censored Data under α-Mixing Condition

We derive rates of uniform strong convergence for kernel density estimators and hazard rate estimators in the presence of right censoring. It is assumed that the failure times (survival times) form a stationary -mixing sequence. Moreover, we show that, by an appropriate choice of the bandwidth, both estimators attain the optimal strong convergence rate known from independent complete samples. The results represent an improvement over that of Cai's paper (cf. Cai (1998b, J. Multivariate Anal., 67, 23–34)).     

Estimation of two normal means which may be common

Consider the problem of estimating the mean of a normal population when independent samples from this as well as a second normal population are available. Pre-test estimators which combine the two sample means if a test of the hypothesis of equal population means accepts but otherwise use only the first sample mean, are compared to limited translation estimators which are derived in the spirit of Bickel (1984, Ann. Statist., 12, 864–879) (we also cover the cases of unknown variances). Our conclusion is that if the accuracy with which the second population mean can be estimated is of the same or better order of magnitude as teh accuracy with which the first can be estimated, then the limited translation estimators largely dominate the pre-test estimators in terms of mean square error loss.This research was supported by grants from the FRD of the CSIR of South Africa.     

Estimation of parameters in a two-parameter exponential distribution using ranked set sample

In situations where the experimental or sampling units in a study can be easily ranked than quantified, McIntyre (1952,Aust. J. Agric. Res.,3, 385–390) proposed that the mean ofn units based on aranked set sample (RSS) be used to estimate the population mean, and observed that it provides an unbiased estimator with a smaller variance compared to a simple random sample (SRS) of the same sizen. McIntyre's concept ofRSS is essentially nonparametric in nature in that the underlying population distribution is assumed to be completely unknown. In this paper we further explore the concept ofRSS when the population is partially known and the parameter of interest is not necessarily the mean. To be specific, we address the problem of estimation of the parameters of a two-parameter exponential distribution. It turns out that the use ofRSS and its suitable modifications results in much improved estimators compared to the use of aSRS.     

Identifiability of mixtures of power-series distributions and related characterizations

The concept of the identifiability of mixtures of distributions is discussed and a sufficient condition for the identifiability of the mixture of a large class of discrete distributions, namely that of the power-series distributions, is given. Specifically, by using probabilistic arguments, an elementary and shorter proof of the Lüxmann-Ellinghaus's (1987,Statist. Probab. Lett.,5, 375–378) result is obtained. Moreover, it is shown that this result is a special case of a stronger result connected with the Stieltjes moment problem. Some recent observations due to Singh and Vasudeva (1984,J. Indian Statist. Assoc.,22, 93–96) and Johnson and Kotz (1989,Ann. Inst. Statist. Math.,41, 13–17) concerning characterizations based on conditional distributions are also revealed as special cases of this latter result. Exploiting the notion of the identifiability of power-series mixtures, characterizations based on regression functions (posterior expectations) are obtained. Finally, multivariate generalizations of the preceding results have also been addressed.     

On Adaptive Combination of Regression Estimators

Consider the problem of choosing between two estimators of the regression function, where one estimator is based on stronger assumptions than the other and thus the rates of convergence are different. We propose a linear combination of the estimators where the weights are estimated by Mallows' C _L . The adaptive estimator retains the optimal rates of convergence and is an extension of Stein-type estimators considered by Li and Hwang (1984, Ann. Statist., 12, 887-897) and related to an estimator in Burman and Chaudhuri (1999, Ann. Inst. Statist. Math. (to appear)).