Next Article in Journal
Acknowledgment to the Reviewers of Epidemiologia in 2022
Previous Article in Journal
The Prevalence and Indications of Intravenous Rehydration Therapy in Hospital Settings: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Trendline and Predictive Analysis of the First-Wave COVID-19 Infections in Malta

1
Department of Mechanical Engineering, Faculty of Engineering, University of Malta, MSD 2080 Msida, Malta
2
Department of Naval Architecture, Ocean, and Marine Engineering, University of Strathclyde, Glasgow G4 0LZ, UK
3
Infection Control Department, Mater Dei Hospital, MSD 2090 Msida, Malta
4
Department of Microbiology, Faculty of Medicine and Surgery, University of Malta, MSD 2080 Msida, Malta
*
Author to whom correspondence should be addressed.
Epidemiologia 2023, 4(1), 33-50; https://doi.org/10.3390/epidemiologia4010003
Submission received: 18 October 2022 / Revised: 26 December 2022 / Accepted: 4 January 2023 / Published: 10 January 2023

Abstract

:
Following the first COVID-19 infected cases, Malta rapidly imposed strict lockdown measures, including restrictions on international travel, together with national social distancing measures, such as prohibition of public gatherings and closure of workplaces. The study aimed to elucidate the effect of the intervention and relaxation of the social distancing measures upon the infection rate by means of a trendline analysis of the daily case data. In addition, the study derived a predictive model by fitting historical data of the SARS-CoV-2 positive cases within a two-parameter Weibull distribution, whilst incorporating swab-testing rates, to forecast the infection rate at minute computational expense. The trendline analysis portrayed the wave of infection to fit within a tri-phasic pattern, where the primary phase was imposed with social measure interventions. Following the relaxation of public measures, the two latter phases transpired, where the two peaks resolved without further escalation of national measures. The derived forecasting model attained accurate predictions of the daily infected cases, attaining a high goodness-of-fit, utilising uncensored government-official infection-rate and swabbing-rate data within the first COVID-19 wave in Malta.

1. Introduction

Upon the identification of a novel coronavirus in China in January 2020, infections caused by the acute respiratory syndrome coronavirus 2 (SARS-CoV-2) rapidly spread worldwide [1]. Coronavirus disease 2019 (COVID-19) cases have been reported to be highly infectious and carry substantial morbidity and mortality [2]. Despite precautionary measures promoted by societal institutions, the impact of COVID-19 on healthcare systems has been significant, overwhelming hospital bed capacity and intensive care resources [3,4,5]. In an effort to ‘flatten the curve’, governments have implemented varying public health measures, which, at the most extreme, have involved lockdowns of large cities and countries [6,7]. These measures, however, have incurred substantial costs, and required large-scale logistical efforts to achieve an effective balance between the protection of human health and the national economy [8]. Studies have indicated that control measures have been effective [9,10], with a full lockdown imposition being the most effectual manner to mitigate the spread of the virus [11]. Yet, as the categories of control measures imposed by a government may vary substantially, very distinct cohorts of the population can truly be affected. As a result, analysing the virus spread in relation to discretised imposed measures has been a constraint in epidemiological analyses [12,13].
In order to constitute data-driven decisions, forecasting models by means of distribution fitting are commonly utilised to implement interventions within speculated timespans. These models have provided quantitative infection case projections that allow policymakers to plan appropriate interventions [14]. Of particular importance is short-term hospital resource allocation, such as emergency expansion and adoption of patient beds, hospital wards, or Nightingale hospitals, patient distributions within critical and intensive units, and diversion/suspension of routine hospital activities to permit functionality within rapidly changing trends. In turn, this endeavour requires predictive modelling that is ideally achievable without the need for specialised expertise or computational infrastructure, such that medical professionals may directly attain the data required for decision-making [15]. Specific enforcement measures must be historically justified to be effective on a national scale, yet epidemiological studies have not been able to isolate such measures due to external imported cases.
Dahal et al. [16] investigated transmission dynamics of COVID-19 at national and regional levels by using laboratory-confirmed RT-PCR data. Utilising dynamic phenomenological growth models, 10- and 20-day sequential forecasts of reproduction numbers were performed at weekly intervals over a timespan of two months (March to May 2021). The forecast indicated a declining trend of COVID-19 cases in Nepal, with a linear trajectory of COVID-19 incidence during the first wave and a sub-exponential growth pattern in the second. The study highlighted how computationally light mathematical models can provide reliable short-term projections and public health planning in such epidemic situations. Using similar techniques, Pajuelo-Reyes et al. [17] described the transmission of COVID-19 in the Amazonas region of Peru between March and July 2020.
Olumoyin et al. [18] argued that mitigation measures may impact epidemiological models with constant parameters. As a result, an algorithm was introduced to estimate the time-varying transmission rate for the COVID-19 pandemic in a background of different mitigation measures. Utilising cumulative and daily reported symptomatic infection numbers, the effect on the basic reproduction number of non-pharmaceutical mitigation measures, such as early detection of positive cases, subsequent contact tracing, and implementation of social distancing, was simulated. Using error metrics, the accuracy of the proposed algorithm was demonstrated for Italy, South Korea, the United Kingdom, and the United States.
The epidemiology of COVID-19 was also studied by Fonseca et al. [19] in the state of Sergipe, Brazil. The authors showed a positive correlation between the number of cases and deaths with demographic density. Using social parameters affecting the COVID-19 pandemic in Brazil, Marinho et al. [20] developed models that could be utilised by public managers for effective decision-making.
Ogwara et al. [21] observed the time-varying reproduction number of SARS-CoV-2 in Georgia, USA as state and countywide policies were implemented that enforced and relaxed COVID-19 public health measures. Policy changes were associated with increases or decreases at different time points. The reproduction number was seen to increase following the reopening of schools for in-person instruction. Additionally, Fain and Dobrovolny [22] developed a correlation model for SARS-CoV-2 infections using a complex partial differential equation mathematical model. The model generated dose–response curves for peak viral load, time of viral peak, viral growth rate, infection duration, and area under the viral curve.
Albeit the various studies conducted to establish and forecast the spread of COVID-19 have been impactful, the analyses have either had limitations in establishing the virus spread causation within a geographic area due to imported cases or incomplete data gathering, or had computationally intensive methodologies in forecasting the virus spread, hence requiring high-end computing architecture and highly trained analysts. In contrast, the primary aim of this work, in recognition of the limited imported cases as a result of the international travel ban, was to elucidate the effect of the distinct enforcement and relaxation social distancing measures upon the virus infection rate within the Maltese population by means of a historical trendline-fitting analysis. Additionally, as short-term hospital resource allocation and adaptation of routine hospital activities was critical throughout the pandemic, the secondary aim was developing a computationally efficient, yet accurate, forecasting model for the COVID-19 infection spread, using SARS-CoV-2 positive-case and swab-testing data in Malta, to be solved on commonly utilised office machines and spreadsheet software for data-driven logistics insights and decision-making.
Malta is a Mediterranean island country with a population of approximately 500,000 inhabitants. Its first SARS-CoV-2 isolate was reported on 7 March 2020. Following the initial case, the country experienced an epidemic until the end of July, with a peak in mid-April, in what was the first wave of COVID-19 in the country. Various interventions were introduced by the government to address the challenge, with international and national travel bans being the most substantial. A fundamental component of the national strategy was an exhaustive track-and-trace strategy, which reached a frequency of almost 2000 swabs per 100,000 residents per week by the end of the first wave.

