Beyond Exploratory Data: Using Cluster Detection Tests to Pinpoint Malaria Hotspots in Nigeria

S. O. Oyamakin, Owolabi Mayowa
doi: https://doi.org/10.1101/2024.06.04.24308429

Abstract

Malaria is not communicable, but it can deteriorate and even kill a person within a day if treatment is not taken. Plasmodium falciparum is the malaria parasite that is most prevalent and deadly in Africa. This study attempts to identify malaria hotspots and potentially high-risk areas in Oyo state, Nigeria. Data on malaria cases from 2013 to 2016 were provided by the Department of Research, Planning, and Statistics (DPRS) of the Oyo State Ministry of Health. Spatial autocorrelation testing, spatiotemporal analysis, and descriptive statistics were performed using SatScan V9.11 and ArcGIS V10.8. The four-year malaria incidence study was conducted in Oyo West LGA and is geographically focused. Four clusters were found using spatial analysis, with the majority of the clusters located in the Oyo South region. Spatial-temporal research revealed three noteworthy groups. The cluster that happened between January 2013 and December 2014 is located in the Oyo South and Oyo Central region. High-high malaria clusters were identified by spatial autocorrelation analysis, particularly in Oyo Central; the yearly clustered pattern of malaria incidence was not coincidental. In order to better achieve effective disease surveillance, this study suggests employing more cluster detection tests (or CDTs) rather than just exploratory data analysis prior to the implementation of health treatments and regulations.

Introduction

The parasites that cause malaria include Plasmodium falciparum and Plasmodium vivax, which bite humans and create a potentially fatal sickness. The deadliest malaria parasite and most common in Africa is Plasmodium falciparum, which poses the most threat. Outside of sub-Saharan Africa, Plasmodium vivax is the predominant malaria parasite. Within 24 hours, P. falciparum malaria can progress to severe illness and death if left untreated. The early symptoms of malaria, which include fever, headache, and chills, usually appear 10 to 15 days after the infective mosquito bite and can be mild or hard to diagnose. Malaria posed a threat to about half of the world's population in 2020. Infants, children under five, pregnant women, people living with HIV/AIDS, and people with low immunity who are traveling to areas where malaria transmission is intense, such as migrant workers, mobile populations, and travelers, are among the population groups that have a significantly higher risk of contracting malaria and developing severe disease (World Malaria report, 2021).

Figure 1:
Figure 1: Transmission cycle of Malaria

Note. The image in figure 1 was sourced from Tulane University website. https://www2.tulane.edu/~wiser/protozoology/notes/images/mal_lc.gif. The image in figure 2 was sourced from Krekora et al. (2017); the course of pregnancy and delivery in a patient with malaria. Ginekologia Polska. 88. 574-575. 10.5603/GP.a2017.0103.

Figure 2:
Figure 2:

Signs and Symptoms of Malaria

Oyo State, popularly referred to as the “Pace Setter,” is the fifth most populous state in Nigeria and one of the 36 states that comprise the Federal Republic of Nigeria. With a total area of roughly 28,454 square kilometers, it is the fourteenth largest state in Nigeria in terms of size. Oyo State shares borders with Ogun State to the south, Kwara State to the north, Ogun State and the Republic of Benin to the west, and Osun State to the east. Oyo State has an equatorial climate with two distinct seasons: the dry season, which runs from November to March, and the wet season, which begins in April and finishes in October. There is also a fair amount of humidity in the state. The Oyos, Ogbomoshos, Oke-Oguns, Ibadans, and Ibarapas comprise the majority of its local population. They are all members of the Yoruba family and are indigenous to the region in southern Africa. Nigeria's Oyo State is home to 29 Local Council Development Areas and 33 Local Governments. The state's climate is ideal for the growth of crops such cashews, cocoa trees, palm trees, plantains, yam, cassava, rice, and plantains (oyostate.gov.ng).

Figure 3:
Figure 3: Geographical map of Oyo state

Note. A geographical map showing the thirty-three local governments in Oyo state, Nigeria.

Estimates indicate that 76% of Nigerians reside in regions where malaria transmission is high. Nigeria had the greatest percentage of malaria-related deaths (24%), accounting for 27% of all cases reported worldwide. Four African countries—Nigeria (31.9%), the Democratic Republic of the Congo (13.2%), Tanzania (4.1%), and Mozambique (3.8%)—accounted for more than half of all malaria-related deaths globally, according to the World Malaria Report 2020. The U.S. Malaria Initiative for States (PMI-S) was originally adopted by eight states in 2021: Akwa Ibom, Benue, Cross River, Ebonyi, Nasarawa, Oyo, Plateau, and Zamfara. The goal of its implementation was to reduce the mortality rates of mothers and children under five, as well as to improve treatment efficacy and accessibility to malaria services at all levels (USAID, 2021). Oyo State's membership in the previously mentioned project indicates that the state still has a high fatality rate from malaria. As a result, Nigeria needs more illness surveillance in order to provide evidence-based treatments. This comprehensive data would facilitate the distribution of health resources based on risk regions and enhance the availability of malaria treatment and prevention. The global community has paid close attention to the confluence of health crises, including the Novel Coronavirus Disease 2019 (COVID 19), cancer, tuberculosis, HIV/AIDS, and malaria, among others. Numerous tests have been developed to ascertain whether clustering is real and meaningful for such purposes.

Several statistical techniques have been used in disease surveillance to identify spatial clusters. Among these techniques, a subclass known as the generic test searches for clusters without speculating on their potential locations. But whether the statistical significance data for every cluster is accessible will depend on the method employed (Takahashi and Shimadzu, 2021). Developed by Moran (1950), Whitemore et al. (1987), Oden (1995), Tango (1995), Rogerson (1999), and Bonetti and Pagano (2005), global clustering tests (GCTs) are methods that do not ascertain the statistical significance of clusters. On the other hand, some techniques—such as those published by Turnbull et al. (1990), Besag and Newell (1991), Kulldorf and Nagarwalla (1995), Kulldorf (1997), and Tango (2000)—are referred to as cluster detection tests (or CDTs). Thus, there's a need to go over some of the study's supporting literature. Dead bird clusters were evaluated by Farzad et al. (2003) as a potential early warning system for West Nile Virus activity. The prospective detection of infectious disease outbreaks on dead birds reported from New York City in 2000 was accomplished using the spatial scan statistic. Its efficacy in offering a preemptive alert regarding West Nile Virus activity in 2001 was also examined. The study suggests that the development of an early warning system for West Nile virus (WNV) outbreaks could provide the basis for more rapid suppression of adult and larval mosquito populations, as well as more focused public awareness campaigns and monitoring activities. Furthermore, a prospective geographic cluster analysis of bird death reports may alert authorities to increased viral activity in birds and mosquitoes early on, allowing them to better allocate their limited laboratory and mosquito-collection resources and prevent human illness caused by the virus. They believed that additional infectious disease surveillance systems, such as those for bioterrorism, could benefit from the adaptation of the scan statistic.

