Estimation of the precision matrix of multivariate Pearson type II model

In this paper, the problem of estimating the precision matrix of a multivariate Pearson type II-model is considered. A new class of estimators is proposed. Moreover, the risk functions of the usual and the proposed estimators are explicitly derived. It is shown that the proposed estimator dominates the MLE and the unbiased estimator, under the quadratic loss function. A simulation study is carried out and confirms these results. Improved estimator of tr (Σ ⁻¹) is also obtained.     

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.     

Gamma-admissibility of estimators in the one-parameter exponential family

Summary Admissibility of estimators under vague prior information on the distribution of the unknown parameter is studied which leads to the notion of gamma-admissibility. A sufficient condition for an estimator of the formδ(x)=(ax+b)/(cx+d) to be gamma-admissible in the one-parameter exponential family under squared error loss is established. As an application of this result two equalizer rules are shown to be unique gamma-minimax estimators by proving their gamma-admissibility.     

Gamma-minimax estimation of a multivariate normal mean

Summary The mean vector of a multivariate normal distribution is to be estimated. A class Γ of priors is considered which consists of all priors whose vector of first moments and matrix of second moments satisfy some given restrictions. The Γ-minimax estimator under arbitrary squared error loss is characterized. The characterization follows from an application of a result of Browder and Karamardian published in Ichiishi (1983) which is a special version of a minimax inequality due to Ky Fan (1972). In particular, it is shown that within the set of all estimators a linear estimator is Γ-minimax. The authors would like to thank the Deutsche Forschungsgemeinschaft for financial support.     

A general class of estimators for estimating population mean using auxiliary information

Summary A general class of estimators for estimating the population mean of the character under study which make use of auxiliary information is proposed. Under simple random sampling without replacement (SRSWOR), the expressions of Bias and Mean Square Error (MSE), up to the first and the second degrees of approximation are derived. General conditions, up to the first order approximation, are also obtained under which any member of this class performs more efficiently than the mean per unit estimator, the ratio estimator and the product estimator. The class of estimators in its optimum case, under the first degree approximation, is discussed. It is shown that it is not possible to obtain optimum values of parameters “a”, “b” and “p”, that are independent of each other. However, the optimum relation among them is given by (b−a)p=ρ C _y/C _x. Under this condition, the expression of MSE of the class is that of the linear regression estimator.     

Minimax estimation under convex loss when the parameter interval is bounded

In the present paper families of truncated distributions with a Lebesgue density forx=(x ₁,...,x _n ) ε ℝ ⁿ are considered, wheref ₀:ℝ → (0, ∞) is a known continuous function andC _n (ϑ) denotes a normalization constant. The unknown truncation parameterϑ which is assumed to belong to a bounded parameter intervalΘ=[0,d] is to be estimated under a convex loss function. It is studied whether a two point prior and a corresponding Bayes estimator form a saddle point when the parameter interval is sufficiently small.     

Minimax invariant estimation of a continuous distribution function under entropy loss

For the invariant decision problem of estimating a continuous distribution function F with two entropy loss functions, it is proved that the best invariant estimators d ₀ exist and are the same as the best invariant estimator of a continuous distribution function under the squared error loss function L (F, d)=∫|F (t) −d (t) |² dF (t). They are minimax for any sample size n≥1.     

On estimating the mean of the selected normal population under the LINEX loss function

Following Parsian and Farsipour (1999), we consider the problem of estimating the mean of the selected normal population, from two normal populations with unknown means and common known variance, under the LINEX loss function. Some admissibility results for a subclass of equivariant estimators are derived and a sufficient condition for the inadmissibility of an arbitrary equivariant estimator is provided. As a consequence, several of the estimators proposed by Parsian and Farsipour (1999) are shown to be inadmissible and better estimators are obtained. Received January 2001/Revised May 2002     

Minimax estimation of a cumulative distribution function by converting to a parametric problem

Let X = (X ₁,...,X _n ) be a sample from an unknown cumulative distribution function F defined on the real line . The problem of estimating the cumulative distribution function F is considered using a decision theoretic approach. No assumptions are imposed on the unknown function F. A general method of finding a minimax estimator d(t;X) of F under the loss function of a general form is presented. The method of solution is based on converting the nonparametric problem of searching for minimax estimators of a distribution function to the parametric problem of searching for minimax estimators of the probability of success for a binomial distribution. The solution uses also the completeness property of the class of monotone decision procedures in a monotone decision problem. Some special cases of the underlying problem are considered in the situation when the loss function in the nonparametric problem is defined by a weighted squared, LINEX or a weighted absolute error.     

Bayesian estimation and prediction based on Rayleigh sample quantiles

Ordered data arise naturally in many fields of statistical practice. Often some sample values are unknown or disregarded due to various reasons. On the basis of some sample quantiles from the Rayleigh distribution, the problems of estimating the Rayleigh parameter, hazard rate and reliability function, and predicting future observations are addressed using a Bayesian perspective. The construction of β-content and β-expectation Bayes tolerance limits is also tackled. Under squared-error loss, Bayes estimators and predictors are deduced analytically. Exact tolerance limits are derived by solving simple nonlinear equations. Highest posterior density estimators and credibility intervals, as well as Bayes estimators and predictors under linear loss, can easily be computed iteratively.     

Stochastic monotonicity of the MLEs of parameters in exponential simple step-stress models under Type-I and Type-II censoring

