text: rechnung von R t herangezogen. Für Fälle, bei denen das Erkrankungsdatum nicht bekannt ist, wird zunächst eine angepasste verallgemeinerte Extremwertverteilung (GEV) auf die Daten, bei denen die Angabe des Meldeverzugs (Differenz zwischen Krankheitsbeginn und Meldedatum) vorhanden ist, angepasst und auf weitere Merkmale wie lokale sowie demografische Differenzen untersucht. Diese GEV-Verteilung wird dann zur Berechnung der Meldeverzüge der unvollständigen Datenpunkte verwendet. Bei der Berechnung des täglichen Wertes von R t zwischen Ende Februar und Ende Oktober konnte ein im Vergleich mit den deutschlandweiten Zahlen ähnlicher Verlauf der Reproduktionszahl festgestellt werden. Erwartbar größeren statistischen Schwankungen wurden im Sommer, vor allem durch geringere Fallzahlen, festgestellt, wobei die Werte für das nördliche Rheinland-Pfalz seit etwa Mitte September konstant über 1 liegen. Die Berechnungen können ebenso auf weitere Regionen und Landkreise übertragen werden.
  --Since the beginning of the corona pandemic in March 2020, various parameters for describing the spread of the disease have been specified for Germany in addition to the daily infection figures (new infections and total infections), which are also used for political decisions. In addition to excess mortality and the weekly incidence, these include the doubling time T 2 and the reproduction number R t . For the latter, various estimates can be found on the website of the Robert-Koch-Institute, see , which are calculated from the case numbers for all of Germany; local differences are not taken into account here. In the present article, the calculations of the RKI on a local level are examined using the example of northern Rhineland-Palatinate and its districts. Here, not the reporting date but the onset of illness is used as a reference for the calculation of R t . For cases where the onset of illness is not known, an adjusted generalized extreme value distribution (GEV) is first fitted to the data for which the reporting delay (difference between the onset of illness and the reporting date) is available and examined for further characteristics such as local as well as demographic differences. This GEV distribution is then used to calculate the reporting delays of incomplete data points. The calculation of the daily value of R t between the end of February and the end of October showed a similar course of the reproductive rate compared to the nationwide figures. Expectably larger statistical fluctuations were observed in the summer, mainly due to lower case numbers. The values for northern Rhineland-Palatinate have been consistently above 1 since about mid-September. The calculations can also be transferred to other regions and administrative districts.
  Datenlage
  Das Robert-Koch-Institut (RKI) veröffentlicht täglich aktualisierte Datensätze zu den gemeldeten COVID-19 Infektionen . Für jede gemeldete Infektion enthält der Datensatz neben demographischen Angaben (Alter, Geschlecht) auch geographische Informationen (Bundesland und Landkreis) sowie Angaben zum Meldedatum (wann wurde der Fall beim jeweiligen Gesundheitsamt registriert?), für Details siehe . Weniger als die Hälfte der Datensätze (Stand 31. Oktober ca. 40% im gesamten nördlichen Rheinland-Pfalz) enthält zusätzliche Informationen zum vermutlichen Datum des Erkrankungsbeginns. In Abbildung 1 wird der Meldeverzug (entspricht der Differenz zwischen Erkrankungsbeginn und Meldedatum) für die im nördlichen Rheinland-Pfalz gemeldeten Fälle dargestellt, ebenso für die Stadt Koblenz sowie für die acht Landkreise (Mayen-Koblenz, Ahrweiler, Altenkirchen, Cochem-Zell, Neuwied, Rhein-Hunsrück-Kreis, Rhein-Lahn-Kreis und Westerwaldkreis).
  Berechnung der lokalen Reproduktionszahl
  Für Datensätze ohne Angabe des Erkrankungsbeginns kann dieser gemäß der angepassten GEV-Verteilung geschätzt werden. Basierend hierauf wird analog zu der in beschriebenen Vorgehensweise die aktuelle Reproduktionszahl R t berechnet; Details hierzu beschreiben wir in Anhang A. Um statistische Schwankungen teilweise zu eliminieren, verwenden wir eine Mittelung der Anzahl der Erkrankungsbeginne über ein Zeitintervall von τ = 4 bzw. τ = 7 Tagen und geben entsprechend die aktuelle, über 4 bzw. 7 Tage gemittelte Reproduktionszahl R 4,t bzw. R 7,t an. Hierbei ist zu beachten, dass die Werte für beide Reproduktionszahlen durch die maximal möglichen Meldeverzüge von 21 Tagen noch in ebendiesem Zeitraum veränderlich bleiben, sodass die jüngsten, aktuellen Resultate gesondert dargestellt werden.
  Bewertung der Resultate
  Betrachtet man die Daten zum Meldeverzug in den einzelnen Landkreisen des nördlichen Rheinland-Pfalz, vgl. Abbildung 1, so fällt auf, dass für den Landkreis Neuwied in lediglich 1% aller Fälle (13 von gesamt 991) der Erkrankungsbeginn verzeichnet ist. Eine Erklärung hierfür, bzw. eine Rückmeldung des zuständigen Gesundheitsamtes fehlen bislang. In den anderen Landkreisen liegt für jeweils zwischen 39% (Koblenz und Mayen-Koblenz) und 72% (Rhein-Lahn-Kreis) der Fälle sowohl das Meldedatum als auch der vermutete Erkrankungsbeginn vor. Der mittlere Meldeverzug als Differenz zwischen Erkrankungsbeginn und Meldedatum schwankt zwischen 4.1 Tagen (Altenkirchen) und 5.7 Tagen (Rhein-Lahn-Kreis). Im Durchschnitt beträgt der Meldeverzug für das gesamte nördliche Rheinland-Pfalz etwa 5.2 Tage. Die jeweilige Verteilung des Meldeverzug verhält sich für alle Landkreise (mit Ausnahme von Neuwied aus den genannten Gründen) ähnlich.
  Bei der Schätzung der 7-Tage Reproduktionszahl für das nördliche Rheinland-Pfalz (Grafik 3) fällt auf, dass zu Beginn der Pandemie (März) die Reproduktionszahl teils sehr große Werte annahm. Dies entspricht dem damaligen raschen Wachstum mit Verdopplungszeiten von minimal 2-3 Tagen. Die Einführung der Kontaktbeschränkungen und Einschränkungen im öffentlichen Leben ab Mitte März reduzierte den R-Wert in der Folge unterhalb der kritischen Marke von 1. Die Oszillationen im Sommer (Juni-Juli) lassen sich teilweise als statistische Schwankungen aufgrund der niedrigen Fallzahlen erklären (es wurden insgesamt nur ca. 220 Krankheitsbeginne im gesamten nördlichen Rheinland-Pfalz in diesem Zeitfenster simuliert). Die niedrigen Fallzahlen während des Sommers verursachen auch die teilweise extrem großen Schwankungen des 4-Tage R-Wertes für die einzelnen Landkreise. In Perioden, in denen es überhaupt keine Neuinfektionen innerhalb eines Landkreises gab, z.b. Mitte Juni-Mitte Juli im Rhein-Hunsrück-Kreis, siehe Abbildungen C.20-C.21, lässt sich kein R-Wert berechnen. Ab September ist ein sukzessives Ansteigen des R-Wertes konstant oberhalb der Marke von 1 (bis auf wenige Ausnahmen) zu beobachten. Dieser Anstieg hat sich im Oktober konsolidiert und wurde für den 10. Oktober auf geschätzt. Ein R-Wert von 1.33 entspricht in etwa einer Verdopplungszeit von 3 Wochen. Die Meldeverzüge für das nördliche Rheinland-Pfalz genügen der Verteilung F". Die Wahrscheinlichkeitsdichte der GEV mit einen scale-Parameter σ, einen shape-Parameter ξ und einen location-Parameter µ ist wie folgt definiert :
  f GEV (x) =        exp − x−µ σ • exp − exp − x−µ σ , falls ξ = 0 (1 + ξ • x−µ σ ) −(1+ 1 ξ ) • exp −(1 + ξ • x−µ σ ) − 1 ξ , falls ξ = 0 und ξ • x−µ σ > −1 0 , sonst . (A.1)
  Die Verteilung der Meldeverzüge der einzelnen Landkreise (außer NR und WW) genügen ebenfalls einer GEV, siehe Die aktuelle Reproduktionszahl R t zum Zeitpunkt t ist definiert als die mittlere Anzahl von Sekundärinfektionen, welche ein am Tag t Infizierter verursacht. Geht man von einem seriellen Intervall σ von 4 Tagen aus, vgl. , so ergibt sich die aktuelle Reproduktionszahl R t als das Verhältnis der heutigen neuen Erkrankungsfälle E t zu denen von vor σ Tagen, also R t = Et Et−σ . Um die statistischen Schwankungen in den täglichen Neuerkrankungen teilweise zu eliminieren, verwenden wir anstelle der tagesaktuellen Werte ein gleitendes Mittel der letzten τ Tage, vergleiche auch :
  R τ,t = 1 τ τ −1 k=0 E t−k 1 τ τ −1 k=0 E t−σ−k = E τ,t E τ,t−σ (A.2) Hierbei bezeichnet E τ,t := 1 τ τ −1 k=0 E t−k
  das gleitende Mittel der Neuerkrankungen der letzten τ Tage. Ergebnisse hierzu für das nördliche Rheinland-Pfalz mit τ = 4 sowie τ = 7 sind in Abbildungen 2 und 3 dargestellt. Der detailliertere Verlauf des R 7,t -Wertes für das nördliche Rheinland-Pfalz seit Anfang August findet sich in Die Reproduktionszahl R t steht in direktem Zusammenhang mit der Verdopplungszeit T 2 der Infektionsfälle. In einer Bevölkerung, die noch größtenteils suszeptibel für die Erkrankung ist, lässt sich die Anzahl der aktuell Infizierten näherungsweise durch die Differentialgleichung
  I ′ = β SI N − γI
  beschreiben, vgl. . Für S ∼ N und R t ∼ β γ erhalten wir I(t) = I 0 e γ(Rt−1)t . Hieraus ergibt sich die Verdopplungszeit zu
  T 2 = log 2 γ(R t − 1) . (A.3)
  Nimmt man eine Genesungsrate γ = 1/10 an -dies entspricht einer mittleren Infektionsdauer von 10 Tagen, vgl. -so ergibt sich ein direkter Zusammenhang zwischen der Reproduktionszahl R t und der Verdopplungszeit T 2 .
  Anhang B. Kohortenanalyse
  Anhang B.1. Verteilungstest Um zu überprüfen, ob ein Datensatz einer gegebenen Verteilung genügt, wird ein Pearson's X 2 -Verteilungstest durchgeführt. Die Nullhypothese lautet H 0 (F ) : " Die Daten genügen der Verteilungsfunktion F" und wird überprüft, indem der Prüfgröße Hierbei ist festzustellen, dass die Fälle bei 5-14-Jährigen sowie über 80-Jährigen signifikant schneller gemeldet werden als der Durchschnitt, bei 60-80-Jährigen jedoch signifikant langsamer. Abbildung B.5 liefert eine dazugehörige grafische Repräsentation. 
  X 2 F := N k=1 (O k − E(F ) k ) 2 E(F ) k die Testgröße X 2 1−α,
  