Before and after adjusting for risk factors, Klassen et al. (2005) looked at the geographic clustering of prostate cancer grade and stage at diagnosis. The study combined projected block group-level disease trends from multilevel models with a spatial scan statistic technique to examine regional variance in prostate cancer grade and stage at diagnosis. According to the study's findings, Maryland prostate cancer cases between 1992 and 1997 fell into four statistically significant clusters with high and low rates of each outcome (later stage at diagnosis and higher histologic grade of tumor). This was accurate even prior to the modifications. After accounting for specific case factors like age, race, and year of diagnosis, the clustering patterns of both outcomes shifted. Before and after adjusting for risk factors, Klassen et al. (2005) looked at the geographic clustering of prostate cancer grade and stage at diagnosis. The study combined projected block group-level disease trends from multilevel models with a spatial scan statistic technique to examine regional variance in prostate cancer grade and stage at diagnosis. According to the study's findings, Maryland prostate cancer cases between 1992 and 1997 fell into four statistically significant clusters with high and low rates of each outcome (later stage at diagnosis and higher histologic grade of tumor). This was accurate even prior to the modifications. After accounting for specific case factors like age, race, and year of diagnosis, the clustering patterns of both outcomes shifted.

Fukuda et al. (2005) examined socioeconomic differences among geographic disease clusters using colon, lung, and breast cancer as case studies from Japan. In order to find spatial clusters of the three common types of cancer (breast, lung, and colon), as well as to identify dimensional socioeconomic factors and minimize the number of variables, the study used the spatial scan statistic and Principal Component Analysis (PCA), which was modified for seven socio-economic indicators. Using Japanese municipal data (N = 3360) from 1993 to 1998, the spatial scan statistic was used to identify spatial clusters of high mortality rates for female breast cancer and male colon and lung cancer. In addition, four sociological clusters (clutter, suburban, rural-poor, and urban-rich) were identified within municipalities based on population density and socioeconomic factors.

METHODOLOGY

Because it works with geographic data, the study of spatial statistics is also known as spatial analysis. Spatial statistics is used in many different fields, and the methods employed depend on the particular research questions being addressed. Examples of these include managing the massive amounts of data generated by GPS (global positioning system) and satellite remote sensing, predicting soil quality in agriculture, and investigating unusually high frequencies of disease events. Spatial data is special because geographic position provides a key that is shared, either exactly or roughly, across data sets from different origins. There are three main types of spatial data: area, lines, and points. Points are features that have a specific location and no further reach in any direction. Cities, industrial zones, and towns are a few instances of point data. Line features are sequences of discrete x, y coordinate pairs with discrete beginning and ending points. Features include things like lines, rivers, and road networks. Shapes known as polygons are created by enclosing an area with a network of connecting lines. Polygons are defined by their area and perimeter. Administrative borders, land use, soil maps, etc. are a few instances of polygon characteristics (Prachi Sahoo, 2012).

DESCRIPTIVE SPATIAL STATISTICS

1 Mean Center

The mean is a critical measure of central tendency for a set of data. When the concept of central tendency is applied to locational point data in two dimensions (X and Y coordinates), the mean center, sometimes referred to as the average location, can be found. It can be found by separately averaging the X and Y

coordinates, as follows: Embedded Image Where; Embedded Image = mean center of X

Embedded Image = mean center of Y

Xi= X coordinate of point i

yi= Y coordinate of point i

n = number of points in the distribution

In many geographical applications, it is appropriate to assign different weights to the points in a spatial distribution. The weights are similar to the frequencies in the analysis of grouped data (weighted mean, etc.). Embedded Image Embedded Image Where; Embedded Image = weighted mean center of X, Embedded Image = weighted mean center of X, and fi= frequency (weight) of point i

2 Euclidean Median

The Euclidean median is an additional useful metric for determining the center. It is more practical to determine the center that minimizes the total unsquared distances rather than the total squared distances. The location that minimizes the total of the Euclidean distances from every other point in a geographical distribution to that central location is called the Euclidean median ( Xe, Ye), or the median center. In this position, the sum is mathematically reduced:also known as the median center, is the place that minimizes the sum of Euclidean distances from all other points in a spatial distribution to that central location. This position reduces the sum mathematically: Embedded Image The simple (unweighted) Euclidean median is logically expanded to become a weighted Euclidean median. Population, sales volume, or any other feature suited to the spatial problem may be represented by the weights or frequencies. The weighted Euclidean median's (Xwe,ywe) coordinates will minimize the expression. Embedded Image

3 Standard Distance

A point pattern's absolute dispersion is measured using the standard distance (SD). The standard distance statistic includes the straight-line or Euclidean distance of each point from the mean center after the geographic coordinates of the mean center have been established. Following are some estimates for standard distance: Embedded Image Weighted standard distance (SWD) is appropriate for those geographic applications requiring a weighted mean center. Embedded Image

4 Spatio-Temporal Analysis - Spatial Scan Statistics

The spatial scan statistic, which is a generalization of the cluster evaluation permutation technique (CEPP) and is a spatial variant of the scan statistic with a variable window size, was proposed by Kulldorff and Nagarwalla (1995) and Kulldorff (1997). Each centroid of a region is given a circular window Z by the spatial scan statistic. The radius of the circle changes from zero to a predetermined upper limit for each of those centroids. The entire region is included in the window if it contains the centroid of an area. In total, a very large number of different but overlapping circular windows are created, each with a different location and size, and each being a potential cluster. Let Zik,K = 1,…, Ki, denote the window composed by the (k − 1) nearest neighbors to region i. Then, all the windows to be scanned by the spatial scan statistic are included in the set Embedded Image

Under the alternative hypothesis, there is an elevated risk within some window Z as compared to outside: Embedded Image Embedded Image

