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 ...
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