2. Methodologies

2.1. Trendline Analysis

2.1.1. Data Collection

The daily number of new COVID-19 cases, together with the respective number of swab-tests taken, were retrieved from the official website of the Ministry of Health in Malta issued from 21 February 2020 onwards. The key interventions introduced by the Maltese government throughout the first wave were also identified, together with the dates when these interventions were implemented and discontinued. These are summarised in Table 1. It should be noted that, as foreign and local travel had ceased within the island, together with an effective track-and-trace practice, the virus spread data were deemed to be solely related to national cases.

2.1.2. Trendline Functions

In an effort to establish the trend by which the cumulative number of infections varied over the first pandemic wave in Malta, a trendline analysis was implemented. Primarily, a linear function (Equation (1)), was utilised upon the case data. Subsequently, logistic functions, based upon the cumulative distribution functions of the Exponential distribution (Equation (2)) and the Weibull distribution (Equations (3)–(5)), were implemented. The Weibull-based logistic function was developed further such that three distinct functions were coupled, establishing a multi-logistic function.
N c = M c 1 n
N c = M c 1 · 1 e n λ c 1
N c = M c 1 · 1 e n λ c 1 k c 1
N c = M c 1 · 1 e n λ c 1 k c 1 + M c 2 · 1 e n λ c 2 k c 2
N c = M c 1 · 1 e n λ c 1 k c 1 + M c 2 · 1 e n λ c 2 k c 2 + M c 3 · 1 e n λ c 3 k c 3
where N c is the cumulative number of positively diagnosed cases, n is the day number,   M c 1 ,     M c 2 ,     M c 3 are the case-dependent magnitude parameters, λ c 1 ,     λ c 2 ,     λ c 3 are the case-dependent scale parameters, and k c 1 ,   k c 2 ,     k c 3 are the case-dependent shape parameters.
In relation to the cumulative infected cases, a daily infected case trendline analysis was undertaken by implementing the arithmetic time-derivative of the cumulative functions (Equations (6)–(10)).
N c ˙ = M c 1
N c ˙ = M c 1 λ c 1 · e n λ c 1
N c ˙ = M c 1 · k c 1 · n k c 1 1 λ c 1 k c 1 · e n λ c 1 k c 1
N c ˙ = M c 1 · k c 1 · n k c 1 1 λ c 1 k c 1 · e n λ c 1 k c 1 + M c 2 · k c 2 · n k c 2 1 λ c 2 k c 2 · e n λ c 2 k c 2
N c ˙ = M c 1 · k c 1 · n k c 1 1 λ c 1 k c 1 · e n λ c 1 k c 1 + M c 2 · k c 2 · n k c 2 1 λ c 2 k c 2 · e n λ c 2 k c 2 + M c 3 · k c 3 · n k c 3 1 λ c 3 k c 3 · e n λ c 3 k c 3
where N c ˙ is the number of positive cases per day.
In a similar manner, the cumulative number of swab-tests over the first wave was analysed by utilising a linear function (Equation (11)) and a Weibull-based logistic function (Equation (12)).
N s = M s 1 · n + 14
N s = M s 1 · 1 e n + 14 λ s 1 k s 1
where N s is the cumulative number of swab-tests,   M s 1 is the swab-test-dependent magnitude parameter, λ s 1 is the swab-test-dependent scale parameter, and k s 1 is the swab-test-dependent shape parameter.
In relation to the cumulative swab-test number, a daily swab-test trendline analysis was undertaken by implementing the arithmetic time-derivative of the cumulative functions (Equations (13) and (14)).
N s ˙ = M s 1
N s ˙ = M s 1 · k s 1 · n + 14 k s 1 1 λ s 1 k s 1 · e n + 14 λ s 1 k s 1
where N s ˙ is the number of swab-tests per day.

