Journal of Statistical Theory and Applications

Volume 17, Issue 2, June 2018, Pages 291 - 306

Performance of the Graybill–Deal Estimator via Pitman Closeness Criterion

Keyu Nie*, Bikas.K. Sinha, A.S. Hedayat
Department of Mathematics, Statistics, and Computer Science, University of Illinois at Chicago
Received 23 September 2016, Accepted 8 May 2017, Available Online 30 June 2018.
10.2991/jsta.2018.17.2.9How to use a DOI?
Graybill-Deal Estimator; Pitman Closeness; Meta Analysis

Pitman closeness criterion is a coverage probability-based criterion to examine the relative performances of estimators. Usually, the performance of the standard Graybill-Deal estimator of the common mean has been examined with respect to the mean squared error (variance). In this study we examine its performance with respect to the Pitman closeness criterion. Specifically, we compare a p-source based Graybill-Deal estimator against its q-source based competitors for q (< p)-source subsets of p-source data. The key references to this paper are [5] and [7].

Copyright © 2018, the Authors. Published by Atlantis Press.
Open Access
This is an open access article under the CC BY-NC license (

1. Introduction

We consider the problem of estimation of the common mean shared by several independent normal populations with unknown and most likely unequal variances. Generally, we have p independent sources with distributions N(μ,σi2) , i = 1,2,…, p. Also, let x¯i and si2 be sample mean and sample variance, respectively; and sx¯i2=si2/ni , σx¯i2=σi2/ni , where ni is the sample size available from the ith source, i = 1,2,…, p.

[1] introduced their estimator μˆGD|2 for p = 2 sources, and claimed that μˆGD|2 was preferable to both sample means with respect to the criteria of mean square error, if and only if sample sizes n1 and n2 were moderate enough (≥ 11), which was corrected by [4] as (n1 ≥ 11,n2 ≥ 11), (n1 = 10,n2 ≥ 19) or (n1 ≥ 19,n2 = 10). Subsequently this result was extended by [3] to include p independent sources. They compared the Graybill-Deal estimator (GDE) of combining p sources

with sample means x¯i , i = 1,2,…, p, and found that μˆGD|p had a smaller variance than each sample mean under the same condition as [1]; that is, either ni > 10 for i = 1,2,…, p, or ni = 10 for some i, and nj > 18 (i, j = 1,2,…, p for each ji). [8] showed that μˆGD|p dominated any μˆGD|q of q (< p) sub-sources, under the same condition as in [3], with respect to mean square criterion.

Here we shall compare μˆGD|p with μˆGD|q by employing Pitman closeness criterion, which was introduced by [6]: We say estimator μˆ1 is better (Pitman-closer) than μˆ2 for the estimation of the parameter µ if and only if P{|μˆ1μ||μˆ2μ|}1/2 . Pitman-closeness criterion is a coverage probability-based criterion and has nothing to do with the loss function. It provides another angle in the sense of ’one-to-one’ comparison measuring the performance efficiency of GDE.

Using the Pitman-closeness criterion, [5] and [7] established that for p = 2, a necessary and sufficient condition for

to hold uniformly in (μ,σx¯12,σx¯22) is that mi = ni − 1 ≥ 4 for each i = 1,2. [7] further established that
for all i = 1,2,…, p and uniformly in (μ,σx¯12,σx¯22,,σx¯p2) , if and only if
holds for all i = 1,2,…, p and uniformly in (σx¯12,σx¯22,,σx¯p2) , where Yj’s are independently distributed as mjχ2(mj) and mj = nj − 1, for j = 1,2,…, p. [7] also gave a sufficient condition by showing that inequality (1.1) holds, for any i = 1,2,…, p, if
for all k in {1,2,…,i − 1,i + 1,…, p}.

In this paper we show that μˆGD|p is a Pitman-closer estimator than μˆGD|q , the GDE of any q (< p) sub-sources, in the sense of Pitman closeness criteria for all p > 2. Without loss of generality, our results are presented by comparing with initial q sub-sources, although for any other q sub-sources, we will obtain similar results.

We provide several sufficient conditions [including the inequality (1.1) (suitably modified) which holds uniformly in (σx¯12,σx¯22,,σx¯p2) for all q sub-sources]. We also provide a necessary condition towards this.

The remainder of this paper is structured as follows. Section 2 describes all the lemmas and preparation for the main results. Section 3 provides a necessary condition when combining all p-sources is preferred. Section 4 presents the corresponding sufficient condition and several corollaries. Section 5 discusses the sample size requirement. Section 6 concludes the paper.

2. Notations and Lemmas

Before introducing our main results, we first introduce the following notations and lemmas.

We borrow the definition from [7]:


where Yi’s are independently distributed and mi = ni − 1 is the degrees of freedom. Obviously we have E(Yi1)=1 .

Let p = q + r and define:

(Uq, Vq) and (Ur, Vr) are mutually independent. As we mentioned before, Uq and Vq represented the initial q subgroups, and hence Ur and Vr represented the remaining r subgroups.

Lemma 2.1.

i=1hσx¯i2Uh2mini=1,2,,h{σx¯i2} holds uniformly in (σx¯12,σx¯22,,σx¯h2) .


Due to Cauchy-Schwarz inequality, we have:


On the other side,


Lemma 2.2.

The inequality E{Uh1}(i=1hσx¯i2)1/2E{Vh} holds uniformly in (σx¯12,σx¯22,,σx¯h2) .


Due to the fact:


The inequality (2.1) is based on Lemma 2.1, and inequality (2.2) is due to Jensen’s inequality.

Lemma 2.3.

The inequality (i=1hσx¯i2)E(Uh1)E(Vh1) holds uniformly in (σx¯12,σx¯22,,σx¯h2) .


See [7].

Lemma 2.4.

The probability P{|μˆGD|pμ||μˆGD|qμ|} is 1/π(E{arctan(UrUq1)}+E{arctan(UrUq1+2VrVq1)}) .


Without loss of generality, it is clear µ may be assumed to be 0. As p = q + r, we need to compute the probability P{|μˆGD|p||μˆGD|q|} .

By applying the fact that:

the probability P{|μˆGD|p||μˆGD|q|} can be written as:

Note that conditionally given sx¯i2 for i = 1,…, p, the ratio

follows a Cauchy distribution.

Therefore, the probability in Eq.(2.3), can be denoted as


3. When is combining preferable ? A Necessary Condition

Theorem 3.1.

P{|μˆGD|pμ||μˆGD|qμ|}1/2 holds uniformly in (μ,σx¯12,σx¯22,,σx¯p2) , only if

which is equivalent to stating:
holds uniformly in (σx¯12,σx¯22,,σx¯p2) .


According to Lemma 2.4, the probability P{|μˆGD|pμ||μˆGD|qμ|} can be denoted as:


The expectations are taken with respect to the independent random variables Yi’s.

Let σx¯i2=σx¯12τi , for i = 1,2,…,q; and σx¯i2=σx¯12τi , for i = q+1,q+2,…,q+r. Here τi’s are in the range of [0,∞), except τ1 = 1. a So Uq=σx¯11Uq* , Vq=σx¯1Vq* , Ur=σx¯1Ur* and Vr=σx¯11Vr* , where


The probability in equation (3.2) will be


The dominated converge theorem implies that:


Next, we will show that γ ≥ 1/2 uniformly in (σx¯12,σx¯22,,σx¯p2) , only if 2E{Ur}E{Uq1)E{Vq}E{Vr1} holds uniformly in (σx¯12,σx¯22,,σx¯p2) .

The derivative of γ with respect to σx¯12 , is given by


To check whether the inequality:

is a necessary condition, we assume that inequality (3.4) is not true for some (τ2,…,τp). Then we notice that
which leads to γ being non-decreasing and contradicts that γ ≥ 1/2 for σx¯12 in a small region near 0.

It is easy to verify that inequality (3.4) is equivalent to inequality (3.1).

4. When is combining preferable ? Several Sufficient Conditions

Theorem 4.1.


then P{|μˆGD|pμ||μˆGD|qμ|}1/2 holds uniformly in (μ,σx¯12,σx¯22,,σx¯p2) .


Continue from equation (3.3) in Theorem 3.1, we obtain the following:


Hence if

then σx¯12πγ0 , which leads to γ(σx¯12,τ2,,τp) being non-decreasing in σx¯12 for any of (τ2,…,τp). The dominated converge theorem implies that γ(σx¯12)1/π(E{arctan0}+E{arctan})=1/2 as σx¯120 . For any other finite (τ2,…,τp) and σx¯12 , we have γ(σx¯12,τ2,,τp)γ(0,τ2,,τp)=1/2 . At the boundary where one or some τi’s go to ∞, we have γ(σx¯12,τ2,,τp)=1or1/2>1/2 . Hence γ ≥ 1/2 holds uniformly in (σx¯12,σx¯22,,σx¯p2) .

Based on Theorem 4.1, we derive several sufficient conditions in the sequel for P{|μˆGD|pμ||μˆGD|qμ|}1/2 to hold uniformly in (μ,σx¯12,σx¯22,,σx¯p2) .

Corollary 4.2.

is a sufficient condition for P{|μˆGD|pμ||μˆGD|qμ|}1/2 to hold uniformly in (μ,σx¯12,σx¯22,,σx¯p2) .


It is easy to state that

is equivalent to

From Theorem 4.1, we only need to show the validity of inequality (4.1).

At the left side of inequality (4.1), we have the following inequality for the second term:


This is because of Ur*2i=q+1q+rτi and (Uq*)2(mini=1,,q(τi))1 following by Lemma 2.1.

We also have:


The inequality (4.3) follows from Jensen’s inequality:

where a=σx¯14(i=q+1q+rτi)mini=1,,q(τi) , A=mini=1,,q(τi)(i=q+1q+rτi)Ur*(Uq*)1 , and B=Vr*(Vq*)1 .

From Lemma 2.3, we know


From Lemma 2.2 and the fact that Uq*mini=1,,q(τi)Uq*1 , we note


So we have the following :


Then at the left side of inequality (4.1), we have the following inequality for the first term:


Combining the above and inequality (4.2), we have




Therefore a sufficient condition for P{|μˆGD|pμ||μˆGD|qμ|}1/2 to hold uniformly in (μ,σx¯12,σx¯22,,σx¯p2) is:


Another corollary is given in the following:

Corollary 4.3.

is a sufficient condition for P{|μˆGD|pμ||μˆGD|qμ|}1/2 to hold uniformly in (μ,σx¯12,σx¯22,,σx¯p2) .


We need to validate inequality (4.1) in Theorem 4.1. On the left side of inequality (4.1), we have the following inequality for the second term:


It is because of the facts Ur*2i=q+1q+rτi and Uq*2i=1qτi following from lemma 2.1.

We also have:


The inequality (4.5) is from Jensen’s inequality.

So at the left side of inequality (4.1), we have the following inequality for the first term:


From Lemma 2.3 we have:


The upper bound of E(Uq*) and (i=1qτi)E(Vq*) is i=1qτi . This is due to Lemma 2.1 and Lemma 2.2, respectively.

Another upper bound for (i=1qτi)E(Vq*) is E(Vq*1) . Since

then we have

Combining inequalities (4.4), (4.6) and (4.7), we have the following:

while combining inequalities (4.4), (4.6) and (4.8), we obtain

So if


Hence it is equivalent to saying that a sufficient condition for P{|μˆGD|pμ||μˆGD|qμ|}1/2 to hold uniformly in (μ,σx¯12,σx¯22,,σx¯p2) is:


5. Sample Size Discussion

In this section we will discuss the sample size requirement to ensure the sufficient condition provided in Corollary 4.3:

Which is the same as:
θi=σx¯i2 , i = 1,2,…, p. We specified the following theorem:

Theorem 5.1.

Let Y1,…,Yp be independent random variables such that Yi~miχ2(mi) , for i = 1,2,…, p,and let p = q + r. Then for any θi > 0 (i = 1,2,…,p),

if there exists 1 < b < 2, such that
for all i = 1,2,…,q, and
for all i = q + 1,q + 2,…,q + r = p.



for i = 1,2,…,q; and define:
for i = q + 1,q + 2,…,q + r = p.

From [7] and [2], we know the following for i = 1,2,…,q:

and similarly for i = q + 1,q + 2,…,q + r = p, we have:

For g1(Yi), we have that:


Based on the above, we see the following


From Cauchy-Schwarz inequality, we can tell:


Applying (5.3) into (5.2), we have the following:


Similarly for g2(Yi), we have that:


Thus, inequality (5.1) is true if there exists 1 < b < 2, such that

for all i = 1,2,…,q, and
for all i = q + 1,q + 2,…,q + r = p.

Remark 1

Observe that inequality (5.4) and inequality (5.5) are respectively equivalent to:


In order to have solutions for inequalities (5.6) and (5.7), we need 1 < b < 2.

Remark 2

One symmetric solution is to let b=2 . Check the Condition (5.6),

  1. (1)

    When q = 1, Condition (5.6) holds if m14227 .

  2. (2)

    When q = 2, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 11, and the other 11 ≤ mj,jmin ≤ 17

    2. (b)

      mmin = 12, and the other 12 ≤ mj,jmin ≤ 37

    3. (c)

      mmin = 13, and the other 13 ≤ mj,jmin ≤ 122

    4. (d)

      mi ≥ 14 for i = 1,2.

  3. (3)

    When q = 3, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 12, and j=13mj,jmin37 . So either all are equal to 12 or some mi = 13 and the other two are each equal to 12.

    2. (b)

      mmin = 13, and j=13mj,jmin122

    3. (c)

      mi ≥ 14 for all i = 1,2,3.

  4. (4)

    When q = 4, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 13, and j=14mj,jmin122

    2. (b)

      mi ≥ 14 for all i = 1,2,3,4.

  5. (5)

    When q = 5, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 13, and j=15mj,jmin122

    2. (b)

      mi ≥ 14 for all i = 1,2,3,4,5.

  6. (6)

    When q = 6, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 13, and j=16mj,jmin122

    2. (b)

      mi ≥ 14 for all i = 1,2,3,4,5,6.

  7. (7)

    When q = 7, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 13, and j=17mj,jmin122

    2. (b)

      mi ≥ 14 for all i = 1,2,3,4,5,6,7.

  8. (8)

    When q = 8, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 13, and j=18mj,jmin122

    2. (b)

      mi ≥ 14 for all i = 1,2,3,4,5,6,7,8.

  9. (9)

    When q = 9, Condition (5.6) holds if one of the following holds:

    1. (a)

      mmin = 13, and j=19mj,jmin122

    2. (b)

      mi ≥ 14 for all i = 1,2,3,4,5,6,7,8,9.

  10. (10)

    When q ≥ 10, Condition (5.6) holds if mi ≥ 14 for all i = 1,2,…,q.

Similar discussion can be found for condition (5.7).

6. Conclusion

In this paper we compared the p-source based GDE with its q-sub-source based competitors under Pitman closeness criterion. We established a necessary condition and several sufficient conditions for the p-source based GDE to be Pitman closer than its q-sub-source based GDE. We further discussed the sample size requirement corresponding to each source, and we found that one sufficient condition is ni ≥ 15 for i = 1,2,…, p.

This sample size requirement is relatively close to the requirement based on mean square errors. Hence in our point of view, the p-source-based GDE dominates any other q-sub-source based GDE not only in terms of mean square error loss function but also in the sense of Pitman closeness criterion of probability coverage, when the sample size of each source is moderately large enough.


Research is supported by the U.S. National Science Foundation Grants DMS-1306394 and DMS-1546674. This work was undertaken during regular visits of Bikas.K. Sinha at University of Illinois, Chicago. He is a retired professor of Indian Statistical Institute, Kolkata. The authors are extremely thankful to two anonymous referees for critically reading the manuscript and offering numerous suggestions towards clearer presentation of our ideas.



In inequality (3.4), we are going to prove, 2E{Ur}E{Uq1}E{Vq}E{Vr1} , it will follow that U’s and V1s are proportional to σ2 .

Journal of Statistical Theory and Applications
17 - 2
291 - 306
Publication Date
ISSN (Online)
ISSN (Print)
10.2991/jsta.2018.17.2.9How to use a DOI?
Copyright © 2018, the Authors. Published by Atlantis Press.
Open Access
This is an open access article under the CC BY-NC license (

Cite this article

AU  - Keyu Nie
AU  - Bikas.K. Sinha
AU  - A.S. Hedayat
PY  - 2018
DA  - 2018/06/30
TI  - Performance of the Graybill–Deal Estimator via Pitman Closeness Criterion
JO  - Journal of Statistical Theory and Applications
SP  - 291
EP  - 306
VL  - 17
IS  - 2
SN  - 2214-1766
UR  -
DO  - 10.2991/jsta.2018.17.2.9
ID  - Nie2018
ER  -