On risk comparisons of some estimators in a linear regression model under inagaki’S loss function and nonnormal error terms

In this paper we consider the risk performances of some estimators for both location and scale parameters in a linear regression model under Inagaki’s loss function We prove that the pre-test estimator for location parameter is dominated by the Stein-rule estimator under Inagaki’s loss function when the distribution of error terms is expressed by the scale mixture of normal distribution and the variance of error terms is unknown.. It is an extension of the results in Nagata (1983) to our situation Also we perform numerical calculations to draw the shapes of the risks.     

Risk of a homoscedasticity pre-test estimator of the regression scale under LINEX loss

In this paper we consider the risk of an estimator of the error variance after a pre-test for homoscedasticity of the variances in the two-sample heteroscedastic linear regression model. This particular pre-test problem has been well investigated but always under the restrictive assumption of a squared error loss function. We consider an asymmetric loss function — the LINEX loss function — and derive the exact risks of various estimators of the error variance.     

The exact distribution and density functions of a pre-test estimator of the error variance in a linear regression model with proxy variables

We consider a linear regression model when some independent variables are unobservable, but proxy variables are available instead of them. We derive the distribution and density functions of a pre-test estimator of the error variance after a pre-test for the null hypothesis that the coefficients for the unobservable variables are zeros. Based on the density function, we show that when the critical value of the pre-test is unity, the coverage probability in the interval estimation of the error variance is maximum.     

Estimating the error variance after a pre-test for an inequality restriction on the coefficients

Estimation of the regression error variance after a preliminary test of an inequality constraint on the coefficient vector is considered. We derive the exact finite sample risk of several inequality restricted and pre-test estimators of σ². These estimators are associated with the maximum likelihood, least squares and minimum mean squared error component estimators. Optimal critical values for the pre-test according to a mini-max regret criterion are numerically calculated. Furthermore, we examine the robustness of the optimal choice of critical values and the risk properties of the estimators of σ² to model mis-specification through the exclusion of relevant regressors.     

On estimating the common mean in two normal distributions after a preliminary test for equality of variances

Given two random samples of equal size from two normal distributions with common mean but possibly different variances, we examine the sampling performance of the pre-test estimator for the common mean after a preliminary test for equality of variances. It is shown that when the alternative in the pretest is one-sided, the Graybill-Deal estimator is dominated by the pre-test estimator if the critical value is chosen appropriately. It is also shown that all estimators, the grand mean, the Graybill-Deal estimator and the pre-test estimator, are admissible when the alternative in the pre-test is two-sided. The optimal critical values in the two-sided pre-test are sought based on the minimax regret and the minimum average risk criteria, and it is shown that the Graybill-Deal estimator is most preferable under the minimum average risk criterion when the alternative in the pre-test is two-sided.     

On the difference in inference and prediction between the joint and independent f-error models for seemingly unrelated regressions

We consider likelihood and Bayesian inferences for seemingly unrelated (linear) regressions for the joint niultivariate terror (e.g. Zellner, 1976) and the independent t-error (e.g. Maronna, 1976) models. For likelihood inference, the scale matrix and the shape parameter for the joint terror model cannot be consistently estimated because of the lack of adequate information to identify the latter. The joint terror model also yields the same MLEs for the regression coefficients and the scale matrix as for the independent normal error model. which are not robust against outliers. Further, linear hypotheses with respect

to the regression coefficients also give rise to the same mill distributions AS for the independent normal error model, though the MLE has a non-normal limiting distribution. In contrast to the striking similarities between the joint t-error and the independent normal error models, the independent f-error model yields AiLEs that are lubust against uuthers. Since the MLE of the shape parameter reflects the tails of the data distributions, this model extends the independent normal error model for modeling data distributions with relatively t hicker tails. These differences are also discussed with respect to the posterior and predictive distributions for Bayesian inference.     

EXACT FINITE-SAMPLE PROPERTIES OF A PRE-TEST ESTIMATOR IN RIDGE REGRESSION

This paper discusses a pre-test regression estimator which uses the least squares estimate when it is “large” and a ridge regression estimate for “small” regression coefficients, where the preliminary test is applied separately to each regression coefficient in turn to determine whether it is “large” or “small.” For orthogonal regressors, the exact finite-sample bias and mean squared error of the pre-test estimator are derived. The latter is less biased than a ridge estimator, and over much of the parameter space the pre-test estimator has smaller mean squared error than least squares. A ridge estimator is found to be inferior to the pre-test estimator in terms of mean squared error in many situations, and at worst the latter estimator is only slightly less efficient than the former at commonly used significance levels.     