N() and μ() stand for the random number of cases and the expected number of cases, respectively, within the given window. The likelihood to observe the observed number of cases inside and outside of each window can be calculated separately. The most likely cluster (MLC) is identified by the window Z* that maximizes the likelihood ratio. With the help of SaTScan Software, Kulldorff's spatial scan statistic has been used in a wide range of epidemiological investigations as well as disease surveillance for the detection of disease clusters (Kulldorff et al., 2007).

Spatial scan statistics can be applied to both aggregated data and data with specific point locations. The latter will be considered in the proposed study. There are m cells in a research area, C sickness cases overall, and N people in total. The null hypothesis states that we assume each cell's population size-dependent expected values are Poisson distributed and that there are no clusters in the map. We define L(Z) as the likelihood under the null hypothesis that there is no cluster in the zone and L_0 as the probability under the alternative hypothesis that there is a cluster in zone Z, per Duczmal and colleagues (2005).

If μ(Z) is the expected number of cases inside the zone Z under the null hypothesis, N is the number of observed cases during the entire study period; and μ(T) the total number of expected cases in the study area across all time periods (Desjardins et al. 2020). A likelihood ratio test is then used to identify clusters using; Embedded Image Where; λ = Test statistic

L(Z) = Likelihood under the alternative hypothesis

L0 = Likelihood under the null hypothesis

nz = The observed number of malaria cases inside the zone Z

μ(Z) = The expected number of malaria cases inside the zone Z

μ(T) = The total number of expected cases of malaria in the zone Z across all time periods

N = The number of observed malaria cases during the entire study period

Once the value of the test statistic has been calculated, it is easy to do the inference. If a zone exhibits a likelihood ratio greater than 1 i.e. Embedded Image, the risk is elevated. The zone that makes up the most likely cluster is identified by this likelihood ratio test, which maximizes over all of the zones (Kulldorf 1997). The likelihood ratio test, whose P-value is produced by using Monte Carlo simulations, is the foundation for the test of the importance of the detected clusters (Dwass 1957). We might wish to take a closer look at secondary clusters with high likelihood values in addition to the most likely cluster. These are defined as clusters that do not overlap with one that is more likely. Reporting these clusters along with the most likely one is frequently of interest. After accounting for demographic factors like age or race, we might also check to see whether any geographic clusters remain. This could point us in the direction of additional spatially associated risk variables that are currently unknown (Kulldorf, 1997).

5 Spatial Autocorrelation

The results of spatial autocorrelation show whether or not the observed value of a variable at one place is independent of its values at other locations. Correlation statistics were created to show interactions between variables, as opposed to autocorrelation statistics, which are intended to show correlations within variables (Artur, 2007). The link between variables across spatial dimensions is shown by spatial autocorrelation. The spatial statistical methods of global spatial autocorrelation and local spatial autocorrelation are used to characterize the relationship between the research regions and quantify the degree of aggregation or dispersion. The Moran's I index is a method for measuring spatial autocorrelation. While the global Moran's I index is used to evaluate the overall spatial autocorrelation and spatial distribution of the study areas, the local Moran's I index may also be used to reflect the local spatial autocorrelation and the particular clustering areas (Tille et al., 2018).

Moran's I index values vary from -1 to + 1. A positive autocorrelation is denoted by an I > 0, and the distribution of cases is aggregated spatially. A negative autocorrelation is shown by an I < 0, and the more examples are dispersed, the closer the value is to − 1. According to Mao et al. (2019), an I = 0 denotes a random distribution of the cases in space.

Classical Measure of Spatial Autocorrelation

  • Moran's I

  • Geary's C

1 Global Moran I

Moran (1950) proposed a measure to estimate the global spatial autocorrelation (β); Embedded Image Where Xi is the observed value at location i, N is the number of locations and, Embedded Image, (i≠j). The weighting function Wij is used to assign weights to every pair of locations in the study area, Wij = 1, if I and j are neighbours, otherwise Wij = 0.

2 Local Morans I

Local morans I (LISA, Anselin, 1995) is used to explore events that are clustered across geographic units. Local morans I is given below, where Embedded Image and S2 is the mean value and variance respectively. Embedded Image

3 Geary's C

Moran's I provides a more general indication, but the Geary coefficient is more sensitive to changes in a small neighborhood. Geary's C does not offer the same inference since it highlights the values that differ between observational pairs rather than their covariation. Embedded Image

DISCUSSION OF RESULTS

This section attempts to investigate the malaria prevalence in the 33 local governments of Oyo State by discovering new malaria clusters. Descriptive spatial statistics, spatio-temporal analysis, and spatial autocorrelation were performed in order to characterize the geographic features, estimate spatial measures of the data, identify and assess malaria clusters, and ultimately determine whether or not the spatial clustering of malaria incidence across study areas and time periods is due to chance. ArcGIS V10.8 software was used to visualize the results, and SatScan V9.11 was employed for spatiotemporal analysis.

DESCRIPTIVE SPATIAL STATISTICS

The map showing the incidence rate of malaria in the years 2013 and 2014 is shown in excellent detail in Figure 1. In LGAs including Iwajowa, Kajola, Ibarapa East, Akinyele, Ido, Ibarapa Central, Ibarapa East, Olorunsogo, and Egbeda, fewer instances of malaria were reported in 2014 as compared to 2013. Malaria incidences increased in areas like as Iseyin, Saki East, Atiba, Oyo West, Ogbomoso North, Ogbomoso South, Surulere, Afijio, Ogo Oluwa, Lagelu, and Ona Ara. Furthermore, malaria occurrences in LGAs such Saki West, Irepo, Atisbo, Ibarapa North, and Itesiwaju did not significantly increase or decrease during either time period, indicating little improvement in those areas. Overall, more people have contracted malaria in 2014 than in 2013. This might be a result of factors like population growth, an increase in malaria testing, or subpar disease management.

Figure 1:
Figure 1:

Spatial distribution of malaria incidence by LGAs in Oyo state, 2013-2014.

Figure 2 below clearly depicts the map displaying the incidence rate of malaria in the years 2015 and 2016. Less cases of malaria were reported in 2016 compared to 2015 in LGAs such Irepo, Oluyole, Orelope, Kajola, Itesiwaju, Olorunsogo, Atisbo, Iwajowa, Iseyin, Afijio, Ido, Akinyele, and Lagelu. Atiba and Ori Ire saw an increase in malaria cases, whilst Saki West, Saki East, Ogbomoso North, and Surulere LGAs saw no improvement. Due to the fact that there were fewer malaria cases reported in 2016, there was an improvement from 2015 to 2016.

