Überdispersion

In der Statistik ist Überdispersion (von lateinisch dispersio = „Zerstreuung“; gelegentlich auch Überstreuung oder Hyperdispersion genannt) ein Phänomen, das oft bei der Modellierung von Zähldaten auftritt. Man spricht von Überdispersion, wenn die empirische Varianz in den Daten größer ist als die vom Modell (z. B. Binomialmodell oder Poisson-Modell) angenommene erwartete „theoretische Varianz“. Die tatsächlich gemessene Variation in den Daten übersteigt also die theoretisch erwartete Variation. Das Diagnostizieren von Überdispersion ist Thematik der Regressionsdiagnostik. Bei Vorliegen eines Poisson-Modells stellt die Anwesenheit von Überdispersion die häufigste in der Praxis auftretende Verletzung der Modellannahmen dieses Modells dar. In der Theorie gilt für eine Poisson-verteilte Zielgröße, dass Varianz und Erwartungswert gleich sind. In der Praxis übersteigt die empirisch beobachtete Varianz jedoch häufig den Erwartungswert. Die Zähldaten streuen also in einem größeren Maße um den Erwartungswert, als durch das Poisson-Modell erwartet wird.

Das Gegenstück der Überdispersion, bei der in der Praxis weniger Variation vorliegt als durch das Modell zu erwarten ist, wird Unterdispersion genannt. Sie tritt allerdings weniger häufig auf als Überdispersion. Die Anwesenheit von Überdispersion hat eine Reihe von negativen Konsequenzen für die Validität eines statistischen Modells und den daraus abgeleiteten Schlussfolgerungen. Beispielsweise ist es möglich, dass das Vorhandensein von Überdispersion die Parameterschätzer verzerrt. Die Hauptgründe für Überdispersion sind unbeobachtete Heterogenität und eine positive Korrelationen zwischen den individuellen Beobachtungen der binären Antwortvariablen.

Überdispersion findet in einer Reihe von biologischen Wissenschaften, wie der Parasitologie und Ökologie und in der Epidemiologie Anwendung.

Ursachen

Überdispersion tritt auf, wenn der Erwartungswert (selbst wenn alle erklärenden Variablen festgesetzt sind) eine gewisse inhärente Variabilität beibehält (sogenannte unbeobachtete Heterogenität). Des Weiteren tritt Überdispersion bei einer positiven Korrelation zwischen den individuellen Beobachtungen der binären Antwortvariablen (eine abhängige Variable, die nur zwei mögliche Werte annehmen kann) auf. Überdispersion tritt auch auf, wenn die Zähldaten in Clustern auftreten oder sich auf irgendeine Weise gemeinsam beeinflussen. Als Konsequenz sind die zugrundeliegenden Ereignisse positiv korreliert und eine Überdispersion der Zähldaten ist die Folge. Auch können Ausreißer oder eine Fehl­spezifikation des Regressionsmodells zu Überdispersion führen.

Auswirkungen

Für verallgemeinerte lineare Modelle hat das Auftreten von Überdispersion schwerwiegende Konsequenzen. Es bedeutet, dass die vom verallgemeinerten linearen Modell ausgegebenen Standardfehler unterschätzt werden und Tests bzgl. der erklärenden Variablen im Allgemeinen signifikanter erscheinen als die Daten nahelegen. Die Konsequenz sind übermäßig komplexe Modelle, was wiederum zu weiteren Problemen führt.

Aufdecken von Überdispersion

Überdispersion kann durch einen Anpassungstest festgestellt werden. Wenn die Residuen­devianz und die Pearson-Chi-Quadrat-Statistik deutlich größer sind als die residualen Freiheitsgrade, ist entweder das angepasste Modell nicht adäquat (eine relevante Variable wurde ausgelassen) oder es liegt Überdispersion vor. Bleibt auch nach dem Anpassen eines Modells mit der stärksten Erklärungskraft (gesättigtes Modell) und nach dem Eliminieren von Ausreißern ein Anpassungsmangel bestehen, wäre die Anwesenheit von Überdispersion eine alternative Erklärung.

