Ähnlichkeitsanalyse

In der Statistik, insbesondere der Multivariaten Statistik, interessiert man sich für die Messung der Ähnlichkeit zwischen verschiedenen Objekten und definiert dazu Ähnlichkeits- und Distanzmaße. Es handelt sich dabei nicht um Maße im mathematischen Sinn, der Begriff bezieht sich ausschließlich auf die Messung einer bestimmten Größe.

In der Regel werden Ähnlichkeitsmaße für nominal oder ordinal skalierte Variablen genutzt und Distanzmaße für metrisch skalierte Variablen (d. h. für Intervall- und Verhältnisskala).

Ähnlichkeitsmaß

Definition

Sei $I=\left\{1,2,\dots ,N\right\}$ eine endliche Menge. Eine Funktion $s\colon I\times I\rightarrow \mathbb {R}$ heißt Ähnlichkeitsmaß oder Ähnlichkeitsfunktion, falls für alle $i,j\in I$ gilt:

• $s(i,j)=s(j,i)$ und
• $s(i,i)\geq s(i,j)$.

Zudem wird oft noch gefordert, dass für alle $i,j\in I$ gilt:

• $s(i,j)\geq 0$ und $s(i,i)=1$.

Die Funktionswerte $s(i,j)$ lassen sich zu einer symmetrischen $N\times N$-Matrix $\left(s(i,j)\right)_{i,j}$ anordnen. Diese Matrix heißt Ähnlichkeitsmatrix. In diesem Kontext wird $s(i,j)$ auch als Ähnlichkeitskoeffizient bezeichnet.

Anwendung in der Bioinformatik

Ähnlichkeitsmatrizen wie PAM oder BLOSUM spielen eine wichtige Rolle beim Sequenzalignment. Ähnliche Proteine, Nukleotide oder Aminosäuren erhalten dabei höhere Scores (d. h. Ähnlichkeitswerte) als unähnliche. Die Ähnlichkeit ist hier durch die chemischen Eigenschaften der Bausteine und ihre Mutationsraten definiert.

Beispiel (AGCT steht für die vier Nukleinbasen Adenin, Guanin, Cytosin und Thymin):

	A	G	C	T
A	10	−1	−3	−4
G	−1	7	−5	−3
C	−3	−5	9	0
T	−4	−3	0	8

Die Moleküle, deren Ähnlichkeit angegeben werden soll, werden in gleicher Reihenfolge spalten- und zeilenweise sortiert. Der Wert $a_{i,j}$ an der Position $(i,j)$ gibt somit an, wie ähnlich das Molekül an der Spaltenposition $i$ dem an der Zeilenposition $j$ ist.

Laut obiger Ähnlichkeitsmatrix sind Cytosin und Tymin (Ähnlichkeits-Score 0) einander ähnlicher als Guanin und Cytosin (Ähnlichkeits-Score -5).

Ähnlichkeitsmaße für binäre Variablen

Für $p$ binäre Variablen und zwei Beobachtungen $i$ und $j$ sei

$n_{00}=\sum _{k=1}^{p}I(x_{ik}=0,x_{jk}=0)$, $n_{01}=\sum _{k=1}^{p}I(x_{ik}=0,x_{jk}=1)$,

$n_{10}=\sum _{k=1}^{p}I(x_{ik}=1,x_{jk}=0)$, $n_{11}=\sum _{k=1}^{p}I(x_{ik}=1,x_{jk}=1)$ und

$p=n_{00}+n_{01}+n_{10}+n_{11}\,$.

Dann kann man folgende Maße definieren:

Ähnlichkeitsmaß	$s(i,j)$
Braun	${\frac {n_{11}}{\max(n_{11}+n_{01},n_{11}+n_{10})}}$
Dice	${\frac {2n_{11}}{n_{01}+n_{10}+2n_{11}}}$
Hamann	${\frac {(n_{00}+n_{11})-(n_{01}+n_{10})}{p}}$
Jaccard (S-Koeffizient)	${\frac {n_{11}}{n_{01}+n_{10}+n_{11}}}$
Kappa	${\frac {1}{1+{\tfrac {p(n_{01}+n_{10})}{2(n_{00}n_{11}-n_{01}n_{10})}}}}$
Kulczynski	${\frac {n_{11}}{n_{01}+n_{10}}}$
Ochiai	${\frac {n_{11}}{\sqrt {(n_{11}+n_{01})(n_{11}+n_{10})}}}$
Phi	${\frac {n_{11}n_{00}-n_{10}n_{01}}{\sqrt {(n_{11}+n_{01})(n_{11}+n_{10})(n_{00}+n_{01})(n_{00}+n_{10})}}}$
Russel Rao	${\frac {n_{11}}{p}}$
Simple Matching (M-Koeffizient)	${\frac {n_{00}+n_{11}}{p}}$
Simpson	${\frac {n_{11}}{\min(n_{11}+n_{01},n_{11}+n_{10})}}$
Sneath	${\frac {n_{11}}{n_{11}+2n_{01}+2n_{10}}}$
Tanimoto (Rogers)	${\frac {n_{00}+n_{11}}{n_{00}+2(n_{01}+n_{10})+n_{11}}}$
Yule	${\frac {n_{00}n_{11}-n_{01}n_{10}}{n_{00}n_{11}+n_{01}n_{10}}}$