Figure 2:
Figure 2:

Spatial distribution of malaria incidence by LGAs in Oyo state, 2015-2016.

Figure 3 below makes it evident that from 2013 to 2016, Oyo South reported the largest number of malaria cases, while Oyo North recorded the lowest. Four LGAs (Ibarapa East, Ibadan North, Ibadan North East, and Ibadan South West) out of the nine LGAs in Oyo South contributed more than 50% of the total malaria cases reported in that region from 2013 to 2016.

Figure 3:
Figure 3:

Spatial distribution of malaria incidence by senatorial districts in Oyo state, 2013-2016.

ESTIMATE MEASURES OF SPATIAL STATISTICS

Table 2 below shows the malaria incidence mean center, Euclidean median, and standard deviation. The mean center and euclidean median are represented by the latitude and longitude 3.792930 and 3.8413320 and 7.89998310 and 7.8012070, respectively. This implies that the malaria incidence centers in Oyo State are concentrated in the LGAs of Oyo West and Afijio. Additionally, latitude 3.792930 and longitude 7.89998310 are the locations of a standard distance of 0.6003910. The standard distance deviation shows that the spread of malaria incidence from Oyo West LGA also covers areas such as Akinyele, Atiba, Afijio, Saki East, Ori Ire, Olorunsogo, Surelere, Atisbo, Iwajowa, Kajola, Iseyin, Itesiwaju, Ogbomoso North, Ogbomoso South, Ibarapa North, Ibarapa Central, Ibarapa East, Ido, Ibadan North, Ibadan North East, Ibadan South East, Ibadan South West, Ibadan North West, Oyo West, Oyo East, Lagelu, Egbeda, Oluyole, Ona Ara and Ogo Oluwa.

View this table:
Table 1: Distribution of the local governments in Oyo state by senatorial districts
View this table:
Table 2: Spatial measures
View this table:
Table 2: Spatial analysis results of malaria incidence in Oyo state

SPATIAL ANALYSIS OF MALARIA INCIDENCE

H0 : Malaria cluster is not statistically significant over space

H1 :Malaria cluster is statistically significant over space

Table 2 below makes it evident that four significant malaria clusters were identified at the 5% level. Ido, Akinyele, Ibadan South West, Ibadan North –West, Ibadan North, and Ibarapa East LGAs have the most probable cluster (RR = 1.55). Oyo West and Oyo East LGA are home to the first secondary cluster (RR = 1.43). Lagelu LGA is home to the second secondary cluster (RR = 1.33), and Kajola is home to the third secondary cluster (RR = 1.12). As a result, the null hypothesis is rejected at the 5% significance level, and it is determined that there is enough data to draw the conclusion that malaria clusters are significant everywhere.

SPATIO-TEMPORAL ANALYSIS OF MALARIA INCIDENCE

H0 : Malaria cluster is not statistically significant over space and time

H1 :Malaria cluster is statistically significant over space and time

Table 3 below makes it evident that three substantial malaria clusters were identified at the 5% level. Between January 2013 and December 2014, the Ido, Akinyele, Ibadan South West, Ibadan North –West, Ibadan North, and Ibarapa East LGAs had the highest probability of clustering (RR = 2.27). The Oyo south region is mostly home to the most likely cluster. The first secondary cluster was discovered between January 2013 and December 2014 at Atisbo, Saki West, Saki East, kajola, Itesiwaju, and Iwajowa (RR = 1.75). The second secondary cluster was discovered between January 2013 and December 2014 at Ogo Oluwa, Ogbomoso South, Ogbomoso North, Atiba, Oyo East, Oyo West, Surulere, Ori Iri, Afijio, and Lagelu LGA (RR = 1.48). Hence, the null hypothesis is rejected at 5% level of significance and it is concluded there is sufficient evidence that malaria cluster are significant across the geographical locations and time periods.

View this table:
Table 3: Spatio-temporal results of malaria incidence in Oyo state, 2013–2016

SPATIAL AUTOCORRELATION (ANSELIN LOCAL MORAN I)

H0 :Spatial pattern of Malaria incidence is random across the study area

H1 : Spatial pattern of Malaria incidence is not random across the study area

Table 4 below shows that at the 5% level, Akinyele, Egbeda, Ibadan North, Ibadan North East, Ibadan North West, Ibadan South East, Ibadan South West, and Oluyole have a considerable High High concentration of malaria incidence. The low to high clustering of malaria incidence in these areas is indicated by the Low High concentration of malaria incidence observed in LGAs like Lagelu and Ido. Another sign of high to low malaria clustering in this area is the high-low concentration in Atiba. There is a strong likelihood that LGAs with high-high and low-high concentrations of malaria incidence will impact their neighbors with low or no clustering. Incidence of malaria is seen in very low concentrations in Ogbomoso North, Ogo Oluwa, and Surulere. Since P values < 0.05 level of significance, the null hypothesis of randomness is therefore rejected and it is concluded that spatial pattern is clustered, which further depicts strong spatial dependence across study area in year 2013. The geographic visual representation of the local Moran I is displayed in figure 6.

View this table:
Table 4: Local Moran I dependence of malaria incidence among Local governments areas (2013)
Figure 4:
Figure 4:

Map showing the Mean center, Euclidean median and Standard distance of malaria incidence in Oyo state, 2013 – 2016.

Figure 5:
Figure 5:

Spatial Analysis of Malaria incidence in Oyo state, 2013-2016.

Figure 6:
Figure 6:

Spatio-Temporal Analysis of Malaria incidence in Oyo state, 2013-2016.

Figure 6:
Figure 6:

Spatial Autocorrelation of Malaria Incidence in Oyo state, 2013-2014 (Anselin Local Moran I)

Table 5 below shows that Akinyele, Ibadan North, Ibadan North East, Ibadan North West, Ibadan South East, Ibadan South West, and Ona Ara have significant High High concentrations of malaria incidence at the 5% threshold. Low to high clustering of malaria incidence in these locations is indicated by the Low High concentration of malaria incidence observed in LGAs such Lagelu, Ido, Oluyole, and Egbeda. There is a strong likelihood that LGAs with high-high and low-high concentrations of malaria incidence will impact their neighbors with low-low or no clustering. Surulere exhibits a very low concentration of malaria incidence. Since P values < 0.05 level of significance, the null hypothesis of randomness is therefore rejected and it is concluded that spatial pattern is clustered, which further depicts strong spatial dependence across study area in year 2014. The geographic visual representation of the local Moran I is displayed in figure 6.