Theoretische Beispiele

Poisson-Modell

Ein Charakteristikum der Poisson-Verteilung ist, dass der Parameter $\lambda$ der Verteilung zugleich der Erwartungswert und Varianz darstellt (siehe Poisson-Verteilung#Eigenschaften). Daher gilt für ein Poisson-Modell, bei dem eine Poisson-Verteilte Antwortvariable (Zielgröße) $Y_{i}\;\sim \;{\mathcal {P}}(\lambda _{i})$ vorliegt:

$\operatorname {E} (Y_{i})=\operatorname {Var} (Y_{i})=\lambda _{i}$.

Aus ähnlichen Gründen wie bei Binomialdaten wird bei Anwendungen der Poisson-Regression häufig eine signifikant höhere empirische Varianz beobachtet. Die einfachste Möglichkeit einer größeren Variabilität als erwartet Rechnung zu tragen, ist die Einführung eines skalaren Faktors $\phi >1$, der Überdispersionsparameter (auch Dispersionsfaktor) genannt wird. In anderen wissenschaftlichen Disziplinen gibt es eine abweichende Notation; beispielsweise wird in der Parasitologie und Epidemiologie der Überdispersionsparameter mit $\kappa$ statt mit $\phi$ bezeichnet. Mithilfe dieses unbekannten Parameters (den es daher zu schätzen gilt) modifiziert man die Varianzformel wie folgt:

$\operatorname {Var} (y_{i})=\phi \lambda _{i}$.

Durch die Wahl $\phi >1$ wird ein Modell erzeugt, bei dem die Varianz größer als der Erwartungswert ist. Man spricht bei dieser Wahl von $\phi$ von einem poissonverteilten verallgemeinerten linearen Modell mit Überdispersion, obwohl genau genommen keine Poisson-Verteilung mehr vorliegt. Für die Poisson-Verteilung in ihrer ursprünglichen Form, bei der die Varianz dem Erwartungswert entspricht, gilt $\phi =1$. Man spricht in diesem Fall auch von Äquidispersion. Wählt man hingegen $\phi <1$ erhält man ein Modell mit Unterdispersion.

Schätzung des Überdispersionsparameters

Mittlere Pearson-Chi-Quadrat-Statistik

Falls der Überdispersionseffekt signifikant ist, ist es sinnvoll den Überdispersionsparameter zu schätzen. Es ist möglich den Überdispersionsparameter mit der mittleren Pearson-Chi-Quadrat-Statistik zu schätzen. Die Pearson-Chi-Quadrat-Statistik, die sich aus der Quadratsumme der Pearson-Residuen ergibt, ist gegeben durch

$X^{2}:=\sum _{i=1}^{n}{\frac {w_{i}(y_{i}-{\hat {\mu }}_{i})^{2}}{\operatorname {V} ({\hat {\mu }}_{i})}}\;{\stackrel {a}{\sim }}\;\phi \chi ^{2}(n-p)$.

Hierbei sind die $w_{i}$ bekannte spezifizierte Gewichte, ${\hat {\mu }}_{i}$ ist der geschätzte Erwartungswert für Beobachtung $i$ und $\operatorname {V} (\cdot )$ ist die Varianzfunktion. Der Pearson-Schätzer für $\phi$ gegeben ist durch die mittlere Pearson-Chi-Quadrat-Statistik ${\hat {\phi }}:=X^{2}/(n-p)$. Diese ist analog zum erwartungstreuen Schätzer der Varianz der Störgrößen definiert, außer dass hier die Pearson-Chi-Quadrat-Statistik die Residuenquadratsumme ersetzt. Sowohl die Residuendevianz als auch die Pearson-Chi-Quadrat-Statistik vergleichen die Anpassung eines Modells relativ zu einem gesättigten Modell und sind asymptotisch Fehler beim Parsen (Konvertierungsfehler. Der Server („https://wikimedia.org/api/rest_“) hat berichtet: „Cannot get mml. Server problem.“): {\displaystyle \phi \cdot \chi ^{2}(n-p)} -verteilt. Hierbei ist $\chi ^{2}(n-p)$ die Chi-Quadrat-Verteilung mit $(n-p)$ Freiheitsgraden.

Mittlere Residuendevianz

Eine weitere Maßzahl, die zur Schätzung des Überdispersionsparameters herangezogen werden kann, ist die mittlere Residuen­devianz. Die Residuendevianz ist eine Maßzahl für die Variabilität der 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} Beobachtungen, nachdem das Modell angepasst wurde. Sie ist definiert durch

$D(y,{\hat {\mu }}):=\sum _{i=1}^{n}w_{i}d(y_{i},{\hat {\mu }}_{i})\;{\stackrel {a}{\sim }}\;\phi \chi ^{2}(n-p)$.

Sie folgt ebenso wie die Pearson-Chi-Quadrat-Statistik einer 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 \phi \cdot \chi^2(n-p)} -Verteilung 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 (n-p)} Freiheitsgraden. Für lineare Regressionsmodelle entspricht sie gerade der Residuenquadratsumme. Wie gewöhnlich ist der Maximum-Likelihood-Schätzer 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 \tilde{\phi} = D(y, \hat{\mu})/ n} nicht erwartungstreue für $\phi$. Daher verwendet man als erwartungstreuen Schätzer für den Überdispersionsparameter 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 \phi} die mittlere Residuendevianz 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 \hat{\phi} = D(y, \hat{\mu}) /(n-p)} .

Unterschiede in der Terminologie zwischen den Disziplinen

Um zu betonen, dass in der Biologie zumeist klumpenartige Anordnungen (sogenannte Cluster) beschrieben werden, wird in der Biologie meist die Bezeichnung „Aggregation“ der Bezeichnung „Überdispersion“ vorgezogen und synonym zu ihr verwendet. Statt davon zu sprechen, dass die „Verteilung Überdispersion aufweist“, spricht man daher in biologischen Anwendungen (wie z. B. der Ökologie und der Parasitologie) oft von einer aggregativen Verteilung (auch geklumpte Verteilung oder gehäufte Verteilung genannt) der untersuchten Zähldaten. Beispielsweise spricht man davon, dass Parasiten eine „aggregative Verteilung“ aufweisen.

Anwendungen in der Biostatistik

Anwendung in der Parasitologie

Eine Ursache dafür, dass eine sehr kleine Anzahl von Wirten eine große Anzahl von Parasiten trägt ist, dass bereits befallene Wirte erfolgreicher parasitiert werden können, da sie bereits geschwächt sind.

Überdispersion charakterisiert in der Parasitologie ein Phänomen der Aggregation einer Mehrheit der Parasiten in einer Minderheit der Wirts­population. Somit hat die Mehrheit der Wirte keine oder nur wenige Parasiten. Eine sehr kleine Anzahl von Wirten trägt jedoch eine große Anzahl von Parasiten. Eine starke Überdispersion bzw. Aggregation lässt sich in den Daten für Nematoden­infektionen bei Teichfröschen feststellen. In diesem Fall wurden 70 % der Parasiten in nur 4 % der Wirte registriert, während 88 % der Wirte nicht infiziert waren und 8 % leichte Infektionen hatten.

Die Parasitenaggregation in Bezug auf Wirte ist ein charakteristisches Merkmal vielzelliger Parasitenpopulationen. Viele Parasiten, die durch direkten Kontakt übertragen werden (solche, die keinen Krankheitsüberträger verwenden), z. B. Ruderfußkrebse, Walläuse, Läuse, Milben, Hakensaugwürmer, viele Fadenwürmer, Pilze und viele Taxa von Protisten, sind nahezu ausnahmslos durch eine aggregative Verteilung charakterisiert, bei der die meisten einzelnen Wirte nur wenige oder keine Parasiten beherbergen.

In der Theorie gibt es viele Erklärungen dafür, warum Parasitenpopulationen meist aggregativ verteilt sind. Zunächst ist es naheliegend, auch dass die Verteilung von Parasitenpopulationen nicht zufällig, sondern aggregativ ist, wenn schon die Übertragungsstadien nicht zufällig sind. Zudem können Wirte, die bereits befallen sind erfolgreicher parasitiert werden, da sie durch eine vorausgehende Parasitierung bereits geschwächt sind. Oft kommt erschwerend hinzu, dass sich Wirtsindividuen je nach Jahreszeit oder Alter erfolgreicher parasitieren lassen. Auch könnte sich die Befallstärke bei einem bereits parasitierten Wirtsindividuum erhöhen, falls sich die Parasiten ungeschlechtlich oder eingeschlechtlich vermehren.

Anwendung in der Ökologie

Auch in der Ökologie werden häufig „klumpenartige“ Verteilungsmuster modelliert. Wenn ein Verteilungsmuster aus zufällig verstreuten „Klumpen“ (Clustern) besteht, liegt ein „klumpenartiges Verteilungsmuster“ vor. Bei einem solchen sind die Untersuchungseinheiten zu Haufen aggregiert. Eine Auszählung der Objekte würde in diesem Fall Untersuchungsgebiete mit sehr vielen Objekten und andere mit recht wenigen oder überhaupt keinem Objekt ergeben. Die Besetzungszahlen streuen also stark von einem Untersuchungsgebiet zum anderen. Ein Beispiel wäre die Verbreitung der Amerikanerkrähe, die eine „aggregative Verteilung“ aufweist. Sie ist in Nordamerika sehr weit verbreitet, weist jedoch eine hohe Abundanz an einzelnen Aggregationspunkten auf und hat somit ein geklumptes Auftreten.

Auch hier gibt es vielfältige Ursachen für das geklumpte Auftreten von Arten. Beispielsweise können bewohnbare Lebensräume oder gewisse Nährstoffe inselartig über ein größeres Gebiet verbreitet sein. Zudem bilden viele Tierarten soziale Verbände, wie z. B. Herden oder Schwärme.

Anwendung in der Epidemiologie

In der Epidemiologie bezeichnet Überdispersion eine hohe individuen-spezifische Variation in der Verteilung der Anzahl der Sekundärübertragungen, die zu „außergewöhnlichen Übertragungsereignissen“ (englisch superspreading events) führen kann. Üblicherweise werden Zähldaten, wie die Anzahl der Sekundärübertragungen, mithilfe der Poisson-Verteilung modelliert. Allerdings bietet es bei Vorliegen von Überdispersion Vorteile, die Anzahl der Sekundärübertragungen statt mit einer Poisson-Verteilung mit einer negativen Binomialverteilung zu modellieren, da sie einen flexiblen Grad an Übertragungsheterogenität abbilden kann und zu Daten aus einer Reihe von Infektionskrankheiten passt. Zudem scheitert die Poisson-Verteilung daran die relevanten Eigenschaften von „außergewöhnlichen Übertragungsereignissen“ zu erfassen, da sich in diesem Fall Erwartungswert und Varianz entsprechen (siehe #Poisson-Modell). Aus oben genannten Gründen lässt sich die Anzahl der Sekundärübertragungen adäquat durch eine negative Binomialverteilung modellieren. Da sowohl die Basisreproduktionszahl 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_{0}} (die „mittlere“ Anzahl an Sekundärübertragungen, die durch ein Individuum in einer anfälligen Population hervorgerufen wurde), als auch der Überdispersionsparameter 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 \kappa} (die individuelle Übertragungsheterogenität) Rückschlüsse auf die Dynamik eines Krankheitsausbruchs zulassen und beide Größen isoliert betrachtet wenig aussagekräftig sind, parametrisiert man oft die Wahrscheinlichkeitsfunktion der negativen Binomialverteilung durch diese beiden Größen gemeinsam (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 \mathcal{NB}(R_{t},\kappa)} ). Anschließend lassen sich beide Größen, die nun die Verteilung charakterisieren, mittels statistischer Verfahren gemeinsam schätzen.

Allerdings ist man oft primär am geschätzten Überdispersionparameter interessiert, da er die Variabilität in der Anzahl der Sekundärfälle quantifiziert und somit als Maß für die Wirkung von Superspreading interpretiert werden kann. Da die Varianz der angenommenen negativen Binomialverteilung gegeben ist durch $R_{0}(1+R_{0}/\kappa )$ liegt ein hoher Grad an individuen-spezifischer Variabilität vor, wenn der Überdispersionsparameter klein, z. B. Fehler beim Parsen (Konvertierungsfehler. Der Server („https://wikimedia.org/api/rest_“) hat berichtet: „Cannot get mml. Server problem.“): {\displaystyle \leq 1} , ist. Für diesen Spezialfall weist die Verteilung der Sekundärübertragungen exponentielle Verteilungsenden auf, d. h. die Eintrittswahrscheinlichkeit von „außergewöhnlichen Übertragungsereignissen“ nimmt steigender Anzahl an Sekundärfällen exponentiell ab.

Allgemein gilt: Je kleiner der Überdispersionsparameter ist, desto stärker ist die Wirkung von Superspreading. Dagegen würde mit steigendem Überdispersionsparameter der Effekt von Superspreading auf die Epidemie abnehmen. Die Interpretation des Überdispersionsparameters wird weiter vereinfacht, indem sich auf den Anteil der Individuen konzentriert wird, der für 80 % der Sekundärübertragungen verantwortlich ist (ein empirisches Muster, bekannt als 80/20-Regel). Ist der Überdispersionsparameter deutlicher kleiner als Eins (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 \kappa << 1} ) bzw. nahe bei Null, so approximiert er den Anteil infizierter Personen, die 80 % der gesamten Sekundärübertragungen verursachen. Beispielsweise würde ein geschätzter Überdispersionsparameter von 0,1 bedeuten, dass die infektiösesten 10 % der Personen etwa 80 % der gesamten Sekundärübertragungen verursachen.

Die Verteilung der Anzahl der Sekundärfälle kann durch eine Fréchet-Verteilung modelliert werden, die „fette Verteilungsenden“ aufweist.

Zusätzlich gibt es für einige Erkrankungen empirische Belege dafür, dass die Verteilung der Anzahl der Sekundärübertragungen „fette Verteilungsenden“ aufweist, d. h. die Verteilung besitzt relativ viel Wahrscheinlichkeitsmasse am auslaufenden Ende der Verteilung. Diese spezielle Eigenschaft bedeutet, dass umso mehr statistische Information in den Extremen und weniger in den Ereignissen steckt, die häufig auftreten. Daher sind „außergewöhnliche Übertragungsereignisse“ zwar extreme, aber dennoch wahrscheinliche Ereignisse, die einen beträchtlichen Beitrag zur Gesamtübertragung leisten. Maßnahmen sollten daher an der endlastigen Natur der Verteilung ansetzen, indem man versucht das Risiko an den Verteilungsenden zu minimieren wie z. B. durch die Immunisierung ausgewählter Personen. Da fette Verteilungsenden jedoch inkonsistent mit der Annahme einer negativen Binomialverteilung sind, werden auch alternative Modellierungsmöglichkeiten zum Negativ-Binomialmodell vorgeschlagen. Die folgende Tabelle gibt die geschätzten Überdispersionsparameter bei ausgewählten Erkrankungen an und ob es Hinweise auf fette Verteilungsenden bei der Verteilung der Sekundärübertragungen gibt:

Krankheit	Geschätzter Überdispersionsparameter	Hinweise auf fette Verteilungsenden
COVID-19	0,1–0,6	ja
SARS	0,16	ja
MERS	0,25	nicht bekannt
Spanische Grippe	1	nicht bekannt
Saisonale Grippe	1	nicht bekannt

Bei COVID-19 geht man bislang von einem geschätzten Überdispersionsparameter von etwa 0,1–0,6 aus. Laut einem Preprint von Akira Endo und Mitautoren, liegt bei COVID-19 – unter Annahme einer Basisreproduktionszahl von 2,5 – der Überdispersionsparameter mit hoher Glaubwürdigkeit etwa bei 0,1 (95 %-Glaubwürdigkeitsintervall: [0,05–0,2]).