- Research article
- Open Access
- Open Peer Review
Municipal health expectancy in Japan: decreased healthy longevity of older people in socioeconomically disadvantaged areas
BMC Public Healthvolume 5, Article number: 65 (2005)
Little is known about small-area variation in healthy longevity of older people and its socioeconomic correlates. This study aimed to estimate health expectancy at 65 years (HE65) at the municipal level in Japan, and to examine its relation to area socio-demographic conditions.
HE65 of municipalities (N = 3361) across Japan was estimated by a linear regression formula with life expectancy at 65 years and the prevalence of those certificated as needing nursing care. The relation between HE65 and area socio-demographic indicators was examined using correlation coefficients.
The estimated HE65 (years) ranged from 13.13 to 17.39 for men and from 14.84 to 20.53 for women. HE65 was significantly positively correlated with the proportion of elderly and per capita income, and negatively correlated with the percentage of households of a single elderly person, divorce rate, and unemployment rate. These relations were stronger in large municipalities (with a population of more than 100,000) than in small and medium-size municipalities.
A decrease in healthy longevity of older people was associated with a higher percentage of households of a single elderly person and divorce rate, and lower socioeconomic conditions. This study suggests that older people in urban areas are susceptible to socio-demographic factors, and a social support network for older people living in socioeconomically disadvantaged conditions should be encouraged.
The ageing of the population has been progressing throughout the world, and Japan is the leading country with population ageing: the proportion of elderly (aged 65 years or over) was 17.4% in 2000 and is estimated to be over 35% in 2050 . The quality of the years lived and well-being in an ageing society with declining mortality and increased life expectancy (LE) have been addressed .
To measure health levels including the aspects of quality of life and well-being, several comprehensive health indicators have been developed, and among them health expectancy (HE) has been commonly used [2–5]. HE is calculated by application of prevalence of morbidity data such as disability prevalence on age-specific person years in a life table [2, 6]. The national health plan in Japan, "Health Japan 21", accompanied by local action plans, aims to prolong HE at the national and local levels through disease prevention and health promotion . However, it is difficult to calculate and monitor HE for health policy planning and evaluation, especially at the local level, because of a lack of data on morbidity. Although previous studies have calculated HE at the prefectural level in Japan, these calculations used inconsistent and various databases of disability prevalence, and could not be applied at the municipal level because of limited availability of data [8–10].
In Japan, long-term care insurance (LTCI) was introduced in 2000 for older people needing nursing care, and citizens aged 40 years or over can receive insurance benefits after application and certification that they need nursing care. Since LTCI unions are managed by municipalities, it is theoretically possible to calculate HE using the municipal data of LTCI as morbidity. However, municipal HE has been calculated in only one prefecture , because the sex- and age-specific prevalence of those certificated as needing nursing care for LTCI by municipalities is not routinely published.
It was confirmed by previous studies that socio-demographic factors influence mortality, morbidity, and other health conditions of older people [12–15]. Socioeconomic status such as education level, occupational social class, and income has been demonstrated to be a critical factor predicting the health status of older people [13, 15]. Family structure and living arrangement have been elicited as other critical factors influencing the health of older people, and conditions such as living alone and being widowed/divorced are related to poorer health status [16, 17].
There was little relation between HE and socio-demographic indicators related to income and family structure at the prefectural level in Japan , while a municipal-level study in one prefecture showed a weak but significant relation between HE and socio-demographic indicators such as the percentage of nuclear households . The relation between health status and regional characteristics might have been obscured due the higher level analysis which included heterogeneous smaller level regions, as a previous study demonstrated that the relation differed according to the area classification by urban-rural difference . Analysis at different levels could lead to an inconsistent relation, and the smaller-area analysis more directly reflects the relation at the individual and household levels [20, 21]. Therefore, studies examining the relation between socio-demographic factors and health including not only death but also morbidity at a smaller level will provide reliable evidence of the influence of socio-demographic factors on health in older people.
The purposes of this study were to estimate HE of older people at the municipal level across Japan, and to elucidate the association between HE and area socio-demographic factors.
Japan consists mainly of four islands (a map is shown in Figures 1 and 2). Hokkaido Island is located in the northern part, and includes only one prefecture (Hokkaido prefecture). The mainland (Honshu) includes several metropolitan areas such as Tokyo, Nagoya, and Osaka. The other islands are Shikoku and Kyushu.
According to Local Autonomy Law, local public entities in Japan are divided into two categories. The first category consists of cities, towns, and villages. The second category consists of prefectures (N = 47). All districts in the country belong to one of the municipalities and at the same time fall within the boundaries of one of the prefectures. As specific cases, the Tokyo prefecture (Tokyo Metropolis) includes 23 special wards ("ku") in addition to cities, towns and villages, and 12 large cities ("cities designated by ordinance") such as Nagoya and Osaka consist of wards ("ku"). In total, there were 3361 units of cities, towns and villages as well as Tokyo special wards and wards of cities designated by ordinance in December 2001.
Estimation of health expectancy
Calculation of HE by the Sullivan method requires life-table data and the sex- and age-specific morbidity . However, municipal data of the sex- and age-specific number of those certificated as needing nursing care for LTCI are not available, and only the total number was available by insurance union. Therefore, this study estimated HE at 65 years (HE65) by a linear regression formula using life expectancy at 65 years (LE65) and the certificated rate in LTCI (P): HE65 = a + b × LE65 + c × P.
The formula was drawn from prefectural data (N = 47) of precise HE65, which was calculated by the Sullivan method with the prefectural life-table in 2000  and the sex- and age-specific certification rate of LTCI in 2002 . The coefficient of determinant (R2) for the regression with LE65 and the crude certification rate was 0.913 for men and 0.906 for women. When the standardized certification ratio (SCR), which was age-standardised by the indirect method, was used as a predictor, R2 was increased to 0.986 for men and 0.982 for women.
The validity of this estimation was confirmed by the municipal data (N = 59) of one prefecture (Shimane prefecture), in which the precise HE65 had been calculated by the Sullivan method . The correlation coefficient between the precise HE65 and estimated HE65 predicted by the formula calculated from prefectural data was 0.976 for men and 0.956 for women.
Municipal LE65 was drawn from the municipal life-table in 2000, using the number of deaths during 1999–2001 and the 2000 census population, and the empirical Bayesian method was applied . The number of those certificated as needing nursing care for LTCI by insurance union in May, 2002, was obtained from the database published by the Welfare and Medical Statistical Agency . SCR was calculated using the municipal population and the age-specific certification rates of the national data from the All-Japan Federation of National Health Insurance Organizations, in which the age categories were 65–69, 70–74, 75–79, and 80+ years . Although most municipalities (N = 2823, 84.0%) had their own insurance unions, there were 61 insurance alliances managed with neighbouring municipalities. In these cases, SCR of the alliance was calculated and used as that of the municipalities for HE estimation.
There was substantial variation in the population size among municipalities, ranging from only a few hundred to one million, and the certification rate in municipalities with a small population size showed extreme fluctuation. To correct the fluctuation of SCR, Bayesian SCR was estimated using hierarchical Poisson regression and Markov chain Monte Carlo (MCMC) method [26–31]. The levels of hierarchy were municipalities for the lower level and secondary medical care zones (SMCZs) for the higher level. The details of calculation of Bayesian SCR are described in the additional file [see Additional file 1]. The range of SCR was decreased from crude SCR of 24.8 to 178.0 (S.D. = 18.3) to Bayesian SCR of 46.8 to 161.3 (S.D. = 16.5).
The regression formula obtained from prefectural data and used for estimation of municipal HE65 was HE65 = 4.772 + 0.706 × LE65 – 1.776 × SCR for men and HE65 = 12.024 + 0.453 × LE65 – 4.576 × SCR for women. Statistical analysis was conducted using SPSS 11.0 for linear regression analysis and MLwiN 1.10 for Bayesian hierarchical Poisson regression. Geographic mapping of HE was conducted using ArcGIS 8.3, with the Japanese map (geographic coordinate system: GRS 1980) obtained from ESRI Japan http://www.esrij.com.
Analysis of relation to socio-demographic indicators
The socio-demographic indicators used in this study were population, population density, proportion of elderly (aged 65 years or over), percentage of households of a single elderly person (aged 65 years or over), percentage of nuclear households, divorce rate, per capita income, and unemployment rate. The data were drawn from the database based on the population census and other national surveys of 1999 and 2000 . Since crude divorce rate showed a large difference in variance within the three categories of municipalities according to the population size mentioned later, Bayesian method with hierarchical Poisson regression was applied to predict the value adjusted for the statistical fluctuation due to small population size, as well as SCR. A summary of the socio-demographic indicators is shown in Table 1.
The relation of HE65 and LE65 to socio-demographic indicators was examined using correlation coefficient (Spearman's) for total municipalities. Then, we examined the correlation according to the population size of municipalities, in which municipalities were divided into three categories according to population size: less than 10,000 (small, N = 1554), 10,000 to 100,000 (medium-size, N = 1469), and more than 100,000 (large, N = 338). Statistical analysis was conducted using SPSS 11.0.
Estimated health expectancy
Table 2 shows LE65 and estimated HE65. HE65 in men (years) ranged from 13.13 to 17.39 with a mean (S.D.) of 15.53 (0.52), and HE65 in women (years) ranged from 14.84 to 20.53 with a mean (S.D.) of 18.01 (0.77). For both men and women, large municipalities showed significantly (p < 0.001) lower HE65 than small and medium-size municipalities. HE65 for the overall Japanese population, which was estimated using the regression formula, was 15.39 for men and 17.62 for women, while HE65 estimated by the Sullivan method with sex- and age-specific long-term care prevalence was 15.40 for men and 17.62 for women, and LE65 in 2000 was 17.56 for men and 22.46 for women .
The geographic mapping of HE65 in men is shown in Figure 1. The central part of the mainland showed accumulation of municipalities with longer HE65, while the northern part of the mainland showed accumulation of municipalities with shorter HE65. Metropolitan areas, especially Osaka, showed relatively shorter HE65. Other regions such as Hokkaido, Kyushu and Shikoku showed a miscellaneous pattern of municipalities with shorter and longer HE65. The geographical pattern in women, shown in Figure 2, was similar to that in men, and the shorter HE65 in the metropolitan areas was more remarkable.
The relation to socio-demographic indicators
Table 3 shows the correlation coefficients of LE65 and HE65 with socio-demographic indicators. As the number of samples was extremely large, most of the correlation coefficients showed statistical significance (p < 0.001). Compared to LE65, HE65 showed a stronger correlation with the percentage of households of a single elderly person and per capita income for men, and with the percentage of nuclear households, divorce rate and unemployment for men and women, while a weaker correlation with the proportion of elderly for men and women. The correlation with the percentage of households of a single elderly person and per capita income for women was reversed between LE65 and HE65.
Concerning the correlation coefficients according to population size of municipalities, the socio-demographic indicators showed a stronger correlation in large municipalities (with a population of more than 100,000) than in small and medium-size municipalities, except for population, population density, and the percentage of nuclear households. Especially, the percentage of households of a single elderly person, divorce rate, per capita income and unemployment rate showed correlation coefficients of more than 0.4 (or less than -0.4) in large municipalities.
This study estimated HE of older people and demonstrated its relation to socio-demographic indicators at the municipal level across Japan. As a result, the estimated HE65 (years) ranged from 13.13 to 17.39 for men and from 14.84 to 20.53 for women, and it was closely related to area socio-demographic conditions such as the percentage of households of a single elderly person, divorce rate, per capita income and unemployment rate. The relation of HE65 to socio-demographic indicators was stronger than that of LE65, and municipalities with a larger population showed a stronger relation between HE65 and socio-demographic factors.
Lower income and education are strongly related to higher morbidity of older people and shorter LE and HE in other countries [13, 15, 33]. In a previous study in Japan, however, area indicators related to income, education and unemployment showed little association with mortality of older people . This finding was consistent with the result of this study that LE65 showed a significant but weak correlation with socio-demographic indicators. HE65 was more strongly correlated with socio-demographic indicators than LE65. These findings indicated that the health status measured by not only mortality but also morbidity is more sensitive for area socioeconomic conditions than that measured by only mortality, and the health level of older people would be substantially decreased by socioeconomic conditions in disadvantaged areas.
A notable finding of this study was the negative relation between HE65 and the percentage of households of a single elderly person and divorce rate. Shorter life expectancy might cause an increase in older people living alone because of an increase in widows/widowers. However, considering that the proportion of a single elderly person did not showed a negative relation to LE65, it is likely that the shorter HE was caused by an increase of older people living alone. It was confirmed by previous studies that living alone including divorce and loss of a spouse has a negative impact on the health of older people . In addition, since most divorced persons were in the young- and middle-aged population, and divorces among older people are more unusual than among the middle-aged population , it seems plausible that the increased divorce rate influence the health status of older people through decreased support in families and communities.
In this study, we found a difference in the relation between HE65 and socio-demographic factors according to the size of municipalities. A previous study, also, showed that urban areas showed a stronger relation between health status and area socioeconomic conditions than rural areas . The first possible reason for this urban-rural difference is the choice of indicator/index for area socioeconomic conditions. A previous study demonstrated that the traditional deprivation index (Townsend index) was related to health only in urban areas, while an index including more multiple domains was related to health in both urban and rural areas . The second possible explanation is the existence of other factors mediating the relations, such as the social support network. The social support network has been demonstrated to influence the health status, and less social support is closely linked to poorer health in older people [12, 36]. Also, social factors are less associated with mortality in a socially cohesive area , and the relation between socioeconomic inequality and mortality is mediated by so-called social capital . Since urban areas show less social support network than rural areas in Japan [39, 40], it is likely that in urban areas, less social support network causes a decline in health status and enhances the influence of other social factors on health status in the elderly population.
A few weaknesses of this study should be mentioned, especially concerning the estimation of disability prevalence and HE. First, the disability prevalence used for HE estimation was based on LTCI data. As a previous study demonstrated, utilization of LTCI services was dependent on not only the level of disability but also socioeconomic conditions such as household income and family structure . However, LTCI in Japan provides consistent insurance across the nation and universally covers all elderly persons aged over 65 years , and thus, the influence of socioeconomic conditions on service utilization and regional variation would seem to be relatively small. In addition, LTCI data are the sole source for routine estimation of HE at present. Second, the disability prevalence (SCR) was estimated by hierarchal regression and MCMC procedure (fully Bayes method) [27–31]. Bayes estimation is commonly used for data smoothing in small area analysis and disease mapping, while alternative methods including empirical Bayes estimation can be applied [29–31, 43–46]. Although SMCZ is generally used as a higher-level regional unit for Bayes estimation of municipal mortality data in Japan [24, 26, 47, 48], an alternative hierarchical setting of regional levels would result in different figures of SCR, and consequently of HE. Finally, HE65 was estimated by a simple linear regression formula based on prefectural data. The validity of the estimated HE was confirmed by the municipal data of one prefecture, and the correlation coefficients between the precise and estimated HE were large enough to use the estimated HE for geographical comparison, instead of the precise HE. However, the application of this HE estimation for individual municipalities should be thought out. To expand the use of HE, the routine publication of detailed data of LTCI including the sex- and age-specific number of those certificated as needing nursing care is expected, which would enable calculation of the precise HE.
The results of this study have several implications in prolonging healthy longevity and supporting older people. The strong association between HE65 and socio-demographic factors indicates the need for efforts in the social context. Especially, there should be a focus on older people living alone. Also, the strong relation between HE65 and socio-demographic factors in large municipalities indicates that older people living in urban areas lack a social support network and are vulnerable to socioeconomically disadvantaged conditions. The social support network including an informal network and community activities will contribute to improving the health status of older people and prolongation of healthy longevity, especially in urban areas.
Although we addressed HE in this study, the results suggest that disability prevalence (SCR) is a useful indicator independently of mortality data including LE. Since HE was estimated using LE and SCR, the difference in the relation with socio-demographic indicators between LE and HE implies that SCR includes quite different dimensions from LE. The use of long-term care prevalence, in combination with mortality indicators (e.g., LE) and comprehensive indicators (e.g., HE), will contribute to concrete health policy making and evaluation for improvement of health, well-being, and quality of life in elderly population.
This study established a method to estimate HE available at the local level using a linear regression formula with life expectancy and long-term care prevalence. Municipal HE65 (years) ranged from 13.13 to 17.39 for men and from 14.84 to 20.53 for women. As a result of correlation analysis, a decreased HE was associated with higher percentage of households of a single elderly person and divorce rate, and lower socioeconomic conditions. This study suggests that older people in urban areas are susceptible to socio-demographic factors, and a social support network for older people living in socioeconomically disadvantaged conditions should be encouraged.
Health and Welfare Statistics, Association: Kokumin Eiseino Doko. 2003, Tokyo: Health and Welfare Statistics Association
Robine JM, Romieu I, Cambois E: Health life expectancy indicators. Bull World Health Organ. 1999, 77: 181-185.
Mathers CD, Sadana R, Salomon JA, Murray CJL, Lopes AD: Healthy life expectancy in 191 countries, 1999. Lancet. 2001, 357: 1685-1691. 10.1016/S0140-6736(00)04824-8.
Bone MR, Bebbington AC, Nicholaas G: Policy application of health life expectancy. J Aging Health. 1998, 10: 136-153.
Pamuk ER, Wagener DK, Molla MT: Achieving national health objectives: the impact on life expectancy and on healthy life expectancy. Am J Public Health. 2004, 94: 378-383.
Sullivan DF: A single index of mortality and morbidity. HSMHA Health Report. 1971, 86: 347-354.
Kawahara K: Health Promotion Law and Community Health Law. 2003, Tokyo: Sunlife Kikaku
Miyashita M, Hashimoto S, Ojima T, Nakamura Y, Hayashi M, Kato M, Fukudome K: Prevalence of the elderly requiring nursing care and health expectancy. Koseino Shihyo. 1999, 46: 25-29.
Segami K: Indicators toward the 21st century. 2000, Tokyo: Health and Welfare Statistics Association
Hasegawa T: Estimation of indicators of premature death and health expectancy for evaluation of Health Japan 21. 2001, Tokyo: Ministry of Health, Labour and Welfare
Shimane Prefectural Institute of Public Health and Environmental Science: Geographical variation of healthy life expectancy in Shimane Prefecture. 2002, Matsue: Shimane Prefectural Institute of Public Health and Environmental Science
Stuck AE, Walthert JM, Nikolaus T, Bula CJ, Hohmann C, Beck JC: Risk factors for functional status decline in community-living elderly people: a systematic literature review. Soc Sci Med. 1999, 48: 445-469. 10.1016/S0277-9536(98)00370-0.
Huisman M, Kunst AE, Mackenbach JP: Socioeconomic inequalities in morbidity among the elderly; a European overview. Soc Sci Med. 2003, 57: 861-873. 10.1016/S0277-9536(02)00454-9.
Bentham G, Eimermann J, Haynes R, Lovett A, Brainard J: Limiting long term illness and its associations with mortality and indicators of social deprivation. J Epidemiol Community Health. 1995, 49: S57-S64.
Bassuk SB, Berkman LF, Amick BC 3rd: Socioeconomic status and mortality among the elderly: findings from four US communities. Am J Epidemiol. 2002, 155: 520-533. 10.1093/aje/155.6.520.
Katz SJ, Kabeto M, Langa KM: Gender disparities in the receipt of home care for elderly people with disability in the United States. JAMA. 2000, 284: 3022-3027. 10.1001/jama.284.23.3022.
Lin G: Regional assessment of elderly disability in the U.S. Soc Sci Med. 2000, 50: 1015-1024. 10.1016/S0277-9536(99)00351-2.
Yamaguchi F, Kakehashi M: An analysis of factors concerned with the active life expectancy and the duration in care of the elderly. J Population Problems. 2001, 57 (4): 51-67.
Jordan H, Roderick P, Martin D: The index of multiple deprivation 2000 and accessibility effects on health. J Epidemiol Community Health. 2004, 58: 250-257. 10.1136/jech.2003.013011.
Krieger N, Chen JT, Waterman PD, Soobader MJ, Subramanian SV, Carson R: Geocoding and monitoring of US socioeconomic inequalities in mortality and cancer incidence: does the choice of area-based measure and geographic level matter?: the Public Health Disparities Geocoding Project. Am J Epidemiol. 2002, 156: 471-482. 10.1093/aje/kwf068.
Pickett KE, Pearl M: Multilevel analysis of neighbourhood socioeconomic context and health outcomes: a critical review. J Epidemiol Community Health. 2001, 55: 111-122. 10.1136/jech.55.2.111.
Ministry of Health, Labour and Welfare: Prefectural life table 2000. 2003, Tokyo: Health and Welfare Statistics Association
All-Japan Federation of National Health Insurance Organizations: Statistics for the long-term care insurance. 2002, Tokyo: All-Japan Federation of National Health Insurance Organizations, [http://www.kokuho.or.jp/shiryou/index.html]
Ministry of Health, Labour and Welfare: Municipal life table 2000. 2003, Tokyo: Health and Welfare Statistics Association
Welfare and Medical Service Agency: Number of the certificated people of long term care insurance. 2003, Tokyo: Welfare and Medical Service Agency, [http://www.wam.go.jp/kaigo/]
Fukuda Y, Nakamura K, Takano T: Municipal socioeconomic status and mortality in Japan: sex and age differences, and trends in 1973–1998. Soc Sci Med. 2004, 59: 2435-2445. 10.1016/j.socscimed.2004.04.012.
Leyland AH, McLeod A: Mortality in England and Wales, 1979–1992. 2000, Glasgow: MRC Social & Public Health Sciences Unit, [http://www.msoc-mrc.gla.ac.uk/Publications/pub/occasional_MAIN.html]
Rashash J, Browne W, Goldstein H, Yang M, Plewis I, Healy M, Woodhouse G, Draper D, Langford I, Lewis : A user's guide to MLwiN. 2002, London: Centre for Multilevel Modelling, Institute of Education, University of London
Nakaya T: Spatial epidemiology. Health GIS. Edited by: Nakaya T, Tanimura S, Nihei N, Horikoshi Y. 2004, Tokyo: Kokon Shoin, 74-121.
Wakefield JC, Best NG, Waller L: Bayesian approaches to disease mapping. Spatial Epidemiology. Edited by: Elliot JC, Wakefield NG, Best NG, Briggs DJ. 2000, New York: Oxford University Press, 104-127.
Bernardinelli L, Montomoli C: Empirical Bayes versus fully Bayesian analysis of geographical variation in disease risk. Stat Med. 1992, 11: 983-1007.
Statistical Information Institute for Consulting and Analysis: Municipal data of social indicators. 2003, Tokyo: Statistical Information Institute for Consulting and Analysis
Gutierrez-Fisac JL, Gispert R, Sola J: Factors explaining the geographical differences in Disability Free Life Expectancy in Spain. J Epidemiol Community Health. 2000, 54: 451-455. 10.1136/jech.54.6.451.
Ministry of Health, Labour and Welfare: Vital statistics of Japan 2001. 2003, Tokyo: Health and Welfare Statistics Association
Barnett S, Roderick P, Martin D, Diamond I, Wrigley H: Interactions between three proxies of health care need at the small area level: an urban/rural comparison. J Epidemiol Community Health. 2002, 56: 754-761. 10.1136/jech.56.10.754.
Kishi R, Horikawa N: Role of the social support network which influences age of death and physical function of elderly people study of trends in and outside of Japan and future problems. Nippon Koshu Eisei Zasshi. 2004, 51: 79-93.
Iwasaki M, Otani T, Ohta A, Sasazawa Y, Kuroiwa M, Suzuki S: Rural-urban differences in sociodemographic, social network and lifestyle factors related to mortality of middle-aged Japanese men from the Komo-Iso cohort study. J Epidemiol. 2002, 12: 93-104.
Kawachi I, Kennedy BP, Lochner K, Prothrow-Stith D: Social capital, income inequality, and mortality. Am J Public Health. 1997, 87: 1491-1498.
Kishi R, Eguchi T, Sasatani H, Yaguchi T: Health status, social networks and support systems of elderly men and women in a large city and in former coal-mining areas – a population-based comparative study of 70-year-old residents of Sapporo and Yubari. Nippon Koshu Eisei Zasshi. 1994, 41: 474-488.
Iwasaki M, Otani T, Sunaga R, Miyazaki H, Xiao L, Wang N, Sasazawa Y, Suzuki S: Social networks and mortality based on the Komo-Ise cohort study in Japan. Int J Epidemiol. 2002, 31: 1208-1218. 10.1093/ije/31.6.1208.
Sugisawa H, Fukaya T, Sugihara Y, Ishikawa H, Nakatani Y, Kim HK: Factors related to under-utilization of in-home services under the long term insurance system. Nippon Koshu Eisei Zasshi. 2002, 49: 425-436.
Ikegami N: Public long-term care insurance in Japan. JAMA. 1997, 278: 1310-1314. 10.1001/jama.278.16.1310.
Best N, Richardson S, Thomson A: A comparison of Bayesian spatial models for disease mapping. Stat Methods Med Res. 2005, 14: 35-59. 10.1191/0962280205sm388oa.
Leyland AH, Davies CA: Empirical Bayes methods for disease mapping. Stat Methods Med Res. 2005, 14: 17-34. 10.1191/0962280205sm387oa.
Lawson AB, Biggeri AB, Boehning D, Lesaffre E, Viel JF, Clark A, Schlattmann P, Divino F: Disease mapping models: an empirical evaluation. Disease Mapping Collaborative Group. Stat Med. 2000, 19: 2217-2241. 10.1002/1097-0258(20000915/30)19:17/18<2217::AID-SIM565>3.0.CO;2-E.
Ghosh M, Rao JNK: Small area estimation an appraisal. Statistical Science. 1994, 9: 55-76.
Nakaya T: An information statistical evaluation of areal units for mortality mapping: toward constructing an analyzing system of disease maps on geographical information system. Theory and Applications of GIS. 1996, 4: 53-60.
Hirako T, Sakeki N, Nakaya T: Application of Bayesian methods on municipal health statistics. Koseino Shihyo. 1999, 46 (10): 3-11.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/5/65/prepub
This study was supported by a Grant-in-Aid for Scientific Research by the Japan Society for the Promotion of Science (Grant No. 14570326 and 16590497).
The authors wish to thank Dr Masahiro Umezaki (Graduate School of Tokyo Medical and Dental University) for technical advice on geographic mapping and Dr Tomoki Nakaya (Ritsumeikan University) for helpful comments on statistics.
The author(s) declare that they have no competing interests.
YF designed the study, analyzed the data, and drafted the article. KN helped to interpret the results and edited the draft. TT supervised the data analysis and writing the article.