Fachbereich Mathematik Prof. Dr. Michael Kohler Dominique Küpper
WS 09/10 02.02.2010
11. Übungsblatt zur
„Mathematik und Statistik für Biologie“
Aufgabe 41
Ward und Quinn sammelten Eikapseln der Raubschneckenart “lepsiella vinosa“ in unterschiedlichen Gebieten einer felsigen Gezeitenküste. Untersucht wurden Unterschiede bei der Fruchtbarkeit der Schneckenart in den verschiedenen Gebieten. In einem bestimmten Gebiet wurden 37 Eikapseln gesammelt. Bei den in diesem Gebiet gesammelten Kapseln lag die durchschnittliche Anzahl der Eier pro Kapsel bei8.07bei einer empirischen Standardabweichung von2.03. Wir nehmen an, dass die ermittelten Werte Realisierungen von unabhängigen und identisch verteilten Zufallsvariablen X1, . . . , X37sind. Außerdem können wir auf Grund weiterer Untersuchungen davon ausgehen, dass die empirische Standardabweichung mit der wirklichen Standardabweichung übereinstimmt.
Wir möchten nun ein zweiseitiges KonfidenzintervallKn(x1, . . . , xn) zum Niveauα= 0.95 für den Erwartungswert der Zufallsvariablen X1 bestimmen. D.h. es soll gelten:
P[EX1 ∈Kn(x1, . . . , xn)] ≥ α= 0.95.
(a) Zeigen Sie mit Hilfe von bekannten Rechenregeln für Erwartungswert und Varianz, dass Pn
i=1Xi−E(Pn i=1Xi) pV(Pn
i=1Xi) =
√n pV(X1)
1 n
Xn
i=1Xi−EX1
.
(b) Zeigen Sie nun, dass (für u >0) aus
√n
√V(X1) 1 n
Pn
i=1Xi−EX1
≤u folgt, dass
EX1 ∈
"
1 n
Xn
i=1Xi−
pV(X1)
√n ·u , 1 n
Xn
i=1Xi+
pV(X1)
√n ·u
# .
(c) Benutzen Sie nun denzentralen Grenzwertsatz und dass für eine standardnormalverteilte (d.h.
N(0,1)-verteilte) Zufallsvariable Z
P(|Z| ≤u0.025) ≥ 0.95 mit u0.025≈1.96 ist, um das gesuchte Konfidenzintervall zu bestimmen.
Lösung:
(a) Nach der Vorlesung istE(Pn
i=1Xi) =Pn
i=1EXi =n·EX1undV(Pn
i=1Xi) =Pn
i=1V(Xi) = n·V(X1), fallsX1, . . . , Xn unabhängig und identisch verteilt sind. Damit ergibt sich:
Pn
i=1Xi−E(Pn i=1Xi) pV(Pn
i=1Xi) =
Pn
i=1Xi−n·EX1
pn·V(X1) = n n1Pn
i=1Xi−EX1
√np V(X1)
=
√n n1Pn
i=1Xi−EX1 pV(X1) =
√n pV(X1)
1 n
Pn
i=1Xi−EX1 .
11. Übung Mathematik und Statistik für Biologie
(b) Auflösen des Betrages ergibt:
√n pV(X1)
1 n
Xn
i=1Xi−EX1
≤u
⇔
√n pV(X1)
1 n
Xn
i=1Xi−EX1
≤u ∧ −
√n pV(X1)
1 n
Xn
i=1Xi−EX1
≤u
⇔ 1 n
Xn
i=1Xi−EX1 ≤
pV(X1)
√n ·u ∧ 1 n
Xn
i=1Xi−EX1≥ −
pV(X1)
√n ·u
⇔ 1 n
Xn
i=1Xi−
pV(X1)
√n ·u≤EX1 ∧ 1 n
Xn
i=1Xi+
pV(X1)
√n ·u≥EX1
Was aber natürlich gleichbedeutend ist zu:
EX1 ∈
"
1 n
Xn
i=1Xi−
pV(X1)
√n ·u , 1 n
Xn
i=1Xi+
pV(X1)
√n ·u
# .
(c) Nach dem zentralen Grenzwertsatz und Teil (a) ist Z := √√n
V(X1) 1 n
Pn
i=1Xi−EX1 eine annähernd N(0,1)-verteilte Zufallsvariable. Damit gilt mit Teil (b):
P(|Z| ≤u0.025) =P
EX1∈
1 n
n
P
i=1
Xi−
√V(X1)
√n ·u0.025 , n1
n
P
i=1
Xi+
√V(X1)
√n ·u0.025
≥ 0.95
Also ist das gesuchte Konfidenzintervall gegeben durch:
8.07− 2.03
√37 ·1.96 , 8.07 + 2.03
√37 ·1.96
≈[7.4159 , 8.7241].
Aufgabe 42
Gegeben seien auf dem Intervall[a,3a]unabhängig identisch gleichverteilte ZufallsvariablenX1, . . . , Xn, n∈Nund a∈R. Wir möchten nun den Parameter
a= 1 2EX1
mit Hilfe von Realisierungen von X1, . . . , Xn schätzen. Geben Sie eine geeignete Schätzung dafür an und zeigen Sie, dass die Schätzung erwartungstreu und konsistent ist.
Hinweis: Modifizieren Sie den Schätzer aus der Vorlesung.
Lösung: Es ist bekannt, dass eine auf dem Intervall [a, b] gleichverteilte Zufallsvariablen den Erwartungswert a+b2 haben, wobei der Schätzer für den Erwartungswert aus der Vorlesung bekannt ist. Damit ist der gesuchte Schätzer
T(x1, . . . , xn) =1
2·x1+. . .+xn
n mit
E(T(X1, . . . , Xn)) =E 1
2n Xn
i=1Xi
= 1 2n
n
X
i=1
E(Xi) = 1
2E(X1) =a
erwartungstreu. Er ist auch konsistent, weil 2·T(X1, . . . , Xn) nach Vorlesung ein konsistenter Schätzer ist.
2
11. Übung Mathematik und Statistik für Biologie
Aufgabe 43
Im Rahmen einer Studie möchte Biologie B. untersuchen, wie weit ein bestimmtes Merkmal bei Honigbienen verbreitet ist. Dazu betrachtet er zufällig ausgewählte Bienen mehrerer Bienenstöcke.
B. überlegt sich folgende Modellannahmen: Jede Biene hat das untersuchte Merkmal mit der selben Wahrscheinlichkeit p unabhängig von allen anderen Tieren. Die Zufallsvariable Xi verwendet er zur mathematischen Modellierung deri-ten Beobachtung.Xi erhält den Wert1, falls diei-te Biene das Merkmal aufweist und 0 sonst. Die Gesamtanzahl aller betrachteten Bienen sein.
Geben Sie einen erwartungstreuen und konsistenten Schätzer für den Parameterp an und begrün- den Sie, dass der von ihnen angegebene Schätzer diese Eigenschaften hat.
Lösung: X1, . . . , Xn sind jeweils binomialverteilt mit Parametern 1 und p. Außerdem sind sie wegen der gemachten Modellannahmen unabhängig. Für die Binomialverteilung mit Parametern 1 und pgilt
p=EX.
Nach der Vorlesung ist
T(x1, . . . , xn) = x1+· · ·+xn
n
ein konsistenter und erwartungstreuer Schätzer für den Erwartungswert, deshalb haben wir mit T(X1, . . . , Xn)einen erwartungstreuen und konsistenten Schätzer fürp.
Aufgabe 44
Wir wollen bei einem Multiple-Choice-Test überprüfen, ob ein Teilnehmer die Antworten durch einfaches Raten ermittelt hat. Dazu beschreiben wir die (richtige) Beantwortung jeder Frage mit einer b(1, p)-verteilten Zufallsvariable und nehmen an, dass die Fragen unabhängig von einander beantwortet wurden. Bei jeder der n= 8 Aufgaben war genau eine der 4angegebenen Antworten richtig. Um nun zwischen den beiden Hypothesen
H0 : p= 0.25 H1 : p >0.25 zu entscheiden, soll der Test
ϕ(x1, . . . , xn) =
(1 falls Pn
i=1xi≥c 0 falls Pn
i=1xi< c
verwendet werden. Bestimmen Sie den kleinsten Parameter c∈ {1, . . . , n} so, dass dies zu einem Test zum Niveau α= 0.05 wird.
Lösung: Für unabhängige b(1, p)-verteilte Zufallsvariablen X1, . . . , Xn ist Pn
i=1Xi eine b(n, p)- verteilte Zufallsvariable. Daher ist die Fehlerwahrscheinlichkeit erster Art dieses Tests gegeben durch:
Pp=0.25[ϕ(X1, . . . , Xn) = 1] =Pp=0.25[
n
X
i=1
Xi ≥c]
=Pp=0.25[
n
X
i=1
Xi=n] +Pp=0.25[
n
X
i=1
Xi=n−1] +. . .+Pp=0.25[
n
X
i=1
Xi =c]
= n
n
·pn·(1−p)0+ n
n−1
·p(n−1)·(1−p)1+. . .+ n
c
·p(n−c)·(1−p)c
3
11. Übung Mathematik und Statistik für Biologie
Für n= 8 und p= 0.25 erhalten wir:
8 8
·0.258·0.750 = 0.000015
8 7
·0.257·0.751 = 0.000366 8
6
·0.256·0.752 = 0.003845
8 5
·0.255·0.753 = 0.023071 8
4
·0.254·0.754 = 0.086517
Also ist c= 5 der kleinste Parameter, so dass Pp=0.25[ϕ(X1, . . . , Xn) = 1]≤α= 0.05.
4