• 検索結果がありません。

Econometrics I: Solutions of Homework #11

N/A
N/A
Protected

Academic year: 2021

シェア "Econometrics I: Solutions of Homework #11"

Copied!
5
0
0

読み込み中.... (全文を見る)

全文

(1)

Econometrics I: Solutions of Homework #11

Hiroki Kato * July 3, 2020

Contents

1 Solutions 1

1.1 (1) . . . . 1

1.2 (2) . . . . 2

1.3 (3) . . . . 3

1.4 (4) . . . . 4

1 Solutions

1.1 (1)

Since we assume u

t

is normally distributed with E(u

t

) = 0 and V (u

t

) = σ

2

, the density function of u

t

is given by

f (u

t

) = 1

(2πσ

2

)

(1/2)

exp (

1 2σ

2

u

2t

) .

By the mutual independent assumption, the joint density function of u

1

, . . . , u

T

is given by

f(u

1

, . . . , u

T

) = 1

(2πσ

2

)

(T/2)

exp (

1 2σ

2

T t=1

u

2t

)

*

e­mail: [email protected]. Room 503. If you find any errors in handouts and materials,

please contact me via e­mail.

(2)

Transforming u

t

, we have the following likelihood function;

1

L(α, β, σ

2

| y

1

, . . . , y

T

) 1

(2πσ

2

)

(T/2)

exp (

1 2σ

2

T t=1

(y

t

α βx

t

)

2

)

. (1)

1.2 (2)

To obtain the log­likelihood function, we take a natural logarithm of (1) as follows:

log L(α, β, σ

2

| y

1

, . . . , y

T

) = T

2 log(2π) T

2 log(σ

2

) 1 2σ

2

T t=1

(y

t

α βx

t

)

2

Taking a first­order derivative with respect to unknown parameters (α, β, σ

2

) yields

∂L(α, β, σ

2

| y

1

, . . . , y

T

)

∂α = 1

σ

2

T t=1

(y

t

α βx

t

) = 0 (2)

∂L(α, β, σ

2

| y

1

, . . . , y

T

)

∂β = 1

σ

2

T t=1

(y

t

α βx

t

)x

t

= 0 (3)

∂L(α, β, σ

2

| y

1

, . . . , y

T

)

∂σ

2

= T

2 1 σ

2

+ 1

4

T t=1

(y

t

α βx

t

)

2

= 0 (4)

Rewriting the equation (2) and the equation (3) yields

T y ¯ αT ˜ βT ˜ x ¯ = 0,

t

y

t

x

t

αT ˜ x ¯ β ˜ ∑

t

x

2t

= 0,

where y ¯ = ∑

t

y

t

/T and x ¯ = ∑

t

x

t

/T .

First, using the first equation, we can obtain the ML estimator of α, denoted by α, as follows; ˜

˜

α = ¯ y β ˜ x. ¯

Substituting this estimator into the second equation, we obtain the ML estimator of β, denoted by

(3)

β, as follows; ˜

β[ ˜ ∑

t

x

2t

T x ¯

2

] + ∑

t

y

t

x

t

T ¯ y = 0 β ˜ =

t

y

t

x

t

T x ¯ y ¯

t

x

2t

T x ¯

2

β ˜ =

t

(y

t

y)(x ¯

t

x) ¯

t

(x

t

x) ¯

2

Solving the equation (4) gives the ML estimator of σ

2

, denoted by σ ˜

2

, as follows;

˜ σ

2

=

T

t=1

(y

t

α ˜ βx ˜

t

)

2

T .

In summary, the ML estimator of θ is

θ ˜ = ( ˜ α, β, ˜ σ ˜

2

)

= (

¯ y β ˜ x, ¯

t

(y

t

y)(x ¯

t

x) ¯

t

(x

t

x) ¯

2

,

T

t=1

(y

t

α ˜ βx ˜

t

)

2

T

)

.

1.3 (3)

By the equation (2), (3), and (4),

2

L(α, β, σ

2

| y

1

, . . . , y

T

)

∂α

2

= T

σ

2

2

L(α, β, σ

2

| y

1

, . . . , y

T

)

