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

12.1 Measurement Error (

N/A
N/A
Protected

Academic year: 2021

シェア "12.1 Measurement Error ("

Copied!
17
0
0

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

全文

(1)

7. Some Formulas:

Let Xnand Ynbe the random variables which satisfy plim Xn =c and plim Yn= d. Then,

(a) plim (Xn+Yn)=c+d (b) plim XnYn =cd

(c) plim Xn/Yn =c/d for d ,0

(d) plim g(Xn)=g(c) for a function g(·)

=⇒ Slutsky’s Theorem (スルツキー定理)

(2)

8. Central Limit Theorem (中心極限定理)

Univariate Case: X1, X2,· · ·, Xn are mutually independently and identically distributed as Xi ∼(µ, σ2).

Then,

XE(X)

V(X)

= X−µ σ/√

n −→ N(0,1), which implies

n(X−µ)= 1

n

n i=1

(Xi−µ) −→ N(0, σ2).

(3)

Multivariate Case: X1, X2,· · ·, Xnare mutually independently and identically distributed as Xi ∼(µ, Σ).

Then,

√1 n

n i=1

(Xi−µ) −→ N(0,Σ) 9. Central Limit Theorem (Generalization)

X1, X2, · · ·, Xn are mutually independently and identically distributed as Xi ∼ (µ, Σi).

Then,

√1 n

n i=1

(Xi−µ) −→ N(0,Σ), where

Σ = lim

n→∞



1 n

n i=1

Σi



.

(4)

10. Definition: Let ˆθnbe a consistent estimator ofθ. Suppose that √

n(ˆθn−θ) converges to N(0,Σ) in distribution.

Then, we say that ˆθnhas an asymptotic distribution (漸近分布): N(θ,Σ/n).

11. X1,X2,· · ·,Xn are random variables with density function f (x;θ).

Let ˆθnbe a maximum likelihood estimator ofθ.

Then, under some regularity conditions. ˆθn is a consistent estimator ofθand the asymptotic distribution of √

n(ˆθ−θ) is given by: N

0,lim (I(θ)

n )1

. 12. Regularity Conditions:

(a) The domain of Xi does not depend onθ.

(b) There exists at least third-order derivative of f (x;θ) with respect toθ, and their derivatives are finite.

(5)

13. Thus, MLE is

(i) consistent

(ii) asymptotically normaland (iii) asymptotically efficient.

(6)

11 Consistency and Asymptotic Normality of OLSE

Regression model: y= Xβ+u, u∼ (0, σ2In).

Consistency:

1. Let ˆβn = (X0X)1X0y be the OLS with sample size n.

Consistency: As n is large, ˆβnconverges toβ. 2. Assume the stationarity assumption for X, i.e.,

1

nX0X −→ Mxx. Then, we have the following result:

1

nX0u −→ 0.

(7)

Proof:

According to Chebyshev’s inequality, for g(Z)≥0, P(g(Z)k)E(g(Z))

k , where k is a positive constant.

Set g(Z)=Z0Z, and Z = 1 nX0u.

Apply Chebyshev’s inequality.

E( (1

nX0u)01 nX0u)

= 1 n2E(

u0XX0u)

= 1 n2E(

tr(u0XX0u))

= 1 n2E(

tr(XX0uu0))

= 1 n2tr(

XX0E(uu0))

= σ2

n2tr(XX0)= σ2

n2tr(X0X)= σ2 n tr(1

nX0X). Therefore,

P( (1

nX0u)01

nX0uk)

≤ σ2 nktr(1

nX0X)−→0×tr(Mxx)=0.

(8)

Note that from the assumption, 1

nX0X −→ Mxx. Therefore, we have:

(1 nX0u)01

nX0u−→0, which implies:

1

nX0u−→0, because (1

nX0u)01

nX0u indicates a quadratic form.

3. Note that 1

nX0X −→ Mxx results in (1

nX0X)1 −→ Mxx1.

=⇒Slutsky’s Theorem

(*) Slutsky’s Theorem g(ˆθ)−→g(θ), when ˆθ−→θ.

(9)

4. OLS is given by:

βˆn=β+(X0X)1X0u= β+(1

nX0X)1(1 nX0u). Therefore,

βˆn−→β+Mxx1×0=β Thus, OLSE is a consitent estimator.

Asymptotic Normality:

1. Asymptotic Normality of OLSE

n( ˆβn−β) −→ N(02M−1xx), when n −→ ∞.

(10)

2. Central Limit Theorem: Greenberg and Webster (1983)

Z1, Z2,· · ·, Znare mutually indelendently distributed with meanµand variance Σi.

Then, we have the following result:

√1 n

n i=1

(Zi−µ) −→ N(0,Σ), where

Σ = lim

n→∞



1 n

n i=1

Σi



. The distribution of Ziis not assumed.

3. Define Zi = xi0ui. Then,Σi =Var(Zi)=σ2x0ixi.

(11)

4. Σis defined as:

Σ = lim

n→∞



1 n

n i=1

σ2x0ixi



= σ2lim

n→∞

(1 nX0X

)

2Mxx,

where

X =









x1 x2

...

xn









5. Applying Central Limit Theorem (Greenberg and Webster (1983), we obtain the following:

√1 n

n i=1

x0iui = 1

nX0u−→ N(0, σ2Mxx). On the other hand, from ˆβn =β+(X0X)1X0u, we can rewrite as:

n( ˆβ−β)=(1

nX0X)1 1

nX0u.

(12)

Var((1

nX0X)−1 1

nX0u )

=E((1

nX0X)−1 1

nX0u((1

nX0X)−1 1

nX0u)0)

=(1

nX0X)−1(1

nX0E(uu0)X)(1

nX0X)−1

2(1

nX0X)1

−→ σ2Mxx1.

Therefore,

n( ˆβ−β) −→ N(0, σ2Mxx1)

=⇒Asymptotic normality (漸近的正規性) of OLSE The distribution of uiis not assumed.

(13)

12 Instrumental Variable ( 操作変数法 )

12.1 Measurement Error (

測定誤差

)

Errors in Variables

1. True regression model:

y= X˜β+u 2. Observed variable:

X = X˜ +V

V: is called the measurement error (測定誤差or観測誤差).

3. For the elements which do not include measurement errors in X, the corre- sponding elements in V are zeros.

(14)

4. Regression using observed variable:

y= Xβ+(uVβ) OLS ofβis:

βˆ =(X0X)1X0y=β+(X0X)1X0(uVβ)

5. Assumptions:

(a) The measurement error in X is uncorrelated with ˜X in the limit. i.e., plim(1

nX˜0V)

=0. Therefore, we obtain the following:

plim(1 nX0X)

=plim(1 nX˜0X˜)

+plim(1 nV0V)

= Σ + Ω

(15)

(b) u is not correlated with V.

u is not correlated with ˜X.

That is,

plim(1 nV0u)

=0, plim(1 nX˜0u)

=0.

6. OLSE ofβis:

βˆ =β+(X0X)1X0(uVβ)=β+(X0X)1( ˜X+V)0(uVβ). Therefore, we obtain the following:

plim ˆβ=β−(Σ + Ω)1Ωβ

(16)

7. Example: The Case of Two Variables:

The regression model is given by:

yt =α+β˜xt+ut, xt = ˜xt+vt. Under the above model,

Σ =plim(1 nX˜0X˜)

= plim



1 1

n

˜xi

1 n

˜xi 1 n

˜x2i



=

(1 µ

µ µ22 )

, whereµandσ2represent the mean and variance of ˜xi.

Ω =plim(1 nV0V)

=plim

(0 0

0 1

n

v2i

)

=

(0 0

0 σ2v

) .

Therefore, plim

(αˆ βˆ )

= (α

β )

((1 µ

µ µ22 )

+

(0 0

0 σ2v

))1(0 0 0 σ2v

) (α β )

= (α

β )

− 1 σ22v

(−µσ2vβ σ2vβ

)

(17)

Now we focus onβ.

βˆ is not consistent. because of:

plim( ˆβ)=β− σ2vβ σ22v

= β

1+σ2v2 < β

参照

関連したドキュメント

In this paper we develop an elementary formula for a family of elements {˜ x[a]} a∈ Z n of the upper cluster algebra for any fixed initial seed Σ.. This family of elements

R.Brown and J-L.Loday [5] noted that if the second dimension G 2 of a simplicial group G, is generated by the degenerate elements, that is, elements coming from lower dimensions,

Analogs of this theorem were proved by Roitberg for nonregular elliptic boundary- value problems and for general elliptic systems of differential equations, the mod- ified scale of

Then it follows immediately from a suitable version of “Hensel’s Lemma” [cf., e.g., the argument of [4], Lemma 2.1] that S may be obtained, as the notation suggests, as the m A

In the case of the Ariki–Koike algebra, that is, the Hecke algebra of the complex reflection group G(l, 1, n), they are Laurent polynomials whose factors determine when Specht

By iterating this procedure, we produce cellular bases for B–M–W algebras on which a large Abelian subalgebra, generated by elements which generalise the Jucys–Murphy elements from

The advection-diffusion equation approximation to the dispersion in the pipe has generated a considera- bly more ill-posed inverse problem than the corre- sponding

In addition, the purpose of this paper is to demonstrate the proposed models and methods with various scenarios for real data analysis for comparing asymmetric distributions for