Fundamental Equations of BLUE and BLUP in the Multivariate Linear Model with Applications

Ordinary least squares estimator (OLSE), best linear unbiased estimator (BLUE), and best linear unbiased predictor (BLUP) in the general linear model with new observations are generalized to the general multivariate linear model. The fundamental equations of BLUE and BLUP in the multivariate linear model are derived by two methods, including the vectorization method and projection method. By using the matrix rank method, some new results of linear BLUE-sufficiency, linear BLUP-sufficiency, and the equality of OLSE, BLUE, and BLUP are given in the multivariate linear model.     

Optimum linear unbiased estimation of scale parameter by absolute values of order statistics in symmetric distributions

A new estimator of the scale parameter by the optimum linear combination of absolute values of order statistics in symmetric location-scale families with known location parameter (without loss of generality assumed to be zero) from complete and Type II censored samples is introduced and is termed as optimum unbiased absolute estimator of the scale parameter. The new estimator of the scale parameter is compared with the corresponding best linear unbiased estimator (BLUE) in the rectangular and normal distributions. Generally it is found that the new estimator is more efficient than the BLUE.     

A projector oriented approach to the best linear unbiased estimator

The paper provides a projector based approach to the best linear unbiased estimator (BLUE). By revisiting the so called generalized projection operator, introduced in Rao (J R Stat Soc Ser B Stat Methodol 36:442–448, 1974), a number of new formulae for BLUE is established. Furthermore, some attention is paid to the coincidence of the BLUE and the ordinary least squares estimator.     

The admissible minimax estimator in Gauss–Markov model under a balanced loss function

In this paper, we investigate the admissible minimax estimator (AME) of regression coefficient in Gauss–Markov model under a balanced loss function. In the class of homogeneous linear estimators, we obtain the AME under two occasions, respectively. We also prove that the AME is a shrinkage estimator of the best linear unbiased estimator (BLUE). Furthermore, we prove that the AME dominates the BLUE under certain conditions.     

On equality of ordinary least squares estimator,best linear unbiased estimator and best linear unbiased predictor in the general linear model

The equality of ordinary least squares estimator (OLSE), best linear unbiased estimator (BLUE) and best linear unbiased predictor (BLUP) in the general linear model with new observations is investigated through matrix rank method, some new necessary and sufficient conditions are given.     

Best unbiased estimation of fixed effects in mixed ANOVA models

In a linear model with an arbitrary variance–covariance matrix, Zyskind (Ann. Math. Statist. 38 (1967) 1092) provided necessary and sufficient conditions for when a given linear function of the fixed-effect parameters has a best linear unbiased estimator (BLUE). If these conditions hold uniformly for all possible variance–covariance parameters (i.e., there is a UBLUE) and if the data are assumed to be normally distributed, these conditions are also necessary and sufficient for the parametric function to have a uniformly minimum variance unbiased estimator (UMVUE). For mixed-effects ANOVA models, we show how these conditions can be translated in terms of the incidence array, which facilitates verification of the UBLUE and UMVUE properties and facilitates construction of designs having such properties.     

Linear approximate ML estimation in scaled Type I generalized logistic distribution based on Type-II censored samples

The scaled (two-parameter) Type I generalized logistic distribution (GLD) is considered with the known shape parameter. The ML method does not yield an explicit estimator for the scale parameter even in complete samples. In this article, we therefore construct a new linear estimator for scale parameter, based on complete and doubly Type-II censored samples, by making linear approximations to the intractable terms of the likelihood equation using least-squares (LS) method, a new approach of linearization. We call this as linear approximate maximum likelihood estimator (LAMLE). We also construct LAMLE based on Taylor series method of linear approximation and found that this estimator is slightly biased than that based on the LS method. A Monte Carlo simulation is used to investigate the performance of LAMLE and found that it is almost as efficient as MLE, though biased than MLE. We also compare unbiased LAMLE with BLUE based on the exact variances of the estimators and interestingly this new unbiased LAMLE is found just as efficient as the BLUE in both complete and Type-II censored samples. Since MLE is known as asymptotically unbiased, in large samples we compare unbiased LAMLE with MLE and found that this estimator is almost as efficient as MLE. We have also discussed interval estimation of the scale parameter from complete and Type-II censored samples. Finally, we present some numerical examples to illustrate the construction of the new estimators developed here.     