2.2. Predictive Analysis

A prediction model was derived by utilising a logarithmic growth rate equation for the daily diagnosed cases (Equation (15)).
K c n = Δ ln N c n Δ n = ln N c n ln N c n 1
where N c n is the cumulative number of positively diagnosed cases on a given day, n is the considered day, and K c n is the daily positive-case logarithmic growth rate.
From the formula, the necessary model output was N c n , and hence, K c n was required to be determined. To statistically predict K c n , a two-parameter Weibull distribution fit was implemented upon the logarithmic growth rate data of the previous days to establish the Weibull scale and shape parameters. The parameters were rolling parameters as the process was done daily throughout the time period. The scale and shape parameters were then utilised within an inverse cumulative distribution in the form of a quantile function for the Weibull distribution to establish a K c n range (Equation (16)).
K c n = Q p ; λ c ; k c = λ c ln 1 F K c n k c = λ c ln 1 p k c
where p is the occurrence probability, λ c is the case-dependent Weibull scale parameter, and k c is the case-dependent Weibull shape parameter.
In addition, to account for the variation in daily swab-testing, and overcome the assumption of a constant swab-testing-rate, a swab-test factor ( c s ) was introduced and coupled with K c n (Equation (17)).
N c n = exp ln N c n 1 + c s K c n
The swab-test coefficient was established to be the ratio between the swab-test logarithmic growth rate on day n and the average swab-test logarithmic growth rate of the considered prior days (Equation (18)).
c s = K s n 1 t · i = n t n 1 K s i = ln N s n ln N s n 1 1 t · i = n t n 1 ln N s i ln N s i 1
where K s n is the swab-test daily logarithmic growth rate, t is the number of considered prior days, and N s n is the cumulative number of swab-tests on a given day.
By incorporating the prior equations, the predictive model was derived to establish the daily number of infected cases over a time period (Equations (19)–(21)).
N c n = exp ln N c n 1 + c s · λ c ln 1 p k c
N c n = exp ln N c n 1 + K s n 1 t · i = n t n 1 K s i · λ c ln 1 p k c
N c n = N c n 1 · exp K s n 1 t · i = n t n 1 K s i · λ c ln 1 p k c

3. Results

3.1. Trendline Analysis

3.1.1. Positive Infected-Case Function

Applying Equations (1)–(5) to the cumulative dataset attained the data-driven variables detailed in Table 2. The linear function attained the least similarity, whereas the triple Weibull-based logistic function attained the highest similarity, with coefficients of determination (R2) of 0.860 and 0.998, respectively. The functions were graphically superimposed upon the dataset, as illustrated in Figure 1.
Implementing Equations (6)–(10) to the daily infected cases dataset attained the data-driven variables detailed in Table 3. The linear-derivative function attained the least similarity, whereas the triple Weibull-based logistic-derivative function attained the highest similarity, with coefficients of determination ( R 2 ) of 0.0 and 0.566, respectively. The functions were graphically superimposed upon the dataset, as illustrated in Figure 2. Furthermore, as the triple Weibull-based logistic-derivative function attained the highest similarity, portraying a high correspondence, the key dates on which the nationally imposed measures were enforced or relaxed (see Table 1) were additionally incorporated, as illustrated in Figure 3.

3.1.2. Swab-Test Function

Applying Equations (11) and (12) to the cumulative swab dataset attained the data-driven variables detailed in Table 4. The linear function attained the least similarity, whereas the Weibull-based logistic function attained the highest similarity, with coefficients of determination ( R 2 ) of 0.901 and 0.999, respectively. The functions were graphically superimposed upon the dataset, as illustrated in Figure 4.
Implementing Equations (13) and (14) to the daily swab count dataset attained the data-driven variables detailed in Table 5. The linear-derivative function attained the least similarity, whereas the Weibull-based logistic-derivative function attained the highest similarity, with coefficients of determination ( R 2 ) of 0.0 and 0.762, respectively. The functions were graphically superimposed upon the dataset, as illustrated in Figure 5. Furthermore, the infection-rate positivity ratio was established by coupling Equations (10) and (14), illustrated in Figure 6. The key dates on which the nationally imposed measures were enforced or relaxed (see Table 1) were additionally incorporated.

