Next Article in Journal
Muscle Activation and Kinematic Analysis during the Inclined Leg Press Exercise in Young Females
Next Article in Special Issue
Assessment of Human Health Risks Posed by Nano-and Microplastics Is Currently Not Feasible
Previous Article in Journal
China’s Provincial Eco-Efficiency and Its Driving Factors—Based on Network DEA and PLS-SEM Method
Previous Article in Special Issue
Fundamental Concepts of Human Thermoregulation and Adaptation to Heat: A Review in the Context of Global Warming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Projection of the Number of Elderly in Different Health States in Thailand in the Next Ten Years, 2020–2030

by
Panupong Tantirat
1,*,
Repeepong Suphanchaimat
1,2,
Thanit Rattanathumsakul
1 and
Thinakorn Noree
2
1
Division of Epidemiology, Disease Control, Ministry of Public Health, Nonthaburi 11000, Thailand
2
International Health Policy Program (IHPP), Ministry of Public Health, Nonthaburi 11000, Thailand
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(22), 8703; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17228703
Submission received: 7 October 2020 / Revised: 12 November 2020 / Accepted: 13 November 2020 / Published: 23 November 2020
(This article belongs to the Special Issue New Challenges and Crucial Topics for 2030 Public Health)

Abstract

:
The objective of this study is to predict the volume of the elderly in different health status categories in Thailand in the next ten years (2020–2030). Multistate modelling was performed. We defined four states of elderly patients (aged ≥ 60 years) according to four different levels of Activities of Daily Living (ADL): social group; home group; bedridden group; and dead group. The volume of newcomers was projected by trend extrapolation methods with exponential growth. The transition probabilities from one state to another was obtained by literature review and model optimization. The mortality rate was obtained by literature review. Sensitivity analysis was conducted. By 2030, the number of social, home, and bedridden groups was 15,593,054, 321,511, and 152,749, respectively. The model prediction error was 1.75%. Sensitivity analysis with the change of transition probabilities by 20% caused the number of bedridden patients to vary from between 150,249 and 155,596. In conclusion, the number of bedridden elders will reach 153,000 in the next decade (3 times larger than the status quo). Policy makers may consider using this finding as an input for future resource planning and allocation. Further studies should be conducted to identify the parameters that better reflect the transition of people from one health state to another.

1. Introduction

The world is now experiencing many aged societies. The definitions of an aged society are multiple and can be categorized into many subgroups, namely, aged, complete-aged, and super-aged societies. An aged society refers to a state where people equal to or more than 60 years constitute over 10% of the total population. Complete-aged and super-aged societies contain a volume of elderly people aged over 60 by more than 20% and 28% of the total population, respectively [1].
Thailand is one among many countries facing aging society challenges. The country has been considered an aged society since 2005, and it may turn to a complete-aged and super-aged society by 2025 and 2030 respectively, if the growth rate of aging population continues at the current pace [2]. A state of dependency of the elderly is one of the critical challenges in the Thai health service system. This challenge is coupled with the high prevalence and incidence of chronic non-communicable diseases (NCDs) such as hypertension, diabetes, stroke, and coronary artery diseases [3,4]. These diseases are likely to undermine the well-being of the elderly and make the elderly prone to a state of to complete dependency, so-called, bedridden.
The Thai Ministry of Public health (MOPH) classified the elderly into three groups, namely, social group, home group, and bedridden group, based on the health status as measured by Activities of Daily Living (ADL). The ADL was an index comprising many health components, such as mobility status, ability to climb up and down the stairs, self-dressing and toileting. The index is scored between 0 and 20. The social group included the elderly with ADL of more than 11, reflecting a group of people who can take care of themselves. The home group included the elderly who can take care of themselves to a certain extent (ADL = 5–11). The bedridden group included the elderly who cannot perform basic living activities by themselves (ADL < 5) [5]. In 2015, MOPH reported that about 79% of the elderly in Thailand (5 million) were in the social group and the rest, 21%, were in home and bedridden groups combined [6]. The National Health Exam Survey also found that 0.2% of the elderly all over the country were in the bedridden group [7]. Projection of health status attracted much attention in many countries, such as studies in Australia (multiple state model) [8], UK (multiple state model) [9], European nations (macro-simulation model) [10], Taiwan (Grey forecasting model) [11], and Sweden (extrapolative model) [12]. In Thailand, there were few studies that attempted to project the volume of elderly in different health status. The latest projection that applied multiple state model in Thailand was in 2014—the situation during that time was far different from the current situation [13]. Thus, updating the current situation of the volume of elderly in different health stages is likely to be helpful for future health policy planning.
One of the main concerns of many health authorities in Thailand, especially the MOPH, is whether and to what extent the Thai healthcare system will have adequate capacity to address the health needs of the elderly in the future. This question cannot be answered unless the projected number of the elderly in the future is first identified. Therefore the objective of this study is to predict the volume of the elderly in different health status categories in Thailand in the next ten years (2020–2030).