Für nicht binäre nominale oder ordinale Variablen definiert man für jede Kategorie der Variablen eine binäre Variable und kann dann die Ähnlichkeitsmaße für binäre Variablen verwenden.

Wahl des Ähnlichkeitsmaßes

Welches Ähnlichkeitsmaß man zur Analyse wählt, hängt von der Problemstellung ab. Es gibt jedoch einige Hinweise, wann sich welches Maß gut eignet in Abhängigkeit von den Eigenschaften der binären Variable:

• Ist die Variable symmetrisch, d. h. beide Kategorien sind gleich wichtig (z. B. Geschlecht), dann ist oft das gleiche Vorhandensein ($n_{11}$) bzw. die gleiche Abwesenheit ($n_{00}$) wichtig für ein Ähnlichkeitsmaß. Dann können Simple Matching, Hamman oder Tanimoto verwendet werden.
• Ist die Variable asymmetrisch, d. h. nur eine Kategorie spielt eine wesentliche Rolle (z. B. Krankheit aufgetreten), dann spielt oft nur das gleiche Auftreten ($n_{11}$) eine Rolle. Dann können Dice, Jaccard, Kulczynskl, Ochiai, Braun, Simpson oder Sneath verwendet werden.
• Kappa, Phi und Yule können sowohl im symmetrischen als auch im asymmetrischen Fall verwendet werden.

Bei der Wahl des Ähnlichkeitmaßes sollten auch Zusammenhänge zwischen den Maßen berücksichtigt werden:

• Dice, Jaccard und Sneath sind monotone Funktionen voneinander:

${\text{Sneath}}\leq {\text{Jaccard}}\leq {\text{Dice}}.$

• Betrachtet man Simpson und Braun, so ist das harmonische Mittel Dice, das arithmetische Mittel Kulczynski und das geometrische Mittel Ochiai. Aus der Ungleichung der Mittelwerte folgt:

${\text{Braun}}\leq {\text{Dice}}\leq {\text{Ochiai}}\leq {\text{Kulczynski}}\leq {\text{Simpson}}.$

• Auch Hamman, Rogers und Simple matching weisen einen Zusammenhang auf.

Distanzmaß

Definition

Sei $I=\left\{1,2,\dots ,N\right\}$ eine endliche Menge. Eine Funktion $d\colon I\times I\rightarrow \mathbb {R}$ heißt Distanzmaß oder Distanzfunktion, falls für alle $i,j\in I$ gilt:

• $d(i,j)=d(j,i)$ sowie
• Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle d(i,j) \geq 0} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle d(i,j) = 0 \Leftrightarrow i=j } .

Die Funktionswerte Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle d(i,j)} lassen sich zu einer symmetrischen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle N \times N} -Matrix Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \left(d(i,j)\right)_{i,j}} anordnen. Diese Matrix heißt Distanzmatrix.

Falls die Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle d} zusätzlich die Dreiecksungleichung erfüllt, ist sie eine Metrik. Häufig wird auch eine Metrik als Distanzfunktion bezeichnet.

Einige Distanzmaße

Für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle p} metrische Variablen und zwei Beobachtungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle i} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle j} kann man folgende Maße definieren:

Distanzmaß	Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle d(i,j)}
Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L_r}	Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \left(\sum_{k=1}^p |x_{ik}-x_{jk}|^r\right)^{1/r}}
Euklidisch Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L_2}	Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sqrt{\sum_{k=1}^p (x_{ik}-x_{jk})^2}}
Pearson	Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sqrt{\sum_{k=1}^p \frac{(x_{ik}-x_{jk})^2}{s_k^2}}} mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle s_k} die Standardabweichung der Variable Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle k}
City-Block Manhattan Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L_1}	Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sum_{k=1}^p |x_{ik}-x_{jk}|}
Gower	Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sum_{k=1}^p \frac{|x_{ik}-x_{jk}|}{r_k}} mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle r_k} die Spannweite der Variable Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle k}
Mahalanobis	Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sqrt{(x_i - x_j)^T S^{-1}(x_i - x_j)}} mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle S} der Stichproben-Kovarianzmatrix der Variablen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle x_i}

Zusammenhang zwischen Ähnlichkeits- und Distanzmaßen

Allgemein kann man ein Distanzmaß aus einem Ähnlichkeitsmaß definieren durch

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle d(i,j)=\sqrt{s(i,i)+s(j,j)-2s(i,j)}} .

Ein so gewonnenes Distanzmaß erfüllt aber im Allgemeinen nicht die Dreiecksungleichung und ist somit keine Metrik.

Siehe auch

• Kosinus-Ähnlichkeit
• Clusteranalyse
• Multidimensionale Skalierung

Literatur

• Joachim Hartung, Bärbel Elpelt: Multivariate Statistik. Lehr- und Handbuch der angewandten Statistik. Oldenbourg Verlag, München 1984, ISBN 3-486-28451-7
• Ludwig Fahrmeir, Alfred Hamerle: Multivariate statistische Verfahren. de Gruyter, Berlin 1984, ISBN 3-11-008509-7