- Open Access
Predictors of compliance with community-directed treatment with ivermectin for onchocerciasis control in Kabo area, southwestern Ethiopia
© Endale et al.; licensee BioMed Central. 2015
- Received: 20 January 2015
- Accepted: 26 January 2015
- Published: 15 February 2015
Compliance with annual ivermectin treatment is a major challenge in community-directed treatment with ivermectin (CDTI) implementation. There are individuals who do not comply with the annual mass treatment, which contributes to the continuity for disease transmission. Hence, ensuring high treatment coverage and sustained compliance should be given due emphasis in the control of onchocerciasis. The aim of this study was to determine CDTI compliance rate and predictors of compliance where the CDTI was in its 9th round in Kabo area, southwestern Ethiopia.
Community-based cross-sectional study was conducted in Kabo area, three weeks after the 9th round of annual ivermectin distribution. Systematic random sampling was used to select head of households and structured, pre-tested questionnaire was used to interview the study participants. Data was analyzed using SPSS version 16. Descriptive statistics was used to compute mean and standard deviation of continuous variables and frequency for categorical variables, while bivariate and multivariate logistic regressions were used to assess the effects of independent variables on the outcome variable. Variables which showed association in multivariate analysis were considered as final predictors of compliance and strength of association was measured through adjusted odds ratio (AOR).
A total of 308 respondents (age range 18-70, mean age ± SD, 32.21 ± 9.64) participated in the study. Of these, 249 (80.8%) reported that they took ivermectin during the 9th round annual treatment. Significantly higher rate of treatment compliance was reported by participants age ≥35 years (AOR = 5.48, 95% CI; 1.97 - 15.23), participants who stayed in the area for more than ten years (AOR = 3.86, 95% CI; 1.83- 8.11), participants who perceive that they are at risk of contracting the disease(AOR = 7.05, 2.70- 18.43), participants who perceive community drug distributors (CDDs) are doing their work well (AOR = 2.35 95% CI; 1.15- 4.83) and participants who know at least one CDD in their village (AOR = 2.83, 95% CI; 1.26- 6.40).
The majority of the study participants in the present study area complied with ivermectin treatment. Nevertheless, intervention packages should consider factors such as age, residence duration and community’s perception of the disease to improve compliance and make drug distribution sustainable.
Onchocerciasis is a debilitating vector-borne disease caused by the parasite, Onchocerca volvulus; transmitted by the bite of black-fly, Simulium damnosum . The disease is public health and socio-economic threat in many African countries [2,3]. Worldwide there are more than 120 million people at risk of contracting the disease, while 18 million people are already infected. More than 99% of the disease burden is in Africa . Onchocerciasis affects the working age population, and it is the second-leading infectious cause of blindness worldwide, being responsible for about 500,000 blindness [5,6].
In Ethiopia, it is estimated that 3 million people are already infected, whereas 7.3 million are at risk of infection. Nine regions surveyed for onchocerciasis were found to be endemic and the endemic areas extend from the northwest part to southwest part of the country [3,7,8]. Onchocerciasis has been targeted for control and eventually elimination, as a disease of public health and socio-economic problems in Ethiopia along with 19 other African countries .
Annual mass treatment with an oral microfilaricidal (ivermectin) is the main control strategy for onchocerciasis in Ethiopia . A community-directed treatment with ivermectin (CDTI) strategy was adopted in line with the African Program for Onchocerciasis Control (APOC) since 1996. CDTI is a cost-effective strategy . Several studies have shown that mass treatment with ivermectin is safe and effective in reducing the transmission of the disease [2,10-14]. The annual mass treatment must continue for about 15 years [15,16], while all eligible members of the community should take the drug in order to ensure an adequate reduction in the transmission of the parasite and protection from further infection . Studies show that not all eligible members of the community receive the annual treatment . Hence, compliance with annual ivermectin treatment has become a major challenge to put CDTI implementation on the ground.
Studies done elsewhere have shown that compliance with ivermectin treatment can be influenced by factors such as age and gender [16,19-21], beliefs about individual susceptibility to onchocerciasis [16,22], programme organization, perceived social influence and support, perceived benefits of ivermectin treatment , characteristics of the client, perceived performance of CDDs/providers [16,22] and methods of dose determination . CDTI is currently on-going in Ethiopia. However, there is limited information regarding treatment compliance rate and factors affecting treatment compliance. Therefore, the aim of this study was to determine the rate of compliance with CDTI and identify factors influencing compliance with CDTI in Kabo area, which would contribute to the success of the on-going program.
Study area and population
Community-based cross-sectional study was conducted in Kabo area, Gambella Regional State, southwestern Ethiopia, between November and December 2012. The study area is located about 667 kilometers south west of Addis Ababa and 340 kilometers from Gambella Town. The major ethnic composition of the area include Amhara, Oromia, Tigre, Kafficho, Kambata, Shakacho, Mezhenger and some others from southern Ethiopia. Most of the inhabitants are engaged in mixed farming like coffee plantation, maize cultivation and rearing livestock. Based on the 2007 Ethiopian National Population and Housing Census, the population of the area is projected to be about 46,583, with 23,429 males and 23,154 females .
Community-based treatment with ivermectin (Mectizan®) for the control of onchocerciasis was initiated in the area in 2004 by World Health Organization (WHO)/APOC in partnership with Ethiopian Federal Ministry of Health (FMOH), The Carter Centre, the local administration and the communities . Since then, ivermectin has been annually distributed to all eligible members of the community in the villages through CDDs , using the height and physical appearance dosing method . According to the information obtained from Gambella Regional Health Bureau (GRHB), in 2011 the annual treatment coverage of ivermectin in the area was around 76% .
Sample size estimation and data collection
To determine the number of participants to be included in the study, single population proportion formula was used with the assumption that: the proportion of ivermectin treated individuals (treatment compliers) in the community was 75.8% , confidence level 95% and degree of precision 5%. To compensate for the non-respondents, 10% was considered and thus, the final sample size was estimated to be 308.
The study area was stratified into four manageable villages based on their proximity. Then, the sample size was distributed proportionally to each village based on their population size. Systematic random sampling technique was used to recruit the households for this study, and heads of households (i.e. husband or wife or any representative), aged 18 years and older, and who stayed in the area for more than one year, apparently healthy, and volunteered to participate were included in this study. Pregnant and lactating women having infants younger than one month of age (during the data collection period) and individuals who were seriously sick were excluded from the study since they had already been excluded from the mass treatment.
A structured, pre-tested questionnaire was used for data collection. The questionnaire was adopted from similar studies conducted previously [16,22]. It has three major components: socio-demographic factors, behavioral factors (such as knowledge, belief and attitude/perception towards the treatment, the disease, and CDDs) and service related factors. The service related factors here refer to activities related to service provision such as timing of service provision (week end, working days, holidays), decision maker of ivermectin distribution (Woreda/District Office, CDDs, Community leaders, HEWs), places of service provision (home to home, fixed village centers and others) and sources of information about the service (family, local leaders, health workers, health extension workers, CDDs, friends, religious leaders, radio and others). The questionnaire was prepared in English and then, translated into Amharic and back to English to check for its consistency. The survey was conducted three weeks after the 9th round of annual ivermectin distribution in the area. The participants were interviewed in local languages by trained data collectors who were selected from the area and the interview made by house-to-house visit.
Data entry and analysis
After the completion of data collection, cleaning, editing and coding was done. Then the data were entered using EpiData software, V.3.1 and analyzed using SPSS version 16.0. Descriptive statistics was used to compute mean and standard deviation of continuous variables and frequency for categorical variables. Bivariate and multivariate logistic regressions were used to assess the effects of independent variables on the outcome variable, while simultaneously controlling for other potential confounding factors. To identify independent predictors of compliance with ivermectin treatment, backward multivariate logistic regression method was used. Variables which showed association in multivariate analysis at 5% level of significance were considered as final predictors of compliance with ivermectin treatment. The strength of association between different exposure variables and the outcome variable was measured through adjusted odds ratio.
The study was carried out after obtaining ethical approval from the Institutional Review Board (IRB) of Aklilu Lemma Institute of Pathobiology, Addis Ababa University. Then, permission was obtained from GRHB, Zonal Health Departments & District Health Offices. Before the interview commenced, informed verbal consent was obtained from all study participants after explaining the objective of the study to the participants.
Socio-demographic characteristics of respondents
Socio-demographic characteristics of 308 study participants recruited from the community, Kabo area, southwestern Ethiopia, 2012
Length of stay
CDTI compliance rate
CDTI compliance rate of the study participants in Kabo area, southwestern Ethiopia, 2012
Variables (N = 308)
Did you receive ivermectin in the recent campaign (this year)?
If no what was the reason? (n = 59)
CDDs not coming to their house
Not being informed
Believing free things are useless for health
Others (fear of side effects, shortage of drug etc)
From where do you collect the tablet? (n = 249)
CDDs come to my house
From village center
From CDDs house
1 (0.4 )
CDTI compliance rate versus socio-demographic factors
CDTI compliance rate versus socio-demographic factors among the study participants in Kabo area, southwestern Ethiopia, 2012
Variables (N = 308)
Compliance to Ivermectin
≤35 Years (213)
>35 Years (95)
Length of stay
<10 Years (74)
≥10 Years (234)
Compliance rate versus behavioral and service related factors
Compliance rate versus behavioral factors among the study participants in Kabo area, southwestern Ethiopia, 2012
Variables (N = 308)
Compliance with Ivermectin
Know (heard about) the disease onchocerciasis
Know the vector (Black fly) for transmission
Know mode of transmission (Black fly bite)
Know sing and symptoms (at least one)
The responsible body for distributing the drug
Health professionals (20)
Know any CDD in their village
Know how CDDs recruited
Attend village meeting on CDTI at least once
Onchocerciasis a serious disease
Onchocerciasis common in their village
Risk perception to onchocerciasis
Yes at risk (252)
Not at risk (34)
Perceived performance of CDDs
Perceived importance of CDTI
As very important (281)
As obligation (27)
Perception towards the program (can control onchocerciasis)
Know a person who stopped treatment
Perceived reason of stopping the treatment
Sever effect of the drug
Lack of good case mgt.
The magnitude of compliance to treatment was significantly higher for respondents who perceive onchocerciasis as a serious disease (91.2%) and common in their village (92.6%) compared to those who perceive onchocerciasis as not serious (30.3%) and not common in their village (46.8%) (p < 0.05). Similarly, those who perceive that they are at risk of contracting the disease had significantly higher treatment compliance rate (97.4%) compared to those who perceive themselves not at risk of getting onchocerciasis (27.6%) (p < 0.05). In addition, respondents who perceived CDDs are performing their work well had significantly higher compliance rate (89.1%) compared to those who deem CDDs are poor in performance (12.1%) (p < 0.05). Furthermore, respondents who perceived that CDTI as very important (87.9%) and the program can control onchocerciasis (89.1%) had significantly higher compliance rate compared to those who look at CDTI as an obligation (7.4%) and those who believe the program can’t control onchocerciasis (33.3%) (p < 0.05) (Table 4).
None of the service-related factors (timing of service provision, decision maker of service periods, places of service provision and sources of information about the service) was associated with compliance to ivermectin treatment.
Independent predictors of compliance with CDTI
Crude analysis of socio-demographic variables (Table 3) showed that age, educational status, ethnicity and length of stay in the study area were significantly associated with compliance to ivermectin treatment (p < 0.05). On the other hand; sex, marital status, religion, occupation, family size and monthly income did not show statistically significant association. While among the behavioral factors (Table 4); awareness of onchocerciasis, perception towards onchocerciasis, perceived performance of CDDs, perception towards the CDTI and the program, knowing a person stopped treatment in their village, awareness regarding a person (professional) responsible for distributing ivermectin and knowing at least one CDD in their village had shown significant association with compliance to ivermectin treatment. Knowledge on the responsible vector for the transmission, mode of transmission and sign and symptoms, previous family history of onchocerciasis, awareness on how CDDs recruited, awareness of other useful effects of ivermectin, methods of ivermectin dose determination, existence of problem with drug distribution, and participation in community meeting regarding CDTI from the behavioral factors did not show statistically significant association with compliance to CDTI in the bivariate analysis.
Independent predictors of CDTI in Kabo area, southwestern Ethiopia, 2012
Compliance with Iv.
Crude OR (95% CI)
Adjusted OR (95% CI)
4.91 (2.03, 11.88)
5.48 (1.97, 15.23)
Length of stay
4.38 (2.40, 8.02 )
3.86 (1.83, 8.11)
Yes at risk
8.33 (3.84, 18.05)
7.05 (2.70, 18.43)
Not at risk
Perceived performance of CDDs
4.18 (2.31, 7.56)
2.35 (1.15, 4.83)
Familiarity with CDDs (at least one)
3.64 (1.92, 6.88)
2.83 (1.26, 6.40)
Respondents with age greater than 35 years old were 5.48 times more likely to comply with CDTI compared to those whose age was less than or equal to 35 years (adjusted odds ratio = 5.48; 95% CI: 1.97 - 15.23). Similarly, individuals who stayed in the study area for more than or equal to ten years were 3.86 times more likely to comply with the treatment compared to those who stayed in the area for less than ten years (adjusted odds ratio = 3.86; 95% CI: 1.83 - 8.11). The respondents who perceived themselves to be at risk of onchocerciasis infection were 7.05 times more likely to comply with the treatment compared to those who did not (adjusted odds ratio = 7.05; 95% CI: 2.70 - 18.43). In addition, respondents who perceived that CDDs are doing their work well were 2.35 times more likely to comply with CDTI compared to their counter parts (adjusted odds ratio = 2.35, 95% CI: 1.15, 4.83). Moreover, respondents who know at least one CDD in their village were 2.83 times more likely to comply with CDTI compared to those who did not know any CDD in their village (adjusted odds ratio = 2.83; 95% CI: 1.26 - 6.40) (Table 5).
The results of this study indicated that 80.8% of the respondents took ivermectin during the 9th round of the annual ivermectin distribution. Although the treatment coverage seems good, there is a need to attain and maintain a high coverage of CDTI; at least 90%, for elimination of onchocerciasis as a public health problem [17,25]. In this study, almost one fifth (19.2%) of eligible respondents in the community were non-compliers to CDTI. This implies that members of the community who do not take the treatment may serve as a reservoir for continued transmission of onchocerciasis in the area. The major reasons for not taking the annual treatment include: being in the field (farming) during the campaign day, homes not being visited by CDDs, not being informed about the campaign day, do not trust the treatment, and other reasoning’s such as fear of side effects. A study conducted in Uganda showed that among those individuals eligible for treatment, 17.4% were reported as non-compliers in the tenth round of ivermectin distribution , although their study did not indicate the reasons. Another study conducted in Nigeria reported that absenteeism during the campaign day was the major reason for missing the treatment .
The need to understand factors contributing to compliance to CDTI is paramount for improving the efforts in the prevention and control of onchocerciasis as an economic and public health problem. Findings of this study showed that age over 35 years, length of stay, perception about risk of contracting onchocerciasis, perception of the performance of CDDs and familiarity with CDDs in their villages were the factors found to be independent predictors of compliance with ivermectin treatment. Older age groups (age >35 years) had a positive association with CDTI compliance and were almost five and half times more likely to comply with the treatment compared to younger ones. The finding is consistent with a study conducted in Nigeria and Cameroon [19,27]. The reason may be attributed to the fact that younger individuals are relatively highly mobile and likely to travel outside the village for work and other opportunities or they could be reluctant and miss the annual treatment. Moreover, a shorter period of ivermectin distribution that lasts for just a few days possibly increases the number of non-compliers. Respondents who stayed in the study area for ten years and more were more likely to comply with the treatment compared to those who stay in the study area less than ten years. This might be due to the fact that the more a person stays in the endemic area, the more he/she acquires relevant information on the benefits of the treatment from the existing intervention package in the area.
The strongest factor associated with compliance to CDTI in this study was perceived personal risk of onchocerciasis. Individuals who perceive that they are at risk of getting the infection were nearly seven times more likely to comply with CDTI compared to those who consider themselves as free from risk of the infection. Other studies also showed that when people believe that personal susceptibility to onchocerciasis is high then, adherence to the treatment is also high [16,18,22,28].On top of this, the finding also demonstrated a pivotal role of CDDs in the success of CDTI. Perceived good performance of CDDs and knowing at least one CDD was associated with increased compliance. Thus, respondents who perceive CDDs are doing their work well were more likely to comply with the treatment and those who know at least one CDD in their village were likely to comply with CDTI compared to their counter parts. This finding is in congruent with the studies conducted in southwestern Ethiopia, Nigeria and Uganda [16,19,22]. The study conducted in southwestern Ethiopia (Bebeka coffee plantation) revealed insufficient knowledge on the side of CDDs resulted in hampered acceptability of CDDs performance by the community . It is important that CDDs are able to carry out their tasks and need regular training .
One of the limitations of this study was that we used a cross-sectional study design, which is a snapshot of a single point in time and misses the seasonal trends of CDTI compliance rate in the study area. Since the documentation and registration system of CDTI implementation was weak in the area, to avoid recall bias we conducted a cross-sectional study just three weeks after the annual ivermectin distribution. The other limitation was that a relatively higher degree of accuracy (5% degree of precision) was used which limits our sample size to 308.
In conclusion, the annual treatment coverage of ivermectin in the study area seems good, but this does not mean that all eligible individuals have been complying with the treatment. There are members of the community that do not comply with the annual treatment particularly younger individuals, recent comers and reluctant. Hence, a special emphasis needs to be given for these segments of population to increase awareness about perceived personal risk of onchocerciasis and so to address the issue of non-compliance. Moreover CDDs need to be continuously motivated and supported to perform their work well and also they should have to develop a rapport with the community members. On top of this, there is a need for further operational research for achieving and sustaining high level of coverage.
The authors thank all the residents of Kabo area, the study participants, data collectors, Gambella Regional Health Bureau and Godere Woreda Health Office for their cooperation during data collection. The study was financially supported by the School of Graduate Studies, Addis Ababa University, through Aklilu Lemma Institute of Pathobiology.
- WHO. Onchocerciasis Control in the WHO African Region: Current Situation and Way Forward. Report of Regional Committee for Africa; 2007. Available: http://www.who.int/iris/handle/10665/1845. [Accessed: July 26, 2014].
- Dadzie Y, Neira M, Hopkins D. Final report of the Conference on the eradicability of Onchocerciasis. Filaria J. 2003;2:2.View ArticlePubMed CentralPubMedGoogle Scholar
- Rasheed M. Onchocerciasis in different regions of Ethiopia. Int J Par Dis. 2007;1:1–15.Google Scholar
- Ken G, Adrian H, Mauricio S. Onchocerciasis from arrival to near elimination. Parasit Vectors. 2011;4:205–18.View ArticleGoogle Scholar
- WHO. Prevention of blindness and visual impairment. Geneva: World Health Organization; 2011.Google Scholar
- WHO. Conceptual and operational framework of onchocerciasis elimination with ivermectin treatment. JAF. 2010;16:5–23.Google Scholar
- Deribe K, Meribo K, Gebre T, Hailu A, Ali A, Aseffa A, et al. The burden of neglected tropical diseases in Ethiopia, and opportunities for integrated control and elimination. Parasit Vectors. 2012;5:240.View ArticlePubMed CentralPubMedGoogle Scholar
- Weldegebreal F, Medhin G, Weldegebriel Z, Legesse M. Assessment of community’s knowledge, attitude and practice about onchocerciasis and community directed treatment with Ivermectin in Quara District, northwestern Ethiopia. Parasit Vectors. 2014;7:98.View ArticlePubMed CentralPubMedGoogle Scholar
- Bulto T, Workineh W, Taticheff S, Kebede A, Tilahun D. Mass treatment of Onchocerchiasis with Ivermectin at Bebeka Coffee Plantation in Ethiopia: Operational aspects for drug distribution and adverse reaction observed. Ethiop J Health Dev. 1998;4:1–5.Google Scholar
- Vieira JC, Cooper PJ, Lovato R, Mancero T, Rivera J, Proaño R, et al. Impact of long-term treatment of onchocerciasis with ivermectin in Ecuador: potential for elimination of infection. BMC Med. 2007;5:9.View ArticlePubMed CentralPubMedGoogle Scholar
- Abiose A, Jones BR, Cousens SN, Murdoch I, Cassels-Brown A. Reduction in incidence of optic nerve disease with annual ivermectin to control onchocerciasis. Lancet. 1993;341:130–4.View ArticlePubMedGoogle Scholar
- Cousens S, Cassels-Brown A, Murdoch I, Babalola O, Jatau D. Impact of annual dosing with ivermectin on progression of onchocercal visual field loss. Bull World Health Organ. 1997;75:229–36.PubMed CentralPubMedGoogle Scholar
- Brieger WR, Awedoba AK, Eneanya CI, Hagan M, Ogbuagu KF. The effects of ivermectin on onchocercal skin disease and severe itching: results of a multicentre trial. Trop Med Int Health. 1998;3:951–61.View ArticlePubMedGoogle Scholar
- Basanez M, Churcher T. River blindness: a success story under threat. PLoS Med. 2006;3:371.View ArticleGoogle Scholar
- Taylor HR, Pacqué M, Muñoz B, Greene BM. Impact of mass treatment of onchocerciasis with ivermectin on the transmission of infection. Science. 1990;250:116–8.View ArticlePubMedGoogle Scholar
- Yirga D, Deribe K, Woldemichael K, Wondafrash M, Kassahun W. Factors associated with compliance with community directed treatment with ivermectin for onchocerciasis control in Southwestern Ethiopia. Parasit Vectors. 2010;3:48.View ArticlePubMed CentralPubMedGoogle Scholar
- WHO. Report of a WHO Expert Committee on onchocerciasis control, WHO Technical Report, vol. 852. Geneva: World Health Organization; 1995.Google Scholar
- Amazigo U, Okeibunor JC, Matovu V, Zouré H, Bump J, Seketeli A. Performance versus predictors: evaluating sustainability in community-directed treatment projects of the African Programme for Onchocerciasis Control. Soc Sci Med. 2007;64:2070–82.View ArticlePubMedGoogle Scholar
- Brieger WR, Otusanya SA, Oke GA, Oshiname FO, Adeniyi JD. Factors associated with coverage in community directed treatment with ivermectin for onchocerciasis control in Oyo State, Nigeria. Trop Med Int Health. 2002;7:11–8.View ArticlePubMedGoogle Scholar
- Semiyaga NB, Lalobo O, Ndyomugyenyi R. Refusal to take ivermectin: the associated ‘risk’ factors in Hoima district, Uganda. Ann Trop Med Parasitol. 2005;99:165–72.View ArticlePubMedGoogle Scholar
- Lakwo TL, Gasarasi DB. Non-adherence to community directed treatment with ivermectin for onchocerciasis control in Rungwe district, southwest Tanzania. East Afr Med J. 2006;83:326–32.View ArticlePubMedGoogle Scholar
- Nuwaha F, Okware J, Ndyomugyenyi R. Predictors of compliance with community-directed ivermectin treatment in Uganda. Trop Med Int Health. 2005;10:659–67.View ArticlePubMedGoogle Scholar
- Gambella Regional Health Bureau (GRHB). Health management information system report. southwest Ethiopia: Gambella; 2012.Google Scholar
- Boatin BA, Hougard JM, Alley ES, Akpoboua LK, Yaméogo L. The impact of Mectizan on the transmission of onchocerciasis. Ann Trop Med Parasitol. 1998;92:46–60.Google Scholar
- Amazigo UV, Obono M, Dadzie KY, Remme J, Jiya J. Monitoring community-directed treatment programmes for sustainability: lessons from the African Programme for Onchocerciasis Control (APOC). Ann Trop Med Parasitol. 2002;96:75–92.View ArticleGoogle Scholar
- Oyibo WA, Fagbenro-Beyiku AF. Evaluation of community compliance with annual ivermectin treatment of onchocerciasis in Patigi, Nigeria. East Afr Med J. 1998;75:237–42.PubMedGoogle Scholar
- Brieger WR, Okeibunor JC, Abiose AO, Wanji S, Elhassan E, Ndyomugyenyi R, et al. Compliance with eight years of annual ivermectin treatment of onchocerciasis in Cameroon and Nigeria. Parasit Vectors. 2011;4:152.View ArticlePubMed CentralPubMedGoogle Scholar
- Becker MH, Maiman LA, Kirscht JP, Haefner DP, Drachman RH, Taylor DW. Patient perceptions and compliance: recent studies of the health belief model. In: Haynes RB, Taylor DW, Sackett DL, editors. Compliance in Health Care, vol. 8. Baltimore: Johns Hopkins University Press; 1979. p. 78–109.Google Scholar
- Richard F, Boatin B, Saurebrey M, Seketeli A. Control of onchocerciasis today: status and challenges. Trends Parasitol. 2001;17:558–63.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0) which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.