2. Method

2.1. Study Design and Development of Conceptual Framework

This study applies a quantitative secondary data analysis. A multi-state model was applied. Two rounds of consultative stakeholder meetings were conducted on 25 February 2020 and 19 March 2020. The participants of the meeting comprised three policy makers, six academics, and two representatives of the health-systems related institutes. The opinions of the attendees were used as the key inputs to construct the model. An additional literature review was performed to identify key parameters and validate the model assumption [8,9,13,14,15,16].
Participants in the stakeholder meetings helped structure the outline of the model. The model focused on the population aged more than 60 years, and consisted of three parts: incoming population, transition state, and dead state. The incoming population represented the transition of population from people aged younger than 60 years to the cohort of 60 years of age and onwards. The transition state represented the probability of people remaining in the same health status or moving to another health status in each year. For example, the people in the social group could move to home, bedridden, or dead groups, or they could remain in the social group. The dead state was the absorbing state of all populations. The overview of the model framework is presented in Figure 1.
The model equation was based on liner difference equations method. The population in the social group at time t + 1 was calculated from the number in the social group population at time t added to people who transferred from the other groups (home or bedridden) to the social group, subtracting the outgoing population (such as those transferring from the social group to other groups included death). All model equations are presented in Table 1.

2.2. Parameter Management

Many parameters were used to estimate the outcome. Λ denoted the volume of population cohort in each year. Δ denoted the proportion of deaths in a specific cohort. Θ denoted the transferal probability between two elder groups, which could be broken down into six scenarios (ΘSH [from social group to home group], ΘSB [from social group to bedridden group], ΘHS [from home group to social group], ΘHB [from home group to bedridden group], ΘBS [from bedridden group to social group], and ΘBH [from bedridden group to home group]). The final outputs were presented in terms of S(t), H(t), B(t), and D(t) which refer to the number of cases in social, home, bedridden, and dead groups at time t, respectively. The time step for the analysis was one year. The analysis commenced when people reached 60 years of age.

2.2.1. Population Cohort in Each Time Frame (Λ)

