Nichtnegative Matrixfaktorisierung
# Nichtnegative Matrixfaktorisierung

Die Nichtnegative Matrixfaktorisierung (NMF) ist ein Verfahren zur Dimensionsreduktion von Daten. Eine Matrix mit Einträgen in den nichtnegativen reellen Zahlen wird dabei linear in Faktoren vom Rang 1 zerlegt. Durch spezielle Algorithmen kann eine Zerlegung gefunden werden, bei der auch die einzelnen Faktoren nichtnegativ sind. Diese Forderung führt in vielen Fällen zu Zerlegungen, die leichter zu interpretieren sind und die Daten als Summe von klar getrennten Komponenten darstellen. Die NMF wird seit ihrer Einführung 1999 in vielen Gebieten der Wissenschaft angewandt.

## Geschichte

In der Chemometrik ist die nichtnegative Matrixfaktorisierung bereits seit den 1970er Jahren bekannt, wobei sie allerdings nicht als Faktorisierung von Matrizen beschrieben wurde. In den 90er Jahren veröffentlichten einige finnische Wissenschaftler Arbeiten zur positiven Matrixfaktorisierung.
Den Durchbruch zur weiteren Verbreitung schaffte die NMF, als die amerikanischen Neurowissenschaftler Daniel D. Lee und Sebastian Seung 1999 in einem Artikel in der Fachzeitschrift Nature ihre Eigenschaften untersuchten und den einfachen multiplikativen Algorithmus zu ihrer Berechnung angaben. Die Faktorisierung erfreut sich seitdem großer Beliebtheit und es sind mathematische Analysen, neue Berechnungsalgorithmen und abgewandelte Faktorisierungen entstanden.

## Definition

Sei $X$ eine beliebige Matrix der Dimension $m\times n$ mit nichtnegativen Einträgen. Weiterhin sei die Anzahl der Komponenten $k\in \mathbb {N}$ gegeben, die kleiner als $m$ und als $n$ sei. Die nichtnegative Matrixfaktorisierung besteht dann aus einer Matrix $W$ der Dimension $m\times k$ und einer Matrix $H$ der Dimension $k\times n$, die beide nichtnegative Einträge besitzen und gemeinsam die Frobeniusnorm der Differenz
    $\lVert X-WH\rVert _{\text{fro}}$
minimieren.
Diese Definition entspricht, abgesehen von der Forderung der Nichtnegativität, genau der der Hauptkomponentenanalyse. Die Faktorisierung ist durch das Minimierungsproblem nicht eindeutig bestimmt. Beispielsweise sind für eine Permutationsmatrix $\Pi \in S_{k}$ die Matrizen $W'=W\Pi$ und $H'=\Pi ^{-1}H$ ebenfalls Minimierer von $\lVert X-W'H'\rVert _{\text{fro}}$, sodass die Ordnung der Faktoren also völlig unbestimmt ist. Auch die Skalierung der Faktoren kann variieren.

## Beispiel

Das folgende Beispiel zeigt, wie man mithilfe der NMF zum Beispiel die Aktivität von Neuronen analysieren kann. Als Ausgangsmatrix $X$ wählen wir eine Matrix, die durch Helligkeitswerte in dem rechts stehenden Video gegeben ist. Das Video zeigt eine optogenetische Aufnahme aus dem Gehirn einer Maus; die Maus wurde also gentechnisch so verändert, dass die Aktivierung der Kalziumkanäle der Neuronen einen Lichtimpuls verursacht.
Die Zeilen der Matrix sollen nun die einzelnen Zeitpunkte darstellen, die Spalten die einzelnen Bildpunkte (Pixel) des Videos. Da das Video aus 400 Einzelbildern besteht, die jeweils $512\times 512$ Pixel groß sind, ergibt sich so eine Matrix der Dimension $400\times 262\,144$.
Als Komponentenanzahl wählen wir $k=6$. Die nichtnegative Matrix $W$ der Dimension $400\times 6$ gibt dann den Zeit-Faktor an, die Matrix $H$ der Dimension $6\times 262\,144$ beschreibt die räumliche Verteilung. Die Zeilen der Matrix $H$ können dann wieder zu Bildern aufgefaltet werden. Es ergeben sich die folgenden Faktoren:
* Ergebnis der NMF sind sechs Komponenten, die jeweils einen zeitlichen und einen räumlichen Faktor haben. Hier sind drei Komponenten gezeigt, die jeweils Gruppen von gemeinsam feuernden Neuronen darzustellen scheinen.
* Die räumlichen Faktoren können auch räumlich überlagert dargestellt werden (hier als Farbkanäle eines RGB-Bildes), um die Neuronen visuell in funktional ähnliche Klassen einzuteilen
Die Matrixfaktorisierung kann somit aufzeigen, dass bestimmte Gruppen von Neuronen gemeinsam feuern. Dass die Gewichtungen nichtnegativ sind, ist hier klar von Vorteil, denn die Faktoren sind leichter interpretierbar, wenn keine Neuronen negativ gewichtet werden.

## Berechnung
### Multiplikative Methode

Lee und Seung gaben folgende multiplikative Update-Regel zur Bestimmung von $W$ und $H$ an:
    $H\leftarrow H\odot {\frac {W^{T}X}{W^{T}WH}}$
    $W\leftarrow W\odot {\frac {XH^{T}}{WHH^{T}}}$
Hierbei bezeichnet $\odot$ das Hadamard-Produkt, also die elementweise Multiplikation, und auch bei den Brüchen sollen Zähler und Nenner in jedem Eintrag dividiert werden. Diese Update-Regeln kann aus dem Gradientenabstieg unter Hinzunahme spezieller Vorfaktoren hergeleitet werden. Der Vorteil eines multiplikativen gegenüber einem additiven Gradientenabstieg ist, dass die Faktoren automatisch das Vorzeichen beibehalten. Einer der Nachteile ist, dass Elemente von $W$ oder $H$, die einmal null sind, nicht mehr positiv werden können.

### Alternierende Least-Sqaures-Näherung

Ein anderes Verfahren zur Bestimmung von $W$ und $H$ ist die Methode der alternierenden Least-Squares-Näherungen. Während das Minimierungsproblem in $W$ und $H$ gemeinsam nicht konvex ist, ist das Minimieren von $W$ für gegebenes $H$ und andersherum ein konvexes Problem und damit leichter zu lösen. In der Tat handelt es sich bei der Minimierung von $\lVert X-WH\rVert _{\text{fro}}$ für festes $H$ einfach um ein least-squares-Regressionsproblem (kleinste-Quadrate-Regression), bis auf die Einschränkung, dass $W$ nichtnegativ sein muss. Für das Regressionsproblem mit nichtnegativen Variablen (NNLS, nonnegative least squares) gibt es zwar keine einfache analytische Darstellung, aber durchdachte Algorithmen, mit deren Hilfe das Optimierungsproblem dann gelöst werden kann. Die Kostenfunktion wird in jedem Schritt garantiert verringert. Die meisten Software-Pakete für die Nichtnegative Matrixfaktorisierung benutzen dieses Verfahren.

### Initialisierung

Da bei den iterativen Optimierungsverfahren nicht garantiert ist, dass sie ein globales Optimum finden, kann die Initialisierung der Faktoren $W$ und $H$ eine wichtige Rolle für das Endergebnis spielen. Statt einer zufälligen Initialisierung hat sich unter anderem die Initialisierung auf Grundlage einer zuvor ausgeführten Singulärwertzerlegung bewährt.