∂β

2

=

t

x

2t

σ

2

2

L(α, β, σ

2

| y

1

, . . . , y

T

)

∂(σ

2

)

2

= T

4

1

σ

6

T t=1

(y

t

α βx

t

)

2

= T

4

1

σ

6

T t=1

u

2t

2

L(α, β, σ

2

| y

1

, . . . , y

T

)

∂α∂β =

t

x

t

σ

2

2

L(α, β, σ

2

| y

1

, . . . , y

T

)

∂α∂σ

2

= 1

σ

4

T t=1

(y

t

α βx

t

) = 1 σ

4

T t=1

u

t

2

L(α, β, σ

2

| y

1

, . . . , y

T

)

∂β∂σ

2

= 1

σ

4

T t=1

(y

t

α βx

t

)x

t

= 1 σ

4

T t=1

u

t

x

t

(4)

By the distributional assumption, we have ∑

t

E(u

t

) = 0, ∑

t

E(u

t

)x

t

= 0, and ∑

t

E(u

2t

) =

t

[E(u

2t

) E(u

t

)] = ∑

t

V (u

t

) = T σ

2

. Thus, the information matrix I(θ) is given by

I(θ) = E

(

2

L(α, β, σ

2

| y

1

, . . . , y

T

)

∂θ∂θ

)

=

 

 

 

T σ2

txt

σ2

0

txt

σ2

tx2t

σ2

0

0 0

T4

 

 

 

1.4 (4)

Step 1: Variance­Covariance Matrix of MLE

By the Cramer­Rao lower bound theorem, the inverse of information matrix I (θ)

1

provides a lower bound of the variance­covariance matrix for unbiased estimators of θ.

2

Then, the variance­

covariance matrix of θ ˜ is

 

 

 

σ2 T

tx2t

t(xt−x)¯2

σ

2x¯

t(xt−x)¯ 2

0

σ

2 x¯

t(xt−x)¯ 2

σ2

t(xt¯x)2

0

0 0

T4

 

 

 

Step 2: Derive the expectation and variance of OLSE

The OLS estimator of β, denoted by β, is equivalent to the ML estimator of ˆ β, denoted by β. That ˜ is,

β ˆ =

t

(y

t

y)(x ¯

t

x) ¯

t

(x

t

x) ¯

2

. Thus, we have E ( ˆ β) = β and V ( ˆ β) = σ

2

/( ∑

t

(x

t

x) ¯

2

). This implies that there is no difference between the MLE and the OLSE of β.

The OLS estimator of σ

2

is given by

ˆ

σ

2

= 1 T 2

T

(y

t

α ˜ βx ˜

t

)

2

.

(5)

This estimator is unbiased: E(ˆ σ

2

) = σ

2

. Since the OLSE is different from the ML estimator of β, β. the ML estimator is biased estimator of ˜ σ

2

, that is,

E(˜ σ

2

) = T 2

T E( ˆ σ

2

) = T 2

T σ

2

.

参照

関連したドキュメント

The torsion free generalized connection is determined and its coefficients are obtained under condition that the metric structure is parallel or recurrent.. The Einstein-Yang

It is suggested by our method that most of the quadratic algebras for all St¨ ackel equivalence classes of 3D second order quantum superintegrable systems on conformally flat

[11] Karsai J., On the asymptotic behaviour of solution of second order linear differential equations with small damping, Acta Math. 61

Viscous profiles for traveling waves of scalar balance laws: The uniformly hyperbolic case ∗..

Reynolds, “Sharp conditions for boundedness in linear discrete Volterra equations,” Journal of Difference Equations and Applications, vol.. Kolmanovskii, “Asymptotic properties of

7, Fan subequation method 8, projective Riccati equation method 9, differential transform method 10, direct algebraic method 11, first integral method 12, Hirota’s bilinear method

Kilbas; Conditions of the existence of a classical solution of a Cauchy type problem for the diffusion equation with the Riemann-Liouville partial derivative, Differential Equations,

Answering a question of de la Harpe and Bridson in the Kourovka Notebook, we build the explicit embeddings of the additive group of rational numbers Q in a finitely generated group