3.2. Predictive Analysis

The forecasting model was applied for the entire first wave (7 March–15 July 2020), and a portion of the second wave (16 July–31 August 2020), to establish the continuity capacity of the model. Employed within a one-day (Figure 7) to a five-day (Figure A1, Figure A2, Figure A3, Figure A4, Figure A5, Figure A6, Figure A7 and Figure A8) prediction framework, the statistical model attained good agreement with the dataset, achieving a global coefficient of determination ( R 2 ) of 0.9995 to 0.9955 between the statistical model median outputs and the actual dataset. Good agreement was also attained for predictions beyond five days. Particularly for the one-day, two-day, and three-day forecasting, solely 4.5%, 7.8%, and 12.3% of the data-points fell outside of the 0th–95th percentile prediction band, respectively. The explicit statistical modelling methodology was therefore deemed to have been validated to a high degree of accuracy.

4. Discussion

Implementing a trendline analysis upon the daily case dataset, as opposed to a sole moving average, permitted superimposing the dates of enforcement and relaxation measures to qualitatively shed light on the effect of the measures along the trend [23]. The infection rate trend increased exponentially throughout the initial days of the pandemic, prior to social measures. Upon the implementation of enforcement measures related to public transport and closure of workplaces, sports facilities, law courts, religious places, and service outlets, the rate-of-change of the daily cases diminished. This acknowledged the effect of the social-distancing enforcement measures, inhibiting the spread of the virus. The subsequent enforcement measures related to the prohibition of public gatherings, and closure of education establishments and non-essential retail/service outlets further decreased the rate-of-change steadily, attaining a peak average infection rate of 11.5 cases per day on 31 March.
The relaxation of public transport measures, re-opening of non-essential retail outlets, and increased size of public gatherings (to a maximum of three persons) on 4 May resulted in the second phase, occurring with an exponential increase in infection rate. The rate, however, peaked rapidly and diminished to a low value within a short timeframe of six days. The third phase was evident following relaxation measures of public gatherings (to a maximum of four persons) and re-opening of service outlets and public places on 22 May. This had a lower rate-of-change than the second phase and rapidly levelled off by 15 July.
The initiation of the second and third phases was consistent with the relaxation of public gathering social distancing measures [24]. As higher numbers of persons were permitted to be in close contact, the possibility of virus transmission increased. Nevertheless, the second and third peaks were much lower and a substantial drop in infection-rate number followed with no new enforcement measures implemented. The reason for this may be attributed to societal diligence in re-orienting itself to effective social distancing [25]. The wearing of masks or visors in shops and on public transport had become obligatory by law on 4 May and may have contributed to the achievement of control in the second and third peaks [26]. A more likely hypothesis, however, may be the effectiveness of test-and-trace that, by this time, had been expanded and consolidated [27]. Approximately 110,000 swab-tests were performed within the first wave, where the peak median swab-tests per day was found to be approximately 1200 on 18 May 2020. Along the second phase time-period, the mean swabbing frequency had increased to over 9000 swabs per week, with peak median positivity ratios of 2.34, 0.95, and 0.34 cases per hundred swabs identified within the first, second, and third phases, respectively. Every positive case was isolated within 24 h of testing with concurrent quarantine of significant contacts. This was possible given the low positivity rate during the second and third phases, which was below 1%. By means of the trendline, it may be argued that, as testing was increasing, had the relaxation measure of 4 May been implemented two weeks later, it would have allowed the infection-rate to level off. Accordingly, the relaxation measure of 1 July would potentially have been implemented earlier.
With regard to the forecasting model, the derivation of an explicit statistical model based upon the logarithmic growth rate was found to be an accurate and computationally feasible methodology, achieving a global R 2 of over 0.99 and a total computation time of less than five seconds on a typical office machine. This methodology attained the logical implementation of solving for the accumulative number of cases on a particular day ( N n ) by determining the logarithmic growth rate ( K n ) via a statistical analysis utilising an inverse cumulative distribution (quantile) function based on a Weibull distribution, together with incorporating a swab-test coefficient ( c n ) to account for the correlation between tests undertaken and positively infected cases.
Incorporating a Weibull distribution fit was advantageous due to its adaptability, permitting the comprehension of both symmetric and non-symmetric distributions whilst interpolating between the exponential distribution and the Rayleigh distribution via a two-parameter implementation. The Weibull cumulative distribution function put forward an explicit function, encompassing two parameters that can be estimated from a dataset, and hence efficiently solvable, in contrast to the application of implicit formulation functionality [28]. This aspect was deemed imperative as high-end statistical approaches tend to lie beyond the statistical knowledge-capacity of medical professionals and the processing power of office machines and commonly utilised spreadsheet software [29,30]. In addition, this model permitted the utilisation of historical data of positive cases and swab-tests, rather than implementing trends that disregard swab-testing correlation. This model may distinctively be utilised for short-to-medium term quantitative risk assessments. Furthermore, week-long facility logistic decisions, such as the number of beds in emergency wards, ventilators, and on-call staff personnel, may be substantiated utilising the model.
By collating and discussing the different statistical modelling and prediction techniques for COVID-19, Yadav and Akhter [31] pinpointed the significance of utilising a single distribution to fit and represent the true virus spread, such that effective data-driven policies may be made. The implemented two-parameter Weibull distribution methodology implemented within this work succeeded in achieving this critical indication. The model, however, is limited by the diminishing capacity of the distribution if left-censored data are utilised [32,33]. Moreover, when furthering the forecasting to a larger timespan, the deviation from the true result may increase substantially as a result of the fundamental approach [34]. Nonetheless, Weibull analyses are typically utilised for medical statistics as the methodology has been acknowledged to sustain accuracy despite an extremely small dataset [35]. In fact, within the context of SARS-CoV-2 statistical analyses, this approach has been notably applied to establish the incubation period of the virus [36,37]. As a result, a rudimentary yet accurate forecasting model was established, encompassing the imperative capacity of accurately predicting positive cases within the termination of one wave and the initiation of another.

