The Behrens-Fisher Problem: An Empirical Likelihood Approach

April 1, 2017 | Author: Frederica Clark | Category: N/A
Share Embed Donate


Short Description

1 Econometrics Working Paper EWP0404 ISSN Department of Economics The Behrens-Fisher Problem: An Empirical Likelihood Ap...

Description

Econometrics Working Paper EWP0404 ISSN 1485-6441

Department of Economics

The Behrens-Fisher Problem: An Empirical Likelihood Approach

Lauren Bin Dong

Department of Economics, University of Victoria Victoria, B.C., Canada V8W 2Y2

Revised, November, 2004

Abstract

A new theoretical solution to the Behrens-Fisher (BF) problem is developed using the maximum empirical likelihood method. The sampling properties of the empirical likelihood ratio (ELR) test for the BF problem are derived using Monte Carlo simulation for a wide range of situations. A comparison of the sizes and powers of the ELR test and the Welch-Aspin test is conducted for a special case of small sample sizes. The empirical results indicate that the ELR test for the BF problem has good power properties.

Keywords:

Behrens-Fisher problem, empirical likelihood, Monte Carlo simulation, testing for normality, size and power

JEL Classifications:

C12, C15

Author Contact: Lauren Dong, Statistics Canada; e-mail: [email protected]; FAX: (613) 951-3292

1

INTRODUCTION

Testing the equality of the means of two normal populations when the variances are both unknown, and not known to be equal, is called the Behrens-Fisher (BF) problem. The Behrens-Fisher problem has been well known since the early 1930’s. One reason for its fame is that it can be proven that there is no exact solution satisfying the classical criteria for good tests. That is, every invariant rejection region of fixed size for the problem must have some unpleasant properties (Zaman, 1996, p. 246). First-best solutions that are uniformly most powerful and invariant either do not exist or have strange properties. We need to look for second-best solutions. In the literature associated with the Behrens-Fisher problem, there have been quite a few “solutions” proposed since the 1930’s. For example, Fisher (1935 and 1941), Welch (1947), Aspin (1948), Cochran and Cox (1950), Qin (1991) and Jing (1995) have all suggested different solutions. Lee and Gurland (1975) presented a detailed comparison of several selected tests and proposed a refined solution to the BF problem. These various solutions can be classified into two categories: approximate tests and asymptotic tests. The purpose of this paper is as follows. First, we derive a new theoretical solution to the Behrens-Fisher problem. The approach that we take is based on the relatively new nonparametric technique, the maximum empirical likelihood (EL) method (Owen, 2001). The way that we exploit the data and information in this paper is distinct from any previous research in both the areas of EL and the BF problem. Our EL approach makes most efficient use of the information available among the studies. Our EL approach ties together nicely the estimation and testing procedures for the equality of two means and/or the distribution of the underlying population. Second, we provide sampling properties for the ELR test using Monte Carlo simulations. The size and the power of the ELR test in finite samples are provided for a range of situations. The results indicate that the EL approach to the BF problem is both efficient and easily applicable. Third, we conduct a power comparison of the empirical likelihood ratio (ELR) test and the Welch-Aspin (WA) (1947) test for the BF problem. The results are interesting, as the ELR test is an asymptotic test, while the WA test is an approximate test in small samples. 1

The EL method makes use of likelihood function and the moment conditions of the data without assuming a specific parametric form for the underlying distribution. It has the flexibility to incorporate various information about the data into the approach. This leads to the efficiency of the method. Applying the EL method to the Behrens-Fisher problem provides us with a new view of the EL method and it demonstrates that the EL method is a useful tool for solving various statistical problems. There have been two empirical likelihood type (EL-type) approaches to the BF problem in the literature; that of Qin (1991) and that of Jing (1995). Details of these two approaches are provided in Section 2.2. The EL approach of this paper is quite distinct from those of Qin and Jing in two respects. First, our way of using the data is different. Based on the knowledge that two samples are independently drawn from two different distributions that are in the same family, normal for the BF problem, we transform the first data set S1 into a data set Sa that has the same theoretical distribution as the second data set S2 . Then we combine the transformed data set Sa and the second data set S2 into a full data set S such that the data have a unique distribution. Second, we exploit the data information in an more efficient way. Based on the knowledge that the first five moments of the data exist, we use these five moment equations as constraints to set up the EL method. Then, we apply the usual EL approach to the full data set S. The ELR test for the BF problem is then constructed. Details of this are provided in Section 2.3. The outline of this paper is as follows. Section 2 provides a brief review of the conventional solutions to the BF problem in the literature. The Welch-Aspin (WA) test and the approaches of Qin (1991) and Jing (1995) are discussed in this section. Section 3 discusses the design of the new EL approach to the BF problem. It provides the ELR test and the EL-type Wald test for the BF problem. Section 4 presents some Monte Carlo experiments and the associated results for the new EL approach. The sizes and the power of the ELR test in finite samples are analyzed in detail across a broad range of situations. The comparison of the ELR test and the WA test is given in this section. Section 5 provides a summary and some conclusions.

2

2

SOLUTIONS TO THE BF PROBLEM

Suppose Si = {xi1 , xi2 , . . . , xini }, i = 1, 2 are two sets of the observed values of two random samples independently drawn from the normal populations N (µi , σi2 ), i = 1, 2, respectively. Let σ2 ρ2 = 22 , (1) σ1 x¯i = s2i

ni 1 X xij , ni j=1

ni 1 X = (xij − x¯i )2 , ni − 1 j=1

C0 =

σ12 n1 σ12 n1

+

=

σ22 n2

1 1 + ρ2 nn12

(2)

(3)

(4)

denote the variance-ratio, sample means, sample variances, and the conventional parameter used in the literature. The problem in hand is to construct a test of the null hypothesis, H0 : µ1 = µ2 , against the alternative hypothesis, Ha : µ1 6= µ2 . Various tests have been developed by Fisher (1935), Welch and Aspin (1947), Cochran and Cox (1950) and others. The critical regions corresponding to each of these conventional tests have the general form of: ˆ |v| > V (C), (5) 2

2

2

2

2

s s s s s ˆ is a function of Cˆ and the where v = (¯ x1 − x¯2 )( n11 + n22 )−1 , Cˆ = ( n11 )( n11 + n22 )−1 . V (C) preassigned nominal significance level α. The distribution of v is no longer Student-t when the two variances are not known to be equal. The critical regions depend on the variance parameters σi2 , i = 1, 2 and sample sizes (n1 , n2 ). These conventional solutions to the BF problem involve utilizing various means to approximate the distribution of the variable v and to control for the size distortion of the tests. The formulars and methods of approximation are too complex to repeat here. Further details and discussion can be found in Lee and Gurland (1975). We have chosen Welch-Aspin test as an example for a brief review.

2.1

Welch-Aspin Test

Welch (1947) and Aspin (1948) independently developed a higher order approximation, ˆ to the distribution of V (C) ˆ in the terms of up to fi−2 and fi−4 , where fi = ni − 1 Vwa (C), 3

and i = 1, 2. Their approach is referred to as the Welch-Aspin (WA) test. The WA test is a highly efficient solution to the BF problem (Weerahandi, 1987). The actual sizes of the Welch-Aspin test in small samples are very close to the nominal significance levels. For example, at a nominal level of 5%, the size of the WA test lies between 4.98% and 5.02% ˆ is lengthy, involving infinite for (n1 , n2 ) = (7, 7). However, the functional form of Vwa (C) series, and it is difficult to work with (Lee and Gurland, 1975). Lee and Gurland (1975) provided detailed comparisons of various tests that were proposed by Fisher (1935), Cochran and Cox (1950), Welch (1937), Welch-Aspin (1947) and others. At the last stage of their comparison, all of the tests were eliminated from their tables, except the Welch-Aspin test because the WA test has very accurate size and good power properties. In addition to this, Lee and Gurland proposed a refined test, which we call the LG test, using two techniques: (i) a simple functional form to approximate the V (C) function; and (ii) a T-transform to accelerate the convergence of the size and the power functions. The method in their research involves solving a minimization problem of the squared difference between the size function of the test and the preassigned nominal level, and the T-transform is not trivial. Lee and Gurland provided refined results that were very close to the results of the WA test. For the reasons mentioned above, we have chosen the WA test and to cite the results of the WA test from Lee and Gurland for our comparison with the ELR test for the case of (n1 , n2 ) = (7, 7) at α = 5%. Of course, we also consider other sample sizes and nominal significance levels when considering the ELR test. A detailed discussion is given in Section 3.

2.2

Approaches of Qin and Jing

Jing (1995) proposed a nonparametric version of the EL approach to the two-sample problem. In his paper, he showed that the nonparametric version of Wilks’ theorem holds in the twosample problem and that the solution is Bartlett correctable. The empirical likelihood ratio statistic has a limiting chi-squared distribution with an error of order O(n−1 ) and with the Bartlett correction, the error is reduced to the order of O(n−2 ), where n is the smaller one of the two sample sizes. A special case of his approach is a particular solution to the Behrens-Fisher problem. The focus of Jing (1995) was primarily on the coverage accuracy and the Bartlett cor-

4

rection of the test. The approach he used works only when the null hypothesis is true. This means that the restriction that two means are equal must be binding, and the Lagrangian multiplier for the constraint must not equal zero. If the restriction is not binding, the solution of his approach reduces to a situation where the probability parameter pi = 1/n, and the estimated empirical mean of the data is just the sample mean. His approach does not appear to offer the potential for dealing with issues relating to power in the context of an ELR test. Qin (1991) generalized Owen’s empirical likelihood to a two-sample problem in which one sample is assumed to come from a distribution that is unknown and the other sample is assumed to come from a known distribution specified upto a parameter, gθ (x2 ). Qin’s approach is a semi-parametric one by combining the EL method and the parametric likelihood method to the two-sample problem. Consider the following assumptions: • µ2 = µ(θ) is differentiable and µ0 (θ) 6= 0 • the known density function gθ (x2 ) is differentiable three times with respect to θ •

R

gθ (x2 )dx2 is twice differentiable

gθ 2 ] , satisfies 0 < I(θ) < ∞ • The Fisher information matrix, I(θ) = E[ ∂ log ∂θ

• |

∂3 ∂θ3

log gθ (x2 ) | is bounded.

Under these assumptions, there exists an EL estimator θˆ of the mean µ1 that is more efficient than the sample mean x¯1 . Qin proved that the empirical likelihood ratio has a limiting distribution of χ2(1) under the null hypothesis µ1 = µ2 = µ(θ). The coverage accuracy is of 1 the order Op (n− 2 ). The connection between the two data sets is brought in only through the single restriction that the two means are equal. The EL approach developed in this paper provides a new theoretical solution to the BF problem and hopefully it can overcome some of the shortcomings mentioned above. The focus of the following sections is on deriving the ELR test for the BF problem and simulation of the sampling properties for the ELR test in the context of solving the Behrens-Fisher problem.

5

3 3.1

THE EL APPROACH

ELR Test

Suppose Si = {xi1 , xi2 , . . . , xini }, i = 1, 2, are the two data sets we have. They are independently drawn from two normal populations N (µi , σi2 ), i = 1, 2. The parameters (µ1 , µ2 , σ12 , σ22 ) are unknown. ni ’s are the sample sizes of the data sets. Without losing generality, we assume that n1 ≥ n2 . Solving the BF problem involves constructing a test for the hypothesis that the two population means equal when the two variances are not known to be equal. The steps of our approach are as follows. First, we transform the data set S1 , into a data set that has the same theoretical distribution as the data set S2 . The transformation has the following form: 1

tj = (x1j − µ1 )(ρ2 ) 2 + µ2 ,

(6)

where x1j ∈ S1 and ρ2 = σ22 /σ12 . Thus, tj ∼ N (µ2 , σ22 ), j = 1, 2, . . . , n1 . The second data set S2 remains unchanged. We denote: tn1 +j = x2j , where x2j ∈ S2 , for j = 1, 2, . . . , n2 . Let n = n1 + n2 . Then, the full data set S = {t1 , t2 , . . . , tn } is of size n and has a distribution of N (µ2 , σ22 ). σ2

The parameter C0 = 1/(1 + σ22 nn21 ) is the one that is used in the literature of the BF 1 problem, (e.g., Welch (1947), Lee and Gurland (1975)), where the σi2 ’s are the population variance parameters, i = 1, 2. Obviously, C0 ∈ (0, 1), and it is a smooth function of the variance ratio parameter, ρ2 , and the ratio of the sample sizes, n1 /n2 . In our study, we keep the ratio of the sample sizes constant and allow the variance ratio to vary such that the corresponding values of C0 cover well the range of (0, 1). The next step is to apply the EL approach to the full data set S. A probability parameter pj is assigned to each data point tj . The empirical likelihood function for the full data set Q is formed as nj=1 pj . The approach of the empirical maximum likelihood method is to maximize the empirical likelihood function subject to the probability constraints, 0 < pj < 1 P and nj=1 pj = 1, and some information constraints. The information that we have is that the data are independent and are normally dis-

6

tributed with a mean µ2 and a variance σ22 . We choose to use the ratio of the variances as one of the parameters, so the parameter vector becomes θ = (µ1 , µ2 , ρ2 , σ22 )0 . In addition, we choose to use the first five unbiased raw moment equations as the information constraints, i.e. m = 5, so that the number of the moment equations, m, is greater than the number of the parameters, p. As we know, if m < p, the system is under-identified; and if m = p, the EL solution becomes exactly the solution from the method of moments. Therefore, five moment equations are necessary and sufficient for the EL method to be effective. We denote the information constraints as Ep h(tj , θ) = 0, which have the following form: n X j=1 n X j=1 n X j=1 n X j=1 n X

pj tj − µ2 = 0

(7)

pj t2j − (µ22 + σ22 ) = 0

(8)

pj t3j − (µ32 + 3σ22 µ2 ) = 0

(9)

pj t4j − (µ42 + 6σ22 µ22 + 3σ24 ) = 0

(10)

pj t5j − (µ52 + 10σ22 µ32 + 15σ24 µ2 ) = 0.

(11)

j=1

The corresponding Lagrangian function is: G=n

−1

n X j=1

log pj − η(

n X

pj − 1) − λ

j=1

0

n X

pj h(tj , θ),

(12)

j=1

where λ is a m × 1 vector, which together with the scalar η are the Lagrangian multipliers. The optimization problem is to maximize the Lagrangian function with respect to pj ’s, λ, and θ. Keeping λ and θ fixed, applying the first order condition with respect to and the probability parameter constraints on pj ’s, we find that η takes the value of unity and the pj ’s can be expressed as functions of λ and θ: pj = n−1 (1 + λ0 h(tj , θ))−1 , j = 1, 2, . . . , n

(13)

Substituting this information back into the Lagrangian function, we get a optimization problem over a reduced number of unknowns, λ and θ.

7

Deriving the first order conditions of the Lagrangian function with respect to the parameter vector θ requires some special attention. We note that the values in the first portion of the data are functions of the unknown parameters. The first derivative of the data with respect to the parameters involves the following terms: ∂tj 1 = (−ρ, 1, (x1j − µ1 ) (ρ2 )−1/2 , 0)0 , j = 1, 2, . . . , n1 , ∂θ 2

(14)

where x1j ∈ S1 . Taking into account this information, the actual first order conditions for the Lagrangian function with respect to the four unknown parameters (µ1 , µ2 , ρ2 , σ22 )0 can be represented as follows: n1 X j=1 n1 X

pj (λ1 + 2λ2 tj + 3λ3 t2j + 4λ4 t3j + 5λ5 t4j )(ρ2 )1/2 = 0 pj (λ1 + 2λ2 tj + 3λ3 t2j + 4λ4 t3j + 5λ5 t4j )

j=1



n X

pj (λ1 + 2λ2 µ2 + 3λ3 (σ22 + µ22 ) + 4λ4 (3σ22 µ2 + µ32 ) + 5λ5 (3σ24 + 6σ22 µ22 + µ42 )) = 0

j=1 n1 X

1 pj (λ1 + 2λ2 tj + 3λ3 t2j + 4λ4 t3j + 5λ5 t4j )(tj − µ1 ) (ρ2 )−1/2 = 0 2 j=1 n X

pj (λ2 + 3λ3 µ2 + 6λ4 (σ22 + µ22 ) + 10λ5 (3σ22 µ2 + µ32 )) = 0.

j=1

Putting these four first order conditions and the five moment equations together, we get a system of nine nonlinear equations. We solve this system using the nonlinear equation solver “Eqsolve” in the Gauss package (Aptech Systems, 2002). ˆ and θˆ as the EL estimators, the solution from the system, for the LaWe denote λ grangian multiplier vector and parameter vector. Then, we obtain the EL estimators for the probability parameters pˆj ’s using the formula (13). The estimated maximum value of the Q likelihood function is obtained as L(Fˆ u ) = nj=1 pˆuj , where u stands for the unconstrained model. In order to construct an ELR test for the BF problem, we also need to have the solutions to the constrained case. With the null hypothesis: H0 : µ1 = µ2 , we simply substitute the restriction µ1 = µ2 into the unconstrained case to get the constrained case, where the parameter vector becomes θ = (µ, ρ2 , σ22 )0 which is of dimension three, with µ = µ1 = µ2 .

8

The five moment equations, Ep hc (tj , θ), are: n X j=1 n X j=1 n X j=1 n X j=1 n X

pj tj − µ = 0

(15)

pj t2j − (µ2 + σ22 ) = 0

(16)

pj t3j − (µ3 + 3σ22 µ) = 0

(17)

pj t4j − (µ4 + 6σ22 µ2 + 3σ24 ) = 0

(18)

pj t5j − (µ5 + 10σ22 µ3 + 15σ24 µ) = 0.

(19)

j=1

where c stands for constrained. Similarly, we have: pcj = n−1 (1 + λ0 hc (tj , θ))−1 , j = 1, 2, . . . , n

(20)

The first derivative of the first portion of the data with respect to the parameters has the following form: ∂tj 1 = (−(ρ2 )1/2 + 1, (x1j − µ) (ρ2 )−1/2 , 0)0 , j = 1, 2, . . . , n1 , ∂θ 2

(21)

where x1j ∈ S1 . The first order conditions of the Lagrangian function with respect to the three parameters are: n1 X

pj (λ1 + 2λ2 tj + 3λ3 t2j + 4λ4 t3j + 5λ5 t4j )((ρ2 )1/2 − 1)

j=1



n X

pj (λ1 + 2λ2 µ + 3λ3 (σ22 + µ2 ) + 4λ4 (3σ22 µ + µ3 ) + 5λ5 (3σ24 + 6σ22 µ2 + µ4 )) = 0

j=1 n1 X

1 pj (λ1 + 2λ2 tj + 3λ3 t2j + 4λ4 t3j + 5λ5 t4j )(tj − µ) (ρ2 )−1/2 = 0 2 j=1 n X

pj (λ2 + 3λ3 µ + 6λ4 (σ22 + µ2 ) + 10λ5 (3σ22 µ + µ3 )) = 0.

j=1

Solving the system of the nonlinear equations of the five moment equations and the three ˆ c and θˆc as the EL estimates, and then, the pˆc ’s using the first order conditions, we get λ j

9

formula (20). The estimated maximum value of the empirical likelihood function under the Q constraint is formed as L(Fˆ c ) = nj=1 pˆcj . The log likelihood ratio statistic has the form: L(Fˆ c ) −2 log R(F ) = −2 log L(Fˆ u ) = 2

n X

0 0 (log(1 + λˆc hc (tj , θˆc )) − log(1 + λˆu h(tj , θˆu )))

j=1

and it has a limiting distribution of χ2(1) under H0 .

3.2

EL-type Wald Test

The EL estimator θˆEL of the parameter vector θ is asymptotically efficient, and it has a limiting distribution of the following form: √

d

n(θˆEL − θ0 ) → N (0, Σ),

where

∂h(y, θ) ∂h(y, θ) |θ0 ]E[h(y, θ)h(y, θ)0 |θ0 ]−1 E[ |θ0 ]]−1 , (22) 0 ∂θ ∂θ and θ0 is the true value of θ. A consistent estimator of the asymptotic covariance matrix Σ can be obtained using the EL estimator θˆEL : Σ = [E[

ˆ = [E[ ∂h(y, θ) | ˆ ]E[h(y, θ)h(y, θ)0 | ˆ ]−1 E[ ∂h(y, θ) | ˆ ]]−1 . Σ θEL θEL θEL ∂θ ∂θ0

(23)

ˆ is a 4 × 4 symmetric matrix. With this, we can The estimated covariance matrix Σ easily obtain an EL-type Wald test for any linear restrictions of the parameters. Suppose cθˆ − r = 0 is a set of j linear restrictions. The EL-type Wald test has the form: ˆ 0 ]−1 (cθˆ − r), W = (cθˆ − r)0 [cΣc

(24)

and it has an asymptotic distribution of χ2(j) , if the restrictions are valid. For the BF problem, the restriction is simply µ1 = µ2 , i.e., c = {1, −1, 0, 0}, then, the 10

EL-type Wald test has the form of:

w = (µ1 − µ2 )2 (a11 + a22 − 2a12 )−1 ,

(25)

ˆ matrix. The EL-type Wald test statistic w where aij , i, j = 1, 2 are the elements of the Σ has an asymptotic distribution of χ2(1) under the null hypothesis. We have explored the EL-type Wald test for the BF problem using Monte Carlo simulations. Intuitively, the Wald test should be computationally easier than the ELR test since the Wald test involves solving the unconstrained case only. However, our study finds that both the computing time and computational difficulty associated with the Wald test are greater than that associated with the ELR test. One possible reason for this is that the Wald test involves computing the consistently estimated covariance matrix for the parameter estimators, and this involves three matrix inversions. For an ill behaved problem, such as the Behrens-Fisher problem, these matrix inversions may cause some difficulties. The results for the Wald test are not as good as those for the ELR test in the BF context. Therefore, we choose to present only the results of the ELR test in this study (The details of the Wald test results are available on request.)

3.3

Testing For Normality

Solving the Behrens-Fisher problem and testing for normality have close connections when using the EL method as described in the previous sections. By definition, the BF problem is the case of testing the equality of two means when two data sets are independently drawn from two normal distributions with general and unknown variances. It was categorized by Tukey (1954) as the fourth level problem of normal sequences of growth in consideration in the area of comparing the typical values of two populations with the aid of a sample drawn from each (Tukey, 1954, p. 713). The EL approach that we described in Section 3.1 provided a new solution to the Behrens-Fisher problem. Further more, suppose these two samples were drawn from two different distributions that are from the same family, and we have interest in making sure that this distribution family is normal. Then, our approach is perfectly suitable for us to conduct an ELR test for normality of the underlying distributions of the two data sets.

11

As described in Dong and Giles (2004), testing for the validity of the moment conditions provides a way of testing for the underlying distribution, normality in this case. Consider two data sets Si = {xi1 , xi2 , . . . , xini }, i = 1, 2. Suppose they are drawn independently from two populations F (µi , σi2 ) of the same family, where F is not known and is not necessarily normal. As described in section 3.1, we apply the EL approach to the full data set S. If the underlying distribution is normal, then the five moment equations Q constraints hold true, and the maximized likelihood function is L(Fˆc ) = nj=1 pˆj . If the hypothesis is not true, then the maximum value of the likelihood function is L(F u ) = n−n . The log empirical likelihood ratio statistic for testing for normality is of the form: −2

X

ˆ 0 h(tj , θ)). ˆ log nˆ pj = 2 log(1 + λ

(26)

It has a limiting null distribution of χ2(1) , where the number of degrees of freedom equals the number of moment equation, five, less the number of parameters, four. Usually, we may be interested in using this technique to test for normality of the two underlying populations of the same distribution family. If we are satisfied with the results and accept the null hypothesis that the underlying populations are normal, then, we can continue the process to solve the Behrens-Fisher problem. However, the first step of testing for normality introduces a pre-testing situation that the size and power of the subsequent test in the procedure may be altered, as described by Giles and Giles (1993). It is well known that sequential testing strategies can result in size and power distortions if the tests in question are not independent of each other. This issue is not explored further in this study. Alternatively, we can avoid the pre-testing issue by conducting an ELR test for normality and an EL-type Wald test for the Behrens-Fisher problem simultaneously. Using the unrestricted model and applying the EL approach described in Section 3.1, we obtain the ˆ and the estimates of the probability parameter, pˆi ’s. EL estimate of the parameter vector, θ, With the pˆi ’s, we can implement the ELR test for normality for the underlying distribution as we described in Dong and Giles (2004). In the meantime, without being influenced by the testing for normality, we can compute the consistent estimate of the covariance matrix of θˆ to perform the EL-type Wald test for the Behrens-Fisher problem as we described at the beginning of this section. If we accept the null hypothesis of the first test that the underlying distribution of the two populations is normal, then the Wald test is well-founded. 12

This alternative approach seems to be able to effectively avoid the pre-testing issue, and it may provide better results. We will explore this alternative in future research.

3.4

Advantages of the EL Approach

The empirical likelihood approach to the BF problem described in this paper utilizes more information from the data sets than do the approaches suggested by other authors. In particular, we exploit the obvious relationship between the two data sets that they both are drawn from normal distributions. We make use of the first five unbiased moment equations of the data sets. In addition, we utilize the empirical likelihood functions. These techniques allow the EL approach of our design to reach a higher asymptotic efficiency for estimation and higher testing power for the new solution to the BF problem. If the ELR test had a well controlled size and good power, then we would say that the ELR test is a good solution to the BF problem and the EL approach is valuable. Monte Carlo simulations that are described in the following section provide an extensive analysis of the sampling properties of the new ELR test.

4

MONTE CARLO EXPERIMENTS

For the Behrens-Fisher problem, the Monte Carlo experimental design is to apply the empirical likelihood method to test for the equality of two means without the knowledge of the variances of the underlying populations. Two random samples are generated independently from two normal distributions: N (µ1 , σ12 ) and N (µ2 , σ22 ) with the true values of the parameters (µ1 , µ1 , σ12 , σ12 )0 = (1, 1, 1, σ22 )0 , where σ22 varies with the variance ratio parameter ρ2 = σ22 /σ12 . The true value of the parameter ρ2 changes according to the values of {0.1, 0.5, 1, 2, 10}. When ρ2 = 1, the two variances equal; when ρ2 is farther apart from one, the two variances are further different from each other where the BF problem becomes severe. The sample size pair ranges from (20, 10) to (250, 125). The ratio of the sample sizes of the two data sets is kept constant, n1 /n2 = 2, for simplicity. These settings lead to the corresponding values of the parameter used in the 13

σ2

literature C0 = 1/(1 + σ22 nn12 ), to vary in accordance of {0.83, 0.5, .33, 0.2, 0.048} which has 1 a good coverage of the region (0, 1). Therefore, our experiments have a good coverage of the parameter space. The number of replications is set to be 5, 000. In conducting the power comparisons, the non-centrality parameter is defined as: δ=

µ1 − µ2 1

(σ12 + σ22 ) 2

.

(27)

The null hypothesis is true when δ = 0. The alternative hypothesis is when δ 6= 0. We have considered are the cases of δ varies according to the values of {1, 1.5, 2, 2.5, 3, 4}. One note we should bear in mind is that, in contrast to the non-centrality parameter used in Lee and Gurland (1975), the parameter δ in our design does not depend on the sample size pair (n1 , n2 ). Therefore, we are able to disentangle the effect upon the power of the ELR test that is arising from the non-centrality parameter δ, or that from sample size pair changes. The log empirical likelihood ratio statistic depends on the parameter ρ2 and the sample size ratio n1 /n2 . Thus, the size and the power of the ELR test are also functions of these parameters. However, the asymptotic distribution of the ELR test statistic does not depend on these nuisance parameters. The ELR test statistic has an asymptotic distribution of χ2(1) under the null hypothesis.

5 5.1

EMPIRICAL RESULTS

Sampling properties

Table 1 in the Appendix presents the computed sizes and the size-adjusted critical values for the ELR test in solving the Behrens-Fisher problem at four nominal significance levels: α = {10%, 5%, 2%, 1%}. These four levels are used so as to have a close look at the behavior of the ELR test at the tail. • The actual size of the ELR test is relatively large comparing to the nominal significance level. The size has an appropriate trend to converge to the correct nominal significance 14

level when the sample sizes increase, given a value of ρ2 . For example, the size changes from 12.14% to 10.66% when the sample pair varies from (20, 10) to (250, 125) at the nominal significance level of 5% when ρ2 = 0.1. We stopped at n = (250, 125) due to computing difficulty. It would be worthwhile to provide further empirical evidence to illustrate the size convergence as the sample size increases in future research. • The size of the ELR test is not very sensitive to the changes in the sample size pair at the nominal significance level of 10%. • The size distortion of the test declines as the value of the parameter ρ2 decreases, holding other factors fixed. For instance, the size changes from 27.36% to 12.14% when ρ2 moves from 10 to 0.1 for small sample sizes (n1 , n2 ) = (20, 10) at α = 5%. The size distortion is the worst when ρ2 = 10 however it is still within the range of our expectations for the ELR test, based on our computational experience. Tables 2 to 4 provide the full analysis of the power of the ELR test across a range of values for the non-centrality parameter δ, the parameter ρ2 , and the sample size pair (n1 , n2 ). The size-adjusted critical values are used to compute the powers of the test in every cases. These values allow us to evaluate the power of the ELR test at the actual significance levels. The findings are as follows. • First, the power increases as the parameter δ increases, given the values of ρ2 , α, and (n1 , n2 ). That is the power of the test if higher when the hypothesis is further away from the truth. For example, the power increases from 25.20% to 95.46% when the parameter δ varies from 1 to 4 at ρ2 = 0.1 and α = 5% for the sample size pair as small as (n1 , n2 ) = (20, 10). These are very encouraging results. They show that the ELR test has good power properties for quite small samples. • Second, the power of the test increases with the sample sizes. That is, the test is more powerful when we have more information in hand. For instance, the power increases from 49.26% to 99.94% when the sample size pair increases from (20, 10) to (250, 125) at ρ2 = 0.1, α = 5%, and δ = 1.5. • Third, the power is higher when the value of ρ2 is farther away from unity, given the values of the parameters of (n1 , n2 ), α, and δ. For example, the power changes from 30.16% to 90.28% when ρ2 changes from 1 to 0.1; the power increases to 66.30% when 15

the parameter ρ2 moves from 1 to 10, holding α = 5%, δ = 2.5, and (n1 , n2 ) = (20, 10) fixed. The rationale behind this result is as follows. (i) When ρ2 = 1, the variances of the two populations are unknown but equal, where the BF problem vanishes. The test for H0 : µ1 = µ2 reduces to the usual t test. (ii) When ρ2 deviates from unity, the sign of the Behrens-Fisher problem shows up, the farther away ρ2 is from unity, the severe the BF problem becomes; in this case, the ELR test becomes more powerful. In another words, under the alternative hypotheses, when the variances of the two populations are unknown and unequal, the ELR test has high power. The result tells us that the ELR test is able to capture the specific information that the samples were drawn from two populations with different variances that are unknown in the BF problem.

Overall, the power results are quite acceptable and the ELR test is recommended to be applied to the BF problem.

5.2

Comparison of the WA and the EL Tests

The WA test is designed for small samples. As discussed in Section 2.1, the solution provided by the Welch-Aspin test results from directly solving the size function equation ˆ =α P (|v| > V (C))

(28)

ˆ function in the terms of, and up to, fi−4 , where fi = ni − 1 is the and approximating V (C) number of degrees of freedom left for each data set, and i = 1, 2. We cite the experimental results of the WA test from Lee and Gurland (1975). The results were specifically for the case of (n1 , n2 ) = (7, 7), α = 5%, and C0 = {0.1, 0.2, 0.3, 0.4, 0.5}. The sample size is very small (n1 , n2 ) = (7, 7), but the test has well controlled size and good power. In this section, we conduct a comparison of the sampling properties between the ELR test and the WA test for this small sample size case and results are provided in Table 5 in the appendix. To make it comparable, the parameter ρ2 in the ELR test takes the values of {9, 4, 2.33, 1.5, 1} in corresponding to the values of the parameter C0 used by Lee and Gurland (1975). The noncentrality parameter δ takes the values of {0.378, 0.7559, 1.1339, 1.5119} in corresponding to

16

the non-centrality parameters of the WA test which changes values according to {1, 2, 3, 4}. √ The two non-centrality parameters have the following relationship: δ = δ0 / 7. From the table, we see that the ELR test, unfortunately, performs poorly for this very small-sample situation. We wish to make more comparisons of these two tests for different samples. However, to our knowledge, there is no other published information available to allow a comparison between the ELR and the WA test for sample size that is greater than seven. For the sample size pair as small as (n1 , n2 ) = (7, 7), the actual size of the ELR test obviously exceeds the nominal significance level but it is still within the usual size range of the ELR test. The size distortion of a test is the difference between the actual size of the test and the nominal significance level. Therefore, the size distortion of the ELR test is large. We compute the power of the ELR test at its actual size level of 5% by using the simulated size-adjusted critical values. The power comparison is made at the same actual size level 5%. The result shows that the power of the ELR test for such a small sample size is low and it is inferior to that of the WA test. It is unfortunate that the ELR test is not able to show its merits for this extremely small sample size pair. The ELR test for the Behrens-Fisher problem is an asymptotic test. The power performance is acceptably good when the sample size pair is as small as (n1 , n2 ) = (20, 10). However, we would not expect it to perform very well when the sample size pair is extremely small, such as (n1 , n2 ) = (7, 7). The power results from Monte Carlo experiment for this specific case (but with various significance levels) are presented in Table 6. Although this particular comparison between the ELR test and the WA test is rather disappointing, it must be kept in perspective. First, we are comparing our asymptotic test with one which is explicitly designed for small samples. Second, this comparison involves a particularly small sample size. The full experimental results for the ELR test solving the Behrens-Fisher problem that are presented in Tables 2 to 5 show that the EL method is able to solve the BF problem, and that the ELR test has good power properties over a wide range of realistic situations.

17

6

COMPUTATIONAL ISSUES

The computing work associated with our EL approach in solving the BF problem is challenging. As Owen (2001) stated, it is computationally challenging to optimize a likelihood function of either parametric or empirical type over some nuisance parameters with other parameters held fixed at test values. The BF problem is well known to be very difficult to solve. This type of difficulty in optimizing empirical likelihood functions is especially clear. (i) The BF problem involves two nuisance parameters; (ii) Along the boundary of the parameter ρ2 , empirical likelihood function is ill-behaved according to our experiment. There are two possible reasons for this difficulty. One reason is historical. Any solution to the Behrens-Fisher problem must have some unpleasant properties (Zaman, 1996, p. 246). The ELR test, like other tests, must have some unfavorable features over certain areas in the parameter space. Our work proves this is true. The second possible reason comes from the design of the empirical likelihood approach. The nature of the EL method is that, in the neighborhood of the solution, the gradient matrix associated with the moment constraints will approach an ill-behaved state of being less than full rank (Mittelhammer et al., 2003). This occurs by design because the basic rationale of the EL method is to modify the sample weights such that the over-identified m empirical moment equations can be satisfied in order to solve for the unique solutions of the p unknowns, where m > p. This creates instability in gradient-based constrained optimization algorithms regarding the representation of the feasible spaces and feasible directions for such problems. Mittelhammer et al. (2003) used a ”concentrating-out” technique that utilized a nonlinear system procedure (NLSYS) and the Nelder-Mead (1965) method to achieve their computational purpose. These techniques worked well for their project. We have explored the possibility of using these techniques in solving the Behrens-Fisher problem. However, the Nelder-Mead method did not work well in our situation. The approach we used in this paper is so-called “direct solve” method. We directly solve the non-linear system of the moment equations and the first order conditions with respect to the parameters. The non-linear equation solving procedure, Eqsolve, in the Gauss package

18

is employed. The numerical solutions are acceptably good. In the process of the Monte Carlo simulations, there are a few samples drawn from the underlying distributions that can not solve the nonlinear system in computing the power of the ELR test. This reflects a typical example of the potential infeasibility that happens in implementing the EL method, in practise. When this happens, we reuse the sample that did not work for the Eqsolve procedure by altering the initial values of the parameters and then implementing the procedure again until the solution is found. In finding a new vector of initial values for the parameters, we apply the essential idea of the Differential Evolution method. Suppose θˆ0 is the vector of unsuccessfully estimated values of θ. The random search direction is formed of the difference between two random vectors θ1 − θ2 . The new initial value for the parameter vector is then in the form: θa = θˆ0 + s(θ1 − θ2 ), where s is a step size taking values of {0.4, 1, 2}. This new initial vector of parameter values is then feeded to the Eqsolve procedure to search for the global maximum. These techniques worked very well in keeping the samples that we have generated. By doing so we effectively avoid throwing away data sets casually and we have practically avoided the problem of selection bias. Therefore, we can claim that we have approximated the exact distribution for the ELR test statistic in finite samples. As a result, the empirical size distortion of the ELR test is effectively improved using this data reuse technique comparing to the ones when we did not use this data recycle technique. However, any gain comes with certain trade off. The computing time is lengthened and it is significantly longer than the EL approach in the testing for normality in Dong and Giles (2004). To give an indication of the extent of the difficulty it involved in solving the BF problem using the ELR method, we provide some examples of the computing times that it took in the Monte Carlo experiments. The computing time is longer when the parameter ρ2 is away from unity. It takes about 10 hours to compute the empirical size for the ELR test when ρ2 = 2, (n1 , n2 ) = (20, 10), and the number of replications is 5,000 using a Pentium 4, 2.4 GHZ PC. It is also very difficult in computation when the null hypothesis is not true, e.g. when the non-centrality

19

parameter δ 6= 0. For example, in computing the power of the ELR test, it takes the same machine around 27 hours for the case of (n1 , n2 ) = (20, 10), ρ2 = 0.5, δ = 1, the number of replications is 5,000. In conclusion, although it is difficult in computing the exact size and power of the ELR test for the BF problem, our EL approach using the Eqsolve algorithm works well for the BF problem. The empirical results are sound.

7

SUMMARY AND CONCLUSIONS

We have developed an new theoretical approach using the EL method to solve the Behrens-Fisher problem in this paper. The fact that the EL method is able to solve the BF problem is important. It shows the flexibility of the EL method in solving various problems in statistics and econometrics. A full range of Monte Carlo experiments are conducted to provide the analysis of the sampling properties of the ELR test. The actual sizes and the size-adjusted critical values in finite samples are simulated. The size-adjusted critical values are used to conduct the analysis of the power properties of the ELR test. The empirical results provide the evidence that the ELR test has good sampling properties across different parameter dimensions: the variance ratio parameter, the sample size pair, and the non-centrality parameter. Generally, the size-adjusted critical values that we have presented in Table 1 are ready to be used by researchers provided that the values of the parameter ρ2 and the sample size pair are conformable with the ones in our study. We have noted that the computing time in solving the BF problem is significantly long. The size of the ELR test statistic, in general, is still larger than the significance levels we consider. In the future study, it would be fruitful to explore some techniques that could reduce both of the computational difficulties and the size distortion and could still maintain the good power properties of the ELR test.

20

Acknowledgements This paper is based on one chapter of the author’s PhD dissertation, completed in the Department of Economics, University of Victoria, in December 2003. The author is very grateful to the thesis supervisor, professor David Giles, for his timely guidance and warm support. Special thanks are also extended to Don Ferguson, Ron Mittelhammer, Min Tsao, Graham Voss and Julie Zhou for their many helpful suggestions and contributions.

21

Appendix: Monte Carlo Results

Table 1: Size and Size-adjusted Critical Values of the ELR Test (n1 , n2 ) :

(20, 10)

(60, 30)

(100, 50)

(250, 125)

(20, 10)

ρ2 = 0.1 10% 0.1838 0.1868 5% 0.1214 0.1124 2% 0.0712 0.0636 1% 0.0488 0.0428 Size-adjusted Critical Values: 10% 4.3067 4.1787 5% 6.5523 6.1960 2% 10.7818 10.1350 1% 14.2266 13.2669

(60, 30)

0.1832 0.1084 0.0580 0.0404

0.1878 0.1066 0.0554 0.0342

0.2202 0.1482 0.0924 0.0686

0.2202 0.1468 0.0852 0.0584

0.2266 0.1448 0.0896 0.0626

0.1964 0.1218 0.0668 0.0424

4.0088 5.8025 9.0218 12.8602

3.9780 5.6186 8.2500 10.3353

5.1494 7.7853 12.4479 16.5913

4.8877 7.1960 10.4005 13.6268

4.9112 7.4334 11.1566 14.1973

4.3564 6.1392 9.3809 11.9212

ρ2 = 2

0.2536 0.1748 0.1068 0.0710

0.2332 0.1466 0.0828 0.0530

0.2926 0.2136 0.1404 0.1046

0.3014 0.2110 0.1422 0.1044

0.2780 0.1966 0.1298 0.0966

0.2448 0.1592 0.0888 0.0570

5.5645 7.5631 10.5112 13.6295

4.7971 6.7822 9.8197 12.0383

6.8276 10.8456 17.7402 22.9714

6.7529 9.3973 13.1147 15.9200

6.4580 9.3711 13.3604 16.7853

5.1142 7.1188 10.0495 12.1949

ρ2 = 10 10% 0.3564 0.3306 5% 0.2736 0.2574 2% 0.1892 0.1826 1% 0.1516 0.1428 Size-adjusted Critical Values: 10% 9.25930 8.4782 5% 15.0960 12.3244 2% 23.8243 18.3507 1% 31.6749 21.8045

(250, 125)

ρ2 = 0.5

ρ2 = 1 10% 0.2424 0.2538 5% 0.1698 0.1788 2% 0.1100 0.1126 1% 0.0814 0.0776 Size-adjusted Critical Values: 10% 5.7988 5.8435 5% 8.9674 7.9674 2% 14.3630 11.8099 1% 18.6082 14.3278

(100, 50)

0.2872 0.2066 0.1384 0.1044

0.2376 0.1576 0.0948 0.0636

6.7945 10.4183 16.2891 21.2274

5.2401 7.4935 11.1926 14.8243

Notes to table: The number of replications is 5,000. The sample sizes are the pair (n1 , n2 ). ρ2 = σ22 /σ12 .

22

Table 2: Power of the ELR Test for the Behrens-Fisher Problem (n1 , n2 ) ρ2 = 0.1 δ:

1

1.5

2

2.5

ρ2 = 0.5 3

4

1

1.5

2

2.5

3

4

0.1688 0.0994 0.0542 0.0371

0.2945 0.1816 0.1028 0.0743

0.4942 0.3101 0.1674 0.1226

0.7038 0.4562 0.2344 0.1562

0.8408 0.608 0.3078 0.2104

0.9508 0.8166 0.5022 0.3314

0.5748 0.4256 0.2710 0.1812

0.9148 0.8350 0.6926 0.5422

0.9832 0.9692 0.9376 0.8764

0.9906 0.9848 0.9798 0.9688

0.7362 0.5988 0.3978 0.2792

0.9768 0.9494 0.8850 0.8032

0.9962 0.9924 0.9840 0.9742

0.9450 0.9052 0.8040 0.7026

0.9992 0.9986 0.9962 0.9928

(20, 10) 10% 5% 2% 1%

0.3826 0.2520 0.1396 0.0906

0.6786 0.4926 0.2262 0.1226

0.9032 0.7764 0.4122 0.2122

0.9684 0.9028 0.5364 0.2836

0.9750 0.9388 0.5910 0.3298

0.9768 0.9546 0.6896 0.4316

(60, 30) 10% 5% 2% 1%

0.6086 0.4800 0.2934 0.2088

0.9314 0.8748 0.7292 0.5964

0.9926 0.9870 0.9724 0.9530

0.9924 0.9846 0.9764 0.9740

0.2870 0.1800 0.1034 0.0682 (100, 50)

10% 5% 2% 1%

0.6976 0.5992 0.4452 0.2972

0.9812 0.9608 0.9074 0.8092

0.9978 0.9970 0.9930 0.9848

0.3180 0.1956 0.0950 0.0584 (250, 125)

10% 5% 2% 1%

0.8254 0.7512 0.6308 0.5474

0.9996 0.9994 0.9972 0.9952

0.4216 0.3130 0.1828 0.1112

Notes to table: The number of replications is 5,000. The sample sizes are the pair (n1 , n2 ). ρ2 = σ22 /σ12 . δ is the non-centrality parameter.

23

Table 3: Power of the ELR Test for the Behrens-Fisher Problem (n1 , n2 ) ρ2 = 1 δ:

1

1.5

2

ρ2 = 2 2.5

3

4

1

1.5

2

2.5

3

4

0.1593 0.1044 0.0658 0.0440

0.1956 0.1303 0.0950 0.0800

0.2672 0.1750 0.1378 0.1240

0.4400 0.3016 0.2354 0.2188

0.6366 0.4612 0.3731 0.3522

0.8352 0.6564 0.5456 0.5272

0.289 0.1804 0.0995 0.0716

0.6084 0.4504 0.2876 0.2170

0.8770 0.7666 0.6190 0.5232

0.9772 0.9462 0.8794 0.8082

0.9968* 0.9942 0.9871 0.9683

0.4190 0.2693 0.1380 0.0808

0.8054 0.6634 0.4756 0.3459

0.9774 0.9472 0.8756 0.7996

0.9970 0.9950 0.9886 0.9776

0.7362 0.6116 0.4572 0.3624

0.9878 0.9746 0.9432 0.9084

0.9996 0.9988 0.9982 0.9976

(20, 10) 10% 5% 2% 1%

0.1564 0.0972 0.0606 0.0420

0.2170 0.1306 0.0898 0.0700

0.3456 0.2078 0.1342 0.1112

0.5116 0.3016 0.1872 0.1584

0.6886 0.4634 0.2806 0.2312

0.9050 0.7318 0.4624 0.3878

(60, 30) 10% 5% 2% 1%

0.2076 0.1440 0.0870 0.0664

0.4050 0.2878 0.1584 0.1144

0.7574 0.6310 0.4192 0.3082

0.9458 0.8962 0.7648 0.6614

0.9846 0.9748 0.9374 0.8930

0.9941 0.9914 0.9850 0.9798

0.1926 0.1342 0.0869 0.0648

(100, 50) 10% 5% 2% 1%

0.2364 0.1572 0.0980 0.0698

0.5486 0.4284 0.2770 0.1766

0.9086 0.8452 0.7322 0.5894

0.9918 0.9828 0.9612 0.9290

0.1988 0.1303 0.0814 0.0620 (250, 125)

10% 5% 2% 1%

0.3104 0.2082 0.1152 0.0776

0.8450 0.7514 0.6086 0.5048

0.9964 0.9914 0.9738 0.9572

0.2744 0.1844 0.1024 0.0682

Notes to table: The number of replications is 5,000. The sample sizes are the pair (n1 , n2 ). ρ2 = σ22 /σ12 . δ is the non-centrality parameter. Results with a star (*) are the ones that the number of replications is less than 5,000 due to computing difficulties.

24

Table 4: Power of the ELR tTest for the Behrens-Fisher Problem (n1 , n2 ) ρ2 = 10 δ:

1

1.5

2

2.5

3

4

0.2056 0.1512 0.1046 0.0764

0.2586 0.2084 0.1782 0.1572

0.4852 0.4182 0.3714 0.3422

0.7378 0.6630 0.6250 0.5928

0.8340 0.7815 0.7482 0.7240

0.8556 0.8088 0.7864 0.7712

0.1900 0.1386 0.0798 0.0600

0.2572 0.2060 0.1594 0.1270

0.7734 0.7144 0.6150 0.5270

0.9894 0.9762 0.8982 0.8072

0.9970 0.9856 0.9277 0.8543

0.9945* 0.9872 0.9287 0.8793

0.1522 0.0812 0.0516 0.0450

0.3674 0.2214 0.1446 0.1296

0.9114 0.8544 0.7945 0.7734

0.9995* 0.9987 0.9960 0.9922

0.9997* 0.9995 0.9973 0.9940

1.0000* 1.0000 1.0000 0.9968

0.2388 0.1276 0.0544 0.0300

0.6006 0.4794 0.3242 0.2248

0.9977* 0.9931 0.9851 0.9747

(20, 10) 10% 5% 2% 1% (60, 30) 10% 5% 2% 1% (100, 50) 10% 5% 2% 1% (250, 125) 10% 5% 2% 1%

Notes to table: The number of replications is 5,000. The sample sizes are the pair (n1 , n2 ). ρ2 = σ22 /σ12 . δ is the non-centrality parameter. Results with a star (*) are the ones that the number of replications is less than 5,000 due to computing difficulties.

25

Table 5: Size and Power Comparison of the Welch-Aspin and the ELR Tests for (n1 , n2 ) = (7, 7), α = .05 δ0

0

C0 0.1 0.2 0.3 0.4 0.5

Vwa 0.0501 0.0500 0.0500 0.0498 0.0498

1 EL 0.2846 0.2306 0.1864 0.1608 0.1510

Vwa 0.2301 0.2349 0.2385 0.2406 0.2413

2 EL 0.0522 0.0554 0.0640 0.0608 0.0762

Vwa 0.5628 0.5753 0.5855 0.5920 0.5942

3 EL 0.0494 0.0582 0.0764 0.0794 0.1048

Vwa 0.8521 0.8631 0.8722 0.8782 0.8803

4 EL 0.0594 0.0604 0.0886 0.1008 0.1380

Vwa 0.9729 0.9767 0.9799 0.9819 0.9826

EL 0.0660 0.0874 0.1164 0.1486 0.2000

Notes to table: The number of replications is 5,000. σ2

σ2

2 σ2 )−1 , C0 n2 2 σ σ2 µ2 )( n1 + n2 )−1 1 2

C0 = ( n1 )( n1 + 1

δ0 = (µ1 −

1

= {.1, .2, .3, .4, .5} corresponds to ρ2 = {9, 4, 2.33, 1.5, 1}. √ = δ 7. δ0 = {1, 2, 3, 4}, which is equivalent to δ = {.378, .7559, 1.1339, 1.5119}.

26

Table 6: Size and Power of the ELR Test for n1 = n2 = 7 ρ2

α

δ=0

δ = 0.3780

δ = 0.7559

δ = 1.1334

δ = 1.5119

1

10% 5% 2% 1%

0.2014 0.1510 0.1072 0.0850

0.1030 0.0522 0.0184 0.0082

0.0992 0.0494 0.0194 0.0098

0.1266 0.0594 0.0210 0.0082

0.1612 0.0660 0.0218 0.0070

1.5

10% 5% 2% 1%

0.2062 0.1608 0.1182 0.0966

0.1134 0.0554 0.0204 0.0106

0.1232 0.0582 0.0224 0.0112

0.1274 0.0604 0.0208 0.0088

0.1770 0.0874 0.0332 0.0102

2.33

10% 5% 2% 1%

0.2338 0.1864 0.1400 0.1110

0.1184 0.0640 0.0206 0.0094

0.1398 0.0764 0.0282 0.0158

0.1562 0.0886 0.0276 0.0106

0.2058 0.1164 0.0374 0.0150

4

10% 5% 2% 1%

0.2838 0.2306 0.1686 0.1366

0.1216 0.0608 0.0278 0.0110

0.1500 0.0794 0.0370 0.0168

0.1756 0.1008 0.0434 0.0172

0.2672 0.1486 0.0566 0.0244

9

10% 5% 2% 1%

0.3384 0.2846 0.2140 0.1760

0.1438 0.0762 0.0242 0.0126

0.1974 0.1048 0.0386 0.0174

0.2652 0.1380 0.0362 0.0180

0.3754 0.2000 0.0508 0.0240

Notes to table: The number of replications is 5,000. 2 C0 = {.1, √ .2, .3, .4, .5} corresponds to ρ = {9, 4, 2.33, 1.5, 1}. δ0 = δ 7. δ0 = {1, 2, 3, 4}, which is equivalent to δ = {.378, .7559, 1.1339, 1.5119}

27

References Aptech Systems, 2002. Gauss 5.0 for Windows NT, (Aptech Systems, Inc., Maple Valley WA). Aspin, A. A., 1948. “An Examination and Further Development of a formula Arising in the Problem of Comparing Two Mean Values”, Biometrika 35, 88-96. Bera, A., Bilias, Y., 2002. “The MM, ME, ML, EL, EF and GMM Approaches to Estimation: A Synthesis”, Journal of Econometrics 107, 51-86. Cochran, W. G., Cox, G. M., 1950. Experiment Designs (John Wiley and Sons, New York). Dong, L.B., Giles, D.E.A., 2004. “An empirical likelihood ratio test for normality”, Working paper EWP0401, Department of Economics, University of Victoria. Fisher, R. A., 1935. “The Fiducial Argument in Statistical Inference”, Annals of Eugenics 6, 391-198. Fisher, R. A., 1941.“The Asymptotic Approach to Behrens’ Integral with further Tables for the d Test of Significance”, Annals of Eugenics 11, 141-172. Giles, J. A., Giles, D. E. A., 1993. “Pre-Testing Estimation and Testing in Econometrics: Recent Developments”, Journal of Economic Surveys 7, 145-197. Jing, B. Y., 1995. “Two-Sample Empirical Likelihood Method”, Statistics and Probability Letters 24, 315-319. Lee, A. F. S., Gurland, J., 1975. “Size and Power of Tests for Equality of Means of Two Normal Populations”, Journal of the American Statistical Association 70, 933-944. Mittelhammer, R., Judge, G., Miller, D., 2000. Econometric Foundations (Cambridge Uni28

versity Press, Cambridge). Mittelhammer, R., G. Judge, and R. Schoenberg (2003), “Empirical Evidence Concerning the Finite Sample Performance of EL-Type Structural Equation Estimation and Inference Methods”, Working Paper, Washington State University, University of California, Berkeley, and Aptech Systems, Inc. Nelder, J. A., Mead, R., 1965. “A Simplex Method for Function Minimization”, Computer Journal 7, 308-313. Owen, A. B., 2001. Empirical Likelihood (Chapman & Hall/CRC, New York). Qin, J., 1991. “Likelihood and Empirical Likelihood Ratio Confidence Intervals in Two Sample Semi-parametric Models”, Technical Report Series University of Waterloo Stat-916. Qin, J., 1993. “Empirical Likelihood in Biased Sample Problems”, Annals of Statistics 21, 1182-1196. Tukey, J. W., 1954. “Unsolved Problems of Experimental Statistics”, Journal of the American Statistical Association 49, 706-731. Weerahandi, S., 1987. “Testing Regression Equality with Unequal Variances”, Econometrica 55, 1211-1215. Welch, B. L., 1947. “The Generalization of ‘Student’s’ Problem When Several Different Population Variances Are Involved”, Biometrika 34, 28-35. Zaman, A., 1996. Statistical Foundations for Econometric Techniques (Academic Press, New York).

29

View more...

Comments

Copyright � 2017 SILO Inc.