The volume of population cohort moving from a particular time point to the following time point (that is, number of people started at time (t) and remained in that health status in time (t + 1) was estimated by trend extrapolation methods with exponential growth. Two parameters were used in this step. First, ΛD denoted the mortality rate of new people who died before transferring to the following year; and second, π denoted the prevalence of people in each health status at any time (t − 1). This meant for the first cohort included in the model, π referred to the percentage of people in a particular health status when they were 59 years of age. The mortality rate of people who died before becoming 60 years old was calculated from the following formula: πdie before become 60-year-old = (n. of people aged 60 years—n. of people in the same cohort at aged 59 years)/(n. of people in the same cohort at aged 59 years). We obtained the number of people in each specific age group from the Civil Registry website [17]; and we used use the mean mortality rate for each cohort year from the Civil Registry data between 1994 and 2019 [17]. The prevalence of the population in each health status (social, home, and bedridden) was based on the study by Charnduwit in 2017, which calculated the prevalence based on population aged between 60 and 64 years [14]. The equation used to estimate the number of population cohort was shown in Table 2.

2.2.2. Transition Probabilities (Θ)

Transition probability is the likelihood (percentage) that an elder in one particular group (social or home or bedridden group) would change to another in the following year. For example, elderlies in home group in year X faced a chance to change their status to bedridden group in year X + 1; the chance would be labeled as transition probability from home to bedridden group. The transition probabilities across health statuses was acquired by the literature review [9]. We adjusted these parameters against the actual data in 2014 [14,17]. The population in 2014 in social, home, bedridden, and dead groups was 8.8 million, 181,000, 107,000, and 284,000, respectively [14,17]. We used ‘Solve Nonlinear Equations Function’ (generalized reduced gradient [GRG] method) in Microsoft Excel 365 (Microsoft (Thailand) Ltd., Bangkok, Thailand), which intended to identify the most fitted parameters (producing the least errors) in the model. GRG method is one of the widely used optimization methods in many prior studies [18,19].

2.2.3. Starting Population (Pop(t = 0))

The volume of population at the start of the analysis (initial reservoirs) was explored. Pop(0) denoted the population aged ≥ 60 years) in 2019 (time 0). S(0), H(0), and B(0) denoted the volume of people aged ≥ 60 years in the social group, the home group and the bedridden group in 2019, respectively. D(0) denoted the number of deaths in 2019. Π denoted the prevalence in specific groups, which was obtained from the prior study by Charnduwit [14]. We assumed that at the beginning of the analysis, D(0) equated 0. The following formula was used to determine initial reservoirs, Table 3.

2.2.4. Mortality Rate (Δ)

A group-specific mortality rate was determined by the following formula, see Table 4. γ denoted a factor of each health specific group that was used to adjust the crude mortality rate. R denoted relative risk of mortality. The numerators were obtained from the Report of Statistics in the Public Health of Thailand, 2018 [20]. The denominators were retrieved from the Civil Registry. We obtained group-specific relative risk (R) from a Danish nationwide population-based cohort study by Ryg et al. [21]. We used published data (not raw data) of Ryg et al.’s study to calculate the relative risk to serve as an input for our model.

2.3. Model Prediction Error

The error of social, home, bedridden, and dead group prediction was calculated by mean absolute percentage error (MAPE). The overall model error was a mean average of the errors in all groups combined (Table 5).

2.4. Model Validation

We validated the model by calculating the numbers in social, home, bedridden, and dead groups in 2014 (based on the initial data in 2013) and adjusted this result with the actual data in 2014.

2.5. Projection of the Amount of Bedridden Patients

We used the optimized parameters to calculate the numbers of elderly in social, home, bedridden, and dead groups based on the equations shown in Table 5. We set the starting point (t = 1) at 2020 and the end point at 2030.

2.6. Sensitivity Analysis

We performed sensitivity analysis by repeating the calculation for 10,000 cycles. The key parameters, such as the prevalence of a specific group for calculated incoming population (π), crude mortality rate (Δ), relative risk of specific group mortality (R), and prevalence of specific group in the elderly (Π), were changed in each cycle, based on their distribution. The transit probability was the main parameter of interest. We presumed four different scenarios with different degree of changes in transit probability (changes by a maximum of 5%, 10%, 15%, and 20%, respectively). Then we assessed the change in the number of bedridden populations by 2030 in these four different scenarios from the base model.

2.7. Statistical Software

All analyses were performed by R 4.0 (The R Foundation for Statistical Computing, Vienna, Austria), RStudio® (© 2020 RStudio, Boston, MA, USA), and Microsoft Excel (Microsoft (Thailand) Ltd., Bangkok, Thailand).

3. Results

3.1. Parameter Identification

Summary of parameters used for projects are presented in Table 6 include population cohort in each time frame, transition probabilities, Starting population, and mortality rate.