5. Conclusions

This study presented a novel statistical model, incorporating swab-testing rates coupled with Weibull-distributed historical data of the SARS-CoV-2 positive-case logarithmic growth rate, to predict the virus infection rate and establish an accurate projection through a numerically explicit framework. The model was validated utilising infection rate data within Malta. Furthermore, an epidemiological elaboration of infection trends was established utilising trendline analyses for the purpose of evaluating social distancing enforcement and relaxation measures upon the virus spread within the population.

Author Contributions

Conceptualization, M.G.B. and M.A.B.; methodology, M.G.B.; software, M.G.B.; validation, M.G.B.; formal analysis, M.G.B. and M.A.B.; investigation, M.G.B. and M.A.B.; resources, M.G.B. and M.A.B.; data curation, M.A.B.; writing—original draft preparation, M.G.B. and M.A.B.; writing—review and editing, M.G.B. and M.A.B.; visualization, M.G.B..; project administration, M.G.B. and M.A.B.; funding acquisition, M.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Figure A1. Two-day infected-case predictive output (first wave).
Figure A1. Two-day infected-case predictive output (first wave).
Epidemiologia 04 00003 g0a1
Figure A2. Two-day infected-case predictive output (second wave segment).
Figure A2. Two-day infected-case predictive output (second wave segment).
Epidemiologia 04 00003 g0a2
Figure A3. Three-day infected-case predictive output (first wave).
Figure A3. Three-day infected-case predictive output (first wave).
Epidemiologia 04 00003 g0a3
Figure A4. Three-day infected-case predictive output (second wave segment).
Figure A4. Three-day infected-case predictive output (second wave segment).
Epidemiologia 04 00003 g0a4
Figure A5. Four-day infected-case predictive output (first wave).
Figure A5. Four-day infected-case predictive output (first wave).
Epidemiologia 04 00003 g0a5
Figure A6. Four-day infected-case predictive output (second wave segment).
Figure A6. Four-day infected-case predictive output (second wave segment).
Epidemiologia 04 00003 g0a6
Figure A7. Five-day infected-case predictive output (first wave).
Figure A7. Five-day infected-case predictive output (first wave).
Epidemiologia 04 00003 g0a7
Figure A8. Five-day infected-case predictive output (second wave segment).
Figure A8. Five-day infected-case predictive output (second wave segment).
Epidemiologia 04 00003 g0a8

