Article Text
Abstract
Objectives: To simultaneously examine the effects of area-level and individual-level socioeconomic position on fatal injuries in children <5 years of age.
Methods: A retrospective cohort study based on the national birth and death registers of Korea. 2 667 060 children born during 1995–8 were followed up from birth to the 5th birthday. Cumulative incidences of fatal injuries were calculated, and through multilevel Poisson regression models, relative risks (RRs) of incidence rate were estimated according to children’s sex, father’s occupation and mother’s education at individual level, and deprivation and degree of urbanity at area level.
Results: Girls had lower risk for fatal injuries than boys (RR 0.81; 95% confidence interval (CI) 0.75 to 0.87). Compared with children with fathers in non-manual occupations, those with fathers in manual (RR 1.45; 95% CI 1.34 to 1.58) or other occupations (RR 1.35; 95% CI 1.13 to 1.62) had higher risk. Children with mothers who were high school graduates (RR 1.23; 95% CI 1.12 to 1.36) or junior school graduates (RR 1.91; 95% CI 1.66 to 2.19) had higher risk than those whose mothers were college graduates. After controlling for individual-level variables, residence in more deprived districts (RR 1.13; 95% CI 1.05 to 1.21) or non-metropolitan regions (urban RR 1.34; 95% CI 1.22 to 1.47 and rural RR 1.61; 95% CI 1.40 to 1.86) was significantly associated with increased risk.
Conclusions: Both individual-level and area-level socioeconomic position influenced the risk for childhood fatal injuries. To reduce the socioeconomic inequalities and the absolute burden in Korea, universal strategies should receive priority, and special efforts in implementation should be directed towards both disadvantaged households and areas.
- NBR, National Birth Registration database
- NDR, National Death Registration database
- PINs, personal identification numbers
- SEP, socioeconomic position
Statistics from Altmetric.com
- NBR, National Birth Registration database
- NDR, National Death Registration database
- PINs, personal identification numbers
- SEP, socioeconomic position
Injury is a leading cause of morbidity and mortality among children in developing as well as in developed countries.1 As a result of effective social interventions and increased awareness, however, mortality in children from injury has been continuously declining for decades in many developed countries.2–4 Injury is not a random event, and like many other health problems, is more concentrated among socioeconomically disadvantaged groups. On the basis of their thorough review, Cubbin and Smith 5 reported that the consistent inverse relationship between socioeconomic position (SEP) and fatal unintentional injuries was found across all age groups, regardless of the SEP measures used.5 In particular, they pointed out that area-level characteristics (eg, neighbourhood deprivation) as well as individual-level SEP may influence the risk for morbidity and mortality from injury.
Why do places matter for injury risk? One reason is that different socioeconomic groups face different hazards within their social contexts. For young children, environmental factors are especially important as their ability to control their surroundings is quite limited. Recent studies using multilevel analysis found that both individual-level and area-level factors influence childhood injury.6,7,8,9,10 However, those multilevel studies have been limited in several respects. Firstly, although various individual-level variables, including children’s age or sex or maternal age or education,6,8,9,10 or a survey respondent’s information,7 were used across studies, they could not properly reflect the SEP of household, which could cause insufficient control for compositional confounding. Secondly, only morbidity from injury was examined through collection of hospital records or self-reports, so bias due to healthcare access or subjective reporting could not be ruled out. Thirdly, all the above-mentioned studies were conducted in Western industrialised societies, and the evidence for socioeconomic inequalities in childhood injury is scarce in emerging economies or developing countries.
Our study has three distinctive features, firstly, we sought to simultaneously identify the effects of area-level as well as individual-level SEP on fatal childhood injury through detailed measurement at both levels. Secondly, mortality from injury is used as an outcome from the nationally representative sample. Finally, it is the first multilevel study about childhood injuries conducted in a non-Western industrialised setting—namely, in South Korea (hereafter Korea), an example of an emerging economy with the highest rank in childhood fatal injuries during 1991–5 among Organization for Economic Co-operation and Development countries.11
METHODS
Sources of data
Data were obtained from three different sources: (1) the National Birth Registration database (NBR) between 1995 and 1998; (2) the National Death Registration database (NDR) between 1995 and 2002; and (3) the National Census in 1995.
From the NBR, we obtained the records of 2 791 962 live births between 1995 and 1998 nationally, and through the personal identification numbers (PINs), they were linked to the NDR for the subsequent 5 years (from birth to the 5th birthday) by the staff of the National Statistical Office. A total of 4698 duplicated records were found, with 256 redundant PINs, which may have resulted from incomplete PINs due to administrative problems in the birth registration process. These duplicates were excluded from analysis. Next, a series of information was abstracted: birth date, sex of newborn, residential address (district), occupation and educational attainment of father and mother. Among children born between 1995 and 1998, a total of 12 255 deaths were identified during the 5-year follow-up period. For the causes of death, the NDR adopted the 10th revision of the International Classification of Disease. A recent Korean study reported that the overall proportion of 1990–2002 deaths <10 years of age certified by doctors was >90%, with continuous increases across successive years, and the cases of external cause were more likely to be diagnosed by doctors.12 We defined fatal injuries when the child died from external cause with one of the following codes: “V”, “W”, “X”, or “Y”.13 Excluding the records with missing values in socioeconomic variables and cause of death, we analysed the records for 2 667 060 children (95.7% of the initial data) and 2926 deaths from injuries (23.9% of total deaths), including 1214 deaths from transportation-related causes, 611 from suffocation, 367 from fall and 324 from drowning.
For individual-level SEP measures, we used mother’s education and father’s occupation from the NBR.3,9,14,15 Father’s occupation was categorised into three groups: (1) non-manual (legislator or administrator, professional, engineer and clerk); (2) manual (service or sales worker, farmer or fisher, skilled labourer, operator and unskilled labourer); and (3) others (unemployed, economically inactive or unknown). Mother’s education was classified as (1) college graduate or above, (2) high school graduate, or (3) junior high school graduate or below. In a preliminary analysis, we examined the effect of father’s education instead of father’s occupation, and found a gradient as observed for mother’s education. However, the models including father’s occupation were found to explain more of the area-level variation than models with father’s education, and even when both variables were simultaneously included in a model, the results changed little. This may be owing to the fact that the study subjects were relatively well educated (91.8% graduated from high school or above) and so the education level was less sensitive to differences in social status.
Considering that both material and informational resources are important for childhood safety, informational resources could be represented by the education level of the mother who is usually a caregiver in Korea, whereas material resources could be better captured by the father’s occupation.
We used districts as the spatial unit of analysis, which is administratively defined in Korea and has a relatively large and heterogeneous population, ranging from 12 000 to 630 000. The district is the smallest municipal unit capable of autonomous policy implementation in Korea and is also the lowest statistical unit for which official data are available. For area-level SEP measures, degree of urbanity,16–19 which is designated by the government (metropolitan, urban and rural), and a deprivation index6,8,20,21 were assigned to 247 districts across the country. The deprivation index was made through modification of the Townsend and Carstairs index.22 On the basis of the National Census data in 1995, we calculated the district-specific proportions of households: living in apartments; without a car; in a crowded condition (>1.5 persons/room); with a female head; in tenancy; and under substandard living conditions (without hot-water supply, a flush toilet or a modern kitchen). In addition, the following district-specific proportions of individuals were calculated: high school graduate or above among adults >25 years of age; unemployed among economically active men; in manual occupations among employed men; >65 years of age, a higher proportion of which could indicate a high level of poverty in Korea where senior pensions are not sufficient for the elderly who lack other sources of income.23 All proportions were z standardised, and a factor analysis with varimax rotation was carried out, which showed that home ownership and unemployment rate belonged to a factor different from the others. This might come from the characteristics of rural region, where the residents are likely to own houses with low market value, and the unemployment rate is relatively low owing to the existence of subsistence farming. Finally, we constructed the district-specific deprivation index by averaging the z-standardised scores of eight proportions (household not living in apartments; without a car; in a crowded condition; with a female head; and under substandard living conditions, and individuals with below high school graduation among adults >25 years of age, >65 years of age and in manual occupations among employed men), and linked it to the NBR based on the residential address. A positive, larger score means higher levels of deprivation. Deprivation indices developed in Western countries may not be transferable to the Asian context. A Japanese study found differences in associations with health between an index that included municipal-level unemployment and overcrowding, compared with an index that included education and per capita income.24 In this study, we used one composite index of deprivation, excluding unemployment rate and home ownership, but adding the variable “degree of urbanity” instead of creating another composite index.
Statistical analysis
To improve mathematical efficiency while securing individual-level information, we constructed a new dataset in an aggregated format, with a hierarchical multilevel structure.25,26 Through cross tabulation by individual-level predictors (child’s sex, father’s occupation and mother’s education), from about 2.6 million records, we obtained a total of 4438 cells at level 1 within 247 districts at level 2.
The response variable, the number of deaths in each cell, was modelled with the offset, the expected death count derived from person-time accounts. Estimates were obtained from multilevel Poisson regression models with log link function, using the penalised quasi-likelihood, second-order approximation procedure.27 We fitted the models through the statistical software MLwiN V.2 (Institute of Education, London, UK), which uses the iterative generalised least squares algorithm.28
We examined a series of models from simple to complex by adding individual-level and area-level predictors, and intra-level and cross-level interaction terms, through which we could evaluate the relationship between fatal injuries and individual-level SEP (fixed parameters) conditional on area-level variations (random parameters). In addition, the degree to which district deprivation and urbanisation could explain the mortality variation across districts (fixed parameters) was examined.
RESULTS
For a total of 2 667 060 births and 13 104 907 person-years, 2926 fatal injuries were identified (cumulative incidence 1.1/1000 children; incidence rate 22/100 000 person-years). As for injury mechanism, transport-related cause (0.46/1000) was the most common, followed by suffocation (0.23), falls (0.14) and drowning (0.12; data not shown). The size of birth cohorts decreased as time passed, from 698 889 to 628 224, and the rates decreased from 1.3 to 0.9/1000 across successive birth cohorts. Within every cohort, mortality was higher in boys than in girls (table 1).
Comparing the cumulative incidence of fatal injuries between socioeconomic groups, it was higher in children born from fathers in manual occupations (1.4/1000) or others (1.2/1000) than in those born from fathers in non-manual occupations (0.8/1000). Mother’s education also showed similar effects; children born from mothers with more education had lower fatal injury rates. As for area level, metropolitan regions had the lowest rate (0.8/1000), less than a half of the rates in the rural region (1.8/1000). In addition, the more deprived the districts, the higher the fatal injury rates. These findings were observed in both boys and girls (table 2).
Figure 1 shows the association between fatal injury rates and area characteristics, in which the circle size is proportional to the population size. Overall, we observed that non-metropolitan regions were more likely to be deprived, and more deprived districts had higher mortality with more variability.
Table 3 shows fixed and random parameters estimated from multilevel Poisson regression models. Significant variation was seen in mortality across districts (model 1). When individual-level variables were added (model 2), all of them (child’s sex, father’s occupation and mother’s education) were statistically significant and random variation between districts decreased by 36%, which meant that about one third of the area-level variation came from differences in the population composition in each district. When the district deprivation variable was added (model 3), the between-districts variation decreased by an additional 25.4%, and the effects of individual-level variables were attenuated. By adding the urbanity variable (model 4), we observed further reduction in the between-districts variation (additional 26.3%) and risk attenuation for district deprivation. As for intralevel interaction, there was no significant interaction between father’s occupation and mother’s education, or between deprivation and urbanity (data not shown). However, we found a significant cross-level interaction; district deprivation had a greater effect on fatal injuries among the children whose parents were from lower SEP strata (model 5 and 6). The between-districts variation became insignificant when both the degree of urbanity and deprivation index were included in the models (table 3).
For common types of injury, including transportation-related causes, suffocation, falls and drowning, we estimated the effects of urbanisation and deprivation, controlling for individual-level variables. Compared with metropolitan regions, rural areas had higher risk for all types of injury. Deprivation showed a clear positive relationship with mortality by transportation-related causes and drowning (fig 2).
DISCUSSION
From multilevel Poisson regression models, we found that male sex of child, lower occupational class of father and lower education of mother were significantly associated with increased risk for fatal injuries among children <5 years. In addition, residence in more deprived districts or non-metropolitan regions was associated with higher risk even after controlling for individual-level variables.
Methodological issues
Our study has several strengths from a methodological perspective. Firstly, this analysis was based on nationwide birth cohorts, ensuring a large sample size and representativeness. Secondly, our outcome was identified through retrospective follow-up. As the exposure had been measured before the outcome, we could avoid the possibility of reverse causation that childhood fatal injuries resulted in the downward mobility of the affected households. Thirdly, although previous multilevel studies used demographic variables such as maternal age or marital status as a proxy for individual-level SEP,6–8 we used direct measures, both father’s occupation and mother’s education. Finally, we adopted multilevel models, through which we could simultaneously examine the effects of individual-level and area-level SEP and explore to what extent the fixed parameters explained the area-level variation in fatal injuries.29 In addition, on the basis of subanalysis, we could examine types of injury in relation to area.
However, our study has some limitations. Firstly, we included only data without missing values on SEP variables. However, when comparing the fatal injury rates between the excluded group and the study subjects, we did not find significant differences (fatal injury rates 1.0 and 1.1/1000 children, respectively, p = 0.75). Secondly, although we found socioeconomic inequalities in injury mortalities, we could not identify whether they resulted from higher incidence of injuries, or lower survival due to worse severity or inappropriate healthcare access. Previous studies showed that incidence in itself was high, and injuries by fatal mechanism were more common in lower SEP groups,20,30 but further studies analysing incident cases and subsequent management are warranted. Thirdly, although we more thoroughly measured the individual-level SEP compared with previous studies, the possibility of residual confounding arising from another aspects such as income or assets, could not be completely ruled out. Fourthly, we used the area information for the districts where the child was born, and therefore the place of injury was not necessarily the same as that of birth. However, the consistent gaps observed across areas strongly suggested the influence of the SEP at area-level. In addition, our study differs from previous studies of contextual variations in injury,6–8,10,31 which have tended to use smaller spatial units of aggregation (eg, neighbourhoods and wards). The causal processes underlying area variations in injury are likely to differ according to the spatial scale of analysis. Thus, neighbourhood variations in injury have been linked to local differences in traffic patterns, as well as supervision of children and play activity. As our study examined variations in injury across larger spatial scales (districts), we were unable to consider these local causal processes. Nonetheless, our findings have value for planning purposes by pointing out area-based socioeconomic disparities in injury mortality rates across larger geographical scales.
Finally, although we identified the area differentials in some common types of injury, there are points to be explored further; for instance, among transportation-related injuries, traffic design and physical environment at area level would be more important for pedestrian, whereas safety device or regulations throughout the country would be so for passengers. However, we could not obtain such detailed information on subtypes of injury. For more elaborate interventions at individual and area levels, further studies on specific types of injury should be followed up, based on updated mortality data.
Public health implications
We found that boys had a higher risk for fatal injuries across all subgroups (tables 1–3), which was consistent with previous studies.6,8,9,10
Children born from lower social groups, such as fathers in manual occupations or mothers with less education, had higher risk than the others, regardless of their sex (tables 2 and 3). This could be explained by lack of supervision, poor housing environment (eg, steep stairway and old building) or inadequate access to information and preventive behaviour.3,5,32 That is, obtaining quality information about child safety, sparing time and efforts for supervision, and living in a safer house is constrained by socioeconomic circumstances. In this regard, universally applicable social measures regardless of socioeconomic resources of individual households, such as legislation, enforcement of product standard or financial outlay, seem desirable rather than educating high-risk groups.33
As for area effects, living in non-metropolitan regions (especially rural) or more deprived districts was independently associated with increased risk for fatal childhood injuries, even after controlling for individual-level SEP. A previous multilevel study of area variation in childhood injury rates in the UK found that family factors could explain most of the area-level variation.31 The difference between that study and our findings could be attributable to: (1) the smaller spatial units (wards) examined in the British study; (2) the control for measures of family safety practices (in addition to socioeconomic factors) in the British study; and (3) the fact that the study areas in the British study included only deprived wards within one region, so that there was insufficient variation in socioeconomic conditions between the wards. We identified that a portion of variation across the districts came from different population compositions (models 1 and 2). However, the greater portion of variation was explained by deprivation and degree of urbanity (models 3 and 4). In addition, as the administrative designation of “rural/urban/metropolitan” on districts was based on their population size and share of primary industries (eg, agriculture and fishing) in their local economies, rural districts were likely to show more variability in mortalities due to small population size and to be more deprived (fig 1). However, these contextual variables were significant after simultaneous control for each other, suggesting that one trait could not be substituted with the other.
This contextual effect might come from a series of factors, ranging from the general resources and investment level of the places to the specific hazards or institutions relevant to injuries, including (1) environmental hazards such as dangerous streets or unprotected industrial and building sites;6,8 (2) municipal goods and services such as police, parks and healthcare facilities;20,34 (3) specific safety policies;9 or (4) socially and culturally determined attitudes or practice.35 Looking at area differences according to types of injury, transportation and drowning showed a clear gradient across deprivation, which might be mediated by unprotected exposure to traffic and lack of safe playgrounds in more deprived districts.6–8,10,20,21 On the other hand, non-metropolitan residence was associated with higher risk for most types of injury even at the same level of deprivation.16,17,19 Especially, injuries by drowning, falls, and transportation showed the largest gaps, which may be related to the physical environment unique to rural areas (eg, neglected water reservoirs),5,17 unsafe housing conditions (eg, lack of window guards)7,10 or insufficient neighbourhood resources (eg, lack of safe playgrounds).5
What is already known
-
Many studies have reported that childhood injury, like other health problems, is more common among socioeconomically disadvantaged groups.
-
However, until now, the effects of individual-level and area-level SEP have not been examined simultaneously, and evidence from non-Western societies is lacking.
What this paper adds
-
Both individual-level and area-level SEP influenced the risk for fatal childhood injuries, and the area differentials were persistent across common types of injury.
-
Injury prevention efforts should be directed towards both disadvantaged households as well as areas.
In addition, when examining cross-level interactions, synergistic effects were observed (models 5 and 6); children born from fathers in lower occupational class or mothers with less education and living in more deprived districts had an even higher risk for fatal injuries than those who were placed in only one group. Special attention should be paid to these groups in monitoring and implementing social interventions.
How can these inequalities in fatal childhood injuries be ameliorated? Some studies advocate targeted interventions for specific types of injury or vulnerable groups.8,20,30,36,37 Other studies emphasise universal strategies to improve the general safety level and societal measures to remedy the underlying structure of social inequalities.3,6,17,33,38 Our findings support the second view; in Korea, the absolute risk for childhood injury is higher compared with other developed countries,11 and we found that most types of injury occur more commonly in lower social groups at individual and area levels.
To reduce the gap between social groups as well as the absolute risk, universal strategies including safety regulation and balanced regional development should receive priority, and special efforts in implementation should be directed towards both disadvantaged households as well as areas. In addition, societal measures to mitigate underlying socioeconomic inequalities are required.
REFERENCES
Footnotes
-
Competing interests: None declared.
-
M-H K and C-Y K were recipients of the 2004–2005 David E Bell Fellowship program at the Harvard Center for Population and Development Studies, and Takemi Fellowship program at the Harvard School of Public Health, respectively. A substantial part of this paper was written during their fellowship.
Linked Articles
- In this issue