- Open Access
Rural-urban differentials of premature mortality burden in south-west China
International Journal for Equity in Health volume 5, Article number: 13 (2006)
Yunnan province is located in south western China and is one of the poorest provinces of the country. This study examines the premature mortality burden from common causes of deaths among an urban region, suburban region and rural region of Kunming, the capital of Yunnan.
Years of life lost (YLL) rate per 1,000 and mortality rate per 100,000 were calculated from medical death certificates in 2003 and broken down by cause of death, age and gender among urban, suburban and rural regions. YLL was calculated without age-weighting and discounting rate. Rates were age-adjusted to the combined population of three regions. However, 3% discounting rate and a standard age-weighting function were included in the sensitivity analysis.
Non-communicable diseases contributed the most YLL in all three regions. The rural region had about 50% higher premature mortality burden compared to the other two regions. YLL from infectious diseases and perinatal problems was still a major problem in the rural region. Among non-communicable diseases, YLL from stroke was the highest in the urban/suburban regions; COPD followed as the second and was the highest in the rural region. Mortality burden from injuries was however higher in the rural region than the other two regions, especially for men. Self-inflicted injuries were between 2–8 times more serious among women. The use of either mortality rate or YLL gives a similar conclusion regarding the order of priority. Reanalysis with age-weighting and 3% discounting rate gave similar results.
Urban south western China has already engaged in epidemiological pattern of developed countries. The rural region is additionally burdened by diseases of poverty and injury on top of the non-communicable diseases.
The allocation of health care resources has become an important issue in many countries experiencing resource limitation, and the assessment of the burden of disease (BOD) is useful for sensible allocation of limited health-related resources . China is a developing country where health resources are very scarce and inequality in health resource allocation between rural and urban regions is apparent [2, 3]. Urban populations, which account for only 30% of whole population, use 80% of total health resources, whereas rural populations, which account for 70% of total populations, only use 20% of total health resources .
Yunnan province is located in the south-west of China, and is one of the poorest provinces in the country. In 2004 the per capita GDP of Yunnan province was US$680, which was 37.6% lower than that of the whole country . Kunming prefecture, the capital, contains two urban, two suburban and eight rural regions, and in 2000 had a population of approximately 5.8 million.
There are three well established systems for death registry in Kunming prefecture, located in an urban region, a suburban region and a rural region respectively. The first two are the surveillance points of death registry of China, whereas the other is one of the surveillance points of the country birth registration. In order to maintain accuracy and completeness of birth and death registration, an independent resurvey on a sample of about 5000 households once every three years is conducted in each of these surveillance points. Compared to the surveillance data, the 2000 census in China actually had more underreported deaths, and was less complete due to the problem of recall bias. These well defined regions and established registration systems provide a good ground for comparison of mortality burden across the three types of regions along with economic transition.
Years of life lost (YLL) is one of the methods of estimating the duration of time lost due to premature death, and is the mortality component of disability adjusted life years (DALY). The conceptual and computational details of years of life lost have been presented elsewhere . The YLL measure not only considers the number of deaths, but also takes into account the age at which death occurred. It is therefore a better tool for quantifying the burden of premature mortality compared to mortality rate per se.
The purpose of this study is to compare the premature mortality burden from each common cause among an urban region, suburban region and rural region of Kunming, broken down by age and gender. This assessment should provide more insight into the details of health inequity across the rural/urban gradient.
Pan Long district, Guan Du district and Shi Lin county were selected as the study regions. Pan Long district is an urban region, with a population of 378,885 (195,590 males and 183,295 females); Guan Du district is a suburban region, with a population of 335,622 (178,099 males and 157,523 females), and Shi Lin county a rural region with a population of 179,746 (93,327 males and 86,419 females). All individuals residing and dying in these three regions in 2003 were included in the analysis.
Causes of death were based on medical death certificate information, maintained by Pan Long Center for Diseases Control, Guan Du Center for Diseases Control and Shi Lin maternal and children hospital. All death reports were grouped by underlying cause of death as defined in the Global Burden of Disease (GBD) study  and coded using the International Classification of Diseases, 9th revision (ICD-9) coding system.
In order to avoid some misreporting of age at death, and misclassification of cause of death, all medical death certificates were verified on the underlying cause of death by a team of two independent physicians. Any discrepancies were reviewed to obtain a consensus. All deaths assigned to ill-defined conditions among the three study regions were redistributed to other more defined causes according to the age and gender distribution of specific conditions, following the conceptual approach in the GBD study .
Overall mortality in Kunming was divided into three broad groups of causes: Group I, communicable, maternal, perinatal and nutritional deficiencies; Group II, non-communicable; and Group III, all injuries. These were then further subdivided into several more specific causes . Age was divided into five groups: 0–4, 5–14, 15–44, 45–59 and 60 years and over. Overall age-specific mortality for each sex was plotted for visual comparison.
Calculation of Years of Life lost (YLL) and mortality rate
In global burden of disease studies, YLL incorporates an age-weighting factor so that life lived in young adults is given a higher social value than that in childhood or older ages in most societies [9, 10]. However, the use of age-weighting has been criticized by some researchers . Traditionally, Chinese culture values living years in the elderly and children, so age-weighting was not used in this study. This study calculated YLL without discounting rate, because deaths occurred in each calendar year and the study has a cross-sectional nature. YLL was also calculated with a 3% discounting rate and a standard age-weighting function. The results of the final ranking were the same so they were omitted from this report and only the YLL calculated without discounting and age-weighting presented. Nevertheless, 3% discounting rate and a standard age-weighting function were included in the sensitivity analysis.
The simplified formula of YLL without discounting and age weighting is:
where d x is the number of deaths at age x, and e x is the life expectancy at that age.
Life expectancy was calculated using the life tables provided in the GBD study- the model life-table west level 26 , which is based on the life expectancy at birth of 82.5 years for females. West level 25 for females is used as the standard for males (in the absence of a Coale and Demeny life table with life expectancy of 80 for males), to ensure comparability with other studies.
Mortality rates were also calculated to compare with the use of YLL as an indicator of mortality burden.
Abridged life tables was used to calculate the average life expectancy at birth of three regions. Ninety-five percent confidence interval of the life expectancy was also computed using the Chiang methodology – Chiang (II) .
Calculations for adjusted rate
The direct standardization method was used to calculate both the age-adjusted YLL rate and the age-sex-adjusted mortality rate. The total population of all 3 study regions was used as the standard population. For comparison of age-specific mortality and YLL rates, the sex distribution was adjusted using the same standard population.
The average life expectancy at birth was 77.1 years (95% CI: 76.5–77.8) for the urban region, 74.3 years (95% CI: 73.5–75.1) for the suburban region and 70.8 years (95% CI: 70.0–71.5) for the rural region.
There were 2,032 deaths registered in the urban region (56.1% males, 43.9% females), 1,927 deaths in the suburban region (53.1% males, 46.9% females) and 1,065 deaths (53.4% males, 46.6% females) in the rural region of Kunming for 2003. Thirty-two deaths in all three regions combined, which were coded as ill-defined conditions, were redistributed proportionally to other more defined causes. All of them were aged over 70 years.
The overall age-specific mortality rates by region are illustrated in Figure 1 (males) and Figure 2 (females). Mortality rates were higher in the rural region throughout the whole age range in both sexes. The urban population had the lowest mortality rate among the three regions in childhood and early adulthood. The rates were similar for all three regions after 60 years of age for both sexes.
Table 1 shows a comparison of YLL by broad cause groups among the three regions. Group I conditions were responsible for the highest mortality burden in the rural region compared to the other two regions, as were group III conditions. Group II conditions accounted for over 50% of the total YLL in each of the three regions, indicating that non-communicable diseases were the major premature mortality burden in the study regions as a whole.
Table 2 presents age breakdown of sex-adjusted YLL/1,000 population by broad cause groups among the three regions. Group I accounted for the highest YLL/1,000 population in the 0–4 age group. Non-communicable diseases were responsible for the highest YLL/1,000 population in the 60+ age group among the three regions. The premature mortality from injuries was the highest in young adults aged 15–44 years in all three regions. The rural region had the highest premature mortality burden in all three death cause groups when the rate per thousand population was used.
Table 3 compares YLL/1,000 population by sex and age among the three regions of Kunming. People aged 60 years and over were responsible for the highest years of life lost compared to the other age groups in each of the three regions. Children aged less than 5 years also had relatively high premature mortality burden in all three regions. In the age group of 60 years and over, the rural region had the lowest value of years of life lost, whereas in the 0–4 years age group, the urban region was responsible for the lowest premature mortality burden.
Table 4 presents expanded information of group of causes of death in terms of adjusted mortality rate/100,000 and YLL/1,000 by region. Although the figures of these two indicators are on a different scale, the ranking within the region are similar. However, there are some differences in ranking with regions. Birth trauma and asphyxia and low birth weight rank low by both indicators in urban and suburban regions, but rank relatively higher by both indicators in the rural region. For group I diseases, birth trauma and asphyxia and low birth weight were important leading causes of years of life lost in the rural region, and had a serious impact on premature health. For group II diseases, stroke was the top leading cause of years of life lost in both the urban and suburban regions, whereas it ranked third in the rural region. COPD was also a leading cause of premature death among the three regions, being the second ranked cause of years of life lost in the urban and suburban region, and the top cause of years of life lost in the rural region. Diabetes was a common problem in the urban region, while ischemic heart disease was predominant in the suburban region. Lung and liver cancer were major causes of premature death in the urban and suburban regions, but were less common in the rural region. For group III diseases, road traffic accidents were among the top 10 leading causes of premature death among all three regions, particularly the rural region, where it ranked fourth in both mortality and YLL. Self-inflicted injury was a major premature mortality concern in the rural region.
Table 5 and 6 are a breakdown of Table 4 into male and female. For group II conditions, stroke, COPD, diabetes, lung cancer and liver cancer in both sexes had similar rankings as listed in Table 4. Geographic discrepancy of group I conditions was also not different between males and females. However, in the urban region, females accounted for higher years of life lost due to diabetes than males, while in the rural region ischemic heart disease was a more common cause of premature death among females than males. In the suburban region, females had higher years of life lost due to self-inflicted injuries than males. Road traffic accidents had a substantial impact on the health of males in each of the three regions, and resulted in a 2 to 3 times higher burden of years of life lost among males compared to females. For all causes, the number of years of life lost was higher in males than females in all three regions.
These findings indicate that the epidemiologic transition is well under way in the study communities. For all three regions, the leading causes of premature death are non-communicable diseases. There are great health disparities across the rural/urban gradient with the rural region having about a 50% higher premature mortality burden compared to the other two regions. Using either mortality rate or YLL gives similar conclusion in order of priority. The rural region is still suffering from infectious diseases and prenatal problems, with five times the burden rates of the other two regions. Among the top non-communicable diseases, stroke is the most common mortality burden in the urban and suburban regions. COPD follows as the second highest burden and is the top burden in the rural region. Ischemic heart disease, although having a lower rate, still poses an important mortality burden in the suburban region. Mortality burden from injuries is however higher in the rural region than in the other two regions. For men living in rural regions, death from road traffic accidents is a major concern. However, self-inflicted injuries are two to eight times higher among women in this setting.
In the rural region, chronic diseases have not displaced but added to the mortality burden from infectious and perinatal problems, and this double burden is a major challenge for health systems in this setting. These results are consistent with other studies in China [14, 15], and in other parts of the developing world, such as Pakistan . Our results differed from a study in sub-Saharan Africa , where childhood and infectious diseases are still the major mortality burden. In urban Kunming, ischemic heart disease, stroke, diabetes, COPD and lung cancer are among the top 10 leading causes of YLL. The situation is similar to that in developed countries such as USA and Switzerland [18, 19]. This may be partly explained by the higher socio-economic status, more hygienic life style and better access to health care. All of these suppress infectious diseases and prenatal problems but can lead to a higher life expectancy and therefore more chronic diseases. On the other hand, comparison among studies has to take into account the difference in discounting rate and age-weighting scheme. Our study calculated YLL without discounting rate and age-weighting, whereas those from developed countries such as US and Switzerland have used 3% discounting rate and age-weighting.
In the rural region, birth trauma and asphyxia and low birth weight are listed in the top 10 causes of YLL. This may be partly explained by the effect of the age distribution. In this study, the 0–1 age group occupies 1.7% of the whole population in the rural region, compared to 0.8% of the urban and 0.7% of the suburban population. The discrepancy in the proportions is due to a higher birth rate and lower life expectancy in the rural region . This two-fold difference in percentage of the infant population however does not explain the four to five time differences in mortality burden of this age group across the rural-urban gradient. China has achieved fairly steady economic growth over the past 20 years. However, the greater disparities in economic development among different regions have led to health disparities among the urban and rural regions, especially for maternal and infant mortality [21, 22].
The leading mortality burden for non-communicable diseases in our study regions include stroke, ischemic heart disease, diabetes, COPD, cancer and injuries. Stroke was a much higher mortality burden compared to ischemic heart disease. The dominance of stroke over ischemic heart disease has been repeatedly found in most parts of mainland China [23, 24]. However, the situation is different from that among ethnic Chinese in Hong Kong and Singapore , and Taiwan , where ethnic Chinese die more from ischemic heart disease than from stroke. Heart disease ranks third among the leading causes of mortality in Taiwan, and is a more common cause of death than stroke. These differences are less likely to be explained by genetic backgrounds than environmental or behavioural factors.
Our findings show that COPD is a major cause of death particularly in the non-urban region. Results of the 1996 National Prevalence Survey in China indicated that more than 60% of men in China currently smoked, whereas this figure was less than 3% among women . Tobacco is the biggest industry in Yunnan, accounting for 30% of China's tobacco production. In 2002 it was responsible for 17.6% of the national total output . More than 2.3 million farmers in the province are engaged in tobacco cultivation. Our results indicate that COPD is a substantial cause of premature death for both genders among all three regions. This is partly explained by the high prevalence of smoking and the fact that COPD accounts for 35% of all smoking related deaths . However, the females in our study had a slightly higher premature mortality burden due to COPD than males in the suburban and rural regions, pointing to the existence of risk factors other than smoking on COPD in women. In 2002 the World Health Organization estimated that 80% of Chinese households use biomass fuel, and coal used in 31% of households . Solid-fuel use has been shown to be an important risk factor for COPD . From household work, women are usually more exposed to this risk factor, which may contribute to the female premature mortality burden. The nature of this additional etiologic factor needs to be further explored.
This study demonstrates a higher premature mortality burden due to injuries in the rural region than in the urban region, especially fatal road traffic accidents among males. This could not be explained by the difference in rate of transportations, which is usually higher in urban regions. Several studies have indicated that the largest proportions of road traffic victims in developing countries are pedestrians, passengers, and cyclists as opposed to drivers, in whom most of the deaths and disabilities in the developed world occur [32, 33]. Sharing of the road by high speed vehicles and walking villagers in addition to other factors of road traffic accidents may be an explanation for the rural predominance .
Self-inflicted injuries are a serious premature mortality burden, especially among females in both the suburban and rural region. In contrast, the rate of suicide is almost the same for men and women in the urban region. The situation is different from those in Micronesia, Hong Kong China and Singapore, where women die less often by suicide than men . High suicide rates are due to multiple factors such as religious beliefs, social isolation, family problems, stressful life events, chronic painful diseases, mental disorders and substance abuse . Psychiatric services have been shown to reduce the rate . Further studies to identify the causes of these injuries and to improve the mental services are needed in our study region.
There are many methodological lessons from this study. Firstly, a good mortality registration system can be accomplished in a developing country such as China. Secondly, mortality rate gives similar results to the standard YLL computation for Group II and group III conditions in terms of priority ranking once the country is overwhelmed by a group of diseases such as diabetes, ischemic heart disease and stroke, which affect a relatively homogenous age group. However, in the longer term and over a broader area, YLL would still be more appropriate when the patterns of disease are less homogeneous. Thirdly, this kind of study can be useful in uncovering less known mortality burdens, such as those from COPD, self-inflicted injuries and road traffic accidents. These can invoke the interest of non-health sector policy makers to get involved.
The strength of this study hinges on the complete vital registration systems, since the study regions are surveillance points for vital statistics in China. The problem of underreporting of deaths found in many studies is thus minimized. Underreporting of deaths has been shown to be more common in infant deaths in a previous study, especially in rural regions . Had this been the case in the current study, the mortality disparity in group I causes would have been underestimated. This study is confined to mortality burden. While this type of burden declines with improvements to the social-economic status of a region, disability burden may be a more important problem in other contexts . Such burden will be the focus of future research. Problems identified solely on the basis of mortality data may be underestimated. For example, in a study in Pakistan , injuries ranked eleventh according to YLL but second according to YLD, resulting in their ranking fifth based on DALY.
The findings suggest that a strong health advocacy should be applied to rural regions of Yunnan province, especially on group I diseases and group III injuries. A continual and consistent effort in prevention and measures to reduce the burden from metabolic syndrome in all regions should be strengthened.
World Development Report 1993: investing in health. 1993, New York: Oxford University Press
Fan SG, Zhang LX, Zhang XB: Growth, Ineqyality, and Povety in Rural China: the Role of Public Investments. IFPRI Research Report 125. 2002
Yu DZ: Changes in Health Care Financing and Health Status: The Case of China in the 1980s. 1992, Economic Policy Series, No 34, International Child Development Centre, UNICEF
Wu SL, Zhu AR: The direction and method to allocate health resources for Chinese government. Journal of Chinese Rural Health Management. 1999, 19 (7): 33-34.
Yunnan Statistical Yearbook: 2004: Compiled by Yunnan Province Statistical Bureau. 2005, Beijing: China Statistics Press
Murray CJL, Lopez AD: The Global burden of Disease: A comprehensive assessment of mortality and disability from diseases, injuries, and risk factors in 1990 projected to 2020. 1996, Cambridge, Mass: Harvard School of Public Health
Murray CJL, Lopez AD: Alternative projections of mortality and disability by cause 1990–2020: Global Burden of Disease study. Lancet. 1997, 349 (9064): 1498-1504. 10.1016/S0140-6736(96)07492-2.
Murray CJL, Lopez AD: The global burden of disease. Cambridge (MA): Harvard school of Public Health on behalf of the World Health Organization and the World Bank; Global Burden of Disease and Injury Series. 1996, 1:
Murray CJL, Lopez A: The incremental effect of age-weighting on YLLs, YLDs, and DALYs: a response. Bulletin of the World Health Organization. 1996, 74 (4): 445-46.
Barendregt JJ, Bonneux L, Van der Mass PJ: DALYs: the age-weights on balance. Bulletin of the World Health Organization. 1996, 74 (4): 439-443.
Anand S, Hanson K: Disability adjusted life years: a critical review. Journal of Health Economics. 1997, 16: 685-702. 10.1016/S0167-6296(97)00005-2.
Coale A, Guo G: Revised regional model life tables at very low levels of mortality. Population Index. 1989, 55 (5): 613-643.
Chiang CL: The life table and its construction. From: Introduction to stochastic processes in Biostatistics. 1968, New York: John Willey & Sons, Chapter 9: 189-241.
Xia Y, Gong YL, Gu XY, Shu BG: The burden of disease study in Shanghai city. Chinese Journal of Health Statistics. 2000, 17 (4): 211-214.
Liu MN, Ji WY, Liu QM, Sui CL: Years of life lost due to premature death in Ha Erbing city. Chinese Journal of Health Statistics. 2002, 19 (4): 225-227.
Hyder AA, Morrow RH: Applying burden of disease methods in developing countries: A case study from Pakistan. American Journal of Public Health. 2000, 90 (8): 1235-1240.
Wurthwein R, Gbangou A, Sauerborn R, Schmidt CM: Measuring the local burden of disease. A study of years of life lost in sub-Saharan Africa. International Journal of Epidemiology. 2001, 30 (3): 501-508. 10.1093/ije/30.3.501.
Schopper D, Pereira J, Torres A, Cuende N, Alonso M, Baylin A, Ammon C, Rougemont A: Estimating the burden of disease in one Swiss canton: what do disability adjusted life years (DALY) tell us?. International Journal of Epidemiology. 2000, 29 (5): 871-877. 10.1093/ije/29.5.871.
Kominski GF, Simon PA, Ho A, Luck J, Lim YW, Fielding JE: Assessing the burden of disease and injury in Los Angeles county using disability-adjusted life years. Public Health Reports. 2002, 117 (2): 185-191.
Zhu WX: The one child family policy. Archives of Disease in Childhood. 2003, 88: 463-464. 10.1136/adc.88.6.463.
Brand I: Medical services in China. Australian and New Zealand Journal of Medicine. 2000, 30: 272-274.
Mason W: Recent trends in infant mortality in China. Mathematical Social Sciences. 1993, 25 (3): 310-312. 10.1016/0165-4896(93)90044-J.
Reed DM: The paradox of high risk of stroke in populations with low risk of coronary heart disease. American Journal of Epidemiology. 1990, 131 (4): 579-588.
Tunstall-Pedoe H, Kuulasmaa K, Amouyel P, Arveiler D, Rajakangas AM, Pajak A: Myocardial infarction and coronary deaths in the World Health Organization MONICA Project: registration procedures, event rates, and case-fatality rates in 38 populations from 21 countries in four continents. Circulation. 1994, 90: 583-612.
Zhang J, Kesteloot H: Differences in all-cause, cardiovascular and cancer mortality between HongKong and Singapore: Role of nutrition. European Journal of Epidemiology. 2001, 17 (5): 469-477. 10.1023/A:1013773517052.
Liau CS: Cardiovascular disease in later life in Taiwan. Geriatrics and Gerontology International. 2004, 4: S304-S306. 10.1111/j.1447-0594.2004.00236.x.
Yang G, Fan L, Tan J, Qi G, Zhang Y, Samet JM, Taylor CE, Becker K, Xu J: Smoking in China: Findings of the 1996 National Prevalence Survey. JAMA. 1999, 282 (3): 1247-1253. 10.1001/jama.282.13.1247.
Yunnan statistical yearbook: 2002: Compiled by Yunnan Province Statistical Bureau. 2003, Beijing: China Statistics Press
Ezzati M, Lopez AD: Regional, disease specific patterns of smoking-attributable mortality in 2000. Tobacco Control. 2004, 13: 388-395. 10.1136/tc.2003.005215.
Bruce N, Perez-Padilla R, Albalak R: The health effects of indoor air pollution exposure in developing countries. 2002, Geneva: World Health Organization
Desai MA, Mehta S, Smith KR: Indoor smoke from solid fuels: assessing the environmental burden of disease at national and local levels. WHO environmental burden of disease series. No. 4. 2004, Geneva: World Health Organization
Odero W, Garner P, Zwi A: Road traffic injuries in developing countries: a comprehensive review of epidemiological studies. Tropical Medicine & International Health. 1997, 2 (5): 445-460. 10.1046/j.1365-3156.1997.d01-296.x.
Suriyawongpaisal P, Kanchanasut S: Road-traffic injuries in Thailand: trends, selected underlying determinants and status of intervention. International Journal of Injury Control and Safety Promotion. 2003, 10: 95-104. 10.1076/icsp.10.1.95.14110.
Wang SY, Hu YL: Road injuries, epidemiological features and risk factors of traffic injury in China. Chinese Journal of Epidemiology. 1997, 18: 142-144.
Kpk LP, Wen-Shing T, (Eds): Suicidal Behavior in the Asia-Pacific Region. 1992, Singapore: Singapore University Press
Goldsmith SK, Pellmar TC, Kleinman AM, Bunney WE, (Eds): Reducing Suicide: A National Imperative. 2002, Washington, DC: National Academies Press
Rotheram-Borus MJ, Piacentini J, Cantwell C, Belin TR, Song J: The 18-month impact of an emergency room intervention for adolescent female suicide attempters. Journal of Consulting and Clinical Psychology. 2000, 68: 1081-93. 10.1037/0022-006X.68.6.1081.
Yang G, Hu J, Rao KQ, Ma J, Rao C, Lopez AD: Mortality registration and surveillance in China: History, current situation and challenges. Population Health Metrics. 2005, 3: 3-10.1186/1478-7954-3-3.
Murray CJL, Lopez AD: Global mortality, disability and contribution of risk factors: Global Burden of Disease study. Lancet. 1997, 349 (9063): 1436-1442. 10.1016/S0140-6736(96)07495-8.
This study is a part of the first author's thesis to fulfill the requirement for a PhD in epidemiology at the Prince of Songkla University. We thank Dr. Alan Geater and Mr Edward McNeil for reviewing the manuscript. They are staff members of Epidemiology Unit, Prince of Songkla University. This research was supported by Pan Long Center for Diseases Control, Guan Du Center for Diseases Control, Shi Lin Health Bureau and Shi Lin Maternal and Children Hospital; we thank all of these sectors for providing the mortality data for this study. The study was supported by a grant from Yunnan Provincial Natural Science Funds (Grant number: 2003C0022Q). Virasakdi Chongsuvivatwong receives his salary from Prince of Songkla University. The funding source has influence neither on the technical part of the research nor on submission of the manuscript.
Cai Le carried out the study and drafted the manuscript. Virasakdi Chongsuvivatwong conceptualized the research idea, participated in the design of the study, interpreted the results and helped to draft the manuscript. Both authors read and approved the final manuscript.