References

  1. Zhou, P.; Yang, X.L.; Wang, X.G.; Hu, B.; Zhang, L.; Zhang, W. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Chang, D.; Chang, X.; Tan, K.J.K. The determinants of COVID-19 morbidity and mortality across countries. Sci. Rep. 2022, 12, 5888. [Google Scholar] [CrossRef]
  3. Khan, S.; Siddique, R.; Ali, A.; Bai, Q.; Li, Z.; Li, H. The spread of novel coronavirus has created an alarming situation worldwide. J. Infect. Public Health 2020, 13, 469–471. [Google Scholar] [CrossRef] [PubMed]
  4. Onder, G.; Rezza, G.; Brusaferro, S. Case-Fatality Rate and Characteristics of Patients Dying in Relation to COVID-19 in Italy. JAMA 2020, 323, 1775–1776. [Google Scholar] [CrossRef] [PubMed]
  5. Mareiniss, D.P. The impending storm: COVID-19, pandemics and our overwhelmed emergency departments. Am. J. Emerg. Med. 2020, 38, 1293–1294. [Google Scholar] [CrossRef] [PubMed]
  6. Laupacis, A. Working together to contain and manage COVID-19. CMAJ 2020, 192, 340–341. [Google Scholar] [CrossRef] [Green Version]
  7. Lau, H.; Khosrawipour, V.; Kocbach, P.; Mikolajczyk, A.; Schubert, J.; Bania, J. The positive impact of lockdown in Wuhan on containing the COVID-19 outbreak in China. J. Travel Med. 2020, 27, 1–7. [Google Scholar] [CrossRef] [Green Version]
  8. Nicola, M.; Alsafi, Z.; Sohrabi, C. The socio-economic implications of the coronavirus pandemic (COVID-19): A review. Int. J. Surg. 2020, 78, 185–193. [Google Scholar] [CrossRef]
  9. Brauner, M.; Mindermann, S.; Sharma, M.; Johnston, D.; Salvatier, J.; Gavenčiak, T.; Stephenson, A.B.; Leech, G.; Altman, G.; Mikulik, V.; et al. Inferring the effectiveness of government interventions against COVID-19. Science 2020, 371, eabd9338. [Google Scholar] [CrossRef]
  10. Li, Y.; Campbell, H.; Kulkarni, D.; Harpur, A.; Nundy, M.; Wang, X.; Nair, H. The temporal association of introducing and lifting non-pharmaceutical interventions with the time-varying reproduction number (R) of SARS-CoV-2: A modelling study across 131 countries. Lancet Infect. Dis. 2021, 21, 193–202. [Google Scholar] [CrossRef]
  11. Davies, N.G.; Barnard, R.C.; Jarvis, C.I.; Russell, T.W.; Semple, M.G.; Jit, M.; Edmunds, W.J. Association of tiered restrictions and a second lockdown with COVID-19 deaths and hospital admissions in England: A modelling study. Lancet Infect. Dis. 2020, 21, 482–492. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Z.; Whittington, J.; Yuan, H.-Y.; Miao, H.; Tian, H.; Stenseth, N.C. Evaluating the effectiveness of control measures in multiple regions during the early phase of the COVID-19 pandemic in 2020. Biosaf. Health. 2021, 5, 264–275. [Google Scholar] [CrossRef] [PubMed]
  13. Hsiang, S.; Allen, D.; Annan-Phan, S.; Bell, K.; Bollinger, I.; Chong, T.; Drunkenmiller, H.; Huang, L.Y.; Hultgren, A.; Krasovich, E.; et al. The effect of large-scale anti-contagion policies on the COVID-19 pandemic. Nature 2020, 371, 262–267. [Google Scholar] [CrossRef] [PubMed]
  14. Giordano, G.; Blanchini, F.; Bruno, R.; Colaneri, P.; Di Filippo, A.; Di Matteo, A. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nat. Med. 2020, 26, 855–860. [Google Scholar] [CrossRef] [PubMed]
  15. Holmdahl, I.; Buckee, C. Wrong but Useful—What Covid-19 Epidemiologic Models Can and Cannot Tell Us. N. Engl. J. Med. 2020, 383, 303–305. [Google Scholar] [CrossRef] [PubMed]
  16. Dahal, S.; Luo, R.; Subedi, R.K.; Dhimal, M.; Chowell, G. Transmission Dynamics and Short-Term Forecasts of COVID-19: Nepal 2020/2021. Epidemiologia 2021, 2, 639–659. [Google Scholar] [CrossRef]
  17. Pajuelo-Reyes, C.; Valencia, H.J.; Montenegro, C.C.; Quezada, E.; Gonzales, L.; Cruz, N.; Canelo, C.; Ordinola, C.; Maicelo Quintana, J.L.; Tejedo, J.R.; et al. Epidemiological Analysis of COVID-19 Cases in Native Amazonian Communities from Peru. Epidemiologia 2021, 2, 490–501. [Google Scholar] [CrossRef]
  18. Olumoyin, K.D.; Khaliq, A.Q.M.; Furati, K.M. Data-Driven Deep-Learning Algorithm for Asymptomatic COVID-19 Model with Varying Mitigation Measures and Transmission Rate. Epidemiologia 2021, 2, 471–489. [Google Scholar] [CrossRef]
  19. Fonseca, L.M.; Sousa, D.S.D.; Cardoso, J.C.; Severino, P.; Cano, A.; Souto, E.B.; Lima, S.O.; Oliveira, C.C.C.D.; Reis, F.P. Epidemiology of COVID-19 in the State of Sergipe/Brazil and Its Relationship with Social Indicators. Epidemiologia 2021, 2, 262–270. [Google Scholar] [CrossRef]
  20. Marinho, P.R.D.; Cordeiro, G.M.; Coelho, H.F.C.; Cabral, P.C. The COVID-19 Pandemic in Brazil: Some Aspects and Tools. Epidemiologia 2021, 2, 243–255. [Google Scholar] [CrossRef]
  21. Ogwara, C.A.; Mallhi, A.K.; Hua, X.; Muniz-Rodriguez, K.; Schwind, J.S.; Zhou, X.; Jones, J.A.; Chopak-Foss, J.; Chowell, G.; Fung, I.C.-H. Spatially Refined Time-Varying Reproduction Numbers of COVID-19 by Health District in Georgia, USA, March–December 2020. Epidemiologia 2021, 2, 179–197. [Google Scholar] [CrossRef]
  22. Fain, B.; Dobrovolny, H.M. Initial Inoculum and the Severity of COVID-19: A Mathematical Modeling Study of the Dose-Response of SARS-CoV-2 Infections. Epidemiologia 2020, 1, 5–15. [Google Scholar] [CrossRef] [PubMed]
  23. Lee, W.; Hwang, S.S.; Song, I.H.M. COVID-19 in South Korea: Epidemiological and spatiotemporal patterns of the spread and the role of aggressive diagnostic tests in the early phase. Int. J. Epidemiol. 2020, 49, 1106–1116. [Google Scholar] [CrossRef]
  24. Lewnard, J.A.; Lo, N.C. Scientific and ethical basis for social-distancing interventions against COVID-19. Lancet Infect. Dis. 2020, 20, 631–633. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Matrajt, L.; Leung, T. Evaluating the Effectiveness of Social Distancing Interventions to Delay or Flatten the Epidemic Curve of Coronavirus Disease. Emerg. Infect. Dis. 2020, 26, 1740–1748. [Google Scholar] [CrossRef]
  26. Chou, R.; Dana, T.; Jungbauer, R.; Weeks, C.; McDonagh, M.S. Masks for Prevention of Respiratory Virus Infections, Including SARS-CoV-2, in Health Care and Community Settings. Ann. Intern. Med. 2020, 173, 542–555. [Google Scholar] [CrossRef] [PubMed]
  27. Kretzschmar, M.E.; Rozhnova, G.; Bootsma, M.C.J.; van Boven, M.; van de Wijgert, J.H.H.M.; Bonten, M.J.M. Impact of delays on effectiveness of contact tracing strategies for COVID-19: A modelling study. Lancet Public Health 2020, 5, e452–e459. [Google Scholar] [CrossRef]
  28. Abernethy, R.B. The New Weibull Handbook: Reliability and Statistical Analysis for Predicting Life, Safety, Supportability, Risk, Cost and Warranty Claims, 5th ed.; Dr. Robert B. Abernethy: North Miami, FL, USA, 2004. [Google Scholar]
  29. Ioannidis, J.P.A.; Cripps, S.; Tanner, M.A. Forecasting for COVID-19 has failed. Int. J. Forecast. 2022, 38, 423–438. [Google Scholar] [CrossRef]
  30. Ribeiro, M.H.D.M.; da Silva, R.G.; Mariani, V.C.; Coelho, L.S. Short-term forecasting COVID-19 cumulative confirmed cases: Perspectives for Brazil. Chaos Solitons Fractals 2020, 135, 109853. [Google Scholar] [CrossRef]
  31. Yadav, S.K.; Akhter, Y. Statistical Modeling for the Prediction of Infectious Disease Dissemination with Special Reference to COVID-19 Spread. Front. Public Health. 2021, 9, 645405. [Google Scholar] [CrossRef]
  32. Clifford Cohen, A. Maximum Likelihood Estimation in the Weibull Distribution Based on Complete and on Censored Samples. Technometrics 1965, 7, 579–588. [Google Scholar] [CrossRef]
  33. Engelhardt, M.; Bain, L.J.; Blumenthal, S. Statistical Analysis of a Weibull Process with Left-Censored Data. In Survival Analysis: State of the Art, 1st ed.; Klein, J.P., Goel, P.K., Eds.; Nato Science (Series E: Applied Sciences); Springer: Dordrecht, The Netherlands, 1992; Volume 211. [Google Scholar]
  34. Cole, R. Data Errors and Forecasting Accuracy. In Economic Forecasts and Expectations: Analysis of Forecasting Behavior and Performance, 1st ed.; National Bureau of Economic Research (NBER): Cambridge, MA, USA, 1969; pp. 42–87. [Google Scholar]
  35. Collett, D. Parametric proportional hazards models. In Modelling Survival Data in Medical Research, 1st ed.; Chapman & Hall/CRC: London, UK, 2015; pp. 171–220. [Google Scholar]
  36. Linton, N.M.; Kobayashi, T.; Yang, Y.; Hayashi, K.; Akhmetzhanov, A.R.; Jung, S.-M.; Yuan, B.; Kinoshita, R.; Nishiura, H. Incubation Period and Other Epidemiological Characteristics of 2019 Novel Coronavirus Infections with Right Truncation: A Statistical Analysis of Publicly Available Case Data. J. Clin. Med. 2020, 9, 538. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lai, C.K.C.; Ng, R.W.Y.; Wong, M.C.S. Epidemiological characteristics of the first 100 cases of coronavirus disease 2019 (COVID-19) in Hong Kong Special Administrative Region, China, a city with a stringent containment policy. Int. J. Epidemiol. 2020, 49, 1096–1105. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Trendline analysis of the cumulative infected cases.
