\(\newcommand{\footnotename}{footnote}\)
\(\def \LWRfootnote {1}\)
\(\newcommand {\footnote }[2][\LWRfootnote ]{{}^{\mathrm {#1}}}\)
\(\newcommand {\footnotemark }[1][\LWRfootnote ]{{}^{\mathrm {#1}}}\)
\(\newcommand \ensuremath [1]{#1}\)
\(\newcommand {\LWRframebox }[2][]{\fbox {#2}} \newcommand {\framebox }[1][]{\LWRframebox } \)
\(\newcommand {\setlength }[2]{}\)
\(\newcommand {\addtolength }[2]{}\)
\(\newcommand {\setcounter }[2]{}\)
\(\newcommand {\addtocounter }[2]{}\)
\(\newcommand {\cline }[1]{}\)
\(\newcommand {\directlua }[1]{\text {(directlua)}}\)
\(\newcommand {\luatexdirectlua }[1]{\text {(directlua)}}\)
\(\newcommand {\protect }{}\)
\(\def \LWRabsorbnumber #1 {}\)
\(\def \LWRabsorbquotenumber "#1 {}\)
\(\def \mathchar {\ifnextchar "\LWRabsorbquotenumber \LWRabsorbnumber }\)
\(\def \mathcode #1={\mathchar }\)
\(\let \delcode \mathcode \)
\(\let \delimiter \mathchar \)
\(\let \LWRref \ref \)
\(\renewcommand {\ref }{\ifstar \LWRref \LWRref }\)
\(\newcommand {\intertext }[1]{\text {#1}\notag \\}\)
\(\newcommand {\mathllap }[2][]{{#1#2}}\)
\(\newcommand {\mathrlap }[2][]{{#1#2}}\)
\(\newcommand {\mathclap }[2][]{{#1#2}}\)
\(\newcommand {\mathmbox }[1]{#1}\)
\(\newcommand {\clap }[1]{#1}\)
\(\newcommand {\LWRmathmakebox }[2][]{#2}\)
\(\newcommand {\mathmakebox }[1][]{\LWRmathmakebox }\)
\(\newcommand {\cramped }[2][]{{#1#2}}\)
\(\newcommand {\crampedllap }[2][]{{#1#2}}\)
\(\newcommand {\crampedrlap }[2][]{{#1#2}}\)
\(\newcommand {\crampedclap }[2][]{{#1#2}}\)
\(\newenvironment {crampedsubarray}[1]{}{}\)
\(\newcommand {\crampedsubstack }{}\)
\(\newcommand {\smashoperator }[2][]{#2\limits }\)
\(\newcommand {\adjustlimits }{}\)
\(\newcommand {\SwapAboveDisplaySkip }{}\)
\(\require {extpfeil}\)
\(\Newextarrow \xleftrightarrow {10,10}{0x2194}\)
\(\Newextarrow \xLeftarrow {10,10}{0x21d0}\)
\(\Newextarrow \xhookleftarrow {10,10}{0x21a9}\)
\(\Newextarrow \xmapsto {10,10}{0x21a6}\)
\(\Newextarrow \xRightarrow {10,10}{0x21d2}\)
\(\Newextarrow \xLeftrightarrow {10,10}{0x21d4}\)
\(\Newextarrow \xhookrightarrow {10,10}{0x21aa}\)
\(\Newextarrow \xrightharpoondown {10,10}{0x21c1}\)
\(\Newextarrow \xleftharpoondown {10,10}{0x21bd}\)
\(\Newextarrow \xrightleftharpoons {10,10}{0x21cc}\)
\(\Newextarrow \xrightharpoonup {10,10}{0x21c0}\)
\(\Newextarrow \xleftharpoonup {10,10}{0x21bc}\)
\(\Newextarrow \xleftrightharpoons {10,10}{0x21cb}\)
\(\newcommand {\LWRdounderbracket }[3]{\underset {#3}{\underline {#1}}}\)
\(\newcommand {\LWRunderbracket }[2][]{\LWRdounderbracket {#2}}\)
\(\newcommand {\underbracket }[1][]{\LWRunderbracket }\)
\(\newcommand {\LWRdooverbracket }[3]{\overset {#3}{\overline {#1}}}\)
\(\newcommand {\LWRoverbracket }[2][]{\LWRdooverbracket {#2}}\)
\(\newcommand {\overbracket }[1][]{\LWRoverbracket }\)
\(\newcommand {\LATEXunderbrace }[1]{\underbrace {#1}}\)
\(\newcommand {\LATEXoverbrace }[1]{\overbrace {#1}}\)
\(\newenvironment {matrix*}[1][]{\begin {matrix}}{\end {matrix}}\)
\(\newenvironment {pmatrix*}[1][]{\begin {pmatrix}}{\end {pmatrix}}\)
\(\newenvironment {bmatrix*}[1][]{\begin {bmatrix}}{\end {bmatrix}}\)
\(\newenvironment {Bmatrix*}[1][]{\begin {Bmatrix}}{\end {Bmatrix}}\)
\(\newenvironment {vmatrix*}[1][]{\begin {vmatrix}}{\end {vmatrix}}\)
\(\newenvironment {Vmatrix*}[1][]{\begin {Vmatrix}}{\end {Vmatrix}}\)
\(\newenvironment {smallmatrix*}[1][]{\begin {matrix}}{\end {matrix}}\)
\(\newenvironment {psmallmatrix*}[1][]{\begin {pmatrix}}{\end {pmatrix}}\)
\(\newenvironment {bsmallmatrix*}[1][]{\begin {bmatrix}}{\end {bmatrix}}\)
\(\newenvironment {Bsmallmatrix*}[1][]{\begin {Bmatrix}}{\end {Bmatrix}}\)
\(\newenvironment {vsmallmatrix*}[1][]{\begin {vmatrix}}{\end {vmatrix}}\)
\(\newenvironment {Vsmallmatrix*}[1][]{\begin {Vmatrix}}{\end {Vmatrix}}\)
\(\newenvironment {psmallmatrix}[1][]{\begin {pmatrix}}{\end {pmatrix}}\)
\(\newenvironment {bsmallmatrix}[1][]{\begin {bmatrix}}{\end {bmatrix}}\)
\(\newenvironment {Bsmallmatrix}[1][]{\begin {Bmatrix}}{\end {Bmatrix}}\)
\(\newenvironment {vsmallmatrix}[1][]{\begin {vmatrix}}{\end {vmatrix}}\)
\(\newenvironment {Vsmallmatrix}[1][]{\begin {Vmatrix}}{\end {Vmatrix}}\)
\(\newcommand {\LWRmultlined }[1][]{\begin {multline*}}\)
\(\newenvironment {multlined}[1][]{\LWRmultlined }{\end {multline*}}\)
\(\let \LWRorigshoveleft \shoveleft \)
\(\renewcommand {\shoveleft }[1][]{\LWRorigshoveleft }\)
\(\let \LWRorigshoveright \shoveright \)
\(\renewcommand {\shoveright }[1][]{\LWRorigshoveright }\)
\(\newenvironment {dcases}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {dcases*}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {rcases}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {rcases*}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {drcases}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {drcases*}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {cases*}{\begin {cases}}{\end {cases}}\)
\(\newcommand {\MoveEqLeft }[1][]{}\)
\(\def \LWRAboxed #1!|!{\fbox {\(#1\)}&\fbox {\(#2\)}} \newcommand {\Aboxed }[1]{\LWRAboxed #1&&!|!} \)
\( \newcommand {\LWRABLines }[1][\Updownarrow ]{#1 \notag \\}\newcommand {\ArrowBetweenLines }{\ifstar \LWRABLines \LWRABLines } \)
\(\newcommand {\shortintertext }[1]{\text {#1}\notag \\}\)
\(\newcommand {\vdotswithin }[1]{\hspace {.5em}\vdots }\)
\(\newcommand {\LWRshortvdotswithinstar }[1]{\vdots \hspace {.5em} & \\}\)
\(\newcommand {\LWRshortvdotswithinnostar }[1]{& \hspace {.5em}\vdots \\}\)
\(\newcommand {\shortvdotswithin }{\ifstar \LWRshortvdotswithinstar \LWRshortvdotswithinnostar }\)
\(\newcommand {\MTFlushSpaceAbove }{}\)
\(\newcommand {\MTFlushSpaceBelow }{\\}\)
\(\newcommand \lparen {(}\)
\(\newcommand \rparen {)}\)
\(\newcommand {\ordinarycolon }{:}\)
\(\newcommand {\vcentcolon }{\mathrel {\mathop \ordinarycolon }}\)
\(\newcommand \dblcolon {\vcentcolon \vcentcolon }\)
\(\newcommand \coloneqq {\vcentcolon =}\)
\(\newcommand \Coloneqq {\dblcolon =}\)
\(\newcommand \coloneq {\vcentcolon {-}}\)
\(\newcommand \Coloneq {\dblcolon {-}}\)
\(\newcommand \eqqcolon {=\vcentcolon }\)
\(\newcommand \Eqqcolon {=\dblcolon }\)
\(\newcommand \eqcolon {\mathrel {-}\vcentcolon }\)
\(\newcommand \Eqcolon {\mathrel {-}\dblcolon }\)
\(\newcommand \colonapprox {\vcentcolon \approx }\)
\(\newcommand \Colonapprox {\dblcolon \approx }\)
\(\newcommand \colonsim {\vcentcolon \sim }\)
\(\newcommand \Colonsim {\dblcolon \sim }\)
\(\newcommand {\nuparrow }{\mathrel {\cancel {\uparrow }}}\)
\(\newcommand {\ndownarrow }{\mathrel {\cancel {\downarrow }}}\)
\(\newcommand {\bigtimes }{\mathop {\Large \times }\limits }\)
\(\newcommand {\prescript }[3]{{}^{#1}_{#2}#3}\)
\(\newenvironment {lgathered}{\begin {gathered}}{\end {gathered}}\)
\(\newenvironment {rgathered}{\begin {gathered}}{\end {gathered}}\)
\(\newcommand {\splitfrac }[2]{{}^{#1}_{#2}}\)
\(\let \splitdfrac \splitfrac \)
\(\newcommand {\LWRoverlaysymbols }[2]{\mathord {\smash {\mathop {#2\strut }\limits ^{\smash {\lower 3ex{#1}}}}\strut }}\)
\(\newcommand{\alphaup}{\unicode{x03B1}}\)
\(\newcommand{\betaup}{\unicode{x03B2}}\)
\(\newcommand{\gammaup}{\unicode{x03B3}}\)
\(\newcommand{\digammaup}{\unicode{x03DD}}\)
\(\newcommand{\deltaup}{\unicode{x03B4}}\)
\(\newcommand{\epsilonup}{\unicode{x03F5}}\)
\(\newcommand{\varepsilonup}{\unicode{x03B5}}\)
\(\newcommand{\zetaup}{\unicode{x03B6}}\)
\(\newcommand{\etaup}{\unicode{x03B7}}\)
\(\newcommand{\thetaup}{\unicode{x03B8}}\)
\(\newcommand{\varthetaup}{\unicode{x03D1}}\)
\(\newcommand{\iotaup}{\unicode{x03B9}}\)
\(\newcommand{\kappaup}{\unicode{x03BA}}\)
\(\newcommand{\varkappaup}{\unicode{x03F0}}\)
\(\newcommand{\lambdaup}{\unicode{x03BB}}\)
\(\newcommand{\muup}{\unicode{x03BC}}\)
\(\newcommand{\nuup}{\unicode{x03BD}}\)
\(\newcommand{\xiup}{\unicode{x03BE}}\)
\(\newcommand{\omicronup}{\unicode{x03BF}}\)
\(\newcommand{\piup}{\unicode{x03C0}}\)
\(\newcommand{\varpiup}{\unicode{x03D6}}\)
\(\newcommand{\rhoup}{\unicode{x03C1}}\)
\(\newcommand{\varrhoup}{\unicode{x03F1}}\)
\(\newcommand{\sigmaup}{\unicode{x03C3}}\)
\(\newcommand{\varsigmaup}{\unicode{x03C2}}\)
\(\newcommand{\tauup}{\unicode{x03C4}}\)
\(\newcommand{\upsilonup}{\unicode{x03C5}}\)
\(\newcommand{\phiup}{\unicode{x03D5}}\)
\(\newcommand{\varphiup}{\unicode{x03C6}}\)
\(\newcommand{\chiup}{\unicode{x03C7}}\)
\(\newcommand{\psiup}{\unicode{x03C8}}\)
\(\newcommand{\omegaup}{\unicode{x03C9}}\)
\(\newcommand{\Alphaup}{\unicode{x0391}}\)
\(\newcommand{\Betaup}{\unicode{x0392}}\)
\(\newcommand{\Gammaup}{\unicode{x0393}}\)
\(\newcommand{\Digammaup}{\unicode{x03DC}}\)
\(\newcommand{\Deltaup}{\unicode{x0394}}\)
\(\newcommand{\Epsilonup}{\unicode{x0395}}\)
\(\newcommand{\Zetaup}{\unicode{x0396}}\)
\(\newcommand{\Etaup}{\unicode{x0397}}\)
\(\newcommand{\Thetaup}{\unicode{x0398}}\)
\(\newcommand{\Varthetaup}{\unicode{x03F4}}\)
\(\newcommand{\Iotaup}{\unicode{x0399}}\)
\(\newcommand{\Kappaup}{\unicode{x039A}}\)
\(\newcommand{\Lambdaup}{\unicode{x039B}}\)
\(\newcommand{\Muup}{\unicode{x039C}}\)
\(\newcommand{\Nuup}{\unicode{x039D}}\)
\(\newcommand{\Xiup}{\unicode{x039E}}\)
\(\newcommand{\Omicronup}{\unicode{x039F}}\)
\(\newcommand{\Piup}{\unicode{x03A0}}\)
\(\newcommand{\Varpiup}{\unicode{x03D6}}\)
\(\newcommand{\Rhoup}{\unicode{x03A1}}\)
\(\newcommand{\Sigmaup}{\unicode{x03A3}}\)
\(\newcommand{\Tauup}{\unicode{x03A4}}\)
\(\newcommand{\Upsilonup}{\unicode{x03A5}}\)
\(\newcommand{\Phiup}{\unicode{x03A6}}\)
\(\newcommand{\Chiup}{\unicode{x03A7}}\)
\(\newcommand{\Psiup}{\unicode{x03A8}}\)
\(\newcommand{\Omegaup}{\unicode{x03A9}}\)
\(\newcommand{\alphait}{\unicode{x1D6FC}}\)
\(\newcommand{\betait}{\unicode{x1D6FD}}\)
\(\newcommand{\gammait}{\unicode{x1D6FE}}\)
\(\newcommand{\digammait}{\mathit{\unicode{x03DD}}}\)
\(\newcommand{\deltait}{\unicode{x1D6FF}}\)
\(\newcommand{\epsilonit}{\unicode{x1D716}}\)
\(\newcommand{\varepsilonit}{\unicode{x1D700}}\)
\(\newcommand{\zetait}{\unicode{x1D701}}\)
\(\newcommand{\etait}{\unicode{x1D702}}\)
\(\newcommand{\thetait}{\unicode{x1D703}}\)
\(\newcommand{\varthetait}{\unicode{x1D717}}\)
\(\newcommand{\iotait}{\unicode{x1D704}}\)
\(\newcommand{\kappait}{\unicode{x1D705}}\)
\(\newcommand{\varkappait}{\unicode{x1D718}}\)
\(\newcommand{\lambdait}{\unicode{x1D706}}\)
\(\newcommand{\muit}{\unicode{x1D707}}\)
\(\newcommand{\nuit}{\unicode{x1D708}}\)
\(\newcommand{\xiit}{\unicode{x1D709}}\)
\(\newcommand{\omicronit}{\unicode{x1D70A}}\)
\(\newcommand{\piit}{\unicode{x1D70B}}\)
\(\newcommand{\varpiit}{\unicode{x1D71B}}\)
\(\newcommand{\rhoit}{\unicode{x1D70C}}\)
\(\newcommand{\varrhoit}{\unicode{x1D71A}}\)
\(\newcommand{\sigmait}{\unicode{x1D70E}}\)
\(\newcommand{\varsigmait}{\unicode{x1D70D}}\)
\(\newcommand{\tauit}{\unicode{x1D70F}}\)
\(\newcommand{\upsilonit}{\unicode{x1D710}}\)
\(\newcommand{\phiit}{\unicode{x1D719}}\)
\(\newcommand{\varphiit}{\unicode{x1D711}}\)
\(\newcommand{\chiit}{\unicode{x1D712}}\)
\(\newcommand{\psiit}{\unicode{x1D713}}\)
\(\newcommand{\omegait}{\unicode{x1D714}}\)
\(\newcommand{\Alphait}{\unicode{x1D6E2}}\)
\(\newcommand{\Betait}{\unicode{x1D6E3}}\)
\(\newcommand{\Gammait}{\unicode{x1D6E4}}\)
\(\newcommand{\Digammait}{\mathit{\unicode{x03DC}}}\)
\(\newcommand{\Deltait}{\unicode{x1D6E5}}\)
\(\newcommand{\Epsilonit}{\unicode{x1D6E6}}\)
\(\newcommand{\Zetait}{\unicode{x1D6E7}}\)
\(\newcommand{\Etait}{\unicode{x1D6E8}}\)
\(\newcommand{\Thetait}{\unicode{x1D6E9}}\)
\(\newcommand{\Varthetait}{\unicode{x1D6F3}}\)
\(\newcommand{\Iotait}{\unicode{x1D6EA}}\)
\(\newcommand{\Kappait}{\unicode{x1D6EB}}\)
\(\newcommand{\Lambdait}{\unicode{x1D6EC}}\)
\(\newcommand{\Muit}{\unicode{x1D6ED}}\)
\(\newcommand{\Nuit}{\unicode{x1D6EE}}\)
\(\newcommand{\Xiit}{\unicode{x1D6EF}}\)
\(\newcommand{\Omicronit}{\unicode{x1D6F0}}\)
\(\newcommand{\Piit}{\unicode{x1D6F1}}\)
\(\newcommand{\Rhoit}{\unicode{x1D6F2}}\)
\(\newcommand{\Sigmait}{\unicode{x1D6F4}}\)
\(\newcommand{\Tauit}{\unicode{x1D6F5}}\)
\(\newcommand{\Upsilonit}{\unicode{x1D6F6}}\)
\(\newcommand{\Phiit}{\unicode{x1D6F7}}\)
\(\newcommand{\Chiit}{\unicode{x1D6F8}}\)
\(\newcommand{\Psiit}{\unicode{x1D6F9}}\)
\(\newcommand{\Omegait}{\unicode{x1D6FA}}\)
\(\let \digammaup \Digammaup \)
\(\renewcommand {\digammait }{\mathit {\digammaup }}\)
\(\newcommand {\smallin }{\unicode {x220A}}\)
\(\newcommand {\smallowns }{\unicode {x220D}}\)
\(\newcommand {\notsmallin }{\LWRoverlaysymbols {/}{\unicode {x220A}}}\)
\(\newcommand {\notsmallowns }{\LWRoverlaysymbols {/}{\unicode {x220D}}}\)
\(\newcommand {\rightangle }{\unicode {x221F}}\)
\(\newcommand {\intclockwise }{\unicode {x2231}}\)
\(\newcommand {\ointclockwise }{\unicode {x2232}}\)
\(\newcommand {\ointctrclockwise }{\unicode {x2233}}\)
\(\newcommand {\oiint }{\unicode {x222F}}\)
\(\newcommand {\oiiint }{\unicode {x2230}}\)
\(\newcommand {\ddag }{\unicode {x2021}}\)
\(\newcommand {\P }{\unicode {x00B6}}\)
\(\newcommand {\copyright }{\unicode {x00A9}}\)
\(\newcommand {\dag }{\unicode {x2020}}\)
\(\newcommand {\pounds }{\unicode {x00A3}}\)
\(\newcommand {\iddots }{\unicode {x22F0}}\)
\(\newcommand {\utimes }{\overline {\times }}\)
\(\newcommand {\dtimes }{\underline {\times }}\)
\(\newcommand {\udtimes }{\overline {\underline {\times }}}\)
\(\newcommand {\leftwave }{\left \{}\)
\(\newcommand {\rightwave }{\right \}}\)
\(\newcommand {\toprule }[1][]{\hline }\)
\(\let \midrule \toprule \)
\(\let \bottomrule \toprule \)
\(\newcommand {\cmidrule }[2][]{}\)
\(\newcommand {\morecmidrules }{}\)
\(\newcommand {\specialrule }[3]{\hline }\)
\(\newcommand {\addlinespace }[1][]{}\)
\(\newcommand {\LWRsubmultirow }[2][]{#2}\)
\(\newcommand {\LWRmultirow }[2][]{\LWRsubmultirow }\)
\(\newcommand {\multirow }[2][]{\LWRmultirow }\)
\(\newcommand {\mrowcell }{}\)
\(\newcommand {\mcolrowcell }{}\)
\(\newcommand {\STneed }[1]{}\)
\( \newcommand {\multicolumn }[3]{#3}\)
\(\newcommand {\tothe }[1]{^{#1}}\)
\(\newcommand {\raiseto }[2]{{#2}^{#1}}\)
\(\newcommand {\ang }[2][]{(\mathrm {#2})\degree }\)
\(\newcommand {\num }[2][]{\mathrm {#2}}\)
\(\newcommand {\si }[2][]{\mathrm {#2}}\)
\(\newcommand {\LWRSI }[2][]{\mathrm {#1\LWRSInumber \,#2}}\)
\(\newcommand {\SI }[2][]{\def \LWRSInumber {#2}\LWRSI }\)
\(\newcommand {\numlist }[2][]{\mathrm {#2}}\)
\(\newcommand {\numrange }[3][]{\mathrm {#2\,\unicode {x2013}\,#3}}\)
\(\newcommand {\SIlist }[3][]{\mathrm {#2\,#3}}\)
\(\newcommand {\SIrange }[4][]{\mathrm {#2\,#4\,\unicode {x2013}\,#3\,#4}}\)
\(\newcommand {\tablenum }[2][]{\mathrm {#2}}\)
\(\newcommand {\ampere }{\mathrm {A}}\)
\(\newcommand {\candela }{\mathrm {cd}}\)
\(\newcommand {\kelvin }{\mathrm {K}}\)
\(\newcommand {\kilogram }{\mathrm {kg}}\)
\(\newcommand {\metre }{\mathrm {m}}\)
\(\newcommand {\mole }{\mathrm {mol}}\)
\(\newcommand {\second }{\mathrm {s}}\)
\(\newcommand {\becquerel }{\mathrm {Bq}}\)
\(\newcommand {\degreeCelsius }{\unicode {x2103}}\)
\(\newcommand {\coulomb }{\mathrm {C}}\)
\(\newcommand {\farad }{\mathrm {F}}\)
\(\newcommand {\gray }{\mathrm {Gy}}\)
\(\newcommand {\hertz }{\mathrm {Hz}}\)
\(\newcommand {\henry }{\mathrm {H}}\)
\(\newcommand {\joule }{\mathrm {J}}\)
\(\newcommand {\katal }{\mathrm {kat}}\)
\(\newcommand {\lumen }{\mathrm {lm}}\)
\(\newcommand {\lux }{\mathrm {lx}}\)
\(\newcommand {\newton }{\mathrm {N}}\)
\(\newcommand {\ohm }{\mathrm {\Omega }}\)
\(\newcommand {\pascal }{\mathrm {Pa}}\)
\(\newcommand {\radian }{\mathrm {rad}}\)
\(\newcommand {\siemens }{\mathrm {S}}\)
\(\newcommand {\sievert }{\mathrm {Sv}}\)
\(\newcommand {\steradian }{\mathrm {sr}}\)
\(\newcommand {\tesla }{\mathrm {T}}\)
\(\newcommand {\volt }{\mathrm {V}}\)
\(\newcommand {\watt }{\mathrm {W}}\)
\(\newcommand {\weber }{\mathrm {Wb}}\)
\(\newcommand {\day }{\mathrm {d}}\)
\(\newcommand {\degree }{\mathrm {^\circ }}\)
\(\newcommand {\hectare }{\mathrm {ha}}\)
\(\newcommand {\hour }{\mathrm {h}}\)
\(\newcommand {\litre }{\mathrm {l}}\)
\(\newcommand {\liter }{\mathrm {L}}\)
\(\newcommand {\arcminute }{^\prime }\)
\(\newcommand {\minute }{\mathrm {min}}\)
\(\newcommand {\arcsecond }{^{\prime \prime }}\)
\(\newcommand {\tonne }{\mathrm {t}}\)
\(\newcommand {\astronomicalunit }{au}\)
\(\newcommand {\atomicmassunit }{u}\)
\(\newcommand {\bohr }{\mathit {a}_0}\)
\(\newcommand {\clight }{\mathit {c}_0}\)
\(\newcommand {\dalton }{\mathrm {D}_\mathrm {a}}\)
\(\newcommand {\electronmass }{\mathit {m}_{\mathrm {e}}}\)
\(\newcommand {\electronvolt }{\mathrm {eV}}\)
\(\newcommand {\elementarycharge }{\mathit {e}}\)
\(\newcommand {\hartree }{\mathit {E}_{\mathrm {h}}}\)
\(\newcommand {\planckbar }{\mathit {\unicode {x210F}}}\)
\(\newcommand {\angstrom }{\mathrm {\unicode {x212B}}}\)
\(\let \LWRorigbar \bar \)
\(\newcommand {\bar }{\mathrm {bar}}\)
\(\newcommand {\barn }{\mathrm {b}}\)
\(\newcommand {\bel }{\mathrm {B}}\)
\(\newcommand {\decibel }{\mathrm {dB}}\)
\(\newcommand {\knot }{\mathrm {kn}}\)
\(\newcommand {\mmHg }{\mathrm {mmHg}}\)
\(\newcommand {\nauticalmile }{\mathrm {M}}\)
\(\newcommand {\neper }{\mathrm {Np}}\)
\(\newcommand {\yocto }{\mathrm {y}}\)
\(\newcommand {\zepto }{\mathrm {z}}\)
\(\newcommand {\atto }{\mathrm {a}}\)
\(\newcommand {\femto }{\mathrm {f}}\)
\(\newcommand {\pico }{\mathrm {p}}\)
\(\newcommand {\nano }{\mathrm {n}}\)
\(\newcommand {\micro }{\mathrm {\unicode {x00B5}}}\)
\(\newcommand {\milli }{\mathrm {m}}\)
\(\newcommand {\centi }{\mathrm {c}}\)
\(\newcommand {\deci }{\mathrm {d}}\)
\(\newcommand {\deca }{\mathrm {da}}\)
\(\newcommand {\hecto }{\mathrm {h}}\)
\(\newcommand {\kilo }{\mathrm {k}}\)
\(\newcommand {\mega }{\mathrm {M}}\)
\(\newcommand {\giga }{\mathrm {G}}\)
\(\newcommand {\tera }{\mathrm {T}}\)
\(\newcommand {\peta }{\mathrm {P}}\)
\(\newcommand {\exa }{\mathrm {E}}\)
\(\newcommand {\zetta }{\mathrm {Z}}\)
\(\newcommand {\yotta }{\mathrm {Y}}\)
\(\newcommand {\percent }{\mathrm {\%}}\)
\(\newcommand {\meter }{\mathrm {m}}\)
\(\newcommand {\metre }{\mathrm {m}}\)
\(\newcommand {\gram }{\mathrm {g}}\)
\(\newcommand {\kg }{\kilo \gram }\)
\(\newcommand {\of }[1]{_{\mathrm {#1}}}\)
\(\newcommand {\squared }{^2}\)
\(\newcommand {\square }[1]{\mathrm {#1}^2}\)
\(\newcommand {\cubed }{^3}\)
\(\newcommand {\cubic }[1]{\mathrm {#1}^3}\)
\(\newcommand {\per }{/}\)
\(\newcommand {\celsius }{\unicode {x2103}}\)
\(\newcommand {\fg }{\femto \gram }\)
\(\newcommand {\pg }{\pico \gram }\)
\(\newcommand {\ng }{\nano \gram }\)
\(\newcommand {\ug }{\micro \gram }\)
\(\newcommand {\mg }{\milli \gram }\)
\(\newcommand {\g }{\gram }\)
\(\newcommand {\kg }{\kilo \gram }\)
\(\newcommand {\amu }{\mathrm {u}}\)
\(\newcommand {\nm }{\nano \metre }\)
\(\newcommand {\um }{\micro \metre }\)
\(\newcommand {\mm }{\milli \metre }\)
\(\newcommand {\cm }{\centi \metre }\)
\(\newcommand {\dm }{\deci \metre }\)
\(\newcommand {\m }{\metre }\)
\(\newcommand {\km }{\kilo \metre }\)
\(\newcommand {\as }{\atto \second }\)
\(\newcommand {\fs }{\femto \second }\)
\(\newcommand {\ps }{\pico \second }\)
\(\newcommand {\ns }{\nano \second }\)
\(\newcommand {\us }{\micro \second }\)
\(\newcommand {\ms }{\milli \second }\)
\(\newcommand {\s }{\second }\)
\(\newcommand {\fmol }{\femto \mol }\)
\(\newcommand {\pmol }{\pico \mol }\)
\(\newcommand {\nmol }{\nano \mol }\)
\(\newcommand {\umol }{\micro \mol }\)
\(\newcommand {\mmol }{\milli \mol }\)
\(\newcommand {\mol }{\mol }\)
\(\newcommand {\kmol }{\kilo \mol }\)
\(\newcommand {\pA }{\pico \ampere }\)
\(\newcommand {\nA }{\nano \ampere }\)
\(\newcommand {\uA }{\micro \ampere }\)
\(\newcommand {\mA }{\milli \ampere }\)
\(\newcommand {\A }{\ampere }\)
\(\newcommand {\kA }{\kilo \ampere }\)
\(\newcommand {\ul }{\micro \litre }\)
\(\newcommand {\ml }{\milli \litre }\)
\(\newcommand {\l }{\litre }\)
\(\newcommand {\hl }{\hecto \litre }\)
\(\newcommand {\uL }{\micro \liter }\)
\(\newcommand {\mL }{\milli \liter }\)
\(\newcommand {\L }{\liter }\)
\(\newcommand {\hL }{\hecto \liter }\)
\(\newcommand {\mHz }{\milli \hertz }\)
\(\newcommand {\Hz }{\hertz }\)
\(\newcommand {\kHz }{\kilo \hertz }\)
\(\newcommand {\MHz }{\mega \hertz }\)
\(\newcommand {\GHz }{\giga \hertz }\)
\(\newcommand {\THz }{\tera \hertz }\)
\(\newcommand {\mN }{\milli \newton }\)
\(\newcommand {\N }{\newton }\)
\(\newcommand {\kN }{\kilo \newton }\)
\(\newcommand {\MN }{\mega \newton }\)
\(\newcommand {\Pa }{\pascal }\)
\(\newcommand {\kPa }{\kilo \pascal }\)
\(\newcommand {\MPa }{\mega \pascal }\)
\(\newcommand {\GPa }{\giga \pascal }\)
\(\newcommand {\mohm }{\milli \ohm }\)
\(\newcommand {\kohm }{\kilo \ohm }\)
\(\newcommand {\Mohm }{\mega \ohm }\)
\(\newcommand {\pV }{\pico \volt }\)
\(\newcommand {\nV }{\nano \volt }\)
\(\newcommand {\uV }{\micro \volt }\)
\(\newcommand {\mV }{\milli \volt }\)
\(\newcommand {\V }{\volt }\)
\(\newcommand {\kV }{\kilo \volt }\)
\(\newcommand {\W }{\watt }\)
\(\newcommand {\uW }{\micro \watt }\)
\(\newcommand {\mW }{\milli \watt }\)
\(\newcommand {\kW }{\kilo \watt }\)
\(\newcommand {\MW }{\mega \watt }\)
\(\newcommand {\GW }{\giga \watt }\)
\(\newcommand {\J }{\joule }\)
\(\newcommand {\uJ }{\micro \joule }\)
\(\newcommand {\mJ }{\milli \joule }\)
\(\newcommand {\kJ }{\kilo \joule }\)
\(\newcommand {\eV }{\electronvolt }\)
\(\newcommand {\meV }{\milli \electronvolt }\)
\(\newcommand {\keV }{\kilo \electronvolt }\)
\(\newcommand {\MeV }{\mega \electronvolt }\)
\(\newcommand {\GeV }{\giga \electronvolt }\)
\(\newcommand {\TeV }{\tera \electronvolt }\)
\(\newcommand {\kWh }{\kilo \watt \hour }\)
\(\newcommand {\F }{\farad }\)
\(\newcommand {\fF }{\femto \farad }\)
\(\newcommand {\pF }{\pico \farad }\)
\(\newcommand {\K }{\mathrm {K}}\)
\(\newcommand {\dB }{\mathrm {dB}}\)
\(\newcommand {\kibi }{\mathrm {Ki}}\)
\(\newcommand {\mebi }{\mathrm {Mi}}\)
\(\newcommand {\gibi }{\mathrm {Gi}}\)
\(\newcommand {\tebi }{\mathrm {Ti}}\)
\(\newcommand {\pebi }{\mathrm {Pi}}\)
\(\newcommand {\exbi }{\mathrm {Ei}}\)
\(\newcommand {\zebi }{\mathrm {Zi}}\)
\(\newcommand {\yobi }{\mathrm {Yi}}\)
\(\require {mhchem}\)
\(\require {cancel}\)
\(\newcommand {\fint }{âĺŊ}\)
\(\newcommand {\hdots }{\cdots }\)
\(\newcommand {\mathnormal }[1]{#1}\)
\(\newcommand {\vecs }[2]{\vec {#1}_{#2}}\)
\(\renewcommand {\A }{\mathcal {A}}\)
\(\newcommand {\B }{\mathcal {B}}\)
\(\renewcommand {\C }{\mathcal {C}}\)
\(\newcommand {\EE }{\mathbb {E}}\)
\(\renewcommand {\N }{\mathcal {N}}\)
\(\renewcommand {\P }{\mathcal {P}}\)
\(\newcommand {\PP }{\mathbb {P}}\)
\(\newcommand {\T }{\mathcal {T}}\)
\(\renewcommand {\U }{\mathcal {U}}\)
\(\newcommand {\X }{\mathcal {X}}\)
\(\newcommand {\Y }{\mathcal {Y}}\)
\(\newcommand {\Bin }{\operatorname {Bin}}\)
\(\newcommand {\Pois }{\operatorname {Pois}}\)
\(\newcommand {\Exp }{\operatorname {Exp}}\)
\(\newcommand {\BetaV }{\operatorname {Beta}}\)
\(\newcommand {\GammaV }{\operatorname {Gamma}}\)
\(\newcommand {\pot }{\mathfrak {P}}\)
\(\renewcommand {\1}{𝟙}\)
\(\newcommand {\id }{\mathrm {id}}\)
\(\renewcommand {\i }{\mathrm {i}}\)
\(\renewcommand {\theta }{\vartheta }\)
\(\newcommand {\Tu }{\underline {T}}\)
\(\newcommand {\To }{\overline {T}}\)
\(\newcommand {\Var }{\mathrm {Var}}\)
\(\newcommand {\Cov }{\mathrm {Cov}}\)
\(\newcommand {\name }[1]{\textsc {#1}}\)
\(\newcommand {\smallpmatrix }[1]{\left (\begin {smallmatrix}#1\end {smallmatrix}\right )}\)
\(\newcommand {\matlab }{{\fontfamily {bch}\scshape \selectfont {}Matlab}}\)
\(\newcommand {\innerproduct }[1]{\left \langle {#1}\right \rangle }\)
\(\newcommand {\norm }[1]{\left \Vert {#1}\right \Vert }\)
\(\renewcommand {\natural }{\mathbb {N}}\)
\(\newcommand {\integer }{\mathbb {Z}}\)
\(\newcommand {\rational }{\mathbb {Q}}\)
\(\newcommand {\real }{\mathbb {R}}\)
\(\newcommand {\complex }{\mathbb {C}}\)
\(\renewcommand {\d }{\mathop {}\!\mathrm {d}}\)
\(\newcommand {\dr }{\d {}r}\)
\(\newcommand {\ds }{\d {}s}\)
\(\newcommand {\dt }{\d {}t}\)
\(\newcommand {\du }{\d {}u}\)
\(\newcommand {\dv }{\d {}v}\)
\(\newcommand {\dw }{\d {}w}\)
\(\newcommand {\dx }{\d {}x}\)
\(\newcommand {\dy }{\d {}y}\)
\(\newcommand {\dz }{\d {}z}\)
\(\newcommand {\dsigma }{\d {}\sigma }\)
\(\newcommand {\dphi }{\d {}\phi }\)
\(\newcommand {\dvarphi }{\d {}\varphi }\)
\(\newcommand {\dtau }{\d {}\tau }\)
\(\newcommand {\dxi }{\d {}\xi }\)
\(\newcommand {\dtheta }{\d {}\theta }\)
\(\newcommand {\tp }{\mathrm {T}}\)
Schätzkriterien
Bemerkung: Gegeben ist wieder ein reguläres statistisches Modell \(\P = \{\PP _\theta \;|\; \theta \in \Theta \}\). Wie kann die Qualität eines Schätzers \(T = T(X)\)
für den Parameter \(\theta \) beurteilt werden?
Der Schätzfehler \(E = |T(X) - q(\theta )|\) ist hierfür ungeeignet, da \(E\) vom unbekannten Parameter \(\theta \) und vom Zufall abhängt (d. h. \(E\) kann erst nach Vorliegen der
Stichprobe berechnet werden).
mittlerer quadratischer Fehler: Sei \(T = T(X)\) ein Schätzer für \(q(\theta ) \in \real \). Der mittlere quadratische Fehler (MQF) von
\(T\) ist definiert durch \(R(\theta , T) := \EE _\theta ((T(X) - q(\theta ))^2)\).
Verzerrung/Bias: \(b(\theta , T) = \EE _\theta (T(X)) - q(\theta )\) heißt Verzerrung oder Bias von
\(T\).
unverzerrt/erwartungstreu:
Gilt \(b(\theta , T) = 0\) für alle \(\theta \in \Theta \), so heißt \(T\) unverzerrt oder erwartungstreu.
Bemerkung: Für den MQF gilt \(R(\theta , T) = \EE _\theta ((T(X) - \EE _\theta (T(X)) + \EE _\theta (T(X)) - q(\theta ))^2)\)
\(= \EE _\theta ((T(X) - \EE _\theta (T(X)))^2) + \EE _\theta ((\EE _\theta (T(X)) - q(\theta ))^2) = \Var _\theta (T(X)) + b(\theta , T)^2\) unabhängig vom Zufall, da \(\EE _\theta
((T(X) - \EE _\theta (T(X)))^2) = 0\). Diese Zerlegung heißt Varianz-Bias-Zerlegung des MQF.
Lemma (Varianz-Bias-Zerlegung): Es gilt \(R(\theta , T) = \Var _\theta (T(X)) + b(\theta , T)^2\).
Beispiel: Seien \(X_1, \dotsc , X_n \sim \N (\mu , \sigma ^2)\) i.i.d. mit \(\theta = (\mu , \sigma ^2) \in \real \times \real ^+\) unbekannt. Weiter oben wurde erwähnt, dass \(\widehat
{\theta } = (\frac {1}{n} \sum _{i=1}^n X_i, \frac {1}{n} (X_i - \overline {X})^2)\) der MLS für \(\theta \) ist, wobei \(\overline {X} \sim \N (\mu , \frac {\sigma ^2}{n})\) gilt.
Für \(q(\theta ) = \mu \) gilt daher \(b(\theta , \overline {X}) = \EE _\theta (\overline {X}) - q(\theta ) = 0\), d. h. \(\overline {X}\) ist ein erwartungstreuer Schätzer
für \(\mu \). Für den MQF gilt \(R(\theta , \overline {X}) = \Var _\theta (\overline {X}) + b(\theta , T)^2 = \frac {\sigma ^2}{n} \to 0\) (\(n \to \infty \)).
Für \(q(\theta ) = \sigma ^2\) ist \(\widehat {\sigma }^2 = \sigma ^2(X) := \frac {1}{n} \sum _{i=1}^n (X_i - \overline {X})^2\) ein brauchbarer Schätzer für \(\sigma ^2\), da
\(\EE _\theta (\widehat {\sigma }^2) = \frac {n-1}{n} \EE _\theta (\frac {1}{n-1} \sum _{i=1}^n (X_i - \overline {X})^2) = \frac {n-1}{n} \sigma ^2 \to \sigma ^2\) (\(n \to \infty \)),
d. h. \(\widehat {\sigma }^2\) ist asymptotisch unverzerrt. Im Gegensatz zur Stichprobenvarianz \(S^2(X) := \frac {n}{n-1} \sigma ^2\) ist die empirische Varianz
\(\sigma ^2(X)\) also kein erwartungstreuer Schätzer für \(\sigma ^2\). Mit \(S := \frac {n\widehat {\sigma }^2}{\sigma ^2} = \sum _{i=1}^n \left (\frac {X_i - \overline {X}}{\sigma
}\right )^2 \sim \chi _{n-1}^2\) gilt \(\EE (S) = n - 1\) und \(\Var (S) = 2(n - 1)\), damit lässt sich der MQF berechnen als
\(R(\theta , \widehat {\sigma }^2) = \frac {\sigma ^2}{n} \Var (S) + (\frac {\sigma ^2}{n} \EE (S) - \sigma ^2)^2 = \left (\frac {\sigma ^2}{n}\right )^2 (2n-1) \to 0\) (\(n \to \infty \)).
Beispiel: Man kann zwei Mittelwertschätzer für \(X_1, \dotsc , X_n \sim \N (\mu , \sigma ^2)\) i.i.d., \(X = (X_1, \dotsc , X_n)\) und \(\theta = (\mu , \sigma ^2)\) unbekannt
mittels MQF vergleichen. Dazu werden für \(q(\theta ) = \mu \) die Schätzer \(T_1(X) := \overline {X}\) und \(T_2(X) := a \overline {X}\) für ein \(a \in (0, 1)\) betrachtet.
Für \(T_1\) gilt \(b(\mu , T_1) = 0\) und \(R(\mu , T_1) = \frac {\sigma ^2}{n}\) wie eben berechnet.
Für \(T_2\) gilt \(b(\mu , T_2) = \EE _\theta (T_2(X)) - \mu = (a - 1)\mu \) und
\(R(\mu , T_2) = \Var _\theta (T_2(X)) + b(\mu , T_2)^2 = \frac {a^2 \sigma ^2}{n} + (a-1)^2 \mu ^2\).
\(T_2\) ist also nicht mehr erwartungstreu, aber für \(\mu \) in einer hinreichend kleinen Umgebung von \(0\) gilt \(R(\mu , T_2) < R(\mu , T_1)\), d. h. für ein solches \(\mu \)
schätzt \(T_2\) besser. Für \(\mu \) hinreichend groß gilt natürlich \(R(\mu , T_1) < R(\mu , T_2)\).
unzulässig: Ein Schätzer \(S = S(X)\) heißt unzulässig, falls es einen Schätzer \(T = T(X)\) gibt mit \(\forall
_{\theta \in \Theta }\; R(\theta , T) \le R(\theta , S)\) und \(\exists _{\theta \in \Theta }\; R(\theta , T) < R(\theta , S)\).
Bemerkung: Es kann keinen perfekten, „besten“ Schätzer \(T\) geben mit \(\forall _{\theta \in \Theta }\; R(\theta , T) \le R(\theta , S)\) für jeden anderen Schätzer \(S\).
Wählt man nämlich ein festes \(\theta _0 \in \Theta \) und setzt \(S(X) := q(\theta _0)\), so gilt \(R(\theta _0, S) = 0\). Für den „besten“ Schätzer \(T\) müsste also
\(\forall _{\theta \in \Theta }\; R(\theta , T) = 0\) gelten, was offensichtlich nicht geht. Daher müssen wir die Klasse der Vergleichskandidaten für einen guten Schätzer auf die Klasse der
unverzerrten Schätzer einschränken (\(S\) ist nicht unverzerrt).
Unverzerrte Schätzer mit gleichmäßig minimaler Varianz
UMVU-Schätzer: Ein unverzerrter Schätzer \(T(X)\) für \(q(\theta )\) heißt UMVU-Schätzer
(uniformly minimal variance unbiased, UMVUE) für \(q(\theta )\), falls für alle unverzerrten Schätzer \(S(X)\) für \(q(\theta )\) gilt, dass \(\forall _{\theta \in
\Theta }\; \Var _\theta (T(X)) \le \Var _\theta (S(X))\).
Satz (Rao-Blackwell): Seien \(T(X)\) ein suffizienter Schätzer
für \(\theta \) und \(S(X)\) ein beliebiger Schätzer für \(q(\theta )\) mit \(\forall _{\theta \in \Theta }\; \EE _\theta (|S(X)|) < \infty \). Dann ist der Schätzer
\(T^\ast (X) := \EE (S(X)|T(X))\) für \(q(\theta )\) unabhängig von \(\theta \) und es gilt \(\forall _{\theta \in \Theta }\; R(\theta , T^\ast ) \le R(\theta , S)\). Ist
zusätzlich \(\Var _\theta (S(X)) < \infty \), so gilt Gleichheit genau dann, wenn \(\forall _{\theta \in \Theta }\; \PP _\theta (T^\ast (X) = S(X)) = 1\).
Bemerkung: Der Schätzer \(T^\ast (X)\) für \(q(\theta )\) ist also mindestens so gut wie \(S(X)\).
\(T^\ast (X)\) ist unabhängig von \(\theta \), weil \(T(X)\) ein suffizienter Schätzer ist.
Beispiel: Für \(X_1, X_2 \sim \Bin (1, p)\) i.i.d. und \(X := (X_1, X_2)\) ist \(T(X) := X_1 + X_2\) eine suffiziente Statistik für \(p\). Wählt man \(S(X) := X_1\), so gilt
für \(t = 0, 1, 2\), dass
\(\EE (S(X)|T(X)=t) = \EE (X_1|X_1+X_2 = t) = \frac {\PP (X_1 = 1, X_1 + X_2 = t)}{\PP (X_1 + X_2 = t)} = \frac {\PP (X_1 = 1, X_2 = t - 1)}{\PP (X_1 + X_2 = t)} = \frac {\PP (X_1 = 1) \cdot
\PP (X_2 = t - 1)}{\PP (X_1 + X_2 = t)}\)
\(= \frac {p \cdot p^{t-1} (1-p)^{1-(t-1)}}{\binom {2}{t} p^t (1 - p)^{2-t}} = \frac {1}{\binom {2}{t}} = \frac {t}{2}\). Somit ist \(T^\ast (X) = \EE (S(X)|T(X)) = \frac {T(X)}{2} = \overline
{X}\).
Berechnet man den MQF von \(T^\ast \), so erhält man
\(\EE _p((T^\ast (X) - p)^2) = \EE _p((\overline {X}-p)^2) = \Var _p(\overline {X}) = \frac {p(1-p)}{2}\). Dies ist echt kleiner als der MQF von \(S\):
\(\EE _p((S(X) - p)^2) = \EE _p((X_1 - p)^2) = \Var _p(X_1) = p(1-p)\). Nach dem Satz von Rao-Blackwell muss dies auch so sein, denn \(\Var _p(S(X)) = p(1-p) < \infty \) und \(\PP _p(\overline {X} =
X_1) < 1\) für alle \(p \in (0, 1)\).
Bemerkung: Für eine Eindeutigkeitsaussage benötigt man den Begriff der vollständigen Statistik. Er besagt, dass \(\forall _{\theta \in \Theta }\; (\EE _\theta (g_1(T(X))) = \EE
_\theta (g_2(T(X))) \;\Rightarrow \; g_1 = g_2)\). Äquivalent dazu ist folgende Definition. Vollständigkeit ist eigentlich eine Eigenschaft von \(\P = \{\PP _\theta \;|\; \theta \in \Theta \}\)
und verlangt eine gewisse Größe von \(\Theta \), um die Implikation zu erzwingen.
vollständig: Eine Statistik \(T(X)\) heißt vollständig, falls für jede messbare Abbildung \(g\) gilt, dass aus \(\forall _{\theta
\in \Theta }\; \EE _\theta (g(T(X))) = 0\) folgt, dass \(\forall _{\theta \in \Theta }\; \PP _\theta (g(T(X)) = 0) = 1\).
Beispiel: Seien \(X_1, \dotsc , X_n \sim \Pois (\theta )\) i.i.d. mit \(\theta \in \Theta := \real ^+\). \(T(X) = X_1 + \dotsb + X_n\) ist eine suffiziente Statistik für \(\theta \) mit \(T(X)
\sim \Pois (n\theta )\). Sei \(g\colon \real \rightarrow \real \) messbar mit \(\forall _{\theta \in \Theta }\; \EE _\theta (g(T(X))) = 0\). Dann gilt \(\EE _\theta (g(T(X))) = \sum _{i=0}^\infty g(i)
\cdot e^{-n\theta } \frac {(n\theta )!}{\i !}\). Dies ist eine Potenzreihe in \(n\theta \), die in einer gewissen Umgebung von \(0\) gleich null ist. Mit dem Eindeutigkeitssatz für Potenzreihen folgt, dass \(g(i)
= 0\) für alle \(i \in \natural _0\), d. h. \(g \equiv 0\) und \(T(X)\) ist eine vollständige Statistik.
Satz (Lehmann-Scheffé): Seien \(T(X)\)
eine vollständige, suffiziente Statistik für \(\theta \) und \(S(X)\) ein unverzerrter Schätzer für \(q(\theta )\).
Dann ist der Rao-Blackwell-Schätzer \(T^\ast (X) := \EE (S(X)|T(X))\) ein UMVU-Schätzer für \(q(\theta )\).
Ist zusätzlich \(\forall _{\theta \in \Theta }\; \Var _\theta (T^\ast (X)) < \infty \), so ist \(T^\ast (X)\) der eindeutige UMVU-Schätzer für \(q(\theta )\).
Bemerkung: Ist \(h(T(X))\) ein unverzerrter Schätzer für \(q(\theta )\) und \(T(X)\) eine vollständige, suffiziente Statistik für \(\theta \), so ist \(h(T(X))\) ein
UMVU-Schätzer für \(q(\theta )\), da wegen
\(\EE (h(T(X))|T(X)) = h(T(X))\) der Schätzer \(S(X) := h(T(X))\) gewählt werden kann.
Satz (Konstruktion von vollständigen, suffizienten Statistiken): Ist \(\P = \{\PP _\theta \;|\; \theta \in \Theta \}\) eine
\(k\)-parametrige Exponentialfamilie und enthält \(c(\Theta ) \subset \real ^k\) mit \(c = (c_1, \dotsc , c_k)\) ein offenes Rechteck in \(\real ^k\), so ist \(T(X) := (T_1(X), \dotsc , T_k(X))\) eine
vollständige, suffiziente Statistik für \(\theta \in \Theta \subset \real ^k\).
Beispiel: Bei \(X_1, \dotsc , X_n \sim \N (\mu , \sigma ^2)\) i.i.d. mit \(X = (X_1, \dotsc , X_n)\) und \(\theta = (\mu , \sigma ^2)\) unbekannt ist \(T(X) := (\sum _{i=1}^n X_i, \sum _{i=1}^n
X_i^2)\) nach dem letzten Satz eine vollständige, suffiziente Statistik für \(\theta \in \Theta := \real \times \real ^+\). \(\overline {X} = \frac {1}{n} T_1(X) =: h(T(X))\) ist ein unverzerrter
Schätzer für \(\mu = q(\theta )\). Damit ist \(\overline {X}\) nach dem Satz von Lehmann-Scheffé ein UMVU-Schätzer, der sogar eindeutig ist, da \(\forall _{\theta \in \Theta }\; \Var
_\theta (h(T(X))) = \Var _\theta (\overline {X}) = \frac {\sigma ^2}{n} < \infty \).
Da \(S^2(X) = \frac {1}{n-1} \sum _{i=1}^n (X_i - \overline {X})^2\) ein erwartungstreuer Schätzer für \(\sigma ^2\) ist und in der Form \(h(T(X)) := \frac {1}{n-1} \left (T_2(X) - \frac
{1}{n} (T_1(X))^2\right )\) geschrieben werden kann, folgt analog, dass auch \(S^2(X)\) ein (bzw. der eindeutige) UMVU-Schätzer für \(\sigma ^2 = q(\theta )\) ist.
Bemerkung: Im Folgenden sei \(\P = \{\PP _\theta \;|\; \theta \in \Theta \}\) immer ein einparametriges, reguläres statistisches Modell, das folgende Regularitätsbedingungen (CR)
erfülle.
Cramér-Rao-Regularitätsbedingungen (CR):
\(\Theta \subset \real \) sei offen.
\(A := \{x \in \real ^n \;|\; p(x, \theta ) > 0\}\) sei unabhängig von \(\theta \) und es gelte
\(\forall _{x \in A} \forall _{\theta \in \Theta }\; (\frac {\partial }{\partial \theta } \ln p(x, \theta ) \text { existiert und ist endlich})\).
Hat \(X\) eine L.-B.-Dichte und ist \(T = T(X)\) eine Statistik mit \(\forall _{\theta \in \Theta }\; \EE _\theta (|T(X)|) < \infty \), so gelte \(\frac {\partial }{\partial \theta } \int
_{\real ^n} T(x) p(x, \theta ) \dx = \int _{\real ^n} T(x) \frac {\partial }{\partial \theta } p(x, \theta ) \dx \).
Bemerkung: Ist \(\P \) eine \(1\)-parametrige Exp.fam. mit \(p(x, \theta ) = \1_A(x) \exp (c(\theta )T(x) + d(\theta ) + S(x))\), wobei \(\forall _{\theta \in \Theta }\; \frac {\partial
}{\partial \theta } c(\theta ) \not = 0\), \(\Theta \subset \real \) offen und \(c\) stetig ist, dann ist (CR) erfüllt.
Fisher-Information: Die Fisher-Information für einen Parameter \(\theta \) ist
gegeben durch
\(I(\theta ) := \EE _\theta \!\left (\left (\frac {\partial }{\partial \theta } \ln p(X, \theta ) \right )^2\right ) \in [0, \infty ]\). Die partielle Ableitung heißt Score-Funktion.
Bemerkung: Es gilt \(I(\theta ) = \int _{\real ^n} \frac {1}{p(x, \theta )} \cdot \left (\frac {\partial }{\partial \theta } p(x, \theta )\right )^2 \dx \).
Für die Score-Funktion gilt \(\EE _\theta \!\left (\frac {\partial }{\partial \theta } \ln p(X, \theta )\right ) = 0\). Damit ist also \(I(\theta ) = \Var _\theta \!\left (\frac {\partial
}{\partial \theta } \ln p(X, \theta )\right )\).
Bemerkung: Gilt \(X = (X_1, \dotsc , X_n)\) mit \(X_1, \dotsc , X_n\) i.i.d., dann gilt
\(I(\theta ) = n \cdot \EE _\theta \!\left (\left (\frac {\partial }{\partial \theta } \ln p_1(X_1, \theta )\right )^2\right )\). Verdoppelt man also den Stichprobenumfang, so verdoppelt sich die
Fisher-Information.
Beispiel: Für \(X \sim \N (\mu , \sigma ^2)\) mit \(\theta = \mu \) unbekannt und \(\sigma ^2\) bekannt gilt \(I(\theta ) = \frac {1}{\sigma ^2}\).
Satz (Informationsungleichung von Cramér-Rao):
Sei \(T(X)\) eine Statistik mit \(\forall _{\theta \in \Theta }\; \Var _\theta (T(X)) < \infty \). Außerdem sei (CR) erfüllt und es gelte \(\forall _{\theta \in \Theta }\; 0 < I(\theta ) <
\infty \). Dann gilt für \(\psi (\theta ) := \EE _\theta (T(X))\), dass \(\psi (\theta )\) für alle \(\theta \in \Theta \) differenzierbar ist und \(\Var _\theta (T(X)) \ge \frac {(\psi
’(\theta ))^2}{I(\theta )}\).
Folgerung: Ist \(T = T(X)\) ein unverzerrter Schätzer für \(\theta \), so gilt unter den Voraussetzungen von eben die Cramér-Rao-Schranke \(\Var _\theta (T(X)) \ge \frac {1}{I(\theta )}\) für alle \(\theta \in \Theta \).
Folgerung: Sind \(X_1, \dotsc , X_n\) i.i.d. mit \(X = (X_1, \dotsc , X_n)\), so gilt unter den Voraussetzungen von eben \(\Var _\theta (T(X)) \ge \frac {(\psi ’(\theta ))^2}{n \cdot I_1(\theta )}\)
mit \(I_1(\theta ) = \EE _\theta \!\left (\left (\frac {\partial }{\partial \theta } \ln p_1(X_1, \theta ) \right )^2\right )\) der Fisher-Information der einzelnen Beobachtung.
Bemerkung: Dafür benötigt man die C.-S.-Ungleichung \(|\Cov (X,Y)| \le (\Var (X))^{1/2} (\Var (Y))^{1/2}\). Man kann zeigen, dass die Cramér-Rao-Schranke nur für
Exponentialfamilien erfüllt sein kann.
Asymptotische Theorie
(schwach) konsistent: Eine Folge von Schätzern \(T_n := T_n(X_1, \dotsc , X_n)\) für \(q(\theta )\) heißt
(schwach) konsistent, falls \(T_n \xrightarrow {\PP _\theta } q(\theta )\) für \(n \to \infty \) und alle \(\theta \in \Theta \).
stark konsistent: Eine Folge von Schätzern \(T_n := T_n(X_1, \dotsc , X_n)\) für \(q(\theta )\) heißt
stark konsistent, falls \(T_n \xrightarrow {n \to \infty } q(\theta )\) \(\PP _\theta \)-f.s. für alle \(\theta \in \Theta \).
Bemerkung: Aus starker Konsistenz folgt immer schwache Konsistenz. Umgekehrt kann man zeigen, dass ein schwach konsistenter Schätzer stark konsistent ist, wenn die stochastische Konvergenz schnell genug
ist.
UMVU-Schätzer sind immer konsistent und MLS sind in der Regel auch konsistent.
Satz (Log-Likelihood-Funktion zum wahren Parameter besser):
Sei \(\P = \{p(\cdot , \theta ) \;|\; \theta \in \Theta \}\) ein reguläres statistisches Modell mit
\(\forall _{\theta , \theta ’ \in \Theta }\; \EE _\theta (\ln p(X, \theta ’)) < \infty \) und
\(\forall _{\theta , \theta ’ \in \Theta ,\; \theta \not = \theta ’}\; \PP _\theta \not = \PP _{\theta ’}\) (\(\P \) ist identifizierbar).
Dann gilt \(\forall _{\theta , \theta ’ \in \Theta ,\; \theta \not = \theta ’}\; \EE _\theta (\ln p(X, \theta )) > \EE _\theta (\ln p(X, \theta ’))\), d. h. die Log-Likelihood-Funktion zum
wahren Parameter \(\theta \) ist im Mittel strikt besser als die Log-Likelihood-Funktion zu einem anderen Parameter \(\theta ’\).
Satz (Konsistenz des MLS):
Seien \(\Theta \subset \real ^k\) kompakt und \(\P = \{p(\cdot , \theta ) \;|\; \theta \in \Theta \}\) ein reguläres statistisches Modell mit
\(\forall _{\theta , \theta ’ \in \Theta }\; \EE _\theta (\ln p(X, \theta ’)) < \infty \) und
\(\forall _{\varepsilon > 0} \exists _{\delta > 0} \forall _{\theta , \theta ’ \in \Theta ,\; \norm {\theta - \theta ’} < \delta } \forall _{x \in \real ^n}\; |\ln p(x,
\theta ) - \ln p(x, \theta ’)| < \varepsilon \)
(\(\ln p(x, \cdot )\colon \Theta \rightarrow \real \) ist gleichmäßig gleichgradig stetig).
Dann ist jeder MLS \(\widehat {\theta }_n\) mit Likel.-Funktion \(L(\theta , (X_1, \dotsc , X_n)) = \prod _{i=1}^n p(X_i, \theta )\) stark konsistent.
Bemerkung: Die Voraussetzungen \(\Theta \subset \real ^k\), \(\Theta \) kompakt und gleichmäßige gleichgradige Stetigkeit der Log-Likelihood-Funktion werden für die starke Konsistenz des
MLS nicht benötigt.
Bemerkung: Ist ein MLS konsistent, so ist die Konvergenzgeschwindigkeit und die asymptotische Verteilung des MLS häufig von großem Interesse.
Satz (asymptotische Normalität des MLS):
Sei \(\P = \{p(\cdot , \theta ) \;|\; \theta \in \Theta \}\) ein reguläres statistisches Modell mit \(\Theta \subset \real \) offen und
\(B := \{x \in \real \;|\ p(x, \theta ) > 0\}\) unabhängig von \(\theta \),
\(\forall _{\theta , \theta ’ \in \Theta ,\; \theta \not = \theta ’}\; \PP _\theta \not = \PP _{\theta ’}\) (\(\P \) ist identifizierbar),
\(\forall _{x \in \real }\; p(x, \cdot )\colon \Theta \rightarrow \real \) dreifach stetig differenzierbar,
-
(4)\(\frac {\partial ^k}{\partial \theta ^k} \int _B p(x, \theta ) d\mu (x) = \int _B \frac {\partial ^k}{\partial \theta ^k} p(x, \theta ) d\mu (x)\) für \(k = 1, 2\) (mit \(\mu
\) dem Zähl- oder L.-B.-Maß),
\(\forall _{\theta \in \Theta } \exists _{c_\theta > 0} \exists _{g_\theta \colon \real \rightarrow \real _0^+,\; \EE _\theta (g_\theta (X_1)) < \infty } \forall _{\theta
^\ast \in \Theta ,\; |\theta ^\ast - \theta | < c_\theta } \forall _{x \in \real }\; \left |\frac {\partial ^3}{\partial \theta ^3} \left .\ln p(x, \theta )\right |_{\theta = \theta ^\ast
}\right | \le g_\theta (x)\),
\(\forall _{\theta \in \Theta }\; I(\theta ) = \Var _\theta \!\left (\frac {\partial }{\partial \theta } \ln p(X, \theta )\right ) \in (0, \infty )\) und
\(\forall _{\theta \in \Theta }\; \widehat {\theta }_n \text { schwach konsistenter MLS für } \theta \).
Dann gilt \(\sqrt {n I(\theta )} (\widehat {\theta }_n - \theta ) \xrightarrow {\text {(d)}} Z\) mit \(Z \sim \N (0, 1)\) für alle \(\theta \in \Theta \).