Optimale Versuchsplanung
Die Optimale Versuchsplanung, auch Optimal Experimental Design genannt ist ein Teilgebiet der statistischen Versuchsplanung.
Ziel ist es, Versuchspläne so zu konstruieren, dass z. B. bei vorgegebener Versuchsanzahl die Identifizierung der interessierenden unbekannten Größen im Sinne eines geeigneten Optimalitätskriteriums bestmöglich gelingt.
Besonders fortgeschritten und mathematisch ausgereift ist die optimale Versuchsplanung für lineare Modelle, insbesondere lineare Regressionsmodelle. Für nichtlineare Modelle gibt es ebenfalls Ansätze zur optimalen Versuchsplanung, in einem speziellen Fall siehe z. B. unter modellbasierte Versuchsplanung.
Beispiel: Einfache lineare Regression[Bearbeiten | Quelltext bearbeiten]
Die Varianzen der geschätzten Regressionsparameter (Anstieg und Absolutglied der Regressionsgeraden) hängen von den Werten der unabhängigen Variablen an den Beobachtungsstellen ab. Falls die Werte der unabhängigen Variablen in einem gewissen Versuchsbereich einstellbar sind, sollte man diese so wählen, dass die geschätzten Regressionsparameter möglichst kleine Varianzen haben. Falls man z. B. 10 Versuche machen möchte und die -Werte im Intervall einstellen kann, ist es in der Regel z. B. nicht optimal, diese 10 Versuche gleichmäßig über zu verteilen. Optimal ist 5 Versuche bei und 5 Versuche bei durchzuführen. Dann kann man, was auch anschaulich leicht einsehbar ist, die Gerade „am sichersten“ festlegen.
![{\displaystyle [a,b]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9c4b788fc5c637e26ee98b45f89a5c08c85f7935)
Optimale Versuchsplanung im allgemeinen linearen Regressionsmodell[Bearbeiten | Quelltext bearbeiten]
Die unbekannten Parameter in einem klassischen Modell der linearen Mehrfachregression werden durch die Kleinste-Quadrate-Schätzung erwartungstreu und varianzoptimal geschätzt (siehe Satz von Gauß-Markov). Varianzoptimal heißt hier: mit kleinster Kovarianzmatrix im Sinne der Loewner-Halbordnung. Die Inverse dieser Kovarianzmatrix wird als (Fishersche) Informationsmatrix bezeichnet. Da diese Matrizen von den unabhängigen Variablen an den Beobachtungsstellen abhängen, sind alle gängigen Optimalitätskriterien Funktionale der Informationsmatrix bzw. der Kovarianzmatrix des Kleinste-Quadrate-Schätzers.
Beispiele für Optimalitätskriterien[Bearbeiten | Quelltext bearbeiten]
D-Optimalität: Maximierung der Determinante der Informationsmatrix. Dieses Kriterium führt zur Minimierung des Volumens des Konfidenzellipsoides für die unbekannten Parameter des linearen Regressionsmodells.
A-Optimalität: Minimierung der Spur der inversen Informationsmatrix, führt zur Minimierung der mittleren Varianz der geschätzten Parameter.
E-Optimalität: Maximierung des minimalen Eigenwertes der Informationsmatrix, führt zur Minimierung der maximal möglichen Varianz der Komponenten des geschätzten Parametervektors.
Einige Kriterien beziehen sich auf die Varianz einer Vorhersage im linearen Regressionsmodell, z. B. die
G-Optimalität: führt zur Minimierung der maximal möglichen Varianz der Vorhersage in einem festgelegten Prognosebereich.
I-Optimalität: minimiert die durchschnittliche Vorhersagevarianz im ganzen Versuchsraum (englisch design space).
Theorie der optimalen Versuchsplanung[Bearbeiten | Quelltext bearbeiten]
Man verzichtet zunächst auf die Abhängigkeit von der Versuchsanzahl und betrachtet Versuchspläne als ein normiertes Maß über dem Versuchsbereich. Für Optimalitätskriterien, die Funktionale der Informationsmatrix sind, kann man sich auf diskrete Versuchsplanmaße beschränken, siehe[1] und findet optimale Pläne mit Mitteln der konvexen Optimierung. Eine große Rolle spielen dabei sogenannte Äquivalenzsätze, siehe,[2] die für ein gegebenes Kriterium ein äquivalentes Optimalitätskriterium liefern, das häufig Möglichkeiten zur iterativen Bestimmung näherungsweise optimaler Pläne bereithält. Zur konkreten Anwendung eines optimalen Versuchsplanmaßes für gegebenen Stichprobenumfang hat man die Gewichte dieses (diskreten) Planmaßes auf Vielfache von zu runden.
Historische Entwicklung[Bearbeiten | Quelltext bearbeiten]
Kirstine Smith (1918) benutzte bereits das Kriterion der G-Optimalität, siehe.[3] Von einer Theorie der optimalen Versuchsplanung kann man seit Gustav Elfving (1952) sprechen, siehe,[4] insbesondere aber seit Jack Kiefer (1959).
In Deutschland haben insbesondere Hans Walter Bandemer (1973) (siehe[5]), Friedrich Pukelsheim (1980) (siehe[6]) und H. Dette zur Entwicklung der optimalen Versuchsplanung beigetragen, siehe z. B.[7]
Standardwerke[Bearbeiten | Quelltext bearbeiten]
- V. V. Fedorov: Theory of Optimal Experiments. Academic Press, New York 1972
- S. D. Silvey: Optimal Design. Chapman and Hall, London 1980
- A. Pazman: Foundations of Optimum Experimental Design. Reichel, Dordrecht 1986
- A. C. Atkinson, A. N. Donev: Optimum Experimental Design. Clarendon Press, Oxford 1992
- F. Pukelsheim: Optimal Design of Experiments. Wiley, New York 1993
Deutschsprachige Werke
- H. Bandemer, (mit einem Autorenkollektiv): Theorie und Anwendung der optimalen Versuchsplanung. Band 1, Akademieverlag, Berlin 1977.
- O. Krafft: Lineare statistische Modelle und optimale Versuchspläne. Vandenhoeck & Ruprecht, Göttingen 1978.
- H. Bandemer, W. Näther: Theorie und Anwendung der optimalen Versuchsplanung. Band 2, Akademieverlag, Berlin 1980.
Einzelnachweise[Bearbeiten | Quelltext bearbeiten]
- ↑ J. Kiefer: Optimum experimental design. In: Journal of Royal Statistical Society, Ser. B. 21, 1959, S. 272–319.
- ↑ J. Kiefer: General equivalence theory for optimum designs (approximate theory). In: Annals of Statistic. 2, 1974, S. 849–879.
- ↑ K. Smith: On the standard deviations of adjusted and interpolated values of an observed polynomial function and its constants and the guidance they give towards a proper choice of the distribution of observations. In: Biometrika. 12, 1918, S. 1–85.
- ↑ G. Elfving: Optimum allocation in linear regression theory. In: Annals of Mathematical Statistics. 23, 1952, S. 255–262.
- ↑ Hans Bandemer, Andreas Bellmann, Wolfhart Jung, Klaus Richter: Optimale Versuchsplanung. Akademieverlag Berlin 1973.
- ↑ F. Pukelsheim: On linear regression designs which maximize information. In: Journal of Statistical Planning and Inference. 16, 1980, S. 339–364.
- ↑ H. Dette: A generalization of D- and D1-optimal designs in polynomial regression. In: Annals of Statistics. 18, 1990, S. 1784–1804.