View this table:
Table 5: Local Moran I dependence of malaria incidence among Local governments areas (2014)

Table 6 below shows that at the 5% level, Ibadan North, Ibadan North East, Ibadan North West, Ibadan South West, Lagelu, and Ona Ara have a statistically significant High High concentration of malaria incidence. In LGAs like Akinyele, Egbeda, Ibadan South East, Ido, and Oluyole, there is a low to high concentration of malaria incidence, which suggests that malaria prevalence is clustered in these regions. There is a strong likelihood that LGAs with high-high and low-high concentrations of malaria incidence will impact their neighbors with low-low or no clustering. The locations of Ibarapa Central, Ibarapa East, and Ibarapa North exhibit a very low concentration of malaria prevalence. Since P values < 0.05 level of significance, the null hypothesis of randomness is therefore rejected and it is concluded that spatial pattern is clustered, which further depicts strong spatial dependence across study area in year 2015. The geographic visual representation of the local Moran I is displayed in figure 7.

View this table:
Table 6: Local Moran I dependence of malaria incidence among Local governments areas (2015)
Figure 7:
Figure 7:

Spatial Autocorrelation of Malaria Incidence in Oyo state, 2015 - 2016 (Anselin Local Moran I)

Table 7 below shows that at the 5% level, Ibadan North, Ibadan North East, Ibadan North West, Ibadan South West, and Ona Ara have a statistically significant High High concentration of malaria incidence. The low to high clustering of malaria incidence in these locations is indicated by the Low High concentration of malaria incidence observed in LGAs such Akinyele, Egbeda, Ibadan South East, Ido, Lagelu, and Oluyole. There is a strong likelihood that LGAs with high-high and low-high concentrations of malaria incidence will impact their neighbors with low-low or no clustering. Ibarapa East, Irepo, Olorunsogo, and Surulere exhibit a very low concentration of malaria incidence. Since P values < 0.05 level of significance, the null hypothesis of randomness is therefore rejected and it is concluded that spatial pattern is clustered, which further depicts strong spatial dependence across study area in year 2016. The geographic visual representation of the local Moran I is displayed in figure 7.

View this table:
Table 7: Local Moran I dependence of malaria incidence among Local governments areas (2016)

SPATIAL AUTOCORRELATION (GLOBAL MORAN I)

Table 8 below displays the 2013 Moran I index calculation result, which is 0.479. The Z score is 4.95 standard deviations, and the p value is 0.001. According to this result, there is a less than 5% chance that the clustered pattern is not a coincidental coincidence. For 2014, the calculated value of the Moran I index is 0.371. The Z score is 4.02 standard deviations, and the p value is 0.002. According to this result, there is a less than 5% chance that the clustered pattern is not a coincidental coincidence. The outcome of the 2015 Moran I index computation was 0.088. The Z score is 1.23 standard deviations, and the p value is 0.114. This result implies that the likelihood that the clustered pattern is a random chance is more than 5%. The outcome of the 2016 Moran I index computation was 0.129. The Z score is 1.71 standard deviations, and the p value is 0.061. This conclusion implies that the likelihood that the clustered pattern is the product of chance is larger than 5%. From 2013 to 2016, the computed value of the Moran I index is 0.273. The Z score is 11.67 standard deviations, and the p value is 0.001. According to this result, there is a less than 5% chance that the clustered pattern is not a coincidental coincidence.

View this table:
Table 8: Moran I Global Spatial Autocorrelation results

CONCLUSION

The results of the study corroborate the hypothesis that spatial analyses identified four major malaria clusters. The most likely cluster is concentrated in the Oyo South region. Spatial-temporal studies also revealed three main malaria clusters. The most likely cluster, which is thought to have formed between January 2013 and December 2014, is located in the Oyo South and Oyo Central region. Six local government units make up the first secondary cluster in Oyo North. In all the geographical and spatio-temporal studies, the same LGAs—Ido, Akinyele, Ibadan South West, Ibadan North -West, Ibadan North, and Ibarapa East—are located in the most likely cluster. The spatial dependence of malaria incidence distribution is illustrated by the findings of the local moran I spatial autocorrelation, which showed a considerable clustering of malaria cases throughout the research area.

Data Availability

All data produced in the present study are available upon reasonable request to the authors

APPENDIX

Figure
Figure
Figure

GLOBAL MORAN I

Note. The image in figure 1 was sourced from Tulane University website. https://www2.tulane.edu/~wiser/protozoology/notes/images/mal_lc.gif. The image in figure 2 was sourced from Krekora et al. (2017); the course of pregnancy and delivery in a patient with malaria. Ginekologia Polska. 88. 574-575. 10.5603/GP.a2017.0103.