Risk behavior of a pre-test estimator for normal variance with the Stein-type estimator

In this paper, we examine the risk behavior of a pre-test estimator for normal variance with the Stein-type estimator. The one-sided pre-test is conducted for the null hypothesis that the population variance is equal to a specific value, and the Stein-type estimator is used if the null hypothesis is rejected. A sufficient condition for the pre-test estimator to dominate the Stein-type estimator is shown.     

The non-optimality of interval restricted and pre-test estimators under squared error loss

We consider the linear regression model with an interval restriction imposed on the coefficients, and examine the sampling performance of a family of Stein interval restricted and pre-test estimators Tor the coefficient vector. The risk, under squared error loss, of these Stein-like estimators are derived, and the inadmissibility of the maximum likelihood interval restricted and pre-test estimators is demonstrated.     

Strawderman-type estimators for a scale parameter with application to the exponential distribution

It is shown that Strawderman's [1974. Minimax estimation of powers of the variance of a normal population under squared error loss. Ann. Statist. 2, 190–198] technique for estimating the variance of a normal distribution can be extended to estimating a general scale parameter in the presence of a nuisance parameter. Employing standard monotone likelihood ratio-type conditions, a new class of improved estimators for this scale parameter is derived under quadratic loss. By imposing an additional condition, a broader class of improved estimators is obtained. The dominating procedures are in form analogous to those in Strawderman [1974. Minimax estimation of powers of the variance of a normal population under squared error loss. Ann. Statist. 2, 190–198]. Application of the general results to the exponential distribution yields new sufficient conditions, other than those of Brewster and Zidek [1974. Improving on equivariant estimators. Ann. Statist. 2, 21–38] and Kubokawa [1994. A unified approach to improving equivariant estimators. Ann. Statist. 22, 290–299], for improving the best affine equivariant estimator of the scale parameter. A class of estimators satisfying the new conditions is constructed. The results shed new light on Strawderman's [1974. Minimax estimation of powers of the variance of a normal population under squared error loss. Ann. Statist. 2, 190–198] technique.     

Bayesian estimation of the linear regression model with an uncertain interval constraint on coefficients

This article considers Bayesian inference in the interval constrained normal linear regression model. Whereas much of the previous literature has concentrated on the case where the prior constraint is correctly specified, our framework explicitly allows for the possibility of an invalid constraint. We adopt a non-informative prior and uncertainty concerning the interval restriction is represented by two prior odds ratios. The sampling theoretic risk of the resulting Bayesian interval pre-test estimator is derived, illustrated and explored. The authors are grateful to the editor and the referee for their helpful comments.     

The density function and the MSE dominance of the pre-test estimator in a heteroscedastic linear regression model with omitted variables

In this paper, we consider a heteroscedastic linear regression model with omitted variables. We derive the density function of the pre-test estimator consisting of the two-stage Aitken estimator (2SAE) and the ordinary least squares estimator (OLSE) after the pre-test for homoscedasticity. We also derive the first two moments based on the density function and show the sufficient condition for the pre-test estimator to dominate the 2SAE in terms of the MSE. Our numerical evaluations show that when this sufficient condition does not hold and when the magnitude of the specification error is large, the pre-test estimator can be dominated by the 2SAE, and further, the 2SAE can be dominated by the OLSE.     

Bayesian inference in a heteroscedastic replicated measurement error model using heavy-tailed distributions

We introduce a multivariate heteroscedastic measurement error model for replications under scale mixtures of normal distribution. The model can provide a robust analysis and can be viewed as a generalization of multiple linear regression from both model structure and distribution assumption. An efficient method based on Markov Chain Monte Carlo is developed for parameter estimation. The deviance information criterion and the conditional predictive ordinates are used as model selection criteria. Simulation studies show robust inference behaviours of the model against both misspecification of distributions and outliers. We work out an illustrative example with a real data set on measurements of plant root decomposition.     

A risk set calibration method for failure time regression by using a covariate reliability sample

Regression parameter estimation in the Cox failure time model is considered when regression variables are subject to measurement error. Assuming that repeat regression vector measurements adhere to a classical measurement model, we can consider an ordinary regression calibration approach in which the unobserved covariates are replaced by an estimate of their conditional expectation given available covariate measurements. However, since the rate of withdrawal from the risk set across the time axis, due to failure or censoring, will typically depend on covariates, we may improve the regression parameter estimator by recalibrating within each risk set. The asymptotic and small sample properties of such a risk set regression calibration estimator are studied. A simple estimator based on a least squares calibration in each risk set appears able to eliminate much of the bias that attends the ordinary regression calibration estimator under extreme measurement error circumstances. Corresponding asymptotic distribution theory is developed, small sample properties are studied using computer simulations and an illustration is provided.     

