- Research article
- Open Access
- Open Peer Review
This article has Open Peer Review reports available.
Does the pharmacy expenditure of patients always correspond with their morbidity burden? Exploring new approaches in the interpretation of pharmacy expenditure
© Calderón-Larrañaga et al; licensee BioMed Central Ltd. 2010
Received: 7 July 2009
Accepted: 11 May 2010
Published: 11 May 2010
The computerisation of primary health care (PHC) records offers the opportunity to focus on pharmacy expenditure from the perspective of the morbidity of individuals. The objective of the present study was to analyse the behaviour of pharmacy expenditure within different morbidity groups. We paid special attention to the identification of individuals who had higher values of pharmacy expenditure than their morbidity would otherwise suggest (i.e. outliers).
Observational study consisting of 75,574 patients seen at PHC centres in Zaragoza, Spain, at least once in 2005. Demographic and disease variables were analysed (ACG® 8.1), together with a response variable that we termed 'total pharmacy expenditure per patient'. Outlier patients were identified based on boxplot methods, adjusted boxplot for asymmetric distributions, and by analysing standardised residuals of tobit regression models.
The pharmacy expenditure of up to 7% of attendees in the studied PHC centres during one year exceeded expectations given their morbidity burden. This group of patients was responsible for up to 24% of the total annual pharmacy expenditure. There was a significantly higher number of outlier patients within the low-morbidity band which matched up with the higher variation coefficient observed in this group (3.2 vs. 2.0 and 1.3 in the moderate- and high-morbidity bands, respectively).
With appropriate validation, the methodologies of the present study could be incorporated in the routine monitoring of the prescribing profile of general practitioners. This could not only enable evaluation of their performance, but also target groups of outlier patients and foster analyses of the causes of unusually high pharmacy expenditures among them. This interpretation of pharmacy expenditure gives new clues for the efficiency in utilisation of healthcare resources, and could be complementary to management interventions focused on individuals with a high morbidity burden.
Analyses of trends in pharmacy expenditures are frequently carried out from a population perspective, focusing on the gross increase in costs [1, 2]. However, the computerisation of primary health care (PHC) records and the resulting large databases offer an excellent opportunity to address this problem considering the health characteristics of individual patients  such as the type and burden of morbidity .
The development of patient classification systems has enabled measurement and profiling the morbidity of patients and populations according to their complexity and expected utilisation of healthcare resources . That is, patients with similar morbidity profiles are likely to require equal delivery of healthcare resources (e.g. prescriptions). The Adjusted Clinical Groups (ACGs) system developed at Johns Hopkins University (Baltimore, MD, USA) is considered to offer great potential within PHC due to its holistic and longitudinal approach in characterising the morbidity burden of patients as well as its capacity to describe the case-mix of a reference population [6, 7]. The commonest applications of ACGs comprise assessment of provider performance, allocation and utilisation of resources, and outcomes analyses [8, 9].
When analysing pharmacy expenditure, the tools mentioned above can also be used to study the behaviour of expenditure within the different categories of health status. Further, focusing on individuals whose expenditure is higher than the average within one specific category of health status seems to be of particular interest. This group of individuals, statistically named as 'outliers' , can also be understood as those whose pharmacy expenditures exceed their burden of morbidity. This new approach is important: (i) for identifying patients for whom pharmacy care and its potential cost may be optimised; and (ii) because it opens up new research channels related to the characterization and prediction of patients and professionals on whom cost-control measures should be focused.
The objective of the present study was to analyse the behaviour of pharmacy expenditure within different morbidity groups. We paid special attention to the identification of individuals who had higher expenditure values than their morbidity would otherwise suggest.
A retrospective observational study was conducted based on information provided by computerised clinical records from six PHC centres in Zaragoza, Spain, in 2005. In an effort to obtain a homogeneous sample with respect to the quality of the records, we only selected centers with ≥ 2 years of experience in the use of computerised clinical records (as compiled through the OMI-AP system). We only included patients aged > 14 years who had been assigned to these centres and who attended such centres at least once during the study year. The final study population was composed of 75,574 patients.
Data were obtained from administrative registries of the Aragon Health Care System (Zaragoza, Spain) after official request and authorization. Personal information was made anonymous according to the Spanish Organic Law of Personal Data Protection 15/1999. This work was part of a project funded by the Carlos III Health Institute, which was approved by the Ethics Committee for Clinical Investigation of Aragon (CEICA).
Demographic variables (age, sex), pharmacy expenditure, and burden of morbidity were recorded for each patient. For the latter, all diagnoses or reasons for visits assigned to each patient during the study period were coded according to the International PHC Classification (IPHCC) . Subsequent conversion ('mapping') was made from the IPHCC to the International Disease Classification (ICD-9-CM) .
Based on the variables of age, sex and diagnosis registered along the study year, a single ACG category was assigned to each patient . The ACG system (version 8.1®) assigns all ICD-9-CM codes to one of 34 diagnosis clusters known as 'Aggregated Diagnosis Codes' (ADGs). Individual diseases or conditions are placed into a single ADG cluster based on the: duration of the condition (e.g. acute, recurrent, or chronic); severity of the condition (e.g. minor and stable versus major and unstable); diagnostic certainty (symptoms versus documented disease), aetiology of the condition (infectious, injury, or other); and specialty care involvement (e.g. medical, surgical, obstetric, haematologic). Each ADG is a grouping of diagnosis codes that are similar in terms of severity and likelihood of persistence of the health condition treated over a relevant period of time. In a second stage, the ACG methodology uses a branching algorithm to place subjects into one of 106 discrete categories based on their assigned ADGs, their age, and their sex. The result is that individuals within a given ACG have experienced a similar pattern of morbidity and resource utilisation over a given year. For the sake of parsimony, ACGs with similar expected use of resources are aggregated into 'morbidity bands' (i.e. low, moderate or high).
The primary measurement variable, which we termed 'total pharmacy expenditure per patient', included the retail price of medicines, supplies and accessories consumed by patients. This was obtained by combining data between the billing database of the pharmacy offices and the OMI-AP patient database. The tables were cross-referenced via the health card identifier (a sequential numeric tag assigned to each patient). All identifying data were removed to guarantee confidentiality.
Identification of outlier patients
There are many statistical procedures for identifying outlier data. Most are designed to prevent a 'masking' effect by groups of outliers or to be applied with multivariate data . The aim of the present study was to identify subgroups of 'unusual' patients in terms of their pharmacy expenditure considering their morbidity burden given the absence of a well-defined external cutoff, rather than the deletion of suspicious data or a robust estimation of location or scale parameters. The applied methods were chosen on the basis of a balanced compromise between their simplicity and their adequacy to the nature of the dependent variable 'pharmacy expenditure'. The statistical procedures for identifying outlier data are detailed below.
(i) Stratified boxplot by ACG (BXP)
where Q 3i and Q 1i are quartiles three and one, respectively, of the distribution of pharmacy expenditure for each ACG. Patients whose spending during the study period exceeded the trimming point (T i ) respective to their ACG category were considered to be outliers.
The constant 1.5 is the commonest value applied in Boxplot techniques. It is widely acknowledged that 0.35% of normally distributed data exceeds the trimming point (T i ) for moderate sample sizes (rate of false-positives). This percentage can be highly increased in non-normal or asymmetric distributions.
(ii) Adjusted boxplot for asymmetric distributions and stratified by ACG (Adj. BXP)
(iii) Standardised residuals (RESID)
Standardization of these residuals is applied to eliminate the effect of the measurement units on account of the dependent variable . If the normality assumption is fulfilled, the percentage of standardised residuals that fall out of the interval ± 2 is ~5%. Thus, cases with a residual > 2 are considered 'extreme' and therefore identified as outliers .
Regarding the regression model, the logarithm of pharmacy expenditure was used as the dependent variable to standardise the distribution . The independent 'ACG category' variable was introduced into the model as a group of dummy variables. The 'pharmacy expenditure' interest variable had a truncation point in the nil expenditure (7.8% of the population recorded nil expenditure), so truncated tobit models were applied . In these cases, the distribution of the variable in question is modelled as a mixture of a continuous distribution and a discrete distribution .
Demographic, clinical, and pharmacy-expenditure characteristics of outlier patients
A descriptive analysis was first carried out on a global basis and for each group of patients (normal/outliers) whereby the mean, median, and variation coefficient of the continuous variables and the distribution of frequencies of the categorical variables were calculated. The differences between outlier patients and the rest were analysed as well as the differences between outlier groups themselves (dependent upon the applied identification method) via the Mann-Whitney U test if the variable was continuous but did not follow a normal distribution (age, mean annual expenditure per patient, and mean annual number of visits per patient) and by a chi-square test for the qualitative sex and morbidity burden variables. Later, a residual study for these latter variables was undertaken in the event of significant differences.
Demographic, clinical, and pharmacy-expenditure characteristics of outlier patients according to the identification method (n = 75,574).
Patients with normal pharmacy expenditure
Patient proportion (%)
Pharmacy expenditure characteristics
Annual average expenditure/patient (€)
Proportion of total expenditure (%)
The STATA 10 software package was used to calculate the truncated tobit models, R free software to calculate function h r (Medcouple, MC) and SPSS 15.0 for the remainder of our statistical analyses.
Summary of the variation in pharmacy expenditure within morbidity bands.
N° of patients
Std. Dev. (€)
Comparisons between normal patients and outlier patients with regard to the annual pharmacy expenditure
There were statistically significant differences between normal patients and outlier patients (regardless of the method used to identify these groups). This is corroborated by the observed/average ratios which, for this group, suggested pharmacy expenditures 3.8-6.1-times greater than the global mean (Table 1, observed/average ratios). The analysis of the standardised residuals of the chi-square test indicated that there was a significantly higher presence of outlier patients within the low-morbidity band which matched up with the higher variation coefficient observed in this group (3.2 vs. 2.0 and 1.3 in the moderate- and high-morbidity bands, respectively). The expected/average ratio showed that the outlier population was 30-60% less ill than the general reference population (Table 1, expected/average ratios). As a result, the observed/expected ratios were 5-13-times greater for outlier patients than for patients with a normal expenditure (Table 1, observed/expected ratios).
Comparisons between groups of outlier patients dependent upon the identification method
The results of the present study indicated that pharmacy expenditure of up to 7% of the population attending the studied PHC centres during one year exceeded expectations given their morbidity burden, and that this group of patients was responsible for up to 24% of the total annual pharmacy expenditure. This information may be of interest in the context of growing concern about the increase in pharmacy expenditures in PHCs, and the application of possible measures for efficient and rational management of drug prescriptions.
Among the existing pharmacy cost-containment policies in Europe, some are focused on one of the key stakeholders in healthcare: the prescriber. This is achieved through the analysis of the quality and costs of their prescribing . Several countries disseminate reviews of drug utilisation to general practitioners (GPs) on the basis that prompt, detailed feedback of individualised prescribing data may be effective in rationalising the use of drugs . The concepts and methodologies shown in the present study could be incorporated in the routine monitoring of the prescribing profile of GPs. This could be used to evaluate not only their performance, but also to target groups of outlier patients and foster analyses of the causes of unusual pharmacy expenditures among them, thereby furthering the design of specific interventions.
The results of the present study indicate that there was a potential margin to improve drug utilisation in the group of patients with a low morbidity burden. That is, there was a lower correspondence between the morbidity profile of patients and their expenditure in those groups of patients who showed a low morbidity burden. A wider variation of pharmacy expenditure was observed within the group of patients with a low morbidity burden (variation coefficient: 3.2 vs. 2.0 and 1.3 in the moderate- and high-morbidity bands, respectively), leading to a higher number of patients whose figures were far removed from the expenditure of most patients belonging to the same group (i.e. outliers). The expenditure of this group of outliers was up to 13-times higher than expected given their burden of morbidity. This interpretation of pharmacy expenditure gives new clues for the efficiency in utilisation of healthcare resources, and could be complementary to management interventions focused on individuals with a high morbidity burden [25, 26].
The trimming methods of the present study have been widely applied in the hospital setting to avoid incorporating outlier cases when valuating the mean utilisation of resources of a given healthcare resource group [27, 28]. However, analysing the group of outliers itself deserves special attention. All three methods are easy to implement, adjusted by morbidity burden and adapted to the large sample size and the non-parametric distribution of the studied variable pharmacy expenditure, which is why other tests (e.g. Grubbs' test, Dixon's Q test) were not used. The three methodologies identify subgroups of the same population, reinforcing the robustness of the selected techniques. However, which method best identifies the group of outlier patients? The empirical evidence suggests that the most sophisticated statistical methods, such as the Adj. BXP for asymmetric distributions or standardised residuals, offer more specificity in identifying outlier patients. The average expenditure of patients detected by these two methods was 800-1000€ higher per year than that of patients detected by the classic BXP method. Moreover, these expenditures had a lower correspondence with their burden of morbidity as reflected in the ratios between the observed and expected expenditure, which were almost twofold greater than that given by the classic BXP method.
According to the World Health Organisation, the rational use of drugs requires that patients receive medications appropriate to their clinical needs at the lowest cost to them and their community. The evidence indicates that there may be factors associated to the practitioner and the patient which deviate them from the rational use of drugs, resulting in an unjustifiable use of scarce resources and widespread health hazards. Among the variables related to the practitioner, the continuous training, the influence of the pharmaceutical industry, and the 'burnout syndrome' as a consequence of chronic stress at the workplace have been identified most frequently . Regarding patient-related factors, the pressure that they can occasionally exert to obtain medications that 'cure' their health problem has been discussed [30, 31]. This problem is aggravated when practitioners' time to see a patient is limited . The pressure can be very difficult to handle because, in most cases, there is little time to explain the preferred treatment option of GPs to sceptical patients .
These findings must be interpreted with caution due to several reasons that may affect their internal and external validity. The main potential limiting factor of the present study could be related to the quality of the diagnostic information registered in the electronic medical records. Even if a series of inclusion criteria were applied during the health centre selection process to guarantee the quality and reliability of clinical data, a two-year period of experience in the use of electronic records may be insufficient. This could explain why patients have pharmacy expenditures for diseases that are not registered in the medical records of GPs, thus over-estimating the lack of correspondence between expenditure and morbidity burden of the population. Using more recent data is required to confirm the internal validity of the results.
Second, there were limitations derived from the study design owing to the restricted access to information. That is, drug use was evaluated only in terms of the expenditure. This ignored aspects related to the number of drug prescriptions or the number of defined daily doses. This variable would have been of great help when analysing the causes of excess use of drugs within the same morbidity group.
A third factor that could compromise the internal validity of the results obtained in the present study was related to the explanatory capacity of the tool used to measure morbidity: ACGs. It could be argued that the ACGs perform weakly in low-morbidity bands or, further, that updates on the effectiveness and cost-effectiveness of health technologies (e.g. drugs) were not considered in the design process. Nevertheless, the statistical power of the ACG system developed at Johns Hopkins University has been internationally validated [34–37] and its different versions are regularly updated .
Fourth, the results of the present study correspond to a population treated in urban health centres. They are therefore not fully representative of the population treated across an entire health system at the PHC level, thus affecting their external validity. Nevertheless, the characteristics of the studied sample, such as the average number of diagnoses per patient , the intensity of use [5, 39], the proportion of women, the aging rate, and the distribution of patients according to ACGs  fit in with those reported in national and international studies.
If the methodology of the present study is to be used to target patients and prescribers amenable to intervention strategies, they should first be validated. This implies a clinical validation by expert panels as well as carrying out economic studies that shed light on the potential savings of this type of interventions.
We began this contribution by asking if the pharmacy expenditure of patients always corresponds with their morbidity burden. The present study identified a group of patients with low morbidity whose pharmacy expenditure exceeded expectations. Identification of the factors that influence this excess in pharmacy expenditure can improve the adaptation of pharmacy services and increase the efficiency of the PHC system.
The study was funded by a grant PI052752 from the Carlos III Health Institute and the Program for the Incorporation of Research Groups into the Spanish Health System (EMER 07/020). The authors thank the various professionals of the Health Centres for their daily reporting of data. This study could not have been conducted without their contribution.
- Rivers PA, Hall NG, Frimpong J: Prescription drug spending: contribution to health care spending and cost containment strategies. J Health Care Finance. 2006, 32: 8-19.PubMedGoogle Scholar
- Instituto Nacional de la Salud, Subdirección General de Coordinación Administrativa: Memoria 2000. Madrid. 2001Google Scholar
- Majeed A: Sources, uses, strengths and limitations of data collected in primary care in England. Health Stat Q. 2004, 21: 5-14.PubMedGoogle Scholar
- Starfield B: Threads and yarns: weaving the tapestry of comorbidity. Ann Fam Med. 2006, 4: 101-103. 10.1370/afm.524.View ArticlePubMedPubMed CentralGoogle Scholar
- Petersen LA, Pietz K, Woodard LD, Byrne M: Comparison of the predictive validity of diagnosis-based risk adjusters for clinical outcomes. Med Care. 2005, 43: 61-67.PubMedGoogle Scholar
- Weiner JP, Starfield BH, Steinwachs DM, Mumford LM: Development and application of a population-oriented measure of ambulatory care case-mix. Med Care. 1991, 29: 452-472. 10.1097/00005650-199105000-00006.View ArticlePubMedGoogle Scholar
- Starfield B, Weiner J, Mumford L, Steinwachs D: Ambulatory care groups: a categorization of diagnoses for research and management. Health Serv Res. 1991, 26: 53-74.PubMedPubMed CentralGoogle Scholar
- Juncosa S, Carrillo E, Bolibar B, Prados A, Gervas J: [Classification systems in case-mix groups in ambulatory care. Perspectives for our primary health care]. Aten Primaria. 1996, 17: 76-82.PubMedGoogle Scholar
- Kolb DS, Clay SB: Ambulatory care groupings: when, how, and the impact on managed care. J Ambul Care Manage. 1994, 17 (1): 29-38.View ArticlePubMedGoogle Scholar
- Hampel FR: Robust Statistics: The Approach Based on Influence Functions. 1986, New York: WileyGoogle Scholar
- Lamberts H, Wood M: ICPC: International Classification of Primary Care. 1987, Oxford: Oxford University PressGoogle Scholar
- Ministry of Health and Consumption: International Classification of Diseases, Ninth Revision, Clinical Modification. Madrid. 1999Google Scholar
- The Johns Hopkins ACG® System: Version 8.1. Reference Manual. Baltimore. 2006Google Scholar
- Hadi A, Rahmatullah A, Werner M: Detection of outliers. WIREs Comp Stati. 2009, 1: 57-70. 10.1002/wics.6.View ArticleGoogle Scholar
- Tukey JW: Exploratory data analysis. 1977, Reading, MA: Addison-WesleyGoogle Scholar
- Hoaglin D, Mosteller F, Tukey J: Understanding robust and exploratory data analysis. 1983, New York: WileyGoogle Scholar
- Hubert M, Vandervieren E: An adjusted boxplot for skewed distributions. Comput Stat Data Anal. 2008, 12: 5186-5201. 10.1016/j.csda.2007.11.008.View ArticleGoogle Scholar
- Gray G, Woodall W: The maximum size of standardized and internally studentized residuals in regression analysis. Am Stat. 1994, 2: 111-113. 10.2307/2684258.Google Scholar
- Belsley DA, Kuh E, Welsh RE: Regression diagnostics: identifying influential data and sources of collinearity. 1980, New York: WileyView ArticleGoogle Scholar
- Kleinbaum D, Kupper L, Muller K, Nizam A: Applied Regression Analysis and other Multivariable Methods. 1998, Pacific Grove CA: Duxbury Press, 3Google Scholar
- Tobin J: Estimation of relationships for limited dependent variables. Econometrica. 1958, 26-24.Google Scholar
- Clavero A, González ML: Una revisión de modelos econométricos aplicados al análsis de demanda y utilización de servicios sanitarios. Hacienda Pública Española. 2005, 173: 129-162.Google Scholar
- Ess S, Schneeweiss S, Szucs T: European Healthcare Policies for Controlling Drug Expenditure. Pharmacoeconomics. 2003, 2: 89-103. 10.2165/00019053-200321020-00002.View ArticleGoogle Scholar
- Farmer A, Légaré F, Turcot L, Grimshaw J, Arvey E, McGowan J, Wolf F: Printed educational materials: effects on professional practice and health care outcomes. Cochrane Database Syst Rev. 2008, 3-CD004398.Google Scholar
- Boyd CM, Boult C, Shadmi E, Leff B, Brager R, Dunbar L, Wolff JL, Wegener S: Guided care for multimorbid older adults. Gerontologist. 2007, 47: 697-704.View ArticlePubMedGoogle Scholar
- Morisky DE, Kominski GF, Afifi AA, Kotlerman JB: The Effects of a Disease Management Program on Self-Reported Health Behaviors and Health Outcomes: Evidence From the "Florida: A Healthy State (FAHS)" Medicaid Program. Health Educ Behav. 2009, 36 (3): 505-17. 10.1177/1090198107311279.View ArticlePubMedGoogle Scholar
- Pirson M, Dramaix M, Leclercq P, Jackson T: Analysis of cost outliers within APR-DRGs in a Belgian general hospital: two complementary approaches. Health Policy. 2006, 1: 13-25. 10.1016/j.healthpol.2005.04.008.View ArticleGoogle Scholar
- Cots F, Elvira D, Castells X, Saez M: Relevance of outlier cases in case mix systems and evaluation of trimming methods. Health Care Manag Sci. 2003, 6: 27-35. 10.1023/A:1021908220013.View ArticlePubMedGoogle Scholar
- Caamaño F, Figueiras A, Lado E, Gestal-Otero J: [Explanatory variables of avoidable cost produced by not prescribing a cheaper equivalent drug]. Gac Sanit. 2001, 6: 513-518.View ArticleGoogle Scholar
- Prosser HWT: A qualitative study of GPs' and PCO stakeholders' views on the importance and influence of cost on prescribing. Soc Sci Med. 2005, 60: 1335-1346. 10.1016/j.socscimed.2004.07.013.View ArticlePubMedGoogle Scholar
- Carthy P, Harvey I, Brawn R, Watkins C: A study of factors associated with cost and variation in prescribing among GPs. Fam Pract. 2000, 1: 36-41. 10.1093/fampra/17.1.36.View ArticleGoogle Scholar
- Lopez-Picazo J, Sanchez J, Rausell V, Sanz J, Sanchez F, Salas E: [Prescription of drugs in family medicine: on the trail of efficiency]. Aten Primaria. 2004, 4: 178-185. 10.1157/13065827.View ArticleGoogle Scholar
- Azpiazu M, García L: [Factors affecting drug expenditure in the primary care centres in a health area]. Aten Primaria. 2002, 2: 84-89.View ArticleGoogle Scholar
- Sibley LM, Moineddin R, Agha MM, Glazier RH: Risk adjustment using administrative data-based and survey-derived methods for explaining physician utilization. Med Care. 2010, 48 (2): 175-182. 10.1097/MLR.0b013e3181c16102.View ArticlePubMedGoogle Scholar
- Sicras-Mainar A, Navarro-Artieda R: [Validating the Adjusted Clinical Groups [ACG] case-mix system in a Spanish population setting: a multicenter study]. Gac Sanit. 2009, 23 (3): 228-231. 10.1016/j.gaceta.2008.04.005.View ArticlePubMedGoogle Scholar
- Zielinski A, Kronogard M, Lenhoff H, Halling A: Validation of ACG Case-mix for equitable resource allocation in Swedish primary health care. BMC Public Health. 2009, 9: 347-10.1186/1471-2458-9-347.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang HY, Weiner JP: An in-depth assessment of a diagnosis-based risk adjustment model based on national health insurance claims: the application of the Johns Hopkins Adjusted Clinical Group case-mix system in Taiwan. BMC Med. 2010, 8: 7-10.1186/1741-7015-8-7.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu L, Cline RR, Schondelmeyer SW, Schommer JC: Pharmaceutical expenditures as a correlate of population health in industrialized nations. Ann Pharmacother. 2008, 42: 368-374. 10.1345/aph.1K352.View ArticlePubMedGoogle Scholar
- Meenan RT, Goodman MJ, Fishman PA, Hornbrook MC, O'Keeffe-Rosetti MC, Bachman DJ: Using risk-adjustment models to identify high-cost risks. Med Care. 2003, 41: 1301-1312. 10.1097/01.MLR.0000094480.13057.75.View ArticlePubMedGoogle Scholar
- Engstrom SG, Carlsson L, Ostgren CJ, Nilsson GH, Borgquist LA: The importance of comorbidity in analysing patient costs in Swedish primary care. BMC Public Health. 2006, 6: 36-10.1186/1471-2458-6-36.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/10/244/prepub
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.