References

  1. Alexander, H., Eric, D., & Micheal, D. (2020). Rapid detection of COVID 19 clusters in the United states using a prospective space-time scan statistic: An update. Applied Geography, 118:102–202.
    OpenUrl
  2. Amanda, W., Ramona, L., & Don, W. (2017). Monitoring for Local Transmission of Zika Virus using Emergency Department Data. Journal of Public Health Informatics(1947-2579), 9(1):115. Retrieved from http://ojphi.org
    OpenUrl
  3. Ann, C. K., Martin, K., & Frank, C. (2005). Geographical Clustering of Prostate Cancer Grade and Stage at Diagnosis, before and after Adjustment for Risk Factors. International Journal of Health Geographics, 4:1. doi:10.1186/1476-072X-4-1
    OpenUrlCrossRefPubMed
  4. Arthur, G. (2007). Reflections on spatial autocorrelation: Regional Science and Urban Economics. 37:491–496.
    OpenUrl
  5. ↵
    Besag, J., & Newell, J. (1991). The detection of clusters in rare diseases. J R Stat Soc Series A., 154:143–55.
    OpenUrlCrossRefWeb of Science
  6. ↵
    Bonetti, M., & Pagano, M. (2005). The interpoint distance distribution as a descriptor of point patterns, with an application to spatial disease clustering. Stat Med, 24:753–73.
    OpenUrlCrossRefPubMedWeb of Science
  7. Callado, C. A. (2003). Fecundidad adolescente en el gran ‘area metropolitan de Costa Rica, Población y Salud en Mesoamerica,. 1,4.
    OpenUrl
  8. Carlin, A., Martin, K., & Barry, R. B. (2017). Geospatial analysis of nonmedical vaccine exemptions and pertussis outbreaks in the United States. Retrieved from PNAS: www.pnas.org/cgi/doi/10.1073/pnas.1700240114
  9. Ciara, E. D., Martin, V., & Christopher, A. G. (2019). Resource Allocation for Epidemic Control Across Multiple Sub-populations. Bulletin of Mathematical Biology, 81:1731–1759. Retrieved from doi:10.1007/s11538-019-00584-2
    OpenUrlCrossRef
  10. David, C. W. (2007). A comparison of spatial clustering and cluster detection techniques for childhood leukemia incidence in Ohio, 1996 – 2003. International Journal of Health Geographics., 6:13. doi:10.1186/1476-072X-6-13.
    OpenUrlCrossRefPubMed
  11. Diego, F. C., Susanne, F. A., & Laith, J. A.-R. (2013). Mapping HIV clustering: a strategy for identifying populations at high risk of HIV infection in sub-Saharan Africa. Diego F Cuadros, Susanne F Awad and Laith J Abu-Raddad. Mapping HIV clustering:International Journal of Health Geographics, 12:28. Retrieved from http://www.ij-healthgeographics.com/content/12/1/28
    OpenUrl
  12. Dussault, G., & Maria, C. F. (2006). Not enough there, too many here: understanding geographical imbalances in distribution of the health workforce. Human Resources for Health, 4:12. Retrieved from http://www.human-resources-health.com/content/4/1/12
    OpenUrl
  13. Esra, O., Bryan, L. W., Su, Y. K., & Melina, S. M. (2005). Comparison of spatial scan statistic and spatial filtering in estimating low birth weight clusters. International Journal of Health Geographics, 4:19. doi:10.1186/1476-072X-4-19.
    OpenUrlCrossRefPubMed
  14. ↵
    Farzad, M., Martin, K., Jessica, J. H., James, R. M., & Varuni, K. (2003). Dead Bird Clusters as an Early Warning System for West Nile Virus Activity. Emerging Infectious Diseases, Vol. 9, No. 6.
  15. Geary, R. (1954). The Contiguity Ratio and Statistical Mapping. The Incorporated Statistician., 5:115–145.
    OpenUrlCrossRef
  16. González, R., Augusto, O., Munguambe, K., Pierrat, C., Pedro, E., Sacoor, C., & et al. (2015). HIV Incidence and Spatial Clustering in a Rural Area of Southern Mozambique. PLoS ONE, 10(7): e0132053. doi:10.1371/journal.pone.0132053
    OpenUrlCrossRefPubMed
  17. Grace, K., & Olusanya, O. (2016). Characteristics, Causes and Spillover Effects of Road Traffic Crashes: Spatial Effects in Road Traffic Crashes Modeling. Lambert Academic Publishing. Retrieved from https://www.researchgate.net/publication/304498345
  18. Hao, L., Miao, G., & Mingxin, Z. (2022). Spatio-temporal distribution of tuberculosis and the efects of environmental factors in China. International Journal of Health Geographics, 22:565. Retrieved from doi:10.1186/s12879-022-07539-4
    OpenUrlCrossRef
  19. Haque, U., Sunahara, T., Hashizume, M., Shields, T., Yamamoto, T., & et, a. (2011). Prevalence, Risk Factors and Spatial Distribution in a Hilly Forest Area of Bangladesh. PLoS ONE, 6(4):18908. doi:10.1371/journal.pone.0018908
    OpenUrlCrossRef
  20. Huang, S. S., Yokoe, D. S., Stelling, J., Placzek, H., Kulldorff, M., Kleinman, K., … Platt, R. (2010). Automated Detection of Infectious Disease Outbreaks in Hospitals: A Retrospective Cohort Study. PLoS Medicine, 32. Retrieved from doi:10.1371/journal.pmed.1000238
    OpenUrlCrossRefPubMed
  21. Huaxiang, R., Xinyu, S., & Xi, Z. (2017). Using the Kulldorff’s scan statistical analysis to detect spatio-temporal clusters of Tuberculosis in Qinghai Province, China, 2009-2016. BMC Infectious Diseases, 17:578. doi:10.1186/s12879-017-2643-y
    OpenUrlCrossRef
  22. Iaron, L. S., Andrey, O. P., Taymara, B. R., Glenda, R., Oliveira, N. F., Ana, L. d., … Eliã, P. B. (2022). Temporal Trend and Spatial Analysis of the HIV Epidemic in Young Men who have Sex with Men in the Second Largest Brazilian Amazonian Province. Retrieved from BMC Infectious Diseases: doi:10.1186/s12879-022-07177-w
    OpenUrlCrossRef
  23. Inkyung, J., Martin, K., & Otukei, J. R. (2010). A spatial scan statistic for multinomial data. Inkyung Jung, Martin Kulldorff, and Otukei John Richard. Stat Med., 29(18): 1910–1918. doi:10.1002/sim.3951
    OpenUrlCrossRef
  24. Jean, G., Belco, P., Alassane, D., Stéphane, R., Ousmane, T., Issaka, S., … Ogobara, K. D. (2006). Space-time clustering of childhood malaria at the household level: a dynamic cohort in a Mali village. doi:10.1186/1471-2458-6-286.
    OpenUrlCrossRefPubMed
  25. Jessie, P., James, C., Gibson, C. B., Samuel, M., Timothy, S., Frank, C., & William, J. M. (2015). Spatial clustering of measles cases during endemic (1998–2002) and epidemic (2010) periods in Lusaka, Zambia. Jessie Pinchoff, James Chipeta, Gibson Chitundu Banda, Samuel Miti, Timothy Shields, Frank Curriero, and William John Moss BMC Infectious Diseases, 15:121. doi:10.1186/s12879-015-0842-y
    OpenUrlCrossRef
  26. Jim, E. B., Askari, A., & W., L. B. (2015). Spatial patterns of epilepsy-related emergency department visits in California. Journal of Public Health Research, 4:441.
    OpenUrl
  27. Juliana, G., Thomas, R. N., Nong, S., Sandy, P. A., Maryam, B. H., & Steve, J. K. (2013). Using statistical methods and genotyping to detect tuberculosis outbreaks. International Journal of Health Geographics, 12:15. Retrieved from http://www.ij-healthgeographics.com/content/12/1/15
    OpenUrl
  28. Kirstine, W., Kristine, B., Mogens, L. L., Gunnar, G., & Annette, K. E. (2021). Geographical variation and clustering are found in atrial fibrillation beyond socioeconomic differences: a Danish cohort study, 1987–2015. International Journal of Health Geographics, 20:11. Retrieved from doi:10.1186/s12942-021-00264-2
    OpenUrlCrossRef
  29. Krekora et al. (2017). The course of pregnancy and delivery in a patient with malaria. Ginekologia Polska, 88. 574–575. doi:10.5603/GP.a2017.0103.
    OpenUrlCrossRef
  30. Kui, L., Songhua, C., Zhang, Y., Tao, L., Bo, X., Wei, W., … Jianmin, J. (2022). Tuberculosis Burden caused by Migrant Population in Eastern China: Evidence from notification records in Zhejiang Province during 2013–2017. Retrieved from BMC Infectious Diseases: doi:10.1186/s12879-022-07071-5
    OpenUrlCrossRef
  31. ↵
    Kulldorf, M., & Nagarwalla, N. (1995). Spatial disease clusters: detection and inference. Stat Med., 14:799–810.
    OpenUrlCrossRefPubMedWeb of Science
  32. ↵
    Kulldorff, M. (1997). A spatial scan statistic. Communications in Statistics: Theory and Methods, 26:1481–1496.
    OpenUrlCrossRefWeb of Science
  33. Kunihiko, T., & Hideyasu, S. (2020). Detecting multiple spatial disease clusters: information criterion and scan statistic approach. International Journal of Health Geographics, 19:33. Retrieved from doi:10.1186/s12942-020-00228-y
    OpenUrlCrossRef
  34. Lan, H., Martin, K., & David, G. (2007). A Spatial Scan Statistic for Survival Data. Biometrics 63, 109–118. doi:10.1111/j.1541-0420.2006.00661.x
    OpenUrlCrossRefPubMedWeb of Science
  35. Liqun, F., Lei, Y., Song, L., Sake, J. d., Dan, F., Xiaona, H., … Wuchun, C. (2006). Spatial analysis of Hemorrhagic Fever with Renal Syndrome in China. Retrieved from BMC Infectious Diseases: http://www.biomedcentral.com/1471-2334/6/77
  36. Mahmoud, T. (2012). Spatial Disease Cluster Detection: An Application to Childhood Asthma in Manitoba, Canada. Biomet Biostat. doi:10.4172/2155-6180.S7-010
    OpenUrlCrossRef
  37. Martin, K., Lan, H., & Kevin, K. (2009). A scan statistic for continuous data based on the normal probability model. Martin Kulldorff, Lan Huang, and Kevin Konty. International Journal of Health Geographics, 8:58. doi:10.1186/1476-072X-8-58.
    OpenUrlCrossRefPubMed
  38. Melanie, B., Richard, H., Farzad, M., & Don, W. (2005). Evaluation of school absenteeism data for early outbreak detection, New York City. doi:10.1186/1471-2458-5-105.
    OpenUrlCrossRefPubMed
  39. Michael, A. H., & Andrew, S. C. (2022). Observing the Spread of Common Illnesses Through a Community: Using Geographic Information Systems (GIS) for Surveillance. doi:10.3122/jabfm.2010.01.090137
    OpenUrlAbstract/FREE Full Text
  40. Mirghani, S. E., Bakri, Y. N., Sayed, M. B., Ibrahim, M. E., Robert, W. S., & Abdisalan, M. N. (2010). The spatial-temporal clustering of Plasmodium falciparum infection over eleven years in Gezira State, Sudan. Malaria Journal. doi:10.1186/1475-2875-9-172
    OpenUrlCrossRef
  41. ↵
    Moran, P. (1950). Notes on continuous stochastic phenomena. Biometrika., 37:17–23.
    OpenUrlCrossRefPubMedWeb of Science
  42. ↵
    Oden, N. (1995). Adjusting Moran’s I for population density. Stat Med., 14:17–26.
    OpenUrlCrossRefPubMedWeb of Science
  43. Osei, F., & Duker, A. (2008). Spatial dependency of V. cholera prevalence on open space refuse dumps in Kumasi, Ghana:. Spatial dependency of V. cho A spatial statistical modeling. International Journal of Health Geographics, 7:62.
    OpenUrl
  44. Oyostate. (n.d.). About Oyo State. Retrieved from https://oyostate.gov.ng/about-oyo-state/
  45. Phillip, J., David, G., Stephen, P., Jonathan, S., Keith, L., Peter, H., … Michael, D. (2013). Identifying Probable Suicide Clusters in Wales Using National Mortality Data. PLoS ONE, 8(8):e71713. doi:10.1371/journal.pone.0071713
    OpenUrlCrossRef
  46. Rafael, d. S., Emilly, G. C., Lucas, F., Rodrigues, d. S., Luiz, H. L., Raiany, R. d., … Sarah, C. G. (2022). The Assistance Gaps in Combating COVID-19 in Brazil: for whom, where and when Vaccination occurs. Retrieved from BMC Infectious Diseases: doi:10.1186/s12879-022-07449-5
    OpenUrlCrossRef
  47. Raid, A., Terri, H., Jacob, C., Daniela, S., & Martin, K. (2020). Geographical surveillance of COVID-19:. Diagnosed cases and death in the United States. Retrieved from doi:10.1101/2020.05.22.20110155
    OpenUrlAbstract/FREE Full Text
  48. Rao, H., Xinyu, S., & Zhang, X. (2017). Using the Kulldorff’s scan statistical analysis to detect spatio-temporal clusters of Tuberculosis in Qinghai Province, China, 2009–2016. Retrieved from BMC Infectious Diseases.
  49. ↵
    Rogerson, P. (1999). The detection of clusters using a spatial version of the Chi square goodness-of-ft statistic. Geogr Anal., 31:130–47.
    OpenUrlWeb of Science
  50. Saman, D., Cole, H., Odoi, A., Myers, M., Carey, D., & et al. (2012). A Spatial Cluster Analysis of Tractor Overturns in Kentucky from 1960 to 2002. PLoS ONE, 7(1): 30532. doi:10.1371/journal.pone.0030532
    OpenUrlCrossRef
  51. Samia, E. M., Bakri, Y. N., Sayed, M. B., Ibrahim, M. E., Robert, W. S., & Abdisalan, M. N. (2010). The spatial-temporal clustering of Plasmodium falciparum infection over eleven years in Gezira State. The Sudan. Malaria Journal 2010. doi:10.1186/1475-2875-9-172
    OpenUrlCrossRef
  52. Shanley, C., Michael, N., Roy, B., Liz, H., John, E., & Bin, J. (2013). Geospatial analyses to identify clusters of adverse antenatal factors for targeted interventions. International Journal of Health Geographics, 12:46. Retrieved from http://www.ij-healthgeographics.com/content/12/1/46
    OpenUrl
  53. Stéphanie, G.-B., Claire, D., Aurélie, G.-G., Denis, H., & Jacqueline, C. (2011). Detection of clusters of a rare disease over a large territory: performance of cluster detection methods. International Journal of Health Geographics, 10:53. Retrieved from http://www.ij-healthgeographics.com/content/10/1/53
    OpenUrl
  54. Sujee, L., Jisu, M., & Inkyung, J. (2021). Optimizing the maximum reported cluster size in the spatial scan statistic for survival data. International Journal of Health Geographics, 20:33. Retrieved from doi:10.1186/s12942-021-00286-w
    OpenUrlCrossRef
  55. Taal, A., Blok, D., Handito, A., Wibowo, S., Sumarsono, Wardana A., … Prakoeswa, C. (2022). Determining Target Populations for Leprosy Prophylactic Interventions: a Hotspot Analysis in Indonesia. Retrieved from BMC Infectious Diseases: doi:10.1186/s12879-022-07103-0
    OpenUrlCrossRef
  56. ↵
    Tango, T. (2000). A test for spatial disease clustering adjusted for multiple testing. Stat Med., 19:191–204.
    OpenUrlCrossRefPubMedWeb of Science
  57. ↵
    Tango, T. (1995). A class of tests for detecting ‘general’ and ‘focused’ clustering of rare diseases. Stat Med., 14:23–34.
    OpenUrl
  58. Tobler, W. (1979). Cellular Geography. In philosophy in geography,edited by S.Gale and G.Olsson. Dordrecht: Reidel.
  59. Toshie, M., Dung, P., Yasuhiro, N., Dan, K., Thang, H. N., Thanh, V. D., & Koichiro, K. (2021). Spatiotemporal distribution of COVID-19 during the first 7 months of the epidemic in Vietnam. BMC Infect Dis, 21:1124. Retrieved from doi:10.1186/s12879-021-06822-0
    OpenUrlCrossRef
  61. ↵
    Turnbull, B., Iwano, E., Burnett, W., & et, a. (1990). Monitoring for clusters of disease: application to leukemia incidence in upstate New York. Am J Epidemiol., 132:136–43.
    OpenUrlPubMedWeb of Science
  62. ↵
    USAID. (2021). Retrieved from https://www.usaid.gov/news-information/news/eliminating-malaria-nigeria-fivewaysusaid%E2%80%99s%E2%80%98president%E2%80%99smalariainitiative#:~:text=An%20estimated%2076%20percent%20of,World%20Malaria%20Report%2C%202020).
  63. Vinas, M., Tuduri, E., Galar, A., Yih, K., Pichel, M., & et al. (2013). Laboratory-Based Prospective Surveillance for Community Outbreaks of Shigella spp. in Argentina. PLoS Negl Trop Dis, 7(12):e2521. doi:10.1371/journal.pntd.0002521
    OpenUrlCrossRefPubMed
  64. Wang, J., Cao, Z., Zeng, D., Wang, Q., Wang, X., & et al. (2014). Epidemiological Analysis, Detection, and Comparison of Space-Time Patterns of Beijing HandFoot-Mouth Disease. PLoS ONE, 9(3):e92745. doi:10.1371/journal.pone.0092745
    OpenUrlCrossRef
  65. Wang, Y., Chen, H., Liao, L., Xiaolong, L., & Zhang, A. (2022). Time trend prediction and spatial-temporal analysis of multidrug-resistant tuberculosis in Guizhou Province, China, during 2014–2020. BMC Infectious Diseases, 22:525. Retrieved from doi:10.1186/s12879-022-07499-9
    OpenUrlCrossRef
  66. Wen, L.-M., Jill, A. M., Recinda, S., Lora, E. F., Caio, R.-L., Jennifer, J. H., & David, J. L. (2013). Pancreatic cancer clusters and arsenic contaminated drinking water wells in Florida. BMC Cancer 2013, 13:111. Retrieved from http://www.biomedcentral.com/1471-2407/13/111
    OpenUrl
  67. ↵
    Whitemore, A., Friend, N., Brown, B., & et, a. (1987). A test to detect clusters of disease. Biometrika., 74:631–5.
    OpenUrlCrossRef
  69. Wu, C.-C., Chu, Y.-H., Sanjay, S., & Chen, C.-H. (2021). Spatially varying effects of measured confounding variables on disease risk. International Journal of Health Geographics, 20:45. Retrieved from doi:10.1186/s12942-021-00298-6
    OpenUrlCrossRef
  70. Yao, Z., Mengyuan, Z., Dianju, K., Wei, S., Changhong, Y., & Rongjie, W. (2021). Spatio-temporal analysis of bacillary dysentery in Sichuan province, China, 2011–2019. BMC Infect Dis, 21:1033. Retrieved from doi:10.1186/s12879-021-06738-9
    OpenUrlCrossRef
  71. Yi, Z., Zhixiong, S., Chunna, M., Chengsheng, J., Cindy, F., Nivedita, S., … Quanyi, W. (2015). Cluster of Human Infections with Avian Influenza (H7N9) Cases: A Temporal and Spatial Analysis. Int. J. Environ. Res. Public Health, 12, 816–828. doi:10.3390/ijerph120100816
    OpenUrlCrossRef
  72. Yoshiharu, F., Masahiro, U., Keiko, N., & Takehito, T. (2005). Variations in societal characteristics of spatial disease clusters: examples of colon, lung and breast cancer in Japan. International Journal of Health Geographics, 4:16. doi:10.1186/1476-072X-4-16
    OpenUrlCrossRefPubMed
Back to top
Posted June 04, 2024.
Download PDF
Data/Code
Email
Beyond Exploratory Data: Using Cluster Detection Tests to Pinpoint Malaria Hotspots in Nigeria
S. O. Oyamakin, Owolabi Mayowa
medRxiv 2024.06.04.24308429; doi: https://doi.org/10.1101/2024.06.04.24308429
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area