18.06.2014 Views

RW I:Discussion Papers - Rheinisch-Westfälisches Institut für ...

RW I:Discussion Papers - Rheinisch-Westfälisches Institut für ...

RW I:Discussion Papers - Rheinisch-Westfälisches Institut für ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Therefore, the contribution to the likelihood is<br />

P (∆w o i |F 1) =<br />

1<br />

√<br />

σ 2 w + σ 2 m<br />

φ<br />

(<br />

∆w o i − X iα<br />

√<br />

σ 2 w + σ 2 m<br />

)<br />

(23)<br />

given one measurement error, and<br />

P (∆w o i |F 2) =<br />

(<br />

)<br />

1 ∆wi o √ φ √ − X iα<br />

σ 2 w +2σm 2 σ 2 w +2σm<br />

2<br />

. (24)<br />

given two measurement errors.<br />

We have now all ingredients to build up the likelihood function to be maximized.<br />

Using indicator functions I(·) taking the value of unity if the condition in the brackets is<br />

satisfied and zero otherwise, we get<br />

L(∆wi o|Ω,X i)= P (N)P (M0)P (∆wi o ∈ C|i ∈ N0)I(∆wo i =0)+<br />

P (N)P (M0)P (∆wi o ∈ U|i ∈ N0)I(∆wo i > 0)+<br />

[P (N)P (M1)[P (∆wi o ∈ C|i ∈ N1) + P (∆wo i ∈ U|i ∈ N1)]+<br />

P (R)P (M0)[P (∆wi o ∈ C|i ∈ R0) + P (∆wo i ∈ U|i ∈ R0)]+<br />

P (R)P (M1)[P (∆wi o ∈ C|i ∈ R1) + P (∆wo i ∈ U|i ∈ R1)]+ (25)<br />

P (R)P (M2)[P (∆wi o ∈ C|i ∈ R2) + P (∆wo i ∈ U|i ∈ R2)]+<br />

P (F )P (M0)P (∆wi o |i ∈ F 0)+<br />

P (F )P (M1)P (∆wi o |i ∈ F 1)+<br />

P (F )P (M2)P (∆wi o|i ∈ F 2)] I(∆wo i ≠0),<br />

where Ω = (α, β rn ,β rf ,γ,σ 2 m,σ 2 w,σ 2 r,P m ) is the vector of parameters to be estimated.<br />

29

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!