Epidemic Indicators

Lead: Peter Julian Cayton, PhD

Collaborators: Robert Leong, MSc, PhD(c), Dominic Ligot, Jan Gil Sarmiento

Description:
These are time series data on the developments of the COVID-19 in the country. It provides up-to-date information on the trends in infections, movements as reaction to policies, and targets for flattening the curve. More importantly, it can provide valuable insights about the community control of COVID-19.

Value to the LGU:

If provided with data from the LGUs, we can provide data that will give insights on the state of COVID-19 in the local area. In particular, the metrics to be displayed can provide a real-time evaluation of how transmissible is COVID-19 in the locality subject to control measures and testing practices imposed as data are collected. The key metrics currently available at the moment (and their utility):

• Effective (real-time) reproduction number - measures the current transmissibility, on average, of a single infection in a community (i.e., how many people, on average, will be infected by a single infection under the prevailing community conditions). This is technically defined as the average number of secondary cases that each infected individual would infect if the conditions remained as they were instantaneously, and is computed as (Cori et al., 2013):

  $$R_t = {{E[\Delta]} \over {\sum \Delta I_{t-s} \omega_s}}$$

  $R_t = {{E[\Delta]} \over {\sum \Delta I_{t-s} \omega_s}}$

  • $ {E[\Delta]} $ := expected number of new infections in day t
  • $ {\sum \Delta I_{t-s} \omega_s} $ := weighted average of the number of secondary cases caused by previous new infections, weighted by the probability distribution of the serial interval

  The serial interval (SI) is the time between the onset of symptom of the first infected to the onset of symptoms of the secondary cases. Assumed in our calculations is that SI follows a lognormal distribution with mean 4.8 days and standard deviation 2.3 days (from Nishiura et al., 2020).

• Delay-adjusted confirmed case fatality (CFR) rate - measures the risk of death of a confirmed case in real-time adjusted for the expected time-delay between case confirmation date and actual death (i.e. to reflect the average time it takes between reporting of confirmation of cases up to reporting of death if it happens). This is computed as (Nishiura et al., 2010):

  $$ CFR_t = {{D_t} \over {{{\sum \Delta I_u F (t-u) } \over C_t } \times (C_t - R_t - D_t) + R_t + D_t } } $$

  $ CFR_t $ = $ D_t $ / (($ {\sum \Delta I_u F (t-u) } $ / $ C_t $ ) x $ (C_t - R_t - D_t) + R_t + D_t) $

  • $ D_t $ := total number of deaths at time t
  • $ R_t $ := total number of recoveries at time t
  • $ C_t $ := total number of confirmed cases at time t
  • $ {{\sum \Delta I_u F (t-u) } \over C_t } $ := delay-adjustment factor applied to cases not yet reported to be closed (either recovered or died) at time t which provides how many cases are expected to have been closed by t considering they have reported as confirmed in time t - u (as not all closed cases, especially recoveries, may have been completely reported)

  The delay distribution F in our calculations assumed a gamma distribution with mean 10.1 days and standard deviation 5.4 days (Shim, et. al. 2020).

• Crude Recovery Rate - is the ratio of the total cumulative recoveries $ R_t $ and the cumulative number of infected cases $ I_t $.

  $$ RecovRate = {R_t \over I_t} $$

  RecovRate = $ R_t $ / $ I_t $

• Percent Change $ p_t $ - is the daily growth in cases, deaths, and recoveries. For $ Y_t $ time series, the percent change series $ p_{t,Y} $ is

  $$ p_{t,Y} = {{Y_t - Y_{t-1}} \over {Y_{t-1}}} $$

  $ p_{t,Y} $ = $ Y_t - Y_{t-1} $ / $ Y_{t-1} $

