- Open Access
Bayesian modelling of the effect of climate on malaria in Burundi
© Nkurunziza et al; licensee BioMed Central Ltd. 2010
- Received: 19 November 2009
- Accepted: 29 April 2010
- Published: 29 April 2010
In Burundi, malaria is a major public health issue in terms of both morbidity and mortality with around 2.5 million clinical cases and more than 15,000 deaths each year. It is the single main cause of mortality in pregnant women and children below five years of age. Due to the severe health and economic cost of malaria, there is still a growing need for methods that will help to understand the influencing factors. Several studies have been done on the subject yielding different results as which factors are most responsible for the increase in malaria. The purpose of this study has been to undertake a spatial/longitudinal statistical analysis to identify important climatic variables that influence malaria incidences in Burundi.
This paper investigates the effects of climate on malaria in Burundi. For the period 1996-2007, real monthly data on both malaria epidemiology and climate in the area of Burundi are described and analysed. From this analysis, a mathematical model is derived and proposed to assess which variables significantly influence malaria incidences in Burundi. The proposed modelling is based on both generalized linear models (GLM) and generalized additive mixed models (GAMM). The modelling is fully Bayesian and inference is carried out by Markov Chain Monte Carlo (MCMC) techniques.
The results obtained from the proposed models are discussed and it is found that malaria incidence in a given month in Burundi is strongly positively associated with the minimum temperature of the previous month. In contrast, it is found that rainfall and maximum temperature in a given month have a possible negative effect on malaria incidence of the same month.
This study has exploited available real monthly data on malaria and climate over 12 years in Burundi to derive and propose a regression modelling to assess climatic factors that are associated with monthly malaria incidence. The results obtained from the proposed models suggest a strong positive association between malaria incidence in a given month and the minimum temperature (night temperature) of the previous month. An open question is, therefore, how to cope with high temperatures at night.
- Markov Chain Monte Carlo
- Credible Interval
- Malaria Incidence
- Lower Akaike Information Criterion
In Burundi, malaria is a major public health issue in terms of both morbidity and mortality with around 2.5 million clinical cases and more than 15,000 deaths each year. During the last decade, for example, malaria cases increased from 550,000 cases in 1991 to 2.8 million in 2001 in a total population estimated at 7 million . In 2001, Burundi was the country most affected by malaria in the world . Malaria is the single main cause of mortality in pregnant women and children below five years of age. Malaria continues to ravage millions of rural Burundians, despite concerted efforts to reduce malaria mortality [3, 4]. This is often attributed to a number of factors, including poverty, limited access to basic health care and specialized health facilities, the cost-sharing system, and the under-funding of the health sector by the government. Currently, the government allocates only 2% - 4% of its national budget towards supporting the health sector. The direct economic costs of malaria that result from treatment and from time away from work or school are enormous, but the overall economic impact of malaria is likely to be much more substantial than suggested by estimates of direct costs alone .
Due to the severe health and economic cost of malaria epidemics, there is still a growing need for methods that will help to understand the influencing factors, allow forecasting, early warning, so that more effective control measures may be implemented . Several studies have been done on the subject yielding different results as which factors are most responsible for the increase in malaria. Epidemiology studies on malaria in Africa have identified, in general, an association between climate variables and malaria [7, 8]. An association between climate variability and the epidemics of malaria has been identified in seven sites of East African highlands in Ethiopia, Kenya and Uganda , suggesting that climate variability had an important influence in initiating epidemics in the highlands of East Africa. There was a significant and positive influence of interactions between maximum temperature, minimum temperature and rainfall on malaria transmission. The work of Pemola and Jauhari  investigated the relationship between climate and malaria incidence using Pearson's correlation analysis. The authors found a high positive correlation between monthly parasite incidence and climatic variables (temperature, rainfall and humidity). Gallup and Sachs  have suggested that the location and severity of malaria are mostly determined by climate and ecology. A significant correlation between malaria risk and elevation, annual maximum temperature and rainfall was also described . Another study found that variation in the malaria transmission intensity was strongly associated with basic climatic factors, noting that even small differences in climate variation can significantly affect malaria transmission intensities . Rainfall, temperature and altitude were the most plausible predictors of malaria prevalence in Botswana .
Others have suggested that malaria was influenced by factors other than climate. Hay et al suggested that the association between local malaria resurgence and regional changes in climate was overly simplistic. Economic, social and political factors may explain recent resurgence in malaria and other mosquito-borne diseases, with no need to invoke climatic changes. Many variables may affect malaria transmission beside climatic changes, such as environmental modification, population growth, limited access to health care, and lack of or unsuccessful malaria control measures . Although climate can affect the incidence of malaria, human economic activities and malaria control strategies play an important role in the incidence of the disease . Relatively high rates of malaria morbidity could result from poor access to health services, inadequate case management, overwhelmed health services, poor immunological competence because of malnutrition, a general disruption to livelihoods because of often-associated flooding, or a combination of these factors .
The above-mentioned studies have led to various controversial conclusions when investigating which factors influence malaria incidence, depending on the country specificity. This makes it difficult to establish common determinants of malaria in all countries . In Burundi, the situation of malaria is poorly documented. Indeed, very few studies have been published on the malaria prevalence in Burundi. Further, published papers concern only three provinces (Ngozi, Kayanza, Karuzi) [20, 21], where "Médecins-sans-Frontières" from Belgium was active. It is believed that the situation of crisis experienced by Burundi since 1993, followed by a break in cooperation with some countries led researchers and international agencies to lose interest on Burundi. So far, no research on malaria has been carried out throughout the entire country in order to establish a generalized framework to investigate this subject.
In this paper, both generalized linear and generalized additive mixed models are proposed to assess the climatic factors that are the highly associated with monthly malaria incidence in Burundi. It is assumed that malaria incidence in a given month in Burundi can be predicted by rainfall, temperature and humidity. Rainfall has a great influence on the mosquito population by increasing the vegetation density and by providing suitable breeding pools for the production and maturation of larvae. Temperature has a great influence on the transmission of malaria, higher temperatures shortening the extrinsic incubation period  and humidity facilitates adult mosquito life span [9, 22, 23].
Burundi is located in East-central Africa, between 2°20 and 4°27 of latitude south and between 28°50 and 30°53 of longitude east; the altitude varies between 775 metres (Lake Tanganyika) and 2,670 metres (Crest Congo - Nil). Burundi has generally a tropical highland climate, with a considerable daily temperature variation in many areas . Temperature also varies considerably from one region to another, mainly as a result of differences in altitude. The central plateau is cool, with temperature averaging 20°C. The area near Lake Tanganyika is warmer, averaging 23°C; the highest mountain areas are cooler, averaging 16°C. Rain is irregular, falling most heavily in the north-west . Dry season varies in length, and there are sometimes longer periods of drought. Most parts of Burundi receive between 130 and 160 cm of rainfall a year . Bounded on the north by Rwanda, in south-east by Tanzania and in west by the Democratic Republic of Congo, Burundi covers an area of 27,834 km2 (of which 2,634 km2 are occupied by Tanganyika Lake) and has a population estimated at about 8 million. In terms of habitat, it remains essentially rural, with 91.6% of the population living in rural areas. The urban population is 8.4% with an annual growth rate of 5.7%. The Burundi population is young: 46.1% are under 15 years of age, while people aged 60 and above represent only 5.4%. With an average density of 266 inhabitants per km2, a population growth rate of 3.44% and a total fertility rate of six children per woman, Burundi is one of Africa's most densely populated country . Burundi is structured in 17 provinces. The epidemiological profile can be summarized as follows. The health system suffers from a shortage of qualified personnel with one doctor per 34,750 inhabitants and one nurse for 3,500 inhabitants. 17.4% of patients do not have access to health care, while 81.5% of patients are forced to go into debt or sell property to pay the health costs. There is a big disparity between the capital Bujumbura and the remainder of the country, as 80% of doctors and more than 50% of nurses are working in Bujumbura. Malaria, which is responsible for more than 50% of hospital deaths in children under five years of age and 40% of all consultations in health centres, is undoubtedly the main public health problem, the main cause of mortality and morbidity in Burundi .
The goal/aim of this study is to assess which climatic factors affect malaria incidence in Burundi. Data on malaria morbidity in Burundi were collected from EPISTAT (Epidemiology and Statistics) , a department of the Burundi Ministry of health, in charge of collecting and storing data on epidemiology all over the country. The health services collect all monthly notification of malaria consultations. Malaria morbidity data/number from 1996 to 2007 were collected. This is the period where data are available in EPISTAT. The well-known nearest neighbour method is used to fill the (~5%) missing data. The estimated population for each province, for the study period, was obtained from the Institute of Statistics and Economic Studies of Burundi (ISTEEBU) . In this study, the number of malaria cases for each province was divided by the total number of population of the province to obtain the incidence rate. Used are the incidence rates per 1,000 inhabitants.
Monthly cumulative precipitation, monthly average of maximum temperature, monthly average of minimum temperature, monthly average of maximum humidity and monthly average of minimum humidity, for the period 1996-2007 were obtained from the Geographic Institute of Burundi (IGEBU) . The record of these variables from 1996 to 2007 has remained uniform, with the same calibration and the same precision. Only 14 stations had enough data for the desired study period. The missing data (2% - 3%) were filled by the same method as in Malaria data (nearest neighbour and cross-validation). Data for three provinces (Bubanza, Bujumbura rural and Cibitoke) were not available for the study period. They are estimated using ordinary kriging .
Modelling the impact of climatic variables on malaria
Generalized linear and generalized additive mixed models
where h is a known link function and γ is a vector of unknown regression parameters.
In most practical regression applications however, the assumption of a strictly linear effect (of the covariates) on the predictor may not be appropriate for continuous covariates.
Here L(.) denotes the likelihood which, under the assumption of conditional independence, is the product of individual likelihood contributions. Markov Chain Monte Carlo (MCMC) simulation methods that allow to draw random samples from the posterior are used to compute the posterior distribution [35, 36].
The aim is to model the dependence of malaria incidence on covariates including rainfall, minimum and maximum temperature, minimum and maximum humidity in Burundi. In the model, the variation of temperature, rainfall and humidity within one province is assumed not significant. This assumption is dictated by the fact that provinces in Burundi are small (17 provinces on an area of 27,834 km2). Taking into account the life cycle of the parasite (13 days for Plasmodium falciparum)  and the incubation period (seven days to four weeks), it is assumed that malaria incidence in a given month is associated with climatic conditions of the same month and those of the previous month. Most of those who become ill in a given month were bitten by mosquitoes in the previous month. The data are available in different scale and units (malaria and humidity data are unit-free, rainfall is measured in centimetres (cm) and temperature in degrees centigrade (°c)). To avoid the effect of scale in our modelling, the data are first standardized. Throughout this study we adopt the following notation for the variables:R n is the rainfall, T x is the maximum temperature, T n is the minimum temperature, H x is the maximum humidity, H n is the minimum humidity and R np , T xp , T np , H xp , H np are the same variables for the previous month.
Generalized linear model
η it is the predictor of malaria incidence assumed to have a gamma distribution,
Y it is the malaria incidence in province i (i = 1,...,17) and month t(t = 1,...,144),
β is a vector of unknown regression parameters β k ,
The regression coefficients β k (k = 0,...,5) are estimated via Markov Chain Monte Carlo (MCMC) simulation that allows to draw random samples from the posterior given by (6).
12,000 iterations of the MCMC are ran with a burn-in phase of 2,000 iterations. A thinning to the Markov Chain is applied to reduce autocorrelations, by requiring the programme to store only every 10th sampled parameter. This leads to a random sample of length 1,000 for every parameter in the model .
Estimates of the coefficients of the linear effects in the generalized linear model.
95% Confidence Interval
Generalized additive mixed model
Here, as above, η it is the predictor of malaria incidence in month t and province i, assumed to have a gamma distribution, X kit and X kip are the five covariates in province i for month t and the previous month p, respectively, f i are unknown smooth functions of the covariates, β0 is the intercept and ε it is the error. The effects of two-factors interaction are assumed to be smaller and are consequently omitted. The main reason is the wish to preserve the simplicity and easy interpretation of the effects, which are often lost by including interactions [29, 38].
The regression coefficients α k (k = 0,...,3) and the nonlinear effect of the continuous covariates are estimated via Markov Chain Monte Carlo (MCMC) as above. The effects of the continuous covariates are modelled by cubic p-splines  with 20 equidistant knots and second order random walk penalty [38, 41]. Positive hyperparameters a = 0.001 and b = 0.0005 for τ2 are chosen to ensure the propriety of the posterior [35, 40, 42]. The sensitivity of the results with respect to changes in the hyperparameters a and b is then checked by re-estimating the model with different choices for the hyperparameters a and b for each effect in the model by setting
(a = 1, b = 0.005); (a = 0.001, b = 0.001); (a = 0.01, b = 0.005); (a = 0.001, b = 0.005) (a = 0.0001, b = 0.0001); (a = 0.0001, b = 0.0005) to assess the dependence of results on minor changes in the model assumptions. The results showed no significant change.
The models are implemented in BayesX, a public domain software for Bayesian inference in structured Additive Regression Models available at http://www.stat.uni-muenchen.de/~bayesx/bayesx.html.
The purpose of this study has been to undertake a spatial/longitudinal statistical analysis to identify important climatic variables that influence malaria incidences in Burundi. Table 1 presents the estimates of the coefficients β i from model (9).
In Table 1, β1, β2 and β3 have negative means. Moreover β1 and β2 have a negative 95% credible interval (CI), suggesting that rainfall and maximum temperature in a given month have a negative effect on malaria incidence of the same month. β4 has a positive mean with a positive 95% CI. This suggests malaria incidence in a given month is positively and significantly associated with minimum temperature of the previous month.
Estimated coefficients of the linear effects in the generalized additive mixed model.
Overall, both models (GLM and GAMM) suggest that malaria incidence in a given month in Burundi is significantly positively associated with the minimum temperature of the previous month. Conversely, rainfall and maximum temperature in a given month are found to have a negative effect on malaria incidence of the same month. Possible explanations of these results are the following. Due to the development cycle of the parasite into mosquitoes and the incubation period, those who became ill in a given month were bitten by mosquitoes in the previous month. Minimum temperature is observed at night and mosquitoes are active only at night; by day-time they hide themselves in houses or vegetation. This explains why this factor, which is observed at the same time (at night) when there is mosquito activity, has a great influence on malaria transmission. Moreover, when the night temperature is high, people do not cover themselves, increasing the risk of being bitten by mosquitoes. Maximum temperature has a negative effect because mosquito's development is interrupted at higher temperature . Too much rainfall may flush away breeding larvae, reducing the numbers of the disease vectors.
This study has exploited available real monthly data on malaria, rainfall, temperature and humidity over 12 years (1996-2007) in the area of Burundi to derive and propose a regression modelling to assess climatic factors that are associated with monthly malaria incidence. The regression is based on the generalized linear model (GLM) and generalized additive mixed model (GAMM). The results obtained suggest a strong positive association between malaria incidence in a given month and the minimum temperature (night temperature) of the previous month. An open question is therefore how to cope with high temperatures at night. One possible suggestion in the context of Burundi would be to popularize the use of bed nets in order to reduce contact with mosquitoes. In contrast, it is found that rainfall and maximum temperature in a given month have a negative effect on malaria incidence of the same month. It is believed that climatic variables alone cannot explain the spatio-temporal distribution of malaria cases in Burundi. An interesting problem under investigation is extending the results in this work by incorporating in the proposed model the spatial effect by accounting for both structured (correlated) and unstructured (uncorrelated) components.
This study was carried out in the frame of the doctorate studies of the first author (HN). He is extremely grateful to the Austrian agency for international cooperation in education and research (ÖAD) for offering him the opportunity to study in Austria through providing him with the necessary financial support. HN would also like to acknowledge the valuable cooperation and support from the workers in Epistat and IGEBU in Burundi by providing him with the data used in this study.
- Gerstl S, Cohuet S, Edoh K, Brasher C, Lesage A, Guthmann JP, Checchi F: Community coverage of an antimalarial combination of artesunate and amodiaquine in Makamba Province, Burundi, nine months after its introduction. Malar J. 2007, 6: 94-10.1186/1475-2875-6-94.PubMed CentralView ArticlePubMedGoogle Scholar
- Ndayiragije A, Niyungeko D, Karenzo J, Niyungeko E, Barutwanayo M, Ciza A, Bosman A, Moyou-Somo R, Nahimana A, Nyarushatsi JP, Barihuta B, Mizero L, Ndaruhutse J, Delacollette C, Ringwald P, Kamana J: Efficacité de combinaisons thérapeutiques avec des dérivés de l'artémisinine dans le traitement de l'accès palustre non-compliqué au Burundi. Trop Med Int Health. 2004, 9: 673-679. 10.1111/j.1365-3156.2004.01255.x.View ArticlePubMedGoogle Scholar
- WHO: Stratégie de coopération de l'OMS avec les pays. République du Burundi. 2005Google Scholar
- UNDP: Rapport national sur le Développement humain au Burundi. 2003Google Scholar
- Greenwood B, Mutabingwa T: Malaria in 2002. Nature. 2002, 415: 670-672. 10.1038/415670a.View ArticlePubMedGoogle Scholar
- Gomez-Elipe A, Otero A, Van Herp M, Aguirre-Jaime A: Forecasting malaria incidence based on monthly case reports and environmental factors in Karuzi, Burundi, 1997-2003. Malar J. 2007, 6: 129-10.1186/1475-2875-6-129.PubMed CentralView ArticlePubMedGoogle Scholar
- Wandiga SO, Opondo M, Olago D, Githeko A, Githui F, Marshall M, Downs T, Opere A, Yanda PZ, Kangalawe R, Kabumbuli R, Kirumira E, Kathuri J, Apindi E, Olaka L, Ogallo L, Mugambi P, Sigalla R, Nanyunja R, Baguma T, Achola P: Vulnerability to climateinduced highland malaria in East Africa. AIACC Working Paper No. 25. 2006Google Scholar
- Gemperli A: Development of spatial statistical methods for modelling point-referenced spatial data in malaria epidemiology. PhD Thesis. 2003, Basel UniversityGoogle Scholar
- Zhou G, Minakawa N, Githeko A, Yan G: Association between climate variability and malaria epidemics in the East African highlands. Proc Natl Acad Sci USA. 2004, 101: 2375-2380. 10.1073/pnas.0308714100.PubMed CentralView ArticlePubMedGoogle Scholar
- Pemola ND, Jauhari RK: Climatic variables and malaria incidence in Dehradun, Uttaranchal, India. J Vector Borne Dis. 2006, 43: 21-28.Google Scholar
- Gallup JL, Sachs JD: The economic burden of malaria. Am J Trop Med Hyg. 2001, 64: 85-96.PubMedGoogle Scholar
- Kazembe LN, Kleinschmidt I, Holtz TH, Sharp BL: Spatial analysis and mapping of malaria risk in Malawi using point-referenced prevalence of infection data. Int J Health Geographics. 2006, 5: 1-9. 10.1186/1476-072X-5-1.View ArticleGoogle Scholar
- Kleinschmidt I: Spatial statistical analysis, modelling and mapping of malaria in Africa. PhD thesis, South Africa. 2001Google Scholar
- Craig MH, Sharp BL, Mabaso ML, Kleinschmidt I: Developing a spatial-statistical model and map of historical malaria prevalence in Botswana using a staged variable selection procedure. Int J Health Geographics. 2007, 6: 1-15. 10.1186/1476-072X-6-1.View ArticleGoogle Scholar
- Hay SI, Cox J, Rogers DJ, Randolph SI, Stern DI, Shanks GD, Myers MF, Snow RW: Climate change and the resurgence of malaria in the East African highlands. Nature. 2002, 415: 905-909. 10.1038/415905a.PubMed CentralView ArticlePubMedGoogle Scholar
- Patz JA, Lindsay SW: New challenges, new tools: the impact of climate Change on infectious diseases. Current Opinion in Microbiology. 1999, 2: 445-451. 10.1016/S1369-5274(99)80078-2.View ArticlePubMedGoogle Scholar
- Richard T: Malaria and climate change. Working papers series, Julian Simon Centre for Policy Research. 2002Google Scholar
- Cox J, Hay SI, Abeku TA, Checchi F, Snow RW: The uncertain burden of Plasmodium falciparum epidemics in Africa. Trends Parasitol. 2007, 23: 142-148. 10.1016/j.pt.2007.02.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Collein R: Implications of climate change on malaria in Karnataka, India. 2000, Senior Honors Thesis in Environmental Science, IndiaGoogle Scholar
- Checchi F, Cox J, Balkan S, Tamrat A, Priotto G, Alberti KP, Guthmann JP: Malaria epidemics and interventions, Kenya, Burundi, southern Sudan and Ethiopia 1999-2004. Emerg Infect Dis. 2006, 12: 1477-1485.PubMed CentralView ArticlePubMedGoogle Scholar
- Guthmann JP, Bonnet M, Ahoua L, Dantoine F, Balkan S, van Herp M, Tamrat A, Legros D, Brown V, Checchi F: Deaths rates from malaria epidemics, Burundi and Ethiopia. Emerg Infect Dis. 2007, 13: 140-143. 10.3201/eid1301.060546.PubMed CentralView ArticlePubMedGoogle Scholar
- Protopopoff N, Van Herp M, Maes P, Reid T, Baza D, D'Alessandro U, Van Bortel W, Coosemans M: Vector control in a malaria epidemic occurring within a complex emergency situation in Burundi: a case study. Malar J. 2007, 6: 93-10.1186/1475-2875-6-93.PubMed CentralView ArticlePubMedGoogle Scholar
- Grover-Kopecm EK, Blumenthal MB, Ceccato P, Dinku T, Omumbo JA, Cannor SJ: Web-based climate information resources for malaria control in Africa. Malar J. 2006, 5: 38-10.1186/1475-2875-5-38.View ArticleGoogle Scholar
- Encyclopedia of the Nations. [http://www.nationsencyclopedia.com/Africa/Burundi-CLIMATE.html]
- Ministry of Public Health in Burundi, EPISTAT:Google Scholar
- Ministry of Economy, Finance and development co-operation, ISTEEBU:Google Scholar
- Ministry of Planning and Environment in Burundi, IGEBU:Google Scholar
- Matheron G: Principles of geostatistics. Economic and Geology. 1963, 58: 1246-1267. 10.2113/gsecongeo.58.8.1246.View ArticleGoogle Scholar
- Fahrmeir L, Tutz G: Multivariate statistical modelling based on generalized linear models, Springer. 2001View ArticleGoogle Scholar
- Hastie TJ, Tibshirani RJ: Generalized additive models. 1997, Chapman & HallGoogle Scholar
- Diggle JP, Ribeiro JP: Model-based geostatistics, Springer. 2007Google Scholar
- Fahrmeir L, Lang S: Bayesian semiparametricregression analysis of multicategorical time-space data. Annals of the Institute of Statistical Mathematics. 2001, 53: 11-30. 10.1023/A:1017904118167.View ArticleGoogle Scholar
- Dobson AJ: An introduction to generalized linear models. 2002, Chapman & Hall, 2Google Scholar
- Kneib T, Fahrmeir L: Structured additive regression for categorical space-time data: A mixed model approach. Biometrics. 2006, 62: 109-118. 10.1111/j.1541-0420.2005.00392.x.View ArticlePubMedGoogle Scholar
- Brezger A, Kneib T, Lang S: Analysing Bayesian structured additive regression models. J Stat Software. 2005, 14: 1-22.View ArticleGoogle Scholar
- Brezger A, Lang S: Generalized structured additive regression based on Bayesian P-splines. Computational Statistics & Data Analysis. 2006, 50: 967-991. 10.1016/j.csda.2004.10.011.View ArticleGoogle Scholar
- Stone CJ: Additive regression and other nonparametric models. Ann Statistics. 1985, 13: 689-705. 10.1214/aos/1176349548.View ArticleGoogle Scholar
- Huang HC, Chen CS: Optimal geostatistical model selection. J Am Stat Assoc. 2007, 102: 1009-1024. 10.1198/016214507000000491.View ArticleGoogle Scholar
- Belitz C, Lang S: Simultaneous selection of variables and smoothing parameters in structured additive regression models. Computational Statistics & Data Analysis. 2008, 53: 61-81. 10.1016/j.csda.2008.05.032.View ArticleGoogle Scholar
- Fahrmeir L, Kneib Tand , Lang S: Penalized structured additive regression for space-time data. A Bayesian perspective. Statistica. 2004, 14: 731-761.Google Scholar
- Sachs J, Malaney P: The economic and social burden of malaria. Nature. 2002, 415: 680-685. 10.1038/415680a.View ArticlePubMedGoogle Scholar
- Lang S, Brezger A: Bayesian P-splines. Journal of Computational & Graphical Statistics. 2004, 13: 183-212. 10.1198/1061860043010.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.