Martins Menetekel Anfang vom Ende?

Wir alle kennen die Botschaften des RKI: Rekord an dem Zuwachs pro Tag an als infiziert Gemeldeten, über 200.000, und 7-Tages-Inzidenz über 1.000 .

Was bedeutet das?

1. Die Nachverfolgung der Infizierten und deren Kontakte ist nicht mehr möglich
2. Noch in dieser Jahreshälfte sind entweder alle in Deutschland infiziert oder von selbst oder durch Impfung immun
3. Aus der Pandemie wird eine endemisch flackernde Normalseuche

Omikron kann tatsächlich das Ende der Pandemie bewirken, und das bei nicht allzu großen Opferzahlen.

Es könnte also sein, dass China mit seiner 0-COVID-Strategie letztendlich zwar sehr viel weniger Todesfälle hatte, aber der Immunisierungsgrad sehr klein ist. Sie müssen also viel länger wachsam sein.

Großbritannien geht diesen Weg, auch Dänemark. Deutschland geht ihn auch, nur wird es verschwiegen.

Jedenfalls haben wir seit Anfang Januar exponentielles Wachstum, dominiert durch Omikron.

Das möchte ich heute erklären.

Da die Zuwächse pro Tag noch ständig wachsen, ist der erste Ansatz eine Annäherung von N’(t) durch eine reine Exponentialfunktion B*b^t .

N’(t) ist unabhängig vom Sockel. Eine Näherung oder Modellierung ist dann gerechtfertigt, wenn N’(t)/N’(t-1) konstant ist, und zwar dann (B*b^t)/B*b^(t-1)) = b . denn B kürzt sich weg.

Die Grafik seit dem 30.12.21 bis zum 21. Januar, immer von Donnerstag bis Donnerstag, bestätigt das:

Besser geht es nicht.

Wie finde ich das B, den Normierungsfaktor.

Wir erinnern uns: Bei jeder Näherungslösung einer Zeitreihe Z(t) durch f(t) soll gelten:

Die Summe vom 30.12. bis 20.1. aller (Z(t) - f(t))^2 muss so klein wie möglich sein. Diese Quadratsumme hängt von der Zeitreihe Z(t) und von den Parametern in der Näherungsfunktion f(t) ab.

Unsere Zeitreihe ist N’(t) und die Näherungsfunktion f(t) ist Bb^t .

Summe über 22 Werte von (N’(t) - Bb^t)^2 soll minimal sein. Was bedeutet das?

Ich muss ein B finden, für das gilt: Wenn ich B etwas ändere, wird die Summe größer!

Aber das tut gerade die Ableitung nach B , die ich dann zu 0 setze: Das Minimum einer Funktion liegt bei einer Nullstelle der Variablen, hier B:

Die Ableitung von (N’(t) - Bb^t)^2 nach B ergibt 2(N’(t) - Bb^t)*(-b^t) , die als Summe = 0 sein muss, und damit ist B = Summe(N’b^t)/Summe(b^t)*(b^t) . Die 2 und das Minuszeichen kürzen sich weg.

Das macht jedes Tabellenprogramm mit links. Wir haben es berechnet, B = 29.980

Das kleine b als Basis der E-Funktion finde ich aus dem Mittelwert aller N’(t)/N’(t-1) und das war 1,06

Unsere Näherungsfunktion für N’ ist 29.980*(1,06^t) und ergibt die Grafik:

Ganz gut. Man könnte jetzt  noch b optimieren, das lohnt sich nicht, denn das exponentielle Wachstum ist immer zeitlich begrenzt.

Wie erhalte ich die Annäherung an die Fallzahlen selbst. Das waren die kumulierten gemeldeten Zuwächse pro Tag, aber geglättet. Die Tageswerte muss man glätten, wie die folgende Grafik zeigt. Man kennt vielleicht die Werte aus der Tagesschau: 

Jetzt haben wir N’(t) angenähert durch Bb^t . Eine Stammfunktion von Bb^t ist Bb^t/ln(b) , denn wenn ich (1/ln(b)Bb^t ableite, kürzt sich ln(b) weg (implizite Ableitungsregel)

ln(b) = k = 0,0602 . Das ist genau N’/N , wenn ich schon vorweg den Sockel 6.629.036 von den Fallzahlen abgezogen habe:

Der letzte Schritt:

Ich setze A = B/0,0602 = 498.392 mit den Originalwerten von B und b

Ich bilde F(t) - Ab^t als neue Zeitreihe, summiere wieder vom 30.12. bis 20.1 und bilde den Mittelwert, das ist der Sockel I, I wie Integrationskonstante, denn. mit Ab^t ist auch I + Ab^t eine Stammfunktion von N’ , bei uns ab heute:

N(t) = 6.629.036 + 498.392*(1,062^t) und so sieht das aus:

Den Sockel werden wir hoffentlich bald für die Verhulst-Annäherung verwenden dürfen, es ei denn, die Durchseuchung ist gewollt und hört erst auf, wenn die Pandemie beendet ist.

In der Endphase herrscht in jedem Fall Verhulst mit Limitierung, denn:

Wenn die Häschen auf der Wiese mehr Gras fressen, als nachwächst, wird ihre Fortpflanzungsrate garantiert kleiner!

ABER: SARS-CoV-2 kann unter Umständen noch mehr Varianten in petto haben, die alles über den Haufen werfen können. Vielleicht werden wir die Pandemie nie los und die chinesische Strategie war die bessere, weil wirksamer gegen Mutationen. Es braucht nur eine Variante, die infektiöser als Omikron ist und zu schwereren Krankheitsverläufen führt, aber gegen die die heutigen Impfstoffe nichts ausrichten, auch nicht die Immunantworten der Infizierten, die gegen die älteren Varianten gebildet wurden. Dann wäre das heutige Motto widerlegt.

Hoffen wir das Beste, lieber Leser!

Über Martin Lindner