THE SAMPLING PERFORMANCE OF INEQUALITY RESTRICTED AND PRE-TEST ESTIMATORS IN A MIS-SPECIFIED LINEAR MODEL

The literature on the sampling properties of the inequality restricted and pre-test estimators typically assumes a properly specified model and focuses on the estimation of the regression coefficient vector. In this paper, we derive and evaluate the risk functions of these estimators for both the prediction vector and the disturbance variance in a model which is mis-specified through the exclusion of relevant regressors. The results suggest that unrestricted estimation is generally preferable to pre-testing or naively imposing restrictions.     

Estimation of a Gamma Scale Parameter Under Asymmetric Squared Log Error Loss

Abstract

Estimation of scale parameter under the squared log error loss function is considered with restriction to the principle of invariance and risk unbiasedness. An explicit form of minimum risk scale-equivariant estimator under this loss is obtained. The admissibility and inadmissibility of a class of linear estimators of the form (cT + d) are considered, where T follows a gamma distribution with an unknown scale parameter η and a known shape parameter ν. This includes the admissibility of the minimum risk equivariant estimator on η (MRE).     

Estimation of the parameter in the problem of the nile

Estimation of the parameter in the problem of the Nile is treated as a decision problem with squared error loss, It is shown that the minimum risk scale equivariant estimator dominates the incomplete sufficient unbiased estimators considered by Iwase and Seto, Sharper bounds for the equivariant estimator are derived which may be used to obtain the values of the same from the sample with sufficient accuracy.     

Risk performance of a pre-test estimator for normal variance with the Stein-variance estimator under the LINEX loss function

In this paper, we derive the exact formula of the risk function of a pre-test estimator for normal variance with the Stein-variance (PTSV) estimator when the asymmetric LINEX loss function is used. Fixing the critical value of the pre-test to unity which is a suggested critical value in some sense, we examine numerically the risk performance of the PTSV estimator based on the risk function derived. Our numerical results show that although the PTSV estimator does not dominate the usual variance estimator when under-estimation is more severe than over-estimation, the PTSV estimator dominates the usual variance estimator when over-estimation is more severe. It is also shown that the dominance of the PTSV estimator over the original Stein-variance estimator is robust to the extension from the quadratic loss function to the LINEX loss function.     

On second-order improved estimation of a gamma scale parameter

This paper concludes our comprehensive study on point estimation of model parameters of a gamma distribution from a second-order decision theoretic point of view. It should be noted that efficient estimation of gamma model parameters for samples ‘not large’ is a challenging task since the exact sampling distributions of the maximum likelihood estimators and its variants are not known. Estimation of a gamma scale parameter has received less attention from the earlier researchers compared to shape parameter estimation. What we have observed here is that improved estimation of the shape parameter does not necessarily lead to improved scale estimation if a natural moment condition (which is also the maximum likelihood restriction) is satisfied. Therefore, this work deals with the gamma scale parameter estimation as a separate new problem, not as a by-product of the shape parameter estimation, and studies several estimators in terms of second-order risk.     

Bayesian Robustness Modelling of Location and Scale Parameters

Abstract. The modelling process in Bayesian Statistics constitutes the fundamental stage of the analysis, since depending on the chosen probability laws the inferences may vary considerably. This is particularly true when conflicts arise between two or more sources of information. For instance, inference in the presence of an outlier (which conflicts with the information provided by the other observations) can be highly dependent on the assumed sampling distribution. When heavy‐tailed (e.g. t) distributions are used, outliers may be rejected whereas this kind of robust inference is not available when we use light‐tailed (e.g. normal) distributions. A long literature has established sufficient conditions on location‐parameter models to resolve conflict in various ways. In this work, we consider a location–scale parameter structure, which is more complex than the single parameter cases because conflicts can arise between three sources of information, namely the likelihood, the prior distribution for the location parameter and the prior for the scale parameter. We establish sufficient conditions on the distributions in a location–scale model to resolve conflicts in different ways as a single observation tends to infinity. In addition, for each case, we explicitly give the limiting posterior distributions as the conflict becomes more extreme.