gms | German Medical Science

49. Jahrestagung der Deutschen Gesellschaft für Medizinische Informatik, Biometrie und Epidemiologie (gmds)
19. Jahrestagung der Schweizerischen Gesellschaft für Medizinische Informatik (SGMI)
Jahrestagung 2004 des Arbeitskreises Medizinische Informatik (ÖAKMI)

Deutsche Gesellschaft für Medizinische Informatik, Biometrie und Epidemiologie
Schweizerische Gesellschaft für Medizinische Informatik (SGMI)

26. bis 30.09.2004, Innsbruck/Tirol

Ein global konvergierender EM-Algorithmus für semiparametrische Mischverteilungsmodelle mit Kovariaten

Meeting Abstract (gmds2004)

Search Medline for

Kooperative Versorgung - Vernetzte Forschung - Ubiquitäre Information. 49. Jahrestagung der Deutschen Gesellschaft für Medizinische Informatik, Biometrie und Epidemiologie (gmds), 19. Jahrestagung der Schweizerischen Gesellschaft für Medizinische Informatik (SGMI) und Jahrestagung 2004 des Arbeitskreises Medizinische Informatik (ÖAKMI) der Österreichischen Computer Gesellschaft (OCG) und der Österreichischen Gesellschaft für Biomedizinische Technik (ÖGBMT). Innsbruck, 26.-30.09.2004. Düsseldorf, Köln: German Medical Science; 2004. Doc04gmds115

The electronic version of this article is the complete one and can be found online at: http://www.egms.de/en/meetings/gmds2004/04gmds115.shtml

Published: September 14, 2004

© 2004 Schlattmann.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by-nc-nd/3.0/deed.en). You are free: to Share – to copy, distribute and transmit the work, provided the original author and source are credited.


Outline

Text

Einleitung

Semiparametrische Mischverteilungsmodelle sind ein wichtiges statistisches Verfahren zur Analyse unbeobachteter Heterogenität, das durch die Annahme einer diskreten mischenden Verteilung zugleich die Klassifikation der individuellen Beobachtung erlaubt. Vielfach stellt sich die Frage, ob beobachtete Kovariate einen Teil der unbeobachteten Heterogenität aufklären können, wie z.B. in der Meta-Analyse klinischer oder epidemiologischer Studien [1].

Methoden

Das Standardverfahren zur Schätzung der Parameter semiparametrischer Mischverteilungsmodelle ist der EM-Algorithmus [2]. Hier wird im E(wartungswert)-Schritt die a-posteriori Wahrscheinlichkeit der Zugehörigkeit der jeweiligen Beobachtung zur latenten Klasse geschätzt. Im M(aximierungs)-Schritt wird der Anteilswert der jeweiligen Klasse und die Parameter der jeweiligen Subpopulation geschätzt. Dies läßt sich leicht auf Kovariate erweitern, indem im M-Schritt die Parameter des populationsspezifischen Regressionsmodells geschätzt werden.. Eine Schwierigkeit besteht darin, daß der EM-Algorithmus oft nur zu einem lokalen Maximum konvergiert. Für den Fall ohne Kovariate schlägt Böhning [3] einen global konvergenten EM-Algorithmus vor, der auf dem Maximum des Gradienten der Richtungsableitung des Likelihoods beruht und den Parameterwert mit minimalem Gradienten gegen denjenigen mit maximalem Gradienten austauscht Wert. Im Fall ohne Kovariate ist das Maximum leicht zu finden. Im Fall mehrerer Kovariate sind deterministische Optimierungsalgorithmen nur bedingt anwendbar. Deshalb wir hier ein global konvergenter EM-Algorithmus vorgestellt, der das Maximum des Gradienten mit Hilfe von Simulated Annealing sucht. Auf diese Weise wird ein global konvergierender EM-Algorithmus für semiparametrische Mischverteilungen mit Kovariaten konstruiert.

Ergebnisse

Die Ergebnisse des Standard-EM und des global konvergenten Algorithmus werden anhand einer Meta-Analyse von epidemiologischen Studien zur Hormonersatztherapie vorgestellt.

Diskussion

Der vorgestellte, global konvergente Algorithmus bietet in der Auswertung praktischer Daten den Vorteil, daß die Verwendung verschiedener Kombinationen von Startwerten entfällt, um Konvergenz zu lokalen Minima aufzuspüren.


Literatur

1.
Blettner M, Schlattmann Meta-Analysis. In: Pigeot I, Ahrens W Hrsg., Handbook of Epidemiology Berlin, Springer, 2004. (im Druck)
2.
McLachlan G J.; Krishnan T. The EM algorithm and extensions, Wiley Series in Probability and Mathematical Statistics. Applied Probability and Statistics. New York, NY: John Wiley & Sons; 1997
3.
Böhning D. The EM-Algorithm with Gradient Function Update for Discrete Mixtures with Fixed (Known) Number of Components. Statistics and Computing 2003; 13: 257-265