• Implied Doubling Rate $ d_t $ - is the defined length of time for the present value of variable to double given the prevailing rate of increase. It is based on the compound accumulation model. Suppose $ A_0 $ is the starting value, $ p $ is the percent change of increase for every time period, and d is the doubling rate, measured in units of the time period. Then, the compound accumulation model for doubling is:

  $$ 2A_0 = {A_0(1+p)^d} $$

  $ 2A_0 $ = $ A_0(1 + p) $

  Solving for $ d $ in terms of $ p $ is:

  $$ d = {{ln(2)} \over {ln(1+p)}} $$

  $ ln(2) $ / $ ln(1 + p) $

  Since there is a time series $ p_{t,Y} $ from time series $ Y_t $, then the implied doubling rate for $ Y_t $ is:

  $$ d_{t,Y} = {{ln(2)} \over {ln(1+p_{t,Y})}} $$

  $ d_{t,Y} $ = $ ln(2) / ln(1 + p_{t,Y}) $

  The two time series $ p_{t,Y} $ and $ d_{t,Y} $ are similar in terms of information, but both are presented for interested readers. Since doubling time is related to growth rate, calculating for equivalent percentage growth from doubling rate:

  $$ p_{t,Y} = {e^{ln(2) \over d_{t,Y}}-1} $$

  $ p_{t,Y} $ = $ {e^{{ln(2)} / d_{t,Y}}-1} $

---

CO-INFORM Risk Scoring Dashboard

Lead: Michael Promentilla, PhD and Jomar Rabajante, PhD

Collaborators: Geminn Louis Apostol, MD, MBA; Dominic Ligot, April Anne Tigue, MSc

Description:
The COVID-19 risk index is defined as a function of the following dimensions namely: Hazard & Exposure, Vulnerability, and Resilience. This conceptual framework is adapted from the InfoRM model developed by the European Commission’s Joint Research Commission (JRC), which is a global multi-hazard disaster risk assessment tool to identify countries at risk of disaster and humanitarian crisis. The risk concepts used in the model are based on several publications in scientific literature and considers the three dimensions of risk: Hazards & Exposure, Vulnerability, and Lack of Coping Capacity. However, this index is non-specific and some of the indicators may not be available or even relevant to the needs of decision-makers in the context of the COVID-19 epidemic in the country. Thus, we propose a conceptual framework to measure the COVID-19 risk index as shown in Figure 1. This approach is modular and allows for a simple and transparent calculation of epidemic risk using a composite index methodology. Each dimension includes different categories where each category can be broken down by a reliable set of indicators. Geometric aggregation is used to compute the composite index or risk score as shown by the following equation:

• $ {H^{w1}} $ := Hazard and Exposure
• $ {V^{w2}} $ := Vulnerability
• $ {L^{w3}} $ := Lack of Resilience

In theory, there is no risk from the COVID-19 epidemic, if there is no exposure, no matter how severe the hazard event is. If there are physical exposure and physical vulnerability, the “hard” risk can be computed and it is considered hazard dependent. The risk is also very low if the person or community is not vulnerable or the resilience to cope and recover is ideal. The “soft” risk can be computed from this second dimension using the concept of vulnerability due to the fragility of the socio-economic system including the susceptibility associated with the low level of awareness, nutritional and health status. These are the social determinants of health and are hazard independent. Likewise, resilience is operationalized and defined by physical infrastructure, health systems capacity, and also institutional and management capacity. Conceptually, better epidemic management to absorb, recover and adapt means higher resilience capacity, which reduces the level of risk from the vulnerability and exposure from this hazard. Likewise, the lack of resilience translates to higher risk overall.

Value to the LGU:

As a novel risk reduction exercise concurrently developed with the pandemic response, this tool offers numerous utility for decision-makers in the Philippines.

• Risk index comparisons and visualization may potentially accelerate not only national but also district-level actions to the developing pandemic. And when the situation has settled, the risk index could be used as a guide for reducing risks and drafting preparedness plans to mitigate future impacts of such health emergencies.
• Information from the CO-INFoRM index could nudge multi-sectoral actions to aptly respond to the complexity and the nuances of the crisis.
• Proactive communication, such as through a regularly updated visualization, strategically aligns the perceived expectations of the public and the media closer to the actual risks. This communication control mitigates potential negative impressions, ensures accountability, and promotes public trust to preserve an institutions brand equity.

Model Framework:

Notes on the Initial Set of Indicators:

Dimension	Indicator	Description and data processing requirement
$ {H^{w1}} $ Hazard & Exposure	Active Case/Outbreak threshold (Number)	The ratio between active infected cases and a number of cases that will potentially lead to an outbreak in the area will be calculated. A lower number means a lower risk. The index will be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization and inversion. The Min value is 0 and the Max value is set to 2. Data were obtained from Epid. Model of Rabajante (R = 2, c = 0.01) for outbreak threshold and DOH for active cases. Sources: Rabajante (2020); DOH
$ {H^{w1}} $ Hazard & Exposure	Projected Number of Cases (Number per 100,000 population)	The projected number of infected cases (population-density adjusted) throughout the epidemic period in the area. A higher number means a higher risk. The data will be first normalized with the population of the region and presented as number of cases per 100,000 population. The index will then be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization. The Min value is 0 cases per 100,000 and Max value is the largest value in the data set. Data were obtained from Epid. Model of Rabajante (R = 2, peak at 25%, base). Sources: Rabajante (2020)
$ {H^{w1}} $ Hazard & Exposure	Projected Number of potential fatalities (Number per 1,000,000 population)	Projected number of fatalities considering the age structure of the population and its corresponding case fatality rates. A higher number means a higher risk. The data will be first normalized with the population of the region and presented as number of cases per million population. The index will then be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization. The Min value is 0% and Max value is the largest value in the data set. Data were obtained from Epid. Model of Rabajante (R = 2, peak at 25%, base). Sources: Rabajante (2020)
$ {H^{w1}} $ Hazard & Exposure	Time-varying Reproductive Number R (Rt)	The average number of secondary cases that each infected individual would infect if the conditions remained as they were instantaneously. A higher number means a higher risk. The index will then be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization. The Min value is 0 and the Max value is set to 2. Data were obtained from COVID-19 Time Series Dashboard of Cayton. Sources: Cayton (2020)
$ {V^{w2}} $ Vulnerability	Poverty incidence (%)	Proportion of poor in the total number of households in the area as a proxy data for the social vulnerability. A higher number means a higher risk. Raw data were obtained from Rabajante. The index is then computed from rescaling the data to a score ranging from 0 to 10 using max-min normalization. The Min value is 0 and Max value is 50%. Data were obtained from Philippine Statistics Authority (Highlights on Household Population, Number of Households, and Average Household Size of the Philippines, 2015 Census of Population) Sources: PSA (2015)
$ {V^{w2}} $ Vulnerability	Elderly Population (%)	The percentage of senior citizen in the population (60 y/o and above). A higher number means a higher risk. The index will then be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization. The Min value is 0% and Max value is the largest value in the data set. Data were obtained from Philippine Statistics Authority (2015 Census of Population). Sources: PSA (2015)
$ {L^{w3}} $ Lack of Resilience	LGU Health expenditure (%)	The amount of LGU budget allocated to health per capita. This serves as proxy data for effective governance to respond against health-related crises. A higher number translates to higher resilience to reduce risk level. The index will then be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization and inversion. The Min value is Php 0 and Max value is the largest value in the data set. Data were obtained from Department of Finance, Bureau of Local Government Finance, Statement of Receipts and Expenditures. Sources: DOF (2018)
$ {L^{w3}} $ Lack of Resilience	Critical care bed demand (number per a million population)	This considers the number of critical care beds. The projected number of critical cared beds needed by healthcare during peak. A lower demand translates to higher resilience to reduce risk level. The data will be first normalized with the population of the region and presented as number of critical care bed per million population. The index will then be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization and inversion. Data were obtained from Epid. Model of Rabajante (R=2, peak at 25%, base). Sources: Rabajante (2020)
$ {L^{w3}} $ Lack of Resilience	Bayanihan Grant (Php per capita)	The amount of financial assistance given to LGU to address COVID-19. This serves as proxy data for effective governance to respond against health-related crises. A higher number translates to higher resilience to reduce risk level. The index will then be computed from rescaling the data to a score ranging from 0 to 100 using max-min normalization and inversion. The Min value is Php 0 and Max value is set to Php 5,236.00 which is based on the 40% allocation for catastrophic payment for health as recommended by WHO. Data were obtained from Department of Budget and Management, Local Budget Circular No. 125. Sources: DOF (2020)