3.1.1. Population Cohort in Each Time Frame (Λ)

The population cohort in each time frame was calculated using the exponential equation, presented as Λ = 302,880 × e0.0359(t), where t started from the calendar year of 1994. The coefficient of determination (R2) was 0.9442. Parameters that were used to indicate the number of incoming populations were πS = 0.9922, πH = 0.0029 (0.000031), πB = 0.0049 (0.000041) and πD = 0.0131 (0.0054).

3.1.2. Transition Probabilities (Θ)

Transition probabilities of ΘSH, ΘSB, ΘHS, ΘHB, ΘBS, and ΘBH equated 0.0169, 0.0071, 0.1257, 0.0782, 0.0470, and 0.0688, respectively. The value of these parameters after adjusting with the actual data was 0.0049, 0.0002, 0.1259, 0.0782, 0.0470, and 0.0709, in consecutive order. After adjustment, the model error reduced from 38.67% to 1.40%.

3.1.3. Starting Population

The overall population volume in 2019 was 11,136,059. The prevalence of the elderly in social, home, and bedridden groups was 0.9683, 0.0199 (0.000046), and 0.0118 (0.000036), respectively.

3.1.4. Mortality Rate (Δ)

The relative risk of the mortality rate of the home group in comparison with the social group was 1.45 (95% CI = 1.4529 − 1.4470), while the relative risk of mortality rate of the bedridden group when compared with the social group was 2.27 (95% CI = 2.2749–2.2656). The crude death rate for the overall population was 0.0307 (0.0011).

3.2. Model Validation

The prediction error in social, home, and bedridden groups was 0.30%, 0.003%, and 4.95%, respectively. The overall model error, given all group-specific errors combined, was 1.75%.

3.3. Projected Amount of Bedridden Patients

The predicted number of the social group population increased from 10.78 million in 2020 to 15.59 million in 2030. The number in the home group also increased (221,610 to 318,980). The number in the bedridden group showed a parabola pattern in the first 4 years followed by a linear increasing pattern from 131,410 to 153,640 thousand between 2020 and 2030. The figures of population estimation in all three groups are presented in Figure 2, Figure 3 and Figure 4.

3.4. Sensitivity Analysis

We found that the mean and median of five different scenarios remained unchanged (about 153,000) given different assumptions. The wider the degree of change in transition probabilities, the greater the variation in the number of bed-ridden patients. The volume of bed-ridden patients ranged from 139,110 to 171,616 given the assumption of 20% change in the transition probabilities, while this figure reduced to 142,672 to 161,817 in the base scenario (0% change in the transition probabilities), see Figure 5.

4. Discussion

Overall, this study found that the number of bedridden patients would reach 153,640 in 2030; about three-fold larger than the volume of bedridden patients at present. The vast volume of predicted numbers might be due to the large amount of elderly entering the cohort in each year and the large degree of variation of the input parameters.
Our study projected a smaller amount of bedridden patients than the suggested number presented by Chandeovwit et al. which indicated that the number of the patients would be about 240,000 [14]. This difference might be due to the difference in the model assumption and design. The study by Chandeovwit et al. used a constant proportion of dependency rate while our study applied a multistate method. Another study by Srithamrongsawat et al. also projected relatively similar amount of patients to our study; that is, the volume of bed ridden patients in 2024 would reach 110,000 (our study forecasted the figure of 133,000 at the same corresponding time) [13]. The prediction error of this study was small (error = 1.75%). The error appeared to be largest in bedridden group (4.95%).
According to the National Health Statistical Office’s recommendation, a full-time care-manager (CM) should take care of about 35–40 home and bedridden patients and a care-giver (CG) should take care of about 4–8 home and bedridden people [6] Based on these figures, it means that there should be about 10,633–15,039 CM and 53,170–131,597 CG in the health care system to meet the standards set by the NHSO. However, the registered CM and CG in 2020, as reported by Department of Health, numbered only 13,175 and 83,751, respectively [22].
In term of methodology, this study faced some limitations. We used multistate modelling because it better reflects the real situation than the static model [15]. However, the model requires a number of parameters to serve as the model input. The most important input was transit probability. The transit probability in this study derived from a parameter estimation rather than empirical study. Moreover, the empirical study needed to estimate transit probability originating from longitudinal research to better reflect the actual transition of a person from one health state to another [9,15]. However, the longitudinal study of health status of the elderly in Thailand was lacking. Another point of concern is that certain demographic variables, such as economic status and education background, have not been considered in this model [13,14]. Some parameters included in this model were quite outdated but they were the most recent data we could access. In addition, most literature used in this study was from foreign studies which might not reflect the actual healthcare system and population demography in Thailand. This phenomenon reflects gaps in knowledge in the field of mathematical modelling for public health use in Thailand.
Policy makers or relevant academics may use the results in this study as inputs for policy decision-making, particularly in terms of resource planning and mobilization. We recommend that further research should be conducted on specific types of bedridden elderly such as groups of patients with diverse socioeconomic backgrounds. Moreover, there should be additional empirical research on the transition of health status in the population. This will enable researchers and mathematical modelers in the field to access data that better reflect the real situation on the ground, and at the same time benefit the accuracy of model prediction in the future.