Optimum linear unbiased estimation of the scale parameter by absolute values of order statistics in the double exponential and the double weibull distributions

The first two moments and product moments of absolute values of order statistics are obtained for the double exponential and the double Weibull distributions. In both of the distributions an optimum linear unbiased estimator of the scale parameter, by absolute values of the order statistics, is obtained from complete and censored samples of size n=3(1)10. It is found that the new estimator is generally more efficient than the best linear unbiased estimator (BLUE) of the scale parameter by order statistcs in both of the distributions.     

Estimation of regression vectors in linear mixed models with Dirichlet process random effects

The Dirichlet process has been used extensively in Bayesian non parametric modeling, and has proven to be very useful. In particular, mixed models with Dirichlet process random effects have been used in modeling many types of data and can often outperform their normal random effect counterparts. Here we examine the linear mixed model with Dirichlet process random effects from a classical view, and derive the best linear unbiased estimator (BLUE) of the fixed effects. We are also able to calculate the resulting covariance matrix and find that the covariance is directly related to the precision parameter of the Dirichlet process, giving a new interpretation of this parameter. We also characterize the relationship between the BLUE and the ordinary least-squares (OLS) estimator and show how confidence intervals can be approximated.     

Optimal estimation in rotation patterns

The aim of this paper is to examine the setting of surveys repeated over time when the elements in the sample are rotated in a predesigned way. On each occasion the best linear unbiased estimator (BLUE) of the current population mean, built on all past responses, is to be found. The most straightforward approach would be to compute the estimator as a solution of a least squares problem with linear restrictions. However, this method has certain drawbacks related to the fact that the size of the response data set increases over time. We follow a different approach based on finding linear recurrence relationships between optimal estimators obtained on successive occasions. Most of the original disadvantages are then corrected. In this context we present the solution to the BLUE estimation problem for some—sufficiently regular—classes of rotation patterns.     

Nonnegative Unbiased Estimation of Scale Parameters and Associated Quantiles Based on a Ranked Set Sample

In this article, we are interested in estimating the scale parameter in location and scale families. It is well known that the best linear unbiased estimator (BLUE) of scale parameter based on a simple random sample (SRS) is nonnegative. However, the BLUE of scale parameter based on a ranked set sample (RSS) can assume negative values. We suggest various modifications of BLUE of scale parameter based on RSS so that the resulting estimators are unbiased as well as nonnegative. Their performances in terms of relative efficiencies are compared and some recommendations are made for normal, logistic, double exponential, two-parameter exponential and Weibull distributions. We also briefly discuss an application of the proposed nonnegative BLUE of scale parameter for quantile estimation for the above populations.     

Computational comparison of the general weighted moments estimators of the scale parameter of a Pareto distribution with a known shape parameter with other estimators based on a multiply type II-censored sample

Wu et al. [Computational comparison for weighted moments estimators and BLUE of the scale parameter of a Pareto distribution with known shape parameter under type II multiply censored sample, Appl. Math. Comput. 181 (2006), pp. 1462–1470] proposed the weighted moments estimators (WMEs) of the scale parameter of a Pareto distribution with known shape parameter on a multiply type II-censored sample. They claimed that some WMEs are better than the best linear unbiased estimator (BLUE) based on the exact mean-squared error (MSE). In this paper, the general WME (GWME) is proposed and the computational comparison of the proposed estimator with the WMEs and BLUE is done on the basis of the exact MSE for given sample sizes and different censoring schemes. As a result, the GWME is performing better than the best estimator among 12 WMEs and BLUE for all cases. Therefore, GWME is recommended for use. At last, one example is given to demonstrate the proposed GWME.     

A note on the equality of the OLSE and the BLUE of the parametric function in the general Gauss–Markov model

In this note we consider the equality of the ordinary least squares estimator (OLSE) and the best linear unbiased estimator (BLUE) of the estimable parametric function in the general Gauss–Markov model. Especially we consider the structures of the covariance matrix V for which the OLSE equals the BLUE. Our results are based on the properties of a particular reparametrized version of the original Gauss–Markov model.      