Figure 1. Trendline analysis of the cumulative infected cases.
Epidemiologia 04 00003 g001
Figure 2. Trendline analysis of the daily infected cases.
Figure 2. Trendline analysis of the daily infected cases.
Epidemiologia 04 00003 g002aEpidemiologia 04 00003 g002b
Figure 3. Comparison between the derived trendline fit (Equation (10)) and the 5-day average, together with the intervention enforcement/relaxation dates.
Figure 3. Comparison between the derived trendline fit (Equation (10)) and the 5-day average, together with the intervention enforcement/relaxation dates.
Epidemiologia 04 00003 g003
Figure 4. Trendline analysis of the cumulative swab-tests.
Figure 4. Trendline analysis of the cumulative swab-tests.
Epidemiologia 04 00003 g004
Figure 5. Trendline analysis of the daily swab-tests.
Figure 5. Trendline analysis of the daily swab-tests.
Epidemiologia 04 00003 g005
Figure 6. Positivity rate throughout the first wave.
Figure 6. Positivity rate throughout the first wave.
Epidemiologia 04 00003 g006
Figure 7. One-day predictive output of infected cases within the first & second COVID-19 waves.
Figure 7. One-day predictive output of infected cases within the first & second COVID-19 waves.
Epidemiologia 04 00003 g007
Table 1. Intervention enforcement/relaxation measures & dates.
Table 1. Intervention enforcement/relaxation measures & dates.
(A)
Lockdown Measures
Enforcement DateRelaxation Date
International travel ban21 March 20201 July 2020
Lockdown of vulnerable persons28 March 20205 June 2020
National travel ban
(essential travel only
between the Maltese islands)
3 April 20204 May 2020
(B)
Social Distancing Measures
Enforcement DateRelaxation Date
Public transport measures (daily decontamination;
passenger screening;
standing passengers disallowed; windows open;
air-conditioning system off; no monetary change)
12 March 20204 May 2020
Closure of workplaces and
distancing of workers
13 March 20205 June 2020
Closure of sports facilities13 March 20205 June 2020
Closure of law courts and local tribunals13 March 20205 June 2020
Closure of religious places13 March 202013 June 2020
Closure of service outlets and public places16 March 202022 May 2020
Closure of education
establishments
21 March 20205 June 2020
Closure of non-essential
retail outlets
23 March 20204 May 2020
Closure of non-essential
service outlets
23 March 202022 May 2020
Measures to protect elderly
and high-risk groups
28 March 20205 June 2020
Prohibition of public gatherings
(limits of 3 persons, 4 persons, and 6 persons)
30 March 20204 May 2020, 22 May 2020,
5 June 2020
Suspension of visits to homes
for the elderly and the national hospital
8 April 202025 May 2020
Table 2. Cumulative infected cases trendline data.
Table 2. Cumulative infected cases trendline data.
Equation M c 1 λ c 1 k c 1 M c 2 λ c 2 k c 2 M c 3 λ c 3 k c 3 R 2
Equation (1)6.486--------0.806
95% CI6.677--------
6.296
Equation (2)885.678.32-------0.970
95% CI940.387.12-------
830.969.53
Equation (3)680.049.611.581------0.991
95% CI691.251.011.655------
668.948.211.508
Equation (4)434.831.382.628234.379.054.687---0.996
95% CI472.033.122.836273.882.325.888---
397.629.642.419194.775.783.486
Equation (5)484.633.642.48664.03100.212.25123.473.5917.980.998
95% CI495.634.412.60682.83103.720.26145.174.7323.93
473.632.862.36745.2396.694.327101.772.4512.03
Table 3. Daily infected cases trendline data.
Table 3. Daily infected cases trendline data.
M c 1 λ c 1 k c 1 M c 2 λ c 2 k c 2 M c 3 λ c 3 k c 3 R 2
Equation (6)4.542--------0.0
95% CI5.427--------
3.657
Equation (7)818.493.88-------0.183
95% CI1077137.1-------
560.050.71
Equation (8)618.850.321.568------0.376
95% CI716.357.791.809------
521.342.841.327
Equation (9)371.130.772.357222.476.425.118---0.459
95% CI476.936.072.836273.882.325.888---
265.329.642.419194.775.783.486
Equation (10)429.633.942.16851.6199.8618.60124.874.3517.370.556
95% CI495.634.412.60682.83103.720.26145.174.7323.93
473.632.862.36745.2396.694.327101.772.4512.03
Table 4. Cumulative swab-test trendline data.
Table 4. Cumulative swab-test trendline data.
M s 1 λ s 1 k s 1 R 2
Equation (11)673.0--0.901
95% CI650.1--
695.9
Equation (12)118,770.3104.92.6520.999
95% CI117,054.6103.72.612
120,486.0106.12.692
Table 5. Daily swab-tests trendline data.
Table 5. Daily swab-tests trendline data.
M s 1 λ s 1 k s 1 R 2
Equation (13)749.8--0.0
95% CI678.8--
820.9
Equation (14)127,657.8110.02.4990.762
95% CI118,635.2105.22.317
136,680.4114.82.681
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Borg, M.G.; Borg, M.A. A Trendline and Predictive Analysis of the First-Wave COVID-19 Infections in Malta. Epidemiologia 2023, 4, 33-50. https://doi.org/10.3390/epidemiologia4010003

AMA Style

Borg MG, Borg MA. A Trendline and Predictive Analysis of the First-Wave COVID-19 Infections in Malta. Epidemiologia. 2023; 4(1):33-50. https://doi.org/10.3390/epidemiologia4010003

Chicago/Turabian Style

Borg, Mitchell G., and Michael A. Borg. 2023. "A Trendline and Predictive Analysis of the First-Wave COVID-19 Infections in Malta" Epidemiologia 4, no. 1: 33-50. https://doi.org/10.3390/epidemiologia4010003

Article Metrics

Back to TopTop