5. Conclusions

We estimated that, by 2030, the number of bedridden elderly would be about 153,000. This figure is about three times larger than the status quo. Policy makers and academics may consider using this finding as an important input for future resource planning, in light of the coming aged society. Further empirical studies should be conducted to identify the parameters that better reflect the transition of people from one health state to another. Additional research that explores subgroups in populations with different economic backgrounds is recommended.

Author Contributions

Conceptualization, P.T. and T.N.; Methodology, P.T., T.R. and R.S.; Software, P.T.; Validation, P.T. and R.S.; Formal Analysis, P.T.; Investigation, P.T.; Resources, P.T. and T.R.; Data Curation, P.T. and T.N.; Writing—Original Draft Preparation, P.T. and R.S.; Writing—Review and Editing, P.T. and R.S.; Visualization, P.T. and R.S.; Supervision, P.T. and T.N.; Project Administration, T.N.; Funding Acquisition, T.N. and R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the International Health Policy Program (IHPP), Ministry of Public Health, Thailand. IHPP received funding support from the Thailand Science Research and Innovation (TSRI) under the Senior Research Scholar on Health Policy and System Research [Contract No. RTA6280007].

Acknowledgments

The authors thank Field Epidemiology Training Program (FETP), Division of Epidemiology, Department of Disease control, Ministry of Public Health, Thailand in providing academic support.

Conflicts of Interest

The authors declare no conflict of interest.

Ethical Consideration

As this study is part of the long-term care project of the International Health Policy Program, MOPH, which aimed to generate evidence for policy decision making of the MOPH; it was not required to obtain ethics approval from the Institute for the Development of Human Research Protections. However, the researchers strictly followed ethical standards in research, that is, all individual information was strictly kept confidential and was not reported in the paper.

