imported>Mendli: sqrt(64) und sqrt(16) korrigiert

2023-06-12T10:32:07Z

sqrt(64) und sqrt(16) korrigiert

Neue Seite

[[Datei:TypeTwoErrorAsFunctionOfMu de.png|mini|hochkant=2|Abhängigkeit des Risikos 2. Art von der wahren Lage des Gegenparameters µ1 (im nebenstehenden Text θ1 genannt) bei einem ein- sowie zweiseitigen Hypothesentest.]]
In der [[Statistik]] ist die '''Operationscharakteristik''', auch '''OC-Kurve''' (''OC'': {{enS}} für ''operating characteristic'') oder '''OC-Funktion'''<ref name="Roenz1994">Bernd Rönz, Hans G. Strohe (1994), ''Lexikon Statistik'', Gabler Verlag, S. 268</ref> genannt, ein Konzept aus der Theorie [[statistischer Test]]s, mit dem ein funktionaler Zusammenhang zwischen der Wahrscheinlichkeit eines Fehlers 2. Art und der tatsächlichen Lage des unbekannten Parameters <math>\theta</math> einer [[Verteilungsfunktion]] <math>F(x|\theta)</math> hergestellt wird.

== Definition ==
[[Datei:PowerFunctionAndSampleSize-one-sided de.png|mini|hochkant=1.35|Einfluss des Stichprobenumfangs auf die [[Gütefunktion]] bzw. Teststärke eines rechts-[[Einseitiger Test|einseitigen Tests]] ]]
[[Datei:PowerFunctionAndSampleSize-two-sided de.png|mini|hochkant=1.35|Einfluss des Stichprobenumfangs auf die Gütefunktion bzw. Teststärke eines zweiseitigen Tests]]
Gegeben ist eine [[Zufallsvariable]] <math>X</math> mit einer [[Verteilungsfunktion]] <math>F(x|\theta)</math>, die von einem unbekannten [[Parameter (Statistik)|Parameter]] <math>\theta</math> abhängt. Für die Schätzung des Parameters werden <math>n</math> Beobachtungen der Zufallsvariablen gemacht. Der Parameter kann dann durch eine [[Schätzfunktion]]

:<math>\hat \theta = t\left(X_1, X_2, \dotsc, X_n\right)</math>

geschätzt werden. Es soll eine Vermutung bezüglich des wahren, unbekannten Parameters statistisch überprüft werden. Es wird also eine [[Nullhypothese|Hypothese]] bezüglich dieses Parameters aufgestellt, die sogenannte Nullhypothese <math>H_0</math>. Man geht nun davon aus, dass bei Wahrheit der Nullhypothese der Schätzwert <math>\hat \theta</math> in der Nähe des wahren Parameters <math>\theta</math> liegen müsste, und lehnt <math>H_0</math> ab, wenn die Distanz zu groß ist, wenn also <math>\hat \theta</math> in den Ablehnungsbereich des Tests fällt. Der Ablehnungsbereich <math>AB</math> wird so festgelegt, dass von allen Stichproben selbst dann, wenn <math>H_0</math> wahr wäre, ein Anteil von <math>\alpha</math> (häufig wählt man <math>\alpha=5\ \%</math>) abgelehnt würde.

Man kann im Hypothesentest zwei Arten von Fehlern begehen:
* Man lehnt <math>H_0</math> ab, obwohl <math>\theta_0</math> der wahre Parameter ist. Es handelt sich also um einen Fehler, den sogenannten α-Fehler oder [[Fehler 1. und 2. Art#Fehler 1. Art|Fehler 1. Art]].
* Man lehnt <math>H_0</math> nicht ab, obwohl ein anderer Parameter <math>\theta_1 </math> der wahre Parameter ist. Das ist der β-Fehler oder [[Fehler 1. und 2. Art#Fehler 1. Art|Fehler 2. Art]].

<math>\alpha</math> wird vor der Testprozedur festgelegt, <math>\beta</math> dagegen hängt vom wahren Parameter <math>\theta_1</math> ab, der in der Regel unbekannt ist. Man kann für die Risikoabschätzung einer falschen Entscheidung die β-Fehler für verschiedene alternative Parameterwerte <math>\theta_1</math> berechnen. Der β-Fehler für einen alternativen Parameter <math>\theta_1</math> berechnet sich als [[Wahrscheinlichkeit]], dass <math>\hat \theta </math> in den Nichtablehnungsbereich <math>NAB</math> der Nullhypothese <math>H_0</math> fällt, wenn bzw. obwohl in Wahrheit <math>\theta_1</math> die Verteilung von <math>\hat \theta</math> regiert:

:<math>\beta = \operatorname{P}\left(\hat \theta \isin NAB \mid \theta_1\right)</math>.

<math>\beta</math> hängt also von <math>\theta_1</math> ab und kann daher auch als Funktion des alternativen Parameters <math>\theta_1</math> dargestellt werden:

:<math>\beta = f(\theta_1)</math>.

Diese Funktion wird als ''Operationscharakteristik'', häufig auch <math>OC(\theta_1)</math> geschrieben, bezeichnet. Die Gegenwahrscheinlichkeit zu <math>\beta</math> ist die Wahrscheinlichkeit, dass <math>H_0</math> abgelehnt und dafür <math>H_1</math> akzeptiert wird, wenn <math>\theta_1</math> der wahre Parameter ist. Hier ist die Ablehnung von <math>H_0</math> zu Gunsten von <math>H_1</math> also erwünscht, weshalb die entsprechende Funktion <math>\gamma(\theta_1) = 1 - OC(\theta_1)</math> auch ''[[Gütefunktion]]'' (und ihr Funktionswert für gegebenes <math>\theta_1</math> ''[[Trennschärfe eines Tests|Trennschärfe]]'' oder ''Teststärke'') genannt wird.

Gütefunktion und Operationscharakteristik stellen damit beide vollständige Charakterisierungen des zugehörigen Tests dar. Man erkennt an ihnen bspw., ob der Test mit wachsender Beobachtungszahl immer besser wird (''Konsistenz'') und ob die Wahrscheinlichkeit, <math>H_0</math> abzulehnen, größer ist, wenn <math>H_1</math> zutrifft, als wenn <math>H_0</math> zutrifft (''Unverfälschtheit'').

== Beispiel ==
[[Datei:BetaFehler.jpg|300px|mini|β-Fehler: Die rote Normalverteilungskurve gibt an, wie der Stichprobenmittelwert ''X'' verteilt wäre, wenn μ = 260 g ist. Die rote Fläche repräsentiert den [[Risiko 1. Art|α-Fehler]] von 5 %. Die blaue Kurve zeigt die Verteilung von ''X'', wenn μ in Wahrheit 255 g betrüge. Die blaue Fläche ist dann die Wahrscheinlichkeit, dass ''X'' ≥ 256,7 ist und ''H''0 damit nicht abgelehnt wird, obwohl die Forellen im Durchschnitt untergewichtig sind. Entsprechendes gilt für die grüne Kurve, bei der das wahre Durchschnittsgewicht der Forellen sogar nur noch 252 g beträgt – wie zu sehen, ist das [[Risiko 2. Art]], sie dennoch als normalgewichtig einzustufen, nun wesentlich kleiner.]]
[[Datei:Operationscharakteristik.png|300px|mini|Operationscharakteristik: Der Ordinatenwert der Grafik gibt den β-Fehler in Abhängigkeit vom unbekannten Parameter μ1 an. Für μ = 260 ist der Wert 0,95, also gerade 1 - α.]]

Ein [[Forelle]]nzüchter liefert seinem Großabnehmer Forellen, die im Durchschnitt mindestens 260 g wiegen sollen. Bei Lieferung wird getestet, ob das Durchschnittsgewicht mindestens 260 Gramm beträgt. Wird die Hypothese abgelehnt, wird die Lieferung beanstandet. Es sei bekannt, dass das Gewicht <math>X</math> der Forellen [[Normalverteilung|normalverteilt]] ist mit der [[Varianz (Stochastik)|Varianz]] <math>\sigma^2 = 64\ \mathrm g^2</math> und einem unbekannten [[Erwartungswert]] <math>\mu</math>. Es werden in einer [[Stichprobe]] <math>n = 16</math> Forellen gewogen, wobei die <math>i</math>-te Forelle <math>x_i\ \mathrm g</math> wiegt. Das Durchschnittsgewicht

:<math>\overline x = \frac{1}{n} \sum_{i=1}^n x_i</math>

dieser Forellen wird ermittelt. Da der [[Stichprobenmittel|Mittelwert]] bei jedem Versuch anders ausfällt, ist diese Größe ebenfalls eine Zufallsvariable <math>\overline X</math> und normalverteilt mit den Parametern

:<math>\mu \;</math> und <math>\;\operatorname{Var}\left(\overline X\right) = \frac {\sigma^2}{n}</math>.

Die Hypothesen lauten nun <math>H_0\colon\mu \ge \mu_0 = 260\;</math> und <math>\;H_1\colon\mu < 260</math>.

Soll der Fehler erster Art beispielsweise <math>\alpha = 0{,}05</math> betragen, ergibt sich der [[Kritischer Wert (Statistik)|kritische Wert]] für die Prüfgröße <math>\overline X</math> als

:<math>\mu_0 - z(1 - \alpha) \cdot \frac {\sigma}{\sqrt n} = 260 - 1{,}65 \cdot \frac {8}{4} = 256{,}7</math>

mit <math>z(1 - \alpha)</math> als <math>(1-\alpha)</math>-Quantil der Standardnormalverteilung.

<math>H_0</math> wird also abgelehnt, wenn <math>\overline x < 256{,}7</math> ist, der Ablehnungsbereich ist <math>(-\infty; 256{,}7)</math>. Ist jetzt tatsächlich <math>\mu = \mu_0</math> wahr, würde in 5 % aller Stichproben <math>\overline x</math> in den Ablehnungsbereich fallen, es würde die Lieferung zu Unrecht zurückgeschickt werden, was dem α-Fehler entspricht.

Es kann aber beispielsweise auch vorkommen, dass das Durchschnittsgewicht in Wahrheit <math>\mu = 255\ \mathrm g</math> beträgt, dass aber zufällig <math>\overline x > 256{,}7\ \mathrm g</math> ist. Das ist der β-Fehler für <math>\mu = 255\ \mathrm g</math>. Die Prüfgröße <math>\overline X</math> ist nun bei unveränderter Varianz in Wahrheit normalverteilt wie

:<math>\overline X \sim \mathcal{N}(255;2)</math>.

Die Wahrscheinlichkeit, dass die Nullhypothese nicht abgelehnt wird, ist dann

:<math>P(\overline X \ge 256{,}7 \mid \mu_1 = 255)</math>

und berechnet sich mit Hilfe der Normalverteilung als

:<math>1-\Phi\left(256{,}7\mid 255;\frac{\sqrt{64}}{\sqrt{16}}\right) = 1 - \Phi_z \left(\frac{256{,}7 - 255}{2}\right)</math>
:<math> = 1-0{,}8023 = 0{,}1977</math>,

wobei <math>\Phi(256{,}7\mid 255;2)</math> der Wert der Normalverteilungsfunktion mit den Parametern 255 und 2 an der Stelle 256,7 ist und <math>\Phi_z</math> der entsprechende Wert der Standardnormalverteilung. Es würde also in ca. 20 % aller Stichproben die Lieferung akzeptiert werden, obwohl die Forellen im Durchschnitt untergewichtig sind. Beträgt dagegen in Wahrheit <math>\mu = 252</math>, ergibt sich der β-Fehler als

:<math>1-\Phi(256{,}7\mid 252;2) = 0{,}0094</math>;

hier ist die Gefahr einer falschen Entscheidung nur noch sehr gering. Die Grafik der Operationscharakteristik zeigt, wie mit wachsender Entfernung von <math>\mu_0</math> der β-Fehler sinkt. Man ist bestrebt, möglichst schnell in den Bereich eines kleinen β-Fehlers zu kommen. Mit der Erhöhung des Stichprobenumfangs kann man den β-Fehler reduzieren. Einen Test mit kleinem β-Fehler nennt man auch ''trennscharf'', weil hier die Verteilungen stark getrennt sind.

== Siehe auch ==
* [[ROC-Kurve]]

== Literatur ==
* Hartung, Joachim/Elpelt, Bärbel/Klösener, Karl-Heinz: ''Statistik – Lehr- und Handbuch der angewandten Statistik.'' 9., durchges. Aufl., Oldenbourg, München 1993, insbesondere Seite 135ff und 381ff.

== Weblinks ==
{{Commonscat|Operating characteristics}}

== Einzelnachweise ==
<references />

[[Kategorie:Testtheorie]]

Operationscharakteristik - Versionsgeschichte

imported>Mendli: sqrt(64) und sqrt(16) korrigiert