In two recent papers by Balakrishnan et al. (J Qual Technol 39:35–47, 2007; Ann Inst Stat Math 61:251–274, 2009), the maximum likelihood estimators [^(q)]₁{\hat{\theta}_{1}} and [^(q)]₂{\hat{\theta}_{2}} of the parameters θ ₁ and θ ₂ have been derived in the framework of exponential simple step-stress models under Type-II and Type-I censoring, respectively. Here, we prove that these estimators are stochastically monotone with respect to θ ₁ and θ ₂, respectively, which has been conjectured in these papers and then utilized to develop exact conditional inference for the parameters θ ₁ and θ ₂. For proving these results, we have established a multivariate stochastic ordering of a particular family of trinomial distributions under truncation, which is also of independent interest.     

Testing forecast accuracy of expectiles and quantiles with the extremal consistent loss functions

Forecast evaluations aim to choose an accurate forecast for making decisions by using loss functions. However, different loss functions often generate different ranking results for forecasts, which complicates the task of comparisons. In this paper, we develop statistical tests for comparing performances of forecasting expectiles and quantiles of a random variable under consistent loss functions. The test statistics are constructed with the extremal consistent loss functions of Ehm et al. (2016). The null hypothesis of the tests is that a benchmark forecast at least performs equally well as a competing one under all extremal consistent loss functions. It can be shown that if such a null holds, the benchmark will also perform at least equally well as the competitor under all consistent loss functions. Thus under the null, when different consistent loss functions are used, the result that the competitor does not outperform the benchmark will not be altered. We establish asymptotic properties of the proposed test statistics and propose to use the re-centered bootstrap to construct their empirical distributions. Through simulations, we show that the proposed test statistics perform reasonably well. We then apply the proposed method to evaluations of several different forecast methods.     

Smooth estimators for estimating order restricted scale parameters of two gamma distributions

We consider the problem of component-wise estimation of ordered scale parameters of two gamma populations, when it is known apriori which population corresponds to each ordered parameter. Under the scale equivariant squared error loss function, smooth estimators that improve upon the best scale equivariant estimators are derived. These smooth estimators are shown to be generalized Bayes with respect to a non-informative prior. Finally, using Monte Carlo simulations, these improved smooth estimators are compared with the best scale equivariant estimators, their non-smooth improvements obtained in Vijayasree, Misra & Singh (1995), and the restricted maximum likelihood estimators. Acknowledgments. Authors are thankful to a referee for suggestions leading to improved presentation.     

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.     

Decision theoretic estimation using record statistics

We consider the problem of estimating the scale parameter θ of the shifted exponential distribution with unknown shift based on a set of observed records drawn from a sequential sample of independent and identically distributed random variables. Under a large class of bowl-shaped loss functions, the best affine equivariant estimator (BAEE) of θ is shown to be inadmissible. Two dominating procedures are proposed. A numerical study is performed to show the extent of risk reduction that the improved estimators provide over the BAEE.     

Sequential estimation of a linear function of normal means under asymmetric loss function

The problem of estimating a linear function of k normal means with unknown variances is considered under an asymmetric loss function such that the associated risk is bounded from above by a known quantity. In the absence of a fixed sample size rule, sequential stopping rules satisfying a general set of assumptions are considered. Two estimators are proposed and second-order asymptotic expansions of their risk functions are derived. It is shown that the usual estimator, namely the linear function of the sample means, is asymptotically inadmissible, being dominated by a shrinkage-type estimator. An example illustrates the use of different multistage sampling schemes and provides asymptotic expansions of the risk functions. Received: August 1999     

Estimation of the scale matrix of a multivariate t-model under entropy loss

This paper deals with the estimation of the scale matrix of a multivariatet-model with unknown location vector and scale matrix to improve upon the usual estimators based on the sample sum of product matrix. The well-known results of the estimation of the scale matrix of the multivariate normal model under the assumption of entropy loss function have been generalized to that of a multivariatet-model. The paper is based on the first author’s unpublished Ph.D. dissertation ‘Estimation of the Scale Matrix of a Multivariate T-model’, University of Western Ontario, Canada. Present address: School of Mathematics and Statistics, The University of Sydney, NSW 2006, Australia.     

Sequential estimation of normal mean under asymmetric loss function with a shrinkage stopping rule

The problem of estimating a normal mean with unknown variance is considered under an asymmetric loss function such that the associated risk is bounded from above by a known quantity. In the absence of a fixed sample size rule, a sequential stopping rule and two sequential estimators of the mean are proposed and second-order asymptotic expansions of their risk functions are derived. It is demonstrated that the sample mean becomes asymptotically inadmissible, being dominated by a shrinkage-type estimator. Also a shrinkage factor is incorporated in the stopping rule and similar inadmissibility results are established. Received September 1997     

Credibility estimators with dependence structure over risks and time under balanced loss function

In this paper, we study the Bühlmann credibility model with constant interest rate and equal dependence structure over risks and time under balanced loss function. By means of orthogonal projection, the inhomogeneous and homogeneous credibility premium estimators are derived, which extend those for the existing models to slightly more general versions. Finally, we investigate the estimation of the structure parameters and present a numerical example to show the effectiveness of the inhomogeneous estimator.     

Estimators of Regression Parameters for Truncated and Censored Data

Estimators of parameters in semi-parametric left truncated and right censored regression models are proposed. In contrast to the majority of existing estimators, the proposed estimators do not require the error term of the regression model to have a symmetric distribution. In addition the estimators use asymmetric “trimming” of observations. Consistency and asymptotic normality of the estimators are shown. Finite sample properties are considered in a small simulation study. For the left truncated case, an empirical application illustrates the usefulness of the estimator.