References

  1. Foundation of Thai Gerontology Research and Development Institute. Situation of the Thai Elderly 2018; AP: Bangkok, Thailand, 2018.
  2. UN. Department of Economic and Social Affairs Population Dynamics. Available online: https://population.un.org/wpp/ (accessed on 28 January 2019).
  3. Divison of Epidemiology. Annual Epidemiology Surveillance Report 2017; MOPH: Nonthaburi, Thailand, 2017.
  4. Divison of Epidemiology. Annual Epidemiological Surveillance Report 2018; MOPH: Nonthaburi, Thailand, 2018.
  5. Department of Health, Ministry of Public Health. Barthel Activities of Daily Living: ADL. Available online: http://hpc9.anamai.moph.go.th/ewt_dl_link.php?nid=508 (accessed on 28 January 2019).
  6. National Health Security Office (NHSO). Guideline of Long-Term Care Services Management for the Dependent Elderly in the National Health Insurance System; NHSO: Bangkok, Thailand, 2016.
  7. Ekkaphakorn, W. Report of the 3rd Public Health Survey by Physical Examination 2003, 3rd ed.; Prapakkham, Y., Bunyaratpan, P., Eds.; MOPH: Nonthaburi, Thailand, 2006.
  8. Leung, E. Projecting the Needs and Costs of Long Term Care in Australia; The University of Melbourne Victoria: Melbourne, Australia, 2003; pp. 1–31. [Google Scholar]
  9. Rickayzen, B.D.; Walsh, D.E.P. A Multi-State Model of Disability for the United Kingdom: Implications for Future Need for Long-Term Care for the Elderly. Br. Actuar. J. 2002, 1, 341–393. [Google Scholar] [CrossRef]
  10. Costa-Font, J.; Wittenberg, R.; Patxot, C.; Comas-Herrera, A.; Gori, C.; di Maio, A.; Pickard, L.; Pozzi, A.; Rothgang, H. Projecting Long-Term Care Expenditure in Four European Union Member States: The Influence of Demographic Scenarios. Soc. Indic. Res. 2008, 86, 303–321. [Google Scholar]
  11. Hsu, K.-T.; Yan, T.-M. Applying Grey Forecasting Models to Disability Rate and Demand for Long-Term Care Service of the Aged in Taiwan. In Proceedings of the 2008 IEEE International Conference on Systems, Man and Cybernetics, Singapore, 12–15 October 2008; pp. 2860–2865. [Google Scholar]
  12. Lagergren, M. Whither Care of Older Persons in Sweden?—A Prospective Analysis Based upon Simulation Model Calculations, 2000–2030. Health Policy 2005, 74, 325–334. [Google Scholar] [CrossRef] [PubMed]
  13. Srithamrongsawat, S.; Bundhamcharoen, K.; Sasat, S.; Odton, P.; Ratkjaroenkhajorn, S. Projection of Demand and Expenditure for Institutional Long Term Care in Thailand; Health Insurance System Research Office: Bangkok, Thailand, 2014. [Google Scholar]
  14. Charnduwit, W. Long-Term Care Insurance System: A System Suitable for Thailand; TDRI: Bangkok, Thailand, 2017. [Google Scholar]
  15. Meira-Machado, L.; de Uña-Álvarez, J.; Cadarso-Suárez, C.; Andersen, P.K. Multi-state models for the analysis of time-to-event data. Stat. Methods Med. Res. 2009, 18, 195–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. World Health Organisation (WHO). Current and Future Long-Term Care Needs. 2002. Available online: http://www.who.int/chp/knowledge/publications/ltc_needs.pdf (accessed on 1 April 2020).
  17. Bureau of Registration Administration. Offical Statistic Registration System. Available online: https://stat.bora.dopa.go.th/stat/statnew/statTDD/ (accessed on 1 April 2020).
  18. Abdel-Aleem, A.; El-Sharief, M.A.; Hassan, M.A.; El-Sebaie, M.G. Optimization of Reliability Based Model for Production Inventory System. Int. J. Manag. Sci. Eng. Manag. 2018, 13, 54–64. [Google Scholar] [CrossRef]
  19. Hashemi, S.H.; Mousavi Dehghani, S.A.; Samimi, S.E.; Dinmohammad, M.; Hashemi, S.A. Performance Comparison of GRG Algorithm with Evolutionary Algorithms in an Aqueous Electrolyte System. Model. Earth Syst. Environ. 2020, 6, 2103–2110. [Google Scholar] [CrossRef]
  20. Strategy and Planning Division. Public Health Statistics A.D.2018; MOPH: Nonthaburi, Thailand, 2018.
  21. Ryg, J.; Engberg, H.; Mariadas, P.; Pedersen, S.G.; Jorgensen, M.G.; Vinding, K.L.; Andersen-Ranberg, K. Barthel index at hospital admission is associated with mortality in geriatric patients: A danish nationwide population-based cohort study. Clin. Epidemiol. 2018, 10, 1789–1800. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Department of Health, Ministry of Public Health. Long Term Care. Available online: http://ltc.anamai.moph.go.th/ (accessed on 28 June 2020).