Remove unwanted variation retrieves unknown experimental designs

Remove unwanted variation (RUV) is an estimation and normalization system in which the underlying correlation structure of a multivariate dataset is estimated from negative control measurements, typically gene expression values, which are assumed to stay constant across experimental conditions. In this paper we derive the weight matrix which is estimated and incorporated into the generalized least squares estimates of RUV-inverse, and show that this weight matrix estimates the average covariance matrix across negative control measurements. RUV-inverse can thus be viewed as an estimation method adjusting for an unknown experimental design. We show that for a balanced incomplete block design (BIBD), RUV-inverse recovers intra- and interblock estimates of the relevant parameters and combines them as a weighted sum just like the best linear unbiased estimator (BLUE), except that the weights are globally estimated from the negative control measurements instead of being individually optimized to each measurement as in the classical, single measurement BIBD BLUE.     

On the natural restrictions in the singular Gauss–Markov model

A Gauss–Markov model is said to be singular if the covariance matrix of the observable random vector in the model is singular. In such a case, there exist some natural restrictions associated with the observable random vector and the unknown parameter vector in the model. In this paper, we derive through the matrix rank method a necessary and sufficient condition for a vector of parametric functions to be estimable, and necessary and sufficient conditions for a linear estimator to be unbiased in the singular Gauss–Markov model. In addition, we give some necessary and sufficient conditions for the ordinary least-square estimator (OLSE) and the best linear unbiased estimator (BLUE) under the model to satisfy the natural restrictions.      

A gm estimation of the location parameters in a spatial linear model

Universal kriging is a form of interpolation that takes into account the local trends in data when minimizing the error associated with the estimator. Under multivariate normality assumptions, the given predictor is the best linear unbiased predictor. but if the underlying distribution is not normal, the estimator will not be unbiased and will be vulnerable to outliers. With spatial data, it is not only the presence of outliers that may spoil the predictions, but also the boundary sites. usually corners, that tend to have high leverage. As an alternative, a weighted one-step generalized M estimator of the location parameters in a spatial linear model is proposed. It is especially recommended in the case of irregularly spaced data.     

Estimation of quantiles of exponential and double exponential distributions based on two order statistics

Asymptotically best linear unbiased estimators (ABLUE) of quantiles, x^., in the two-parameter (location-scale) exponential and double exponential families are obtained as linear combinations of two suitably chosen order statistics. Exact formulae for the linear combinations are given as functions of £. The derived estimators in both cases compare favorably with the usual nonparametric estimator. Also, in the exponential case the derived estimator compares favorably with the Sarhan-Greenberg BLUE based on a complete sample     

The BLUE's covariance matrix revisited: A review

In this paper we comment on and review some unexpected but interesting features of the BLUE (best linear unbiased estimator) of the expectation vector in the general linear model and in particular, the BLUE's covariance matrix. Most of these features appear in the literature but are rather scattered or hidden.     

Nonlinear unbiased estimation in the linear regression model with nonnormal disturbances

In the application of the linear regression model there continues to be wide-spread use of the Least Squares Estimator (LSE) due to its theoretical optimality. For example, it is well known that the LSE is the best unbiased estimator under normality while it remains best linear unbiased estimator (BLUE) when the normality assumption is dropped. In this paper we extend an approach given in Knautz (1993) that allows improvement of the LSE in the context of nonnormal and nonsymmetric error distributions. It will be shown that there exist linear plus quadratic (LPQ) estimators, consisting of linear and quadratic terms in the dependent variable, which dominate the LS estimator, depending on second, third and fourth moments of the error distribution. A simulation study illustrates that this remains true if the moments have to be estimated from the data. Computation of confidence intervals using bootstrap methods reveal significant improvement compared with inference based on the LS especially for nonsymmetric distributions of the error term.     

Editor's Notes

Two examples are given that illustrate the Gauss–Markov theorem. The examples will help students appreciate the use of the word best in connection with the term BLUE (best linear unbiased estimator). The examples will also provide several instructive exercises.