Model to Describe Fast Shutoff of CoVID-19 Pandemic Spread ========================================================== * Genghmun Eng ## Abstract Early CoVID-19 growth obeys: ![Graphic][1], with *Ko* = [(ln 2)/(*tdbl*)], where *tdbl* is the pandemic growth *doubling time*. Given ![Graphic][2], the daily number of new CoVID-19 cases is ![Graphic][3]. Implementing society-wide *Social Distancing* increases the *tdbl doubling time*, and a linear function of time for *tdbl* was used in our *Initial Model*: ![Formula][4] to describe these changes, with *Go* = [*KA*/*γo*]. However, this equation could not easily model some quickly decreasing *ρ*[*t*] cases, indicating that a second *Social Distancing* process was involved. This second process is most evident in the initial CoVID-19 data from *China, South Korea*, and *Italy*. The *Italy* data is analyzed here in detail as representative of this second process. Modifying *Zo*[*t*] to allow exponential cutoffs: ![Formula][5] provides a new *Enhanced Initial Model (EIM)*, which significantly improves datafits, where ![Graphic][6]. Since large variations are present in *ρdata* [*t*], these models were generalized into an orthogonal function series, to provide additional data fitting parameters: ![Formula][7] Its first term can give *No*[*t*] or *NE* [*t*], for *Z*[*t*] → *Zo*[*t*] or Z[t] → *Ze*[*t*]. The *Lm*(*Z*) are *Laguerre Polynomials*, with *L*(*Z*) = 1, and {*gm*; *m* = *0,MF*} are constants derived from each dataset. When *ρ*[*t*] = *dN*[*t*]/*dt* gradually decreases, using *Zo*[*t*] provided good datafits at small *MF* values, but was inadequate if *ρ*[*t*] decreased faster. For those cases, *ZE*[*t*] was used in the above *N*(*Z*) series to give the most general *Enhanced Orthogonal Function [EOF]* model developed here. Even with *MF* = 0, *qo* = 0, this *EOF* model fit the *Italy* CoVID-19 data for *ρ*[*t*] = *dN*[*t*]/*dt* fairly well. When the *ρ*[*t*] post-peak behavior is not Gaussian, then *ZE*[*t*] with *δo* ≠ 0, *qo* = 0; which we call *ZA*[*t*], is also likely to be a sufficient extension of the *Zo*[*t*] model. The *EOF* model also can model a gradually decreasing *ρ*[*t*] tail using small {*δo*, *qo*} values [with *6 Figures*]. ## 1 Introduction Let ![Graphic][8] be the total number of CoVID-19 cases in any given locality, with ![Graphic][9] being the predicted number of daily new CoVID-19 cases, so that: ![Formula][10] ![Formula][11] Early CoVID-19 growth often obeys ![Graphic][12], with *Ko* = [(ln2)/*tdbl*, where *tdbl* is the pandemic *doubling time*. The start of society-wide *Social Distancing* at ![Graphic][13] can gradually lengthen *tdbl* for ![Graphic][14]. The ![Graphic][15] exponential growth phase is not applicable for estimating *Social Distancing* effects. For ![Graphic][16], an *Initial Model* for CoVID-19 pandemic shutoff was first developed**1** using a linear function of time to describe the *tdbl* changes: ![Formula][17] Given measured ![Graphic][18], the data end-points {*NI*, *NF*} help to set {*Ko*, *αS*}. An *Orthogonal Function Model [OFM]* was developed next**2**, with Eq. [1.2] as the first term of the orthogonal function series. Each new *OFM* term provides another fitting parameter, to progressively better match ![Graphic][19] and ![Graphic][20]. The *OFM* improves on the *Initial Model*, and it works best with gradually decreasing ![Graphic][21] [“Slow *Shutoff*”]. In contrast, when ![Graphic][22] decreased quickly [“Fast *Shutoff*”], the *Initial Model* was not a good datafit, and a few-term *OFM* series only gave small improvements. This result indicates there is an inherent limit to what the gradually changing *tdbl doubling time* of Eq. [1.2] can model. For these cases, typified by CoVID-19 pandemic evolution in Italy, data often showed a stage where ![Graphic][23] or ![Graphic][24], which likely represents a second process, independent of the gradually changing *tdbl doubling time*. An *Enhanced Initial Model (EIM)* is developed here to include this second process. The prior *OFM* methods can then be applied, giving an *Enhanced Orthogonal Function (EOF)* model for this more general case. ### 1.1 Review of Prior Models The *Initial Model* of Eq. [1.2] is still needed as the first part of the *OFM*. The *Initial Model* starts with measured data end-points {*NI*, *NF*}, where ![Graphic][25] is the largest data time interval so that ![Graphic][26]. Usually *αS* in Eq. [1.2] was chosen first, and *Ko* or *tdbl* adjusted to match the *NF* data end-point, using an *ExcelTM_Goal-Seek* or its equivalent. The final {*Ko*, *αS*} values were the pair with the minimum *root-mean-square* (*rms*) error between the given data and the Eq. [1.2] model. The above {*Ko*, *αS*} also provides a *t* = 0 estimate for the pandemic start, and gives {*KA*, *γo*} as new data fitting parameters: ![Formula][27] ![Formula][28] ![Formula][29] ![Formula][30] Using (*tF* − *tI*), with *No*[*t* = *tI*] = *NI* and *No*[*t* = *tF*] = *NF*, sets: ![Formula][31] ![Formula][32] which determines {*KA*, *γo*} in *Zo*[*t*] for Eq. [1.3c]: ![Formula][33] ![Formula][34] The Eq. [1.4a] value for *tI* is what determines the new *t* = 0 point, as an extrapolation for when *No*[*t* = 0] = **1**. In addition: ![Formula][35] ![Formula][36] provides an estimate for the total number of cases ![Graphic][37] at the pandemic end, and determines a function for the *ρo*[*t*] long-time tail. The {0 < *t* < *tI*} period prior to the start of *Social Distancing*, extrapolates what pandemic progress would have looked like, if *Social Distancing* had begun at *t* = 0. An *Orthogonal Function Model [OFM]* was then developed**2** to better model the different observed *ρdata*[*t*] shapes, as an improvement of the *Initial Model*: ![Formula][38] ![Formula][39] ![Formula][40] ![Formula][41] with *Lm*(*Z*) being the *Laguerre Polynomials*, and *Lm*(*Z* = 0) = *Lo*(*Z*) = 1. Using *Z* = *Zo*[*t*] from Eq. [1.3c] gives *N*(*Z*) → *N*(*Zo*) and *R*(*Z*) → *R*(*Zo*). The {*gm*; *m* = (0, *MF*)} constants in Eq. [1.7a] can be arranged in a ![Graphic][42] —vector form, with comparable constants for R(Z) from Eq. [1.7b] arranged in a ![Graphic][43] —vector form. For *MF* = 2, it allows Eq. [1.7d] to be written as: ![Formula][44] Once these *cm* values are determined for *R*(*Z*) in Eqs. [1.7b]-[1.7c], an *OFM* feature is that *c*, by itself, becomes the *OFM* best estimate for the total number of CoVID-19 cases at the pandemic end. For large enough *MF* values, and a monotonic *Z*—function, this *OFM* can provide successively better approximations to almost any given set of *Ndata*[*Z*], with *Z*[*t*] → *Zo*[*t*] of Eq. [1.3c] being a specific case. The *OFM* implicitly uses a *Linear Y-axis*, so its results differ from the *Initial Model* datafit on a *Logarithmic Y-axis*. As an example, compare the *Initial Model* result of *No*(*Zo*) = *G* exp[−*Zo*] with the Eq. [1.7a] *OFM* result of *N*(*Zo*) = *g* exp[−*Zo*] for *MF* = 0. In the *Initial Model, Go* is fixed so that *No*(*Zo*) exactly matches {*NI*, *NF*} at the {*tI*, *tF*} boundaries. In the *OFM*, *g* = *Go* is no longer required, so that the *OFM N*(*Zo*) best datafit is not constrained to exactly match {*NI*, *NF*} at {*tI*, *tF*}. The above *R*(*Z*) and *Z*[*t*] gives *N*[*t*] and *ρ*[*t*] as an explicit functions of time: ![Formula][45] Since Eq. [1.6b] gives *ρ*{*t*} ~ [1/*t*2], the aim here is to model faster decaying functions such as *ρ*{*t*} ~ [exp(−*δot*)] or *ρ*{*t*} ~ [exp(−*qot*2)]. ### 1.2 Updated *Initial Model* Results for Italy The *Zo*[*t*] model of Eq. [1.3a] was applied to *bing.com* data**9** for Italy, starting with ![Graphic][46] CoVID-19 cases as an early pandemic point, up through June 15, 2020. Here, *tI* is when mandatory *Social Distancing* was introduced at ![Graphic][47]; with *tF* being when ![Graphic][48]. Data prior to *Social Distancing* (*t* < *tI*) was excluding from this *Social Distancing* analysis. **Figure 1** compares the *ρdata*[*t*] results with the updated *No*[*t*] and *ρo*[*t*] = *dNo*[*t*]/*dt* predictions using *Zo*[*t*]. The ![Graphic][49] (Day 1) to ![Graphic][50] (Day 17) interval was examined for estimating a *t* = 0 pandemic start where *Ndata*[*t* = 0] → **1**. A best fit value of *toffset* = 9.10055 days was found, giving *tI* = (17 - 9.10055) = 8.89945 days for ![Graphic][51] (Day 17), while ![Graphic][52] (Day 114) gives *tF* = (114 − 9.10055) = 104.89945 days, so that: ![Formula][53] ![Formula][54] ![Formula][55] ![Formula][56] with {*KA*, *γo*} ≈ {4.2405, 0.33078} and (*tF* − *tI*) = 97 days. In **Fig. 1**, the X-axis uses this *t* = 0 point where *No*[*t* = 0] → 1, and it shows what *Social Distancing* effects would have been, if it had been operating throughout the *t* > 0 period. This *ρo*[*t*] prediction still has a much more gradual drop than the data. This discrepancy indicates that a second *Social Distancing* process is operating, besides just the gradual t^ lengthening of the *Initial Model*. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/11/2020.08.07.20169904/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/08/11/2020.08.07.20169904/F1) Figure 1: Updated *Initial Model* for ITALY, CoVID-19 data to 6/15/2020. Number of daily CoVID-19 cases calculated as if *Social Distancing* started at 3/2/2020, but only data from 3/10/20 actual *Social Distancing* start, with *N* = 10, 149 cases, was used in calculations. **Figure 2** compares the *No*[*t*] predictions for this model, to the measured *Ndata*[*t*]. Systematic deviations are evident, with the net *rms* error on a *Logarithmic Y-axis* being *rmserror* = 0.097828. To cure these defects, an enhanced *Z*[*t*] model is developed next. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/11/2020.08.07.20169904/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/08/11/2020.08.07.20169904/F2) Figure 2: Comparison of *Initial Model* Best-Fit to Measured Data. Best-fit done on *Logarithmic Y-axis*, using data from 3/10/20 actual *Social Distancing* start with *N* = 10, 149 cases, through 6/15/20 with *N* = 237, 290. ## 2 Developing Enhanced *Z*[*t*] Models To generalize *Z*[*t*] beyond Eq. [1.3c], it is convenient to use the *t* = {0+, ∞−} domain, and require convergence of *Z*[*t*] → *Zo*[*t*] in some limit, along with: ![Formula][57] ![Formula][58] and that the *MF* = 0 case of Eq. [1.7a] remains as: ![Formula][59] How to choose an appropriate *Z*[*t*], as part of an *Enhanced Initial Model (EIM*), which also allows a *ρ*[*t*] ~ [exp(−*δot*)] or *ρ*[*t*] ~ [exp(−*qot*2)] stage, is examined next. It can be motivated by studying a simple *NT*[*t*] test-case, where *ρT*[*t*] itself is a pure exponential decay, as a function of time: ![Formula][60] ![Formula][61] while also preserving *NT*{*t* → 0} = 1. Comparing Eq. [2.3a] to the Eq. [2.2], sets *ZT*[*t*] for this test-case: ![Formula][62] ![Formula][63] At large times, Eq. [2.4b] gives: ![Formula][64] since (− ln[1 − *x*] ≈ *x*) for small x, which shows that if *ρT*[*t*] has an exponential tail, then *ZT*[*t*] also has an exponential tail. A simple generalization for *Z*[*t*] in Eq. [1.7a]-[1.7d] would be either: ![Formula][65] ![Formula][66] As with Eq. [1.3c], the original CoVID-19 exponential growth factor *KA* remains only as part of the *Go* scaling factor, while the {*δo* → 0} limit of Eq. [2.6b] converges back to the Eq. [1.3c] *Initial Model*. Using Eq. [2.6a] for *ZE*[*t*] in the Eq. [2.2] *N*(*ZE*) example gives: ![Formula][67] which exhibits the following variety of long-time limits: ![Formula][68] ![Formula][69] ![Formula][70] ![Formula][71] Here, any *qo* ≠ 0 Gaussian component in *ZE*[*t*] gives a *ρ*[*t*] tail that is also a pure Gaussian. An exponential component (*qo* ≡ 0) in *ZE*[*t*] gives a time- modified exponential *ρ*[*t*] tail, while having {*q* ≡ 0, *δo* ≡ 0} in *ZE*[*t*] gives the prior *ρ*[*t*]~(1/*t*2) result of Eq. [1.6b]. ## 3 Pandemic Fast vs Slow Shutoffs Each *ZE*[*t*] function modifies *N*[*t*] predictions for the pandemic start, pandemic end, and the mid-range where *ρ*[*t*] has its pandemic peak. The Eq. [2.6a] *ZE*[*t*] function especially alters the calculated CoVID-19 pandemic tail. For either *qo* ≠ 0 or *δo* ≠ 0, Eqs. [2.6a]-[2.6b] gives a pandemic *Fast Shutoff*, compared to the gradually decreasing *Zo*[*t*] of Eq. [1.3c] in the *Initial Model*, which is a pandemic *Slow Shutoff*. However *Zo*[*t*] from the *Initial Model*, and *ZA*[*t*] from Eq. [2.6b] both gave long-term *ρ*[*t*] tails that decay much slower than the *qo* ≠ 0 Eq. [2.6a] Gaussian. If data does not show evidence of a Gaussian pandemic *Fast Shutoff*, assuming the post-peak *ρ*[*t*] data will be Gaussian is likely to provide optimistically inaccurate *N*[*t*] predictions for CoVID-19 pandemic evolution. Apparently, this is exactly what was done by the *University of Washington IHME (Institute of Health Metrics and Evaluation)* in their widely publicized initial preprint**3** of 27 March 2020, with this Gaussian model continuing throughout their subsequent updates**4-6** up through 29 April 2020. *IHME* changed everything in their 4 May 2020**7–8** update. They no longer used *ρ*[*t*] Gaussian tails, and it doubled or tripled their predicted CoVID-19 pandemic death rates. Thus, unless the post-peak *ρ*[*t*] exhibits Gaussian behavior, the *ZA*[*t*] with *δo* =0 is likely the most important modification to *Zo* [t], which is the pandemic *Fast Shutoff* model used here. Using *ZA*[*t*], the pandemic *Fast Shutoff* can be extrapolated to calculate a pandemic start point where *N*[*t* = 0] = 1. We then examine if this *ZA*[*t*] must also carry over to the *ρ*[*t*] long-term tail, Since the long-term low *pdata*[*t*] tail may differ among localities, and is not well known, the *δo* ≠ 0 case of Eq. [2.6b] could end with a *Slow Shutoff*, giving: ![Formula][72] ![Formula][73] with *Go* as in Eq. [2.1a]. This Eq. [3.1a] *ZB*[*t*] function has the {*δo*, *λo*} *Mitigation Measure* operating at the start of *Social Distancing*, but reverting to the *Initial Model* in the long-time limit. Combining Eq. [2.6b] and Eq. [3.1a] cases gives this *Enhanced Initial Model (EIM)* equation: ![Formula][74] where *κ* = 0 is a pure exponential, and *κ* = 1 has a modified tail that includes its own long-term shutoff. Comparing *κ* = {0,1} in Eq. [3.2] provides a simple test for which model matches CoVID-19 data better in any locality. Any other *κ* > 0 value then recovers the more general Ao = *κγo* case. This Eq. [3.2] *ZB*[*t*] replaces *Zo*[*t*] of Eq. [1.3c], and its *EIM* companion *NB*[*t*] is: ![Formula][75] while using *Z*[*t*] → *ZB*[*t*] in Eqs. [1.7a]-[1.7d] gives an *Enhanced Orthogonal Function [EOF]* model. ## 4 Finding {*KA*, *γo*, *δo*} for *ZB*[*t*] from Data If *δo* = 0, the prior Eqs. [1.4a]-[1.5b] for {*KA*, *γo*, *δo* = 0} and {*tI*, *tF*} can be used, with the initial best-fit {*Ko*, *αS* } values determined by minimizing the *rms* error between Eq. [1.2] and the measured data on a *Logarithmic Y-axis*. Unfortunately, Eqs. [1.4a]-[1.5b] cannot be used when *δo* = 0, although finding a good {*KA*, *γo*, *δo*} starting point is still needed for the *EIM*: ![Formula][76] prior to any *EOF* analysis. The *κ* = 1 case also has this special symmetry: ![Formula][77] ![Formula][78] ![Formula][79] ![Formula][80] which can re-written as: ![Formula][81] ![Formula][82] ![Formula][83] ![Formula][84] ![Formula][85] while the *κ* = 0 case is: ![Formula][86] ![Formula][87] For *κ* = {0, 1}, the *t* = 0 point, *Go* from Eq. [2.1a], and the *NI*(*tI*) and *NF*(*tF*) initial and final points, give these equations to help set {*KA*, *γo*, *δo*}: ![Formula][88] ![Formula][89] ![Formula][90] for the *EIM*. Minimizing the *rms* error between the Eq. [3.3] {*KA*, *γo*, *δo*, *tI*} functions and measured data on a *Logarithmic Y-axis* can be done as follows. Start with estimated values for ![Graphic][91] in: ![Formula][92] ![Formula][93] ![Formula][94] ![Formula][95] where *tdata* is the data start time. Set a preliminary value for ![Graphic][96] first, to fix the time scale for the *Ndata* [t] measured values: ![Formula][97] ![Formula][98] ![Formula][99] Next, pick values for ![Graphic][100] for ![Graphic][101] in Eq. [4.6b], allowing direct comparison between *N*[*t*] and ![Graphic][102] at each data point: ![Formula][103] ![Formula][104] The resulting calculated values for both ![Graphic][105] can often be much too high or low, compared to the {*NI*, *NF*} measured data, but those values can be renormalized to: ![Formula][106] ![Formula][107] Here, *SI* is the renormalization coefficient, and ![Graphic][108] allows easy comparison to the measured *Ndata*[*t*] since ![Graphic][109]. Given ![Graphic][110], the ![Graphic][111] value that is needed to obey ![Graphic][112] can be set by using *Excel*™*_Goal-Seek* or its equivalent, which also sets a particular *SI* value. Next, the ![Graphic][113] value is adjusted to find the specific {*KA*, *γo*} parameter pair that gives *SI* = 1. This process is needed because these ![Graphic][114] cases do not allow easy determination of *tI* as in Eq. [1.4a], or for {*KA*, *γo*}, given *tI*, as in Eqs. [1.5a]-[1.5b]. The *rms* error on a *Logarithmic Y-axis*, between this ![Graphic][115] and the *Ndata*[*t*] is one of many ![Graphic][116] choices. The minimum rms error among all these *SI* = 1 cases and the *Ndata*[*t*], when varying ![Graphic][117] gives the best {*KA*, *γo*, *δo*, *toffset*} values for Eqs. [4.3a]-[4.3e]. ## 5 *Enhanced Initial Model [EIM]* Results for Italy The Eq.[1.3a] *Initial Model* results were shown in **Figs. 1-2**. Mandatory *Social Distancing* was introduced at ![Graphic][118] which is the *tI* data point, with ![Graphic][119] being the *tF* data point. The *EIM* was then applied to the same data, to highlight the improvements that can be obtained from using the *EIM* of *ZA*[*t*] and *NA*(*ZA*), in place of *Zo*[*t*] and *No*(*Zo*). **Figures 3-4** show the resulting *EIM* best-fits for *ZA*[*t*] and *NA*(*ZA*), along with using a *ρA*[*t*] tail that is a pure exponential decay. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/11/2020.08.07.20169904/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/08/11/2020.08.07.20169904/F3) Figure 3: *Enhanced Initial Model (EIM)* for ITALY CoVID-19 data to 6/15/20. *EIM* best fit with *NA*[*t*] ~ exp(−*ZA*[*t*]) using enhanced *ZA*[*t*] function having an exponential decay. Adding in exponential decay term gives significantly improved fit, compared to prior *OFM*. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/11/2020.08.07.20169904/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/08/11/2020.08.07.20169904/F4) Figure 4: *Enhanced Initial Model (EIM)* datafit for ITALY CoVID-19 data to 6/15/20. *EIM* best fit with *NA*[*t*] ~ exp(−*ZA*[*t*]) using enhanced *ZA*[*t*] function having an exponential decay. Datafit minimizing *rms* error on *Logarithmic Y-axis* gives significant improvement vs *OFM*. For the *EIM*, a new best estimate of *toffset* = 2.866 days was found within the ![Graphic][120] (Day 1) to ![Graphic][121] (Day 17) data, setting the *EIM* t = 0 point. Then *tI* = (17 — 2.866) = 14.134 days, while ![Graphic][122] (Day 114) gives *tF* = (114 — 2.866) = 111.134 days, along with: ![Formula][123] ![Formula][124] ![Formula][125] ![Formula][126] This X-axis *t* = 0 point is a hypothetical *EIM* pandemic starting point, if *Social Distancing* had been operating throughout the initial CoVID-19 period. **Figure 3** has a predicted CoVID-19 pandemic peak of ~5, 217 / day at *t* = 29.134 days on 3/25/2020, with ~243,100 total cases at the pandemic end. This datafit has > 4*X* error reduction over the *δo* = 0 case, as summarized next: ![Formula][127] The *EIM* with *ZA*[*t*] and *δo* = 0.0336 gives a *ρA*[*t*] curve that is in excellent agreement with the **Fig. 3** measured *ρdata*[*t*] data. Comparing *NA*[*t*] and *Ndata*[*t*] in **Fig. 4** also shows an excellent match over the whole *Logarithmic Y-axis* range, used in the *rms* error minimization. The **Figs. 3-4** *κ* = 0 results were then compared to the *κ* =1 case, using Zb[t] and Nb(Zb) of Eqs. [3.2]-[3.3]. The *κ* = 1 case has a *Social Distancing* factor that gradually turns off the *EIM* exponential decay. The resulting *rms* error best fits converged to *δo* → 0, as follows: ![Formula][128] This *κ* = 1 result shows that the process giving rise to exponential tails in the *ρdata*[*t*] for a pandemic *Fast Shutoff* is independent of the *Social Distancing* process that gradually lengthens the *tdbl doubling time*, as measured by *γo*. Comparing the *κ* = {0,1} cases shows that the pure exponential tail with *ZA*[*t*] = *ZB*[*t*]*κ*=0 matches the *ρdata*[*t*] data best. ## 6 *Enhanced Orthogonal Functions* for Italy Any monotonic *Z*[*t*] can convert measured *Ndata*[*t*] data into *Ndata*(*Z*). Using the Eq. [2.6b] *ZA*[*t*] → *Z*[*t*] in Eqs. [1.7a]-[1.7d] extends the *EIM* into an *EOF* model, where Eq. [1.9] gives: ![Formula][129] ![Formula][130] The Eq. [5.2] *δo* ≠ 0 entries and Eq. [5.1a]-[5.1d] boundary conditions give: ![Formula][131] ![Formula][132] These data-driven *ZA*[*t*] limits are used next, along with *Z*[*t*] → *ZA*[*t*] in: ![Formula][133] ![Formula][134] where *m* = {0, *MF*} sets how many terms are in the Eq. [6.3a] series. Generally *MF* = 2 is used here. The *Lm*(*Z*) are the *Laguerre Polynomials*, with the first few *Lm*(*Z*) being: ![Formula][135] ![Formula][136] ![Formula][137] ![Formula][138] ![Formula][139] ![Formula][140] Some important properties of the *Laguerre Polynomials* are: ![Formula][141] ![Formula][142] ![Formula][143] ![Formula][144] ![Formula][145] Here Eq. [6.5a] defines an *orthogonal function set*. The “**n**!” (**n***-factorial*) in Eq. [6.5d], for **n** an integer, is defined as the product: ![Formula][146] ![Formula][147] where *factorials* with negative integers are not allowed. For *MF* > 2, the following equations developed by Watson**10**, and improved by Gillis and Weiss**11**, helps in evaluating Eq. [6.3b]: ![Formula][148] ![Formula][149] ![Formula][150] ![Formula][151] where ALL terms in the Eq. [6.7c] sum for *n* = {0, (*r* + *s*)} have an implicit requirement that all negative *factorial* arguments are excluded**10–11**. Since the *N*(*Z*) of Eq. [6.3a] has {*gm*; *m* = (0, *MF*)}, and *N*(*Z*) also appears in each *gn*-equation of Eq. [6.3b], how to determine each *gm* by itself, can be done as follows. First define: ![Formula][152] ![Formula][153] When the *Ndata*(*Z*) is comprised of *j* = {1, 2, …*J*} discrete values between ![Graphic][154] with each *Zj* having an ![Graphic][155] value, the Eq. [6.8a] integral needs to be replaced by a sum. Let *Z* = *Z*1 and *Z*J+1 = *Zj*, the *Qn* replacement for Eq. [6.8a] is then: ![Formula][156] ![Formula][157] Eq. [6.3b] can then be re-written as a 3 × 3 matrix **M**3, which relates a data-driven ![Graphic][158]-vector to a resultant ![Graphic][159]-vector: ![Formula][160] ![Formula][161] ![Formula][162] where (**M**3)−1 is the matrix inverse of **M**3. When ![Graphic][163], this **M**3 becomes the Identity Matrix. The following *Km.n*(*Z*) integrals set *Km,n*: ![Formula][164] ![Formula][165] The *Km.n*(*Z*) integrals can be determined using Eq. [6.5c], which gives: ![Formula][166] ![Formula][167] ![Formula][168] ![Formula][169] ![Formula][170] ![Formula][171] The extract {*g*, *g*1, *g*2}, the 3 × 3 symmetric **M**3 matrix needs inversion: ![Formula][172] ![Formula][173] ![Formula][174] which determines {*g*, *g*1, *g*2} from the {*Q*, *Q*1, *Q*2} data. A best-fit *N*(*Z*) for *Z* = {0+,∞−} results, along with an equivalent fit for *R*(*Z*) using Eq. [1.7d]. Instead of having to find the best {*g*, *g*1, *g*2} triplet, one could find the best ![Graphic][175] by just using using {*Q*, *Q*1} and an **M**2 sub-matrix; or one could find the best ![Graphic][176] by itself by just using {*Q*} and an **M**1 sub-matrix: ![Formula][177] ![Formula][178] ![Formula][179] ![Formula][180] ![Formula][181] ![Formula][182] Once the {*gm*; *m* = (0, *MF*)} constants are found and used in Eq. [1.8], its *c* value provides the new *EOF* estimate for the predicted total number of CoVID-19 cases at the pandemic end, refining the initial Eq. [5.2] *N*[*t* → ∞] *EIM* value. ## 7 *EOF* Model Results for Italy The *EOF* model starts with the *EIM* of Eq. [2.6b] using *ZA*[*t*], and the *bing.com* Italy data**9**, which gives: ![Formula][183] ![Formula][184] ![Formula][185] ![Formula][186] ![Formula][187] via Eqs. [5.1a]-[5.1d], [5.2], and [6.2a]-[6.2b]. For these ![Graphic][188] values, with (*tF* − *tI*) = 97 days, the **M**3 matrix of *Km,n* entries, via Eq. [6.8b], is: ![Formula][189] It has a rather small det[**M**3] = 0.04024218 value, with an inverse of: ![Formula][190] A convolution of *Lm*(*ZA*) functions with the measured ![Graphic][191] dataset vector of Eqs. [6.9a]-[6.9b], along with the above (**M**3)−1, gives this final ![Graphic][192] –vector12: ![Formula][193] determining the constants for *N*(*ZA*) in Eq. [1.7a]. The coefficients for *R*(*ZA*), which sets the predicted number of daily new CoVID-19 cases, are: ![Formula][194] determining the constants needed for *R*(*ZA*) in Eq. [1.7b]. Using these {*g*, *g*1, *g*2} values along with Eq. [1.8] gives: ![Formula][195] as a new predicted total number of CoVID-19 cases at the pandemic end for the *EOF* model, which is a ~4.17% or 10, 139 increase in the number of cases, compared to the *EIM* value of Eq. [7.1c]. Using Eq. [6.1b] for *ZA*[*t*], and substituting the Eq. [7.3] ![Graphic][196] values into Eq. [1.7b] gives *R*(*ZA*). The *ρ*[*t*] in Eq. [6.1a] is derived from *R*(*ZA*) using Eq. [1.7b], with the resulting *EOF ρ*[*t*] plotted in **Figure 5**, along with the *t* > *tI* raw data for the daily new CoVID-19 cases. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/11/2020.08.07.20169904/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/08/11/2020.08.07.20169904/F5) Figure 5: *Enhanced Orthogonal Functions (EOF)* for ITALY CoVID-19 data to 6/15/20. Orthogonal functions give additional parameters for further datafit improvement, using a 3-term series (*MF* = 2). Result shows *EIM*, by itself, provides most of the improvement. The **Figure 5** *EOF* model also gives a *t* > *tI* extrapolation, which shows what the combination of processes would look like, if they all had been operating continuously from the CoVID-19 pandemic start. The companion *N*[*t*] analytic result, along with the *t* > *tI* raw data for the total number of CoVID-19 cases is show in **Figure 6**. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/11/2020.08.07.20169904/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/08/11/2020.08.07.20169904/F6) Figure 6: *Enhanced Orthogonal Functions (EOF)* for ITALY CoVID-19 data to 6/15/20. Data for total number of CoVID-19 cases versus time, compared with 3-term *EOF* model for *N*(*ZA*[*t*]) shows excellent match after *Social Distancing* start at *N* = 10, 149. Comparing the size and timing of the *ρ*[*t*] pandemic peak, and its Day 200 value, between the *EIM* (**Figs. 3-4**) and *EOF* model (**Figs. 5-6**), gives: ![Formula][197] showing the *EOF* model predicts more cases total and more daily new CoVID-19 cases at Day 200, as well as modifying the pandemic peak predictions. While the above analysis used *MF* = 2, with the Eq. [7.4] ![Graphic][198] setting the best {*g*, *g*1, *g*2} values, this *EOF* model also provides estimates for the simpler *MF* = {0, 1} cases, as outlined by Eqs. [6.14a]-[6.14f]. For *MF* = 1, the best two ![Graphic][199] values were gotten by only using {*Q*, *Q*1} and an **M**2 sub-matrix of **M**3. For *MF* = 0, the best ![Graphic][200] by itself is derived by using {*Q*} and the **M**1 sub-matrix. These alternative estimates give: ![Formula][201] These additional calculations give the following progression of estimates for *N*[*t* → ∞], which is the final number of CoVID-19 cases at the pandemic end: ![Formula][202] based on increasing the number of data fitting parameters used with the original data. This summary shows the *N*[*t* → ∞] projections are fairly stable, with an average and 1σ standard deviation: ![Formula][203] among these different calculations, where 1σ is ~2.06% of the overall average. Comparing the results among **Figs. 1-6** also highlights these items: 1. All *ρ*[*t*] functions have a sharp rise, and slower decreasing tail. The fastest changing *ρdata*[*t*] tail, as in the Italy CoVID-19 *(Fast Shutoff*) data, was successfully modeled by adding in an exponential term, as in Eq. [6.1b]. 2. The datafits in **Fig. 4** and **Fig. 6** show that the extra parameters in the *EIM* and *EOF* model fits the *ρdata*[*t*] shape progressively better. 3. The *EOF* model shows only relatively small changes of ~2.06% in the *N*[*t* → ∞] limits (Eq. [7.10]), as an estimate of uncertainty in the *EIM*. 4. The *Enhanced Initial Model (EIM)* function captures much of the progression to a pandemic *Fast Shutoff*, as seen in the Italy data. The *ρ*[*t*] tail may still differ from these predictions, due to factors such as: 1. The CoVID-19 dynamics may change in the long-term low *ρ*[*t*] regime; 2. A “second wave” or multiple waves of *ρ*[*t*] resurgence may occur, which are beyond the scope of this CoVID-19 pandemic modeling. ## 8 Summary and Conclusions The early stages of the CoVID-19 coronavirus pandemic starts off with a nearly exponential rise in the number of infections with time. Defining *N*[*t*] as the expected total number of CoVID-19 cases vs time, this basic function: ![Formula][204] ![Formula][205] ![Formula][206] models *Social Distancing* effects as a gradual lengthening of the pandemic growth *doubling time*, which enables pandemic shutoff with only a small population of infected persons. The Eq. [8.1b] *Zo*[*t*] was our *Initial Model***1**, and gives a CoVID-19 *Slow Shutoff* with a long-term *ρo*[*t*]~[1/t2] tail. Previously we showed**1–2** that this *Zo*[*t*] model fits many *Ndata*[*t*] and *ρdata*[*t*] cases. However, some data had a CoVID-19 *Fast Shutoff*, with a *ρdata*[*t*] ~ [exp(−*δot*)] exponential tail, such as in Italy**9**, where a Gaussian tail *ρ*[*t*] ~ [exp(−*qot*2)] would have decreased too quickly. An *Enhanced Initial Model (EIM)* was developed here, using this *ZA*[*t*] function: ![Formula][207] ![Formula][208] ![Formula][209] ![Formula][210] We also examined if the exp(−*δot*) exponential decay could also be subject to a *Slow Shutoff*, giving exp[−*δot* / (1 + *γot*)] instead of exp(−*δot*), but that did not match the Italy data. To allow more data fitting parameters beyond just {*KA*, *γo*, *δo*}, an orthogonal function method was developed**2**: ![Formula][211] ![Formula][212] ![Formula][213] ![Formula][214] which is applicable to a generic *Z*[*t*] function, with *N*[*t*] = *N*(*Z*[*t*]), where *Z*[*t*] → *Zo*[*t*] and *Z*[*t*] → *ZA*[*t*] are special cases. Larger *MF* with more {*Lm*(*Z*); *m* = (0, + *MF*)} terms can match almost any *arbitrary function*, enabling fits to a variety of *N*[*t*] and *ρ*[*t*] shapes. The {*gm*; *m* = (0, + *MF*)} are constants determined from each dataset. The *Lm*(*Z*) are the *Laguerre Polynomials*, with several important properties given in Eqs. [6.4a]-[6.5e]. Using *ZA*[*t*] in Eqs. [8.3a]-[8.3d] results in this *Enhanced Orthogonal Function [EOF]* model, which is applicable to both *Slow* or *Fast Shutoff* CoVID-19 pandemic data. The *ρ*[*t*] expected number of daily new CoVID-19 cases is: ![Formula][215] ![Formula][216] Methods were developed to derive the {*KA*, *γo*, *δo*} values, and to determine the {*gm*; *m* = (0, + *MF*)} and {*cm*; *m* = (0, + *MF*)} constants from data. Whereas our *Initial Model* and *EIM* were *MF* = 0 cases, the *MF* = 2 case was used here to examine the Italy CoVID-19 data, as an *EOF* model example. The *bing.com* data for Italy up to ~6/15/2020 was then analyzed, with **Figures 3-6** giving the new Italy results. Both the *EIM* and the *EOF* model provided good datafits, giving similar *N*[*t* → ∞] results for the final number of CoVID-19 pandemic cases, differing by only ~2% at the 1σ level. The *ρ*[*t*] post-peak behavior best indicates if a *δo* ≠ 0 model (CoVID-19 pandemic *Fast Shutoff*) is applicable. The *δo* ≠ 0 case likely is a second *Social Distancing* process, that operates along with, but is independent of the gradual pandemic *doubling time* changes. That *doubling time* change gives rise to a CoVID-19 pandemic *Slow Shutoff* (*γo* ≠ 0), and that process still operates concurrently with the *δo* ≠ 0 CoVID-19 pandemic *Fast Shutoff*. This analysis shows a wide variety of CoVID-19 data can be modeled using {*KA*, *γo*, *δo*, *toffset*} as parameters, covering: (I) an exponential rise at CoVID- 19 pandemic start; (II) a gradual lengthening of *doubling times* for a pandemic *Slow Shutoff*; and (III) an exponential decay for pandemic *Fast Shutoffs*. ## Data Availability All data used is in the public domain or was on the CoVID-19 website maintained by bing.com. * Received August 7, 2020. * Revision received August 7, 2020. * Accepted August 11, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## 10 References 1. 1.[https://www.MedRxiv.org/content/10.1101/2020.05.04.20091207v1](https://www.MedRxiv.org/content/10.1101/2020.05.04.20091207v1), [https://doi.org/10.1101/2020.05.04.20091207](https://doi.org/10.1101/2020.05.04.20091207), “Initial Model for the Impact of Social Distancing on CoVID-19 Spread”, Genghmun Eng. 2. 2.[https://www.MedRxiv.org/content/10.1101/2020.06.30.20143149v1](https://www.MedRxiv.org/content/10.1101/2020.06.30.20143149v1), [https://doi.org/10.1101/2020.06.30.20143149](https://doi.org/10.1101/2020.06.30.20143149), “ Orthogonal Functions for Evaluating Social Distancing Impact on Covid-19 Spread “, Genghmun Eng 3. 3.[https://www.medrxiv.org/content/10.1101/2020.03.27.20043752v1](https://www.medrxiv.org/content/10.1101/2020.03.27.20043752v1), “Forecasting COVID-19 impact on hospital bed-days, ICU-days, ventilator-days and deaths by US state in the next 4 months”, IHME COVID-19 Health Service Utilization Forecasting Team. 4. 4.[https://www.geekwire.com/2020/univ-washington-epidemiologists-predict-80000-covid-19-deaths-u-s-july/](https://www.geekwire.com/2020/univ-washington-epidemiologists-predict-80000-covid-19-deaths-u-s-july/) “Univ. of Washington researchers predict 80,000 COVID-19 deaths in U.S. by July”, Alan Boyle, GeekWire, March 26, 2020. 5. 5.[https://www.yahoo.com/finance/news/coronavirus-modelers-raise-projected-u-041641553.html](https://www.yahoo.com/finance/news/coronavirus-modelers-raise-projected-u-041641553.html), “ Coronavirus modelers raise projected U.S. death toll and lengthen state-by-state recovery timeline”, Alan Boyle, GeekWire, April 27, 2020. b 6. 6.[https://covid19.healthdata.org](https://covid19.healthdata.org), update of 29 April 2020. 7. 7.[http://www.healthdata.org/covid/updates](http://www.healthdata.org/covid/updates) “COVID-19: What’s New for May 4, 2020: Updated IHME COVID-19 projections: Predicting the Next Phase of the Epidemic”, IHME COVID-19 Health Service Utilization Forecasting Team. 8. 8.[https://finance.yahoo.com/news/pandemic-projection-puts-u-death-220824741.html](https://finance.yahoo.com/news/pandemic-projection-puts-u-death-220824741.html), “New pandemic projection puts U.S. death toll at nearly 135,000, due to less social distancing”, Alan Boyle, GeekWire, May 4, 2020. 9. 9.[www.bing.com/covid](http://www.bing.com/covid): ’Bing COVID-19 Tracker’, and [https://www.bing.com/covid?form=CPVD07](https://www.bing.com/covid?form=CPVD07). 10. 10. G. N. Watson, “A Note of the Polynomials of Hermite and Laguerre”, Journal of the London Mathematical Society, 13(1938), pp. 29-32. 11. 11. J. Gillis and G. Weiss, “Products of Laguerre Polynomials”, Math. Comput., 14(69), Jan. 1960, pp. 60-63 12. 12.For correct results, calculations done by swp55.exe requires no commas in the integer matrix or vector entries, when using: ScientificWorkplaceTM\_Compute\_EvaluateNumerically. It is due to swp55.exe also accepting European notation where 10,149, 10.149. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/graphic-1.gif [5]: /embed/graphic-2.gif [6]: /embed/inline-graphic-4.gif [7]: /embed/graphic-3.gif [8]: /embed/inline-graphic-5.gif [9]: /embed/inline-graphic-6.gif [10]: /embed/graphic-4.gif [11]: /embed/graphic-5.gif [12]: /embed/inline-graphic-7.gif [13]: /embed/inline-graphic-8.gif [14]: /embed/inline-graphic-9.gif [15]: /embed/inline-graphic-10.gif [16]: /embed/inline-graphic-11.gif [17]: /embed/graphic-6.gif [18]: /embed/inline-graphic-12.gif [19]: /embed/inline-graphic-13.gif [20]: /embed/inline-graphic-14.gif [21]: /embed/inline-graphic-15.gif [22]: /embed/inline-graphic-16.gif [23]: /embed/inline-graphic-17.gif [24]: /embed/inline-graphic-18.gif [25]: /embed/inline-graphic-19.gif [26]: /embed/inline-graphic-20.gif [27]: /embed/graphic-7.gif [28]: /embed/graphic-8.gif [29]: /embed/graphic-9.gif [30]: /embed/graphic-10.gif [31]: /embed/graphic-11.gif [32]: /embed/graphic-12.gif [33]: /embed/graphic-13.gif [34]: /embed/graphic-14.gif [35]: /embed/graphic-15.gif [36]: /embed/graphic-16.gif [37]: /embed/inline-graphic-21.gif [38]: /embed/graphic-17.gif [39]: /embed/graphic-18.gif [40]: /embed/graphic-19.gif [41]: /embed/graphic-20.gif [42]: /embed/inline-graphic-22.gif [43]: /embed/inline-graphic-23.gif [44]: /embed/graphic-21.gif [45]: /embed/graphic-22.gif [46]: /embed/inline-graphic-24.gif [47]: /embed/inline-graphic-25.gif [48]: /embed/inline-graphic-26.gif [49]: /embed/inline-graphic-27.gif [50]: /embed/inline-graphic-28.gif [51]: /embed/inline-graphic-29.gif [52]: /embed/inline-graphic-30.gif [53]: /embed/graphic-23.gif [54]: /embed/graphic-24.gif [55]: /embed/graphic-25.gif [56]: /embed/graphic-26.gif [57]: /embed/graphic-29.gif [58]: /embed/graphic-30.gif [59]: /embed/graphic-31.gif [60]: /embed/graphic-32.gif [61]: /embed/graphic-33.gif [62]: /embed/graphic-34.gif [63]: /embed/graphic-35.gif [64]: /embed/graphic-36.gif [65]: /embed/graphic-37.gif [66]: /embed/graphic-38.gif [67]: /embed/graphic-39.gif [68]: /embed/graphic-40.gif [69]: /embed/graphic-41.gif [70]: /embed/graphic-42.gif [71]: /embed/graphic-43.gif [72]: /embed/graphic-44.gif [73]: /embed/graphic-45.gif [74]: /embed/graphic-46.gif [75]: /embed/graphic-47.gif [76]: /embed/graphic-48.gif [77]: /embed/graphic-49.gif [78]: /embed/graphic-50.gif [79]: /embed/graphic-51.gif [80]: /embed/graphic-52.gif [81]: /embed/graphic-53.gif [82]: /embed/graphic-54.gif [83]: /embed/graphic-55.gif [84]: /embed/graphic-56.gif [85]: /embed/graphic-57.gif [86]: /embed/graphic-58.gif [87]: /embed/graphic-59.gif [88]: /embed/graphic-60.gif [89]: /embed/graphic-61.gif [90]: /embed/graphic-62.gif [91]: /embed/inline-graphic-31.gif [92]: /embed/graphic-63.gif [93]: /embed/graphic-64.gif [94]: /embed/graphic-65.gif [95]: /embed/graphic-66.gif [96]: /embed/inline-graphic-32.gif [97]: /embed/graphic-67.gif [98]: /embed/graphic-68.gif [99]: /embed/graphic-69.gif [100]: /embed/inline-graphic-33.gif [101]: /embed/inline-graphic-34.gif [102]: /embed/inline-graphic-35.gif [103]: /embed/graphic-70.gif [104]: /embed/graphic-71.gif [105]: /embed/inline-graphic-36.gif [106]: /embed/graphic-72.gif [107]: /embed/graphic-73.gif [108]: /embed/inline-graphic-37.gif [109]: /embed/inline-graphic-38.gif [110]: /embed/inline-graphic-39.gif [111]: /embed/inline-graphic-40.gif [112]: /embed/inline-graphic-41.gif [113]: /embed/inline-graphic-42.gif [114]: /embed/inline-graphic-43.gif [115]: /embed/inline-graphic-44.gif [116]: /embed/inline-graphic-45.gif [117]: /embed/inline-graphic-46.gif [118]: /embed/inline-graphic-47.gif [119]: /embed/inline-graphic-48.gif [120]: /embed/inline-graphic-49.gif [121]: /embed/inline-graphic-50.gif [122]: /embed/inline-graphic-51.gif [123]: /embed/graphic-76.gif [124]: /embed/graphic-77.gif [125]: /embed/graphic-78.gif [126]: /embed/graphic-79.gif [127]: /embed/graphic-80.gif [128]: /embed/graphic-81.gif [129]: /embed/graphic-82.gif [130]: /embed/graphic-83.gif [131]: /embed/graphic-84.gif [132]: /embed/graphic-85.gif [133]: /embed/graphic-86.gif [134]: /embed/graphic-87.gif [135]: /embed/graphic-88.gif [136]: /embed/graphic-89.gif [137]: /embed/graphic-90.gif [138]: /embed/graphic-91.gif [139]: /embed/graphic-92.gif [140]: /embed/graphic-93.gif [141]: /embed/graphic-94.gif [142]: /embed/graphic-95.gif [143]: /embed/graphic-96.gif [144]: /embed/graphic-97.gif [145]: /embed/graphic-98.gif [146]: /embed/graphic-99.gif [147]: /embed/graphic-100.gif [148]: /embed/graphic-101.gif [149]: /embed/graphic-102.gif [150]: /embed/graphic-103.gif [151]: /embed/graphic-104.gif [152]: /embed/graphic-105.gif [153]: /embed/graphic-106.gif [154]: /embed/inline-graphic-52.gif [155]: /embed/inline-graphic-53.gif [156]: /embed/graphic-107.gif [157]: /embed/graphic-108.gif [158]: /embed/inline-graphic-54.gif [159]: /embed/inline-graphic-55.gif [160]: /embed/graphic-109.gif [161]: /embed/graphic-110.gif [162]: /embed/graphic-111.gif [163]: /embed/inline-graphic-56.gif [164]: /embed/graphic-112.gif [165]: /embed/graphic-113.gif [166]: /embed/graphic-114.gif [167]: /embed/graphic-115.gif [168]: /embed/graphic-116.gif [169]: /embed/graphic-117.gif [170]: /embed/graphic-118.gif [171]: /embed/graphic-119.gif [172]: /embed/graphic-120.gif [173]: /embed/graphic-121.gif [174]: /embed/graphic-122.gif [175]: /embed/inline-graphic-57.gif [176]: /embed/inline-graphic-58.gif [177]: /embed/graphic-123.gif [178]: /embed/graphic-124.gif [179]: /embed/graphic-125.gif [180]: /embed/graphic-126.gif [181]: /embed/graphic-127.gif [182]: /embed/graphic-128.gif [183]: /embed/graphic-129.gif [184]: /embed/graphic-130.gif [185]: /embed/graphic-131.gif [186]: /embed/graphic-132.gif [187]: /embed/graphic-133.gif [188]: /embed/inline-graphic-59.gif [189]: /embed/graphic-134.gif [190]: /embed/graphic-135.gif [191]: /embed/inline-graphic-60.gif [192]: /embed/inline-graphic-61.gif [193]: /embed/graphic-136.gif [194]: /embed/graphic-137.gif [195]: /embed/graphic-138.gif [196]: /embed/inline-graphic-62.gif [197]: /embed/graphic-141.gif [198]: /embed/inline-graphic-63.gif [199]: /embed/inline-graphic-64.gif [200]: /embed/inline-graphic-65.gif [201]: /embed/graphic-142.gif [202]: /embed/graphic-143.gif [203]: /embed/graphic-144.gif [204]: /embed/graphic-145.gif [205]: /embed/graphic-146.gif [206]: /embed/graphic-147.gif [207]: /embed/graphic-148.gif [208]: /embed/graphic-149.gif [209]: /embed/graphic-150.gif [210]: /embed/graphic-151.gif [211]: /embed/graphic-152.gif [212]: /embed/graphic-153.gif [213]: /embed/graphic-154.gif [214]: /embed/graphic-155.gif [215]: /embed/graphic-156.gif [216]: /embed/graphic-157.gif