Figure 1. Framework of the model used in this study.
Figure 1. Framework of the model used in this study.
Ijerph 17 08703 g001
Figure 2. Number of the elderly in social group in Thailand, 2020–2030.
Figure 2. Number of the elderly in social group in Thailand, 2020–2030.
Ijerph 17 08703 g002
Figure 3. Number of the elderly in home group in Thailand, 2020–2030.
Figure 3. Number of the elderly in home group in Thailand, 2020–2030.
Ijerph 17 08703 g003
Figure 4. Number of the elderly in bedridden group in Thailand, 2020–2030.
Figure 4. Number of the elderly in bedridden group in Thailand, 2020–2030.
Ijerph 17 08703 g004
Figure 5. Number of bedridden patients given different degrees of uncertainty of the transition probabilities (scenarios with 5%, 10%, 15%, and 20% changes).
Figure 5. Number of bedridden patients given different degrees of uncertainty of the transition probabilities (scenarios with 5%, 10%, 15%, and 20% changes).
Ijerph 17 08703 g005
Table 1. Summary of the liner difference equations used to estimate volume of patients in social, home, bedridden, and dead groups.
Table 1. Summary of the liner difference equations used to estimate volume of patients in social, home, bedridden, and dead groups.
Model EquationFormula
1S(t + 1) = S(t) + ΛS(t) + H(t) × ΘHS + B(t) × ΘBS − S(t) × ΘSH − S(t) × ΘSB − S(t) × ΔS
2H(t + 1) = H(t) + ΛH(t) + S(t) × ΘSH + B(t) × ΘBH − H(t) × ΘHS − H(t) × ΘHB − H(t) × ΔH
3B(t + 1) = B(t) + ΛB(t) + H(t) × ΘHB + S(t) × ΘSB − B(t) × ΘBS − B(t) × ΘBH − BG(t) × ΔB
4D(t + 1) = D(t) + S(t) × ΔS + H(t) × ΔH + B(t) × ΔB
Table 2. Summary of the equations used to estimate the number of population cohort at time t.
Table 2. Summary of the equations used to estimate the number of population cohort at time t.
EquationFormulaDescription
1ΛD = πdie before becoming 60-year-old × ΛNumber of deaths at time t − 1
2Λ′ = Λ − ΛDNumber of the elders entering time t
3ΛS = πS × Λ′Number of social group entering time t
4ΛH = πH × Λ′Number of home group entering time t
5ΛB = πB × Λ′Number of bedridden group entering time t
Table 3. Summary of the equations used to estimate the volume of population at the start of the analysis (initial reservoirs).
Table 3. Summary of the equations used to estimate the volume of population at the start of the analysis (initial reservoirs).
EquationFormulaDescription
1Pop(0) = S(0) + H(0) + B(0) + D(0)Total volume of population aged ≥60 years consisted of people aged ≥60 years in the social group, the home group, the bedridden group, and the death group at time 0.
2S(0) = Pop(0) × ΠSVolume of people aged ≥60 years in the social group resulted from population aged ≥60 years multiplied by prevalence of social group.
3H(0) = Pop(0) × ΠHVolume of people aged ≥60 years in the home group resulted from population aged ≥60 years multiplied by prevalence of home group.
4B(0) = Pop(0) × ΠBVolume of people aged ≥60 years in the bedridden group resulted from population aged ≥60 years multiplied by prevalence of bedridden group.
5D(0) = 0Volume of dead people in the model. We assumed there was no death at the beginning of the analysis.
Table 4. Summary of the equations used to estimate group-specific mortality rate.
Table 4. Summary of the equations used to estimate group-specific mortality rate.
EquationFormulaDescription
1 Δ = Π s Δ S +   Π H Δ H +   Π B Δ B Π S + Π H + Π B Crude mortality rate was a prevalence weight average of group-specific mortality.
2ΔS = γ s × ΔSocial group mortality rate was social group specific severity factors multiply by crude mortality.
3ΔH = γ H × ΔHome group mortality rate was home group specific severity factors multiply by crude mortality.
4ΔB = γ B × ΔBedridden group mortality rate was bedridden group specific severity factors multiply by crude mortality.
5 R H = Δ H Δ S = Δ × γ H Δ × γ s Relative risk of home mortality was calculated from home group mortality rate over social group mortality rate.
6 R B = Δ B Δ S = Δ × γ B Δ × γ s Relative risk of bedridden mortality was calculated from bedridden group mortality rate over social group mortality rate.
7 γ s = 1 Π S + Π H R H + Π B R B Social group mortality was calculated from group specific prevalence and relative risk. This equation was rewritten form of equations 1–6.
Table 5. Summary of the equations used to calculate model error.
Table 5. Summary of the equations used to calculate model error.
EquationFormulaDescription
1 E r r o r j = 1 N i = 1 N O b s e r v e r d i j P r e d i c t e d i j O b s e r v e r d i j Mean absolute percentage error (MAPE) of group j was the summation of absolute difference divided by observed value.
2 M o d e l   e r r o r = 1 K j = 1 K E r r o r j Model error was mean average of all group-specific errors combined.
Table 6. Summary of parameters used for projecting the amount of bedridden patients.
Table 6. Summary of parameters used for projecting the amount of bedridden patients.
NoGroup of VariablesParametersMeanSDReference (Ref)
1New populationΛΛ = 302,880 × e0.0359(t)Bureau of Registration administration [17]
2Prevalence of specific group who age equal 59-year-oldπS0.9922-Model calibration from Bureau of Registration administration [17]
3πH0.00290.000031Charnduwit [14]
4πB0.00490.000041
5πD0.01310.0054Model calibration from Bureau of Registration administration [17]
6Mortality rateΔ0.03070.0011Bureau of Registration administration and Strategy and Planning Division [17,20]
7Relative mortality rateRH1.450.0010 (SE of ln RR)Ryg [21]
8RB2.270.0010 (SE of ln RR)
9Prevalence of specific group in elderlyΠS0.9683-Model calibration from Charnduwit [14]
10ΠH0.01990.000046Charnduwit [14]
11ΠB0.01180.000036
12Mortality rate in specific groupΔS0.0494-Model calibration from Bureau of Registration administration, Strategy and Planning Division, and Ryg [17,20,21]
13ΔH0.1465-
14ΔB0.2050-
15Transit probability from social groupΘSH0.0169-Model calibration from Rickayzen [9]
16ΘSB0.0071-
17Transit probability from home groupΘHS0.1257-
18ΘHB0.0782-
19Transit probability from bedridden groupΘBS0.0470-
20ΘBH0.0688-
21Initial Total populationPop(t = 0)11,136,059-Bureau of Registration administration [17]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tantirat, P.; Suphanchaimat, R.; Rattanathumsakul, T.; Noree, T. Projection of the Number of Elderly in Different Health States in Thailand in the Next Ten Years, 2020–2030. Int. J. Environ. Res. Public Health 2020, 17, 8703. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17228703

AMA Style

Tantirat P, Suphanchaimat R, Rattanathumsakul T, Noree T. Projection of the Number of Elderly in Different Health States in Thailand in the Next Ten Years, 2020–2030. International Journal of Environmental Research and Public Health. 2020; 17(22):8703. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17228703

Chicago/Turabian Style

Tantirat, Panupong, Repeepong Suphanchaimat, Thanit Rattanathumsakul, and Thinakorn Noree. 2020. "Projection of the Number of Elderly in Different Health States in Thailand in the Next Ten Years, 2020–2030" International Journal of Environmental Research and Public Health 17, no. 22: 8703. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17228703

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop