Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Eviews5 / EViews5 / Docs / EViews 5 Users Guide.pdf
Скачиваний:
152
Добавлен:
23.03.2015
Размер:
11.51 Mб
Скачать

662—Chapter 21. Discrete and Limited Dependent Variable Models

Negative Binomial

If we maximize the negative binomial log likelihood, given above, for fixed η2 , we obtain the QMLE of the conditional mean parameters β . This QML estimator is consistent even if the conditional distribution of y is not negative binomial, provided that mi is correctly specified.

EViews sets η2 = 1 by default, which is a special case known as the geometric distribution. You may specify any other (positive) value by changing the number in the Fixed variance parameter field box. For the negative binomial QMLE, EViews by default reports the robust QMLE standard errors.

Views of Count Models

EViews provides a full complement of views of count models. You can examine the estimation output, compute frequencies for the dependent variable, view the covariance matrix, or perform coefficient tests. Additionally, you can select View/Actual, Fitted, Residual…

and pick from a number of views describing the ordinary residuals oi i− ( i ˆ ) , e = y m x , β

or you can examine the correlogram and histogram of these residuals. For the most part, all of these views are self-explanatory.

Note, however, that the LR test statistics presented in the summary statistics at the bottom of the equation output, or as computed under the View/Coefficient Tests/Redundant Variables - Likelihood Ratio… have a known asymptotic distribution only if the conditional distribution is correctly specified. Under the weaker GLM assumption that the true variance is proportional to the nominal variance, we can form a quasi-likelihood ratio,

QLR =

ˆ 2

ˆ

2

is the estimated proportional variance factor. This QLR sta-

LR σ

, where σ

 

tistic has an asymptotic χ2 distribution under the assumption that the mean is correctly specified and that the variances follow the GLM structure. EViews does not compute the

ˆ

2

based upon the stan-

QLR statistic, but it can be estimated by computing an estimate of σ

 

dardized residuals. We provide an example of the use of the QLR test statistic below.

If the GLM assumption does not hold, then there is no usable QLR test statistic with a known distribution; see Wooldridge (1996).

Procedures for Count Models

Most of the procedures are self-explanatory. Some details are required for the forecasting and residual creation procedures.

Forecast… provides you the option to forecast the dependent variable

y

i

or the pre-

dicted linear index x

ˆ

 

 

 

 

 

are

β . Note that for all of these models the forecasts of y

 

ˆ

 

i ˆ

ˆ

ˆ

 

 

 

i

 

given by yi =

m( xi, β)

where m(xi, β) =

exp ( xiβ) .

 

 

 

 

 

Demonstrations—663

Make Residual Series… provides the following three types of residuals for count models:

Ordinary

eoi

ˆ

 

= yim( xi, β)

Standardized (Pearson)

 

ˆ

 

esi

yim( xi, β)

 

= ------------------------------

 

 

ˆ ˆ

 

 

v( xi, β, γ)

Generalized

 

eg =(varies)

 

 

 

where the γ represents any additional parameters in the variance specification. Note that the specification of the variances may vary significantly between specifications.

For example, the Poisson model has v( x

ˆ

ˆ

, β) =

m( x , β) , while the exponential

ˆ

ˆ

2

i

 

i

has v( xi, β) =

m( xi, β)

 

.

 

 

The generalized residuals can be used to obtain the score vector by multiplying the generalized residuals by each variable in x . These scores can be used in a variety of LM or conditional moment tests for specification testing; see Wooldridge (1996).

Demonstrations

A Specification Test for Overdispersion

Consider the model:

NUMBi = β1 + β2IPi + β3FEBi + i ,

(21.48)

where the dependent variable NUMB is the number of strikes, IP is a measure of industrial production, and FEB is a February dummy variable, as reported in Kennan (1985, Table 1).

The results from Poisson estimation of this model are presented below:

664—Chapter 21. Discrete and Limited Dependent Variable Models

Dependent Variable: NUMB

Method: ML/QML - Poisson Count

Date: 09/14/97 Time: 10:58

Sample: 1 103

Included observations: 103

Convergence achieved after 4 iterations

Covariance matrix computed using second derivatives

Variable

Coefficient

Std. Error

z-Statistic

Prob.

 

 

 

 

 

C

1.725630

0.043656

39.52764

0.0000

IP

2.775334

0.819104

3.388254

0.0007

FEB

-0.377407

0.174520

-2.162539

0.0306

 

 

 

 

R-squared

0.064502

Mean dependent var

5.495146

Adjusted R-squared

0.045792

S.D. dependent var

3.653829

S.E. of regression

3.569190

Akaike info criterion

5.583421

Sum squared resid

1273.912

Schwarz criterion

5.660160

Log likelihood

-284.5462

Hannan-Quinn criter.

5.614503

Restr. log likelihood

-292.9694

Avg. log likelihood

-2.762584

LR statistic (2 df)

16.84645

LR index (Pseudo-R2)

0.028751

Probability(LR stat)

0.000220

 

 

 

 

 

 

 

 

Cameron and Trivedi (1990) propose a regression based test of the Poisson restriction

v( xi, β) = m( xi, β) . To carry out the test, first estimate the Poisson model and obtain the fitted values of the dependent variable. Click Forecast and provide a name for the forecasted dependent variable, say NUMB_F. The test is based on an auxiliary regression of

2

 

ˆ 2

 

 

 

 

 

eoi

yi on yi and testing the significance of the regression coefficient. For this example,

the test regression can be estimated by the command:

 

 

 

 

ls (numb-numb_f)^2-numb numb_f^2

 

 

 

yielding the following results:

 

 

 

 

 

 

 

Dependent Variable: (NUMB-NUMB_F)^2-NUMB

 

 

 

 

 

Method: Least Squares

 

 

 

 

 

 

 

Date: 09/14/97 Time: 11:05

 

 

 

 

 

 

Sample: 1 103

 

 

 

 

 

 

 

Included observations: 103

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Variable

Coefficient

Std. Error

t-Statistic

Prob.

 

 

 

 

 

 

 

 

 

 

NUMB_F^2

0.238874

0.052115

4.583571

0.0000

 

 

 

 

 

 

 

 

 

 

R-squared

0.043930

Mean dependent var

6.872929

 

 

 

Adjusted R-squared

0.043930

S.D. dependent var

17.65726

 

 

 

S.E. of regression

17.26506

Akaike info criterion

8.544908

 

 

 

Sum squared resid

30404.41

Schwarz criterion

8.570488

 

 

 

Log likelihood

-439.0628

Durbin-Watson stat

1.711805

 

 

 

 

 

 

 

 

 

The t-statistic of the coefficient is highly significant, leading us to reject the Poisson restriction. Moreover, the estimated coefficient is significantly positive, indicating overdispersion in the residuals.

Demonstrations—665

An alternative approach, suggested by Wooldridge (1996), is to regress esi − 1

ˆ

, on yi . To

perform this test, select Proc/Make Residual Series… and select Standardized. Save the results in a series, say SRESID. Then estimating the regression specification:

sresid^2-1 numbf yields the results:

Dependent Variable: SRESID^2-1

Method: Least Squares

Date: 10/06/97 Time: 16:05

Sample: 1 103

Included observations: 103

Variable

Coefficient

Std. Error

t-Statistic

Prob.

 

 

 

 

 

NUMBF

0.221238

0.055002

4.022326

0.0001

 

 

 

 

R-squared

0.017556

Mean dependent var

1.161573

Adjusted R-squared

0.017556

S.D. dependent var

3.138974

S.E. of regression

3.111299

Akaike info criterion

5.117619

Sum squared resid

987.3786

Schwarz criterion

5.143199

Log likelihood

-262.5574

Durbin-Watson stat

1.764537

 

 

 

 

 

Both tests suggest the presence of overdispersion, with the variance approximated by v = m( 1 + 0.23m) .

Given the evidence of overdispersion and the rejection of the Poisson restriction, we will re-estimate the model, allowing for mean-variance inequality. Our approach will be to estimate the two-step negative binomial QMLE specification (termed the quasi-generalized pseudo-maximum likelihood estimator by Gourieroux, Monfort, and Trognon (1984a, b))

using the estimate of ˆ 2 derived above. To compute this estimator, simply select Negative

η

Binomial (QML) and enter 0.22124 in the edit field for Fixed variance parameter.

We will use the GLM variance calculations, so you should click on Option in the Equation Specification dialog and mark the Robust Covariance and GLM options. The estimation results are shown below:

666—Chapter 21. Discrete and Limited Dependent Variable Models

Dependent Variable: NUMB

Method: QML - Negative Binomial Count

Date: 10/11/97 Time: 23:53

Sample: 1 103

Included observations: 103

QML parameter used in estimation: 0.22124

Convergence achieved after 3 iterations

GLM Robust Standard Errors & Covariance

Variance factor estimate = 2.465660162

Covariance matrix computed using second derivatives

Variable

Coefficient

Std. Error

z-Statistic

Prob.

 

 

 

 

 

C

1.724906

0.102543

16.82135

0.0000

IP

2.833103

1.919447

1.475999

0.1399

FEB

-0.369558

0.377376

-0.979285

0.3274

 

 

 

 

R-squared

0.064374

Mean dependent var

5.495146

Adjusted R-squared

0.045661

S.D. dependent var

3.653829

S.E. of regression

3.569435

Akaike info criterion

5.174385

Sum squared resid

1274.087

Schwarz criterion

5.251125

Log likelihood

-263.4808

Hannan-Quinn criter.

5.205468

Restr. log likelihood

-522.9973

Avg. log likelihood

-2.558066

LR statistic (2 df)

519.0330

LR index (Pseudo-R2)

0.496210

Probability(LR stat)

0.000000

 

 

 

 

 

 

 

 

The header indicates that the estimated GLM variance factor is 2.4, suggesting that the negative binomial ML would not have been an appropriate specification. Nevertheless, the negative binomial QML should be consistent, and under the GLM assumption, the standard errors should be consistently estimated. It is worth noting that the coefficients on IP and FEB, which were strongly statistically significant in the Poisson specification, are no longer significantly different from zero at conventional significance levels.

Quasi-likelihood Ratio Statistic

As described by Wooldridge (1996), specification testing using likelihood ratio statistics requires some care when based upon QML models. We illustrate here the differences between a standard LR test for significant coefficients and the corresponding QLR statistic.

From the results above, we know that the overall likelihood ratio statistic for the Poisson model is 16.85, with a corresponding p-value of 0.0002. This statistic is valid under the assumption that m(xi, β) is specified correctly and that the mean-variance equality holds.

We can decisively reject the latter hypothesis, suggesting that we should derive the QML estimator with consistently estimated covariance matrix under the GLM variance assumption. While EViews does not automatically adjust the LR statistic to reflect the QML assumption, it is easy enough to compute the adjustment by hand. Following Wooldridge, we construct the QLR statistic by dividing the original LR statistic by the estimated GLM variance factor.

Соседние файлы в папке Docs