RESEARCH
Under-registration of deaths in Thailand in 2005-2006: results of cross-matching data from two sources
Sous-enregistrement des décès en Thaïlande sur la période 2005-2006: résultats des données croisées émanant de deux sources
Subregistro de muertes en Tailandia en el periodo comprendido entre 2005 y 2006: resultados de la obtención de datos de pruebas cruzadas procedentes de dos fuentes
Patama Vapattanawong*; Pramote Prasartkul
Institute for Population and Social Research, Mahidol University, 999 Phuttamonthon 4 Rd. Salaya, Phuttamonthon, Nakhonpathom 73170, Thailand
ABSTRACT
OBJECTIVE: To estimate the magnitude of under-registration of deaths, by age and sex, in Thailand.
METHODS: The data in this study were derived from two sources: the Thai Survey of Population Changes (SPC) 2005-2006, a consecutive multi-round household survey conducted over a 12-month period, and Thailand's vital registration records. SPC death entries for people of all ages were matched to 2005-2006 death records from vital registration. The principles of a dual records system were applied to estimate the magnitude of under-registration of deaths, classified by age and sex, using the Chandrasekaran-Deming formula.
FINDINGS: Overall under-registration of deaths during 2005-2006 was 9.00% (95% confidence interval, CI: 8.95-9.05) for males and 8.36% (95% CI: 8.31-8.41) for females. For both males and females, under-registration decreased as age increased. Under-registration was greatest among people of either sex aged 1-4 years, whereas it was < 10% among people 60 years of age and older, both males and females.
CONCLUSION: These findings provided correction factors that can be used for adjusting mortality data from the registration system.
RÉSUMÉ
OBJECTIF: Estimer l'ampleur du sous-enregistrement des décès par âge et sexe en Thaïlande.
MÉTHODES: Les données de cette étude proviennent de deux sources: l'enquête thaïlandaise sur les changements démographiques (ECD) de 2005-2006, une enquête consécutive à tours multiples conduite auprès des ménages et menée sur une période de 12 mois, ainsi que les registres nationaux d'état civil. Les entrées de l'ECD pour les décès de personnes de tous âges ont été comparées aux mentions de décès faites aux registres d'état civil pour la période 2005-2006. Les principes d'un système à double enregistrement ont été appliqués afin d'évaluer l'ampleur du sous-enregistrement des décès, classés par âge et sexe, à l'aide de la formule de Chandrasekaran-Deming.
RÉSULTATS: Dans l'ensemble, le sous-enregistrement des décès sur la période 2005-2006 a atteint 9,00% (intervalle de confiance de 95%, IC: 8,95-9,05) chez les hommes et 8,36% (IC de 95%: 8,31-8,41) chez les femmes. Pour les deux sexes, le sous-enregistrement des décès diminue au fur et à mesure que l'âge augmente. Le sous-enregistrement était le plus important chez les enfants âgés de 1 à 4 ans, quel que soit leur sexe, alors qu'il était inférieur à 10% chez les hommes et femmes de 60 ans et plus.
CONCLUSION: Ces résultats ont fourni des facteurs de correction qui peuvent être utilisés pour ajuster les données relatives à la mortalité à partir du système d'enregistrement.
RESUMEN
OBJETIVO: Estimar la magnitud del subregistro de muertes, por edad y sexo, en Tailandia.
MÉTODOS: Los datos incluidos en este estudio proceden de dos fuentes: el Estudio tailandés de variaciones de población (EVP) de 2005-2006, un estudio de hogares con visitas múltiples consecutivas realizado en un periodo de 12 meses y los registros civiles de Tailandia. Las notificaciones de muerte del EVP para personas de todas las edades se combinaron con los registros de muertes para el periodo comprendido entre 2005 y 2006 del registro civil. Se aplicaron los principios de un sistema de registros dual para estimar la magnitud del subregistro de muertes, clasificado por edad y sexo, mediante el uso de la fórmula Chandrasekaran-Deming.
RESULTADOS: El subregistro global de muertes durante el periodo comprendido entre 2005 y 2006 fue del 9,00% (Intervalo de confianza del 95%, IC: 8,95-9,05) para hombres y del 8,36% (IC del 95%: 8,31-8,41) para mujeres. Tanto para hombres como para mujeres, el subregistro descendió a medida que incrementaba la edad. El subregistro fue mayor entre personas de ambos sexos con edades comprendidas entre 1 y 4 años, e inferior al 10% entre personas mayores de 60 años, tanto hombres como mujeres.
CONCLUSIÓN: Estos resultados proporcionaron factores de corrección que se pueden utilizar para ajustar los datos de mortalidad del sistema de registro.
Introduction
Mortality is a demographic indicator not only of health outcomes in various population groups, but also of a country's state of development. Death data are fundamentally crucial for identifying health problems and for monitoring health programmes. They are also used to measure and compare mortality rates in subpopulations within the same country or across different countries. Accordingly, the quality of mortality data is the key to reliable indicators and should therefore be carefully assessed before the data are used.
In Thailand, the most important source of mortality data is the vital registration system, which was initiated in 1916. Civil registration law mandates that every vital event (births and deaths) be registered at the offices of the district or municipality registrars, which are under the interior ministry. According to this law, births must be registered within 15 days of delivery, while deaths must be registered within 24 hours. For deaths, the registration process is not straightforward because two steps are involved: (i) notifying an authorized person about the death, and (ii) validating and registering the death.1 After both registration steps are completed, the name of the deceased is deleted from the household roster and the details concerning the death are entered into the registration system. Before the end of each year the interior ministry counts the births, deaths and total population and makes the figures publicly available.
In addition to the vital registration system, Thailand has a second important source of population data: the population census, which is conducted every 10 years. Theoretically, these two sources of demographic data should provide very similar population figures, but in reality the data sometimes conflict. During the first Population Seminar held in Bangkok in March 1963, disagreements arose surrounding the rate of population growth in Thailand because vital registration and population census data were in conflict. To learn the real rate of population growth, the National Statistical Office, which is in charge of the population census, was requested by the cabinet to conduct a national survey known as the Survey of Population Change (SPC). The first SPC was undertaken in 1964-1966. It later became an intercensal population survey conducted every 10 years at about the mid-period between the two censuses. The main purposes of the SPC are to obtain reliable estimates of vital event rates and of the rate of population growth, as well as to measure the under-registration of vital events. Thailand's most recent SPC was conducted in 2005-2006.
Although estimating under-registration is a main objective of each SPC, the methods used in the surveys have not been the same. The differences in such methods are described in detail elsewhere.2 In brief, the estimation methods used in the first three SPCs (1964-1966, 1974-1976, 1984-1986) consisted of a dual records procedure, while those used in the 1995-1996 and 2005-2006 SPCs were based directly on survey questions, whether vital events were registered or not. According to the SPC results, during the 1960s and 1970s approximately 60% of the deaths in the Thai population were registered.3,4 In the mid-1980s the fraction rose to 76%5 and in the two most recent SPCs (1995-1996 and 2005-2006) 95% completeness was reported.6,7 The very large increase in vital event registration observed in the 1995-1996 SPC relative to previous SPCs raised doubts, and the completeness of death registration in Thailand remains a controversial issue requiring further scrutiny.
This study aimed to investigate the quality of mortality data in Thailand's civil registration system. Its objectives were to estimate the magnitude of under-registration of deaths and its associated age and sex patterns in 2005-2006 by cross-matching data from two sources.
Methods
Data sources
The two sources of data used in this study were the SPC 20052006 and vital registration records. The SPC 2005-2006 was a national, longitudinal, multi-round household survey with a large sample: 2050 sample rural villages/urban blocks and 82 000 households.7 This survey was designed to determine changes in household composition. To overcome the errors generated by the respondents' limited ability to correctly recall events that occurred over an extended period, each sample household was visited five times within a year, or about once every three months. The survey period began in July 2005 and ended in July 2006.
All households in the sample rural villages and urban blocks were visited by the enumerators to obtain information on the usual household residents. In the initial interview, all persons residing in each household were listed and their basic demographic characteristics were recorded. A total of 327 735 individuals were counted in this first interview. As an additional check, to obtain as complete a coverage of births and deaths as possible, the enumerator also asked in every household: (i) if any women were known to be pregnant, and, if so, the number of months of pregnancy; and (ii) if any individuals were known to be seriously ill. For the second and subsequent interviews, the enumerator determined and recorded any changes in household composition that had taken place since the previous visit on account of births, deaths or movements in and out. If such changes had taken place, the date when they occurred was recorded. An additional question was asked about any babies born or deceased since the enumerator's last visit.
To obtain vital registration data, individual death records for 2005-2006 from the interior ministry were used. Each record was in electronic format and had a 13-digit unique personal identifier (PID). Death records from before and after the SPC survey period were excluded from the study.
The individual records of the SPC had been linked to death records from vital registration using a 13-digit unique PID never recorded in previous SPCs except the last one. Of a total of 327 735 records in the SPC, 37 110 were excluded because their unique PIDs were missing or incomplete. Thus, only 290 625 records (88.7%) could be used to link to vital registration data. The method allowed for four possible outcomes from the cross-matching: (i) deaths reported in the SPC matched to deaths in the vital event registry; with (ii) deaths reported in the SPC but not matched to deaths in the registry; with (iii) people reported as survivors in the SPC but matched to deaths in the registry; and (iv) people reported as survivors in the SPC but not matched to deaths in the registry. All of these outcomes except the fourth were used to analyse the under-registration of deaths by age and sex.
Data analysis
This study applied the Chandrasekaran-Deming formula to carry out a dual records system estimation. The dual records system is a method used to estimate the completeness of a data source by cross-matching its records case by case against those of another source. We used the Chandrasekaran-Deming formula primarily to estimate the data that were missing from both systems, i.e. the number of people known to have died but not reported in either data source. We used this approach based on four crucial assumptions: (i) that the population of interest was closed; (ii) that events recorded in both systems could be cross-matched; (iii) that the probability of an event occurring in one system was independent of the probability of it occurring in the other, and (iv) that each case had an equal probability of ascertainment in either source. The under-registration of deaths was calculated by constructing a 2 x 2 contingency table (Table 1) containing the results of cross-matching the deaths from the SPC with those from vital registration records.
To estimate total events (N) according to the Chandrasekaran-Deming formula:8
where M (matched) stands for the events found in both systems; NS (not surveyed) stands for the events missed by the survey but found in vital registration records; and NR (not registered) stands for the events missed by registration but found in the survey. An estimate of the variance of N was calculated using the following formula:
where:
and
The 95% confidence intervals (CIs) of N were calculated using the formula:
The formula of p1 gave the estimated completeness of death registration. Thus, the under-registration (Du) could be estimated from:
and the variance of q1 was calculated as:9
Results
Cross-matching of deaths
The diagram of the cross-matching process between deaths listed in the SPC and vital registration data is shown in Fig. 1. Of 290 625 individual records from the SPC, 1882 were deaths that could be matched to vital registration data; 179 were deaths that could not be matched; 556 were people listed as survivors in the SPC who could be matched to vital registration death records; and 288 008 were people listed as survivors in the SPC who could not be matched to vital registration death records (Table 2).
As shown in Table 2, the SPC listed a total of 2061deaths as having occurred during 2005-2006. Over that period, 2438 deaths were included in the vital event registry. The figures in this table also reflect counting errors, since 556 people listed in the death registry were reported as being alive during the survey. The table also shows the magnitude of the under-registration: 179 people reported to be dead during the survey were not registered in the vital event registry.
Estimation of under-registration of deaths
A second 2 x 2 contingency table (Table 3) was created with the death data in Table 2 to determine how many people had been missed by both systems (the cell with '???'). The estimated total number of deaths occurring during the 2005-2006 survey and its 95% CI were:
where the result was rounded to the nearest integer. In addition, the estimated under-registration of deaths, and its 95% CI could be calculated as:
To estimate the under-registration of deaths by age and sex, we repeated the above calculation using data classified by age and sex instead. Table 4 presents the age- and sex-specific numbers of deaths in 2005-2006 as recorded by the SPC and vital registration records, and Table 5 shows the estimated percentage of under-registration of deaths. In 2005-2006 under-registration was 8.69% (95% CI: 8.65-8.72) for both sexes combined and, for all ages overall, it was 9.00% (95% CI: 8.95-9.05) for males and 8.36% (95% CI: 8.31-8.41) for females.
When we stratified by broad age groups, the highest under-registration was found among children 1-4 years old: 54.55% (95% CI: 54.25-54.84) in males and 71.43% (95% CI: 71.11-71.75) in females. Under-registration of male and female deaths decreased as age increased, with the sole exception of deaths among females less than 1 year of age, which showed much lower under-registration than deaths among male infants (8.33%; 95% CI: 7.87-8.80 versus 34.78%, 95% CI: 34.48-35.08, respectively). Furthermore, deaths among older males and females showed less than 10% under-registration.
Discussion
The aim of this study was to estimate the magnitude of under-registration of deaths by making use of two valuable sources of data that already existed in Thailand. Although the method used herein has been known for more than half a century, it is still useful, particularly in countries having several sources of mortality data, such as Egypt10 and the Philippines.11 This method is frequently used in studies and is often referred to in epidemiology and ecology as capture-recapture.9,12-18
In Thailand, the dual records procedure or a similar method based on cross-matching of data from two different sources has been practised since the mid-1960s, and in that early period it was applied in a manner very similar to the one shown in this study in the sense that data from independent sources - in our case the survey and the vital event registration system - were cross-matched. However, some practical differences do exist. Whereas decades ago the deaths reported in the survey were manually cross-matched with the ones included in the registration system on the basis of name and other characteristics common to both systems, such as sex, age of deceased, place of death and place of residence,3-5 in this study manual cross-matching was replaced by computerized cross-matching using a 13-digit unique PID. In addition, to ensure correct cross-matching, we double checked the age and sex of each pair of cross-matching cases. Thus, this study was stronger than earlier studies. Moreover, the dual records procedure made it possible to generate estimates for population subgroups. Thus it is possible to estimate under-registration of deaths at a subnational level or among subgroups of the population if the data are available.
The estimated under-registration of deaths for both sexes combined and for males and females separately was 8.69%, 9.00% and 8.36%, respectively.
These figures are approximately 1.8, 1.7 and 1.9 times higher than the percentages estimated directly by Thailand's National Statistical Office from the questions in the SPC.7 This is because our study captured more deaths, particularly 556 deaths that were in the vital event registration system but that had not been counted in the survey. Two factors could account for the substantial under-reporting of deaths in the SPC (~20% of all deaths in both systems): (i) the assumption that the population of interest was closed and (ii) the quality of the counting of deaths in the survey. The life/death status of the 556 people who were in the death registry but not in the SPC had been carefully reassessed in the fifth (last) visit paid to the households during the survey period (data not shown). Of these people, 78.8% were found to be regular residents, 17.4% had moved away and 3.8% were temporarily absent or transient residents. Thus, the quality of the counting of deaths may have been the main reason that deaths were under-reported in the SPC, followed by the mistaken assumption that the population surveyed was a closed one. Similar findings have not been reported in other countries.11,19
Our figures may be an underestimation of the actual rate of underreporting. One reason is that not all SPC 2005-2006 sample data could be analysed because some information (11.3%) was incomplete; either the PID was missing (10.3%) or it had fewer than 13 digits (0.2%) or was duplicated (0.8%). Another reason is that only one unique identifier was used as a linkage key. If some identifiers were incorrectly recorded during the interview or when entered into the system, the results of the cross-matching could have been affected. Thus, this study arguably provides a minimum estimate of the rate of under-registration of deaths in Thailand during 2005-2006. Its findings lend support to the notion that the use of different estimation methods accounts for the large improvement in the completeness of death registration observed after the mid-1990s.
Conclusion
In Thailand, overall under-registration of deaths in 2005-2006 was around 8.69%, but the percentage varied by age and sex. The age- and sex-specific estimates of under-registration from this study can be used as correction factors to improve the accuracy of the death statistics in the country's registration system. Several studies have calculated the magnitude of under-registration overall, but few have made age- and sex-specific estimates. By using mortality data from two important sources - the SPC 2005-2006 and the vital registration system - this study provides some degree of assurance that the results are reliable. Age- and sex-specific numbers of deaths, after adjustment by the correction factors generated in this study, can be used as numerators to make various mortality estimates. The newly-adjusted age- and sex-specific death rates can even be used, for example, to construct life tables of Thailand's general population and subpopulation groups.
ACKNOWLEDGEMENTS
We thank Thailand's National Statistical Office and the Ministry of Interior for allowing us to use the Survey of Population Change and the mortality data in the vital registration system for this study.
Competing interests: None declared.
References
1. Prasartkul P, Vapattanawong P. The completeness of death registration in Thailand: Evidence from demographic surveillance system of the Kanchanaburi Project. World Health Popul 2006;8:43-51. PMID:18277108
2. Hill K, Vapattanawong P, Prasartkul P, Porapakkham Y, Lim SS, Lopez AD. Epidemiologic transition interrupted: a reassessment of mortality trends in Thailand, 1980-2000. Int J Epidemiol 2007;36:374-84. doi:10.1093/ije/dyl257 PMID:17182635
3. Report on the 1964-1966 Survey of Population Change. Bangkok: National Statistical Office; 1969.
4. Report on the 1974-1976 Survey of Population Change. Bangkok: National Statistical Office; 1979.
5. Report on the 1985-1986 Survey of Population Change. Bangkok: National Statistical Office; 1989.
6. Report on the 1995-1996 Survey of Population Change. Bangkok: National Statistical Office; 1998.
7. Report on the 2005-2006 Survey of Population Change. Bangkok: Statistical Forecasting Bureau; 2007.
8. Sekar CC, Deming WE. On a method of estimating birth and death rates and the extent of registration. J Am Stat Assoc 1949;44:101-15. doi:10.2307/2280353
9. Bernillon P, Lievre L, Pillonel J, Laporte A, Costagliola D. Record-linkage between two anomymous databases for a capture-recapture estimation of underreporting of AIDS cases: France 1990-1993. Int J Epidemiol 2000;29:168-74. doi:10.1093/ije/29.1.168 PMID:10750619
10. El-shalakani M. Estimating the completeness of births and deaths registration in Egypt through dual record systems. Genus 1985;41:119-32. PMID:12267531
11. Carter KL, Williams G, Tallo V, Sanvictores D, Madera H, Riley I. Capture-recapture analysis of all-cause mortality data in Bohol, Philippines. Popul Health Metr 2011;9:9. doi:10.1186/1478-7954-9-9 PMID:21496336
12. Regehr EV, Lunn NJ, Amstrup SC, Stirling I. Supplemental materials for the analysis of capture-recapture data for polar bears in western Hudson Bay, Canada, 1984-2004. Virginia: US Department of the Interior & US Geological Survey; 2007.
13. Corrao G, Bagnardi V, Vittadini G, Favilli S. Capture-recapture methods to size alcohol related problems in a population. J Epidemiol Community Health 2000;54:603-10. doi:10.1136/jech.54.8.603 PMID:10890872
14. Laska EM. The use of capture-recapture methods in public health. Bull World Health Organ 2002;80:845. PMID:12481204
15. Farber JE. Dual system estimates of housing units based on administrative records. In: Proceedings of the Survey Research Methods Section [Internet]. Alexandria: American Statistical Association;2002. Available from: http:// www.amstat.org/sections/srms/proceedings/ [accessed 29 July 2011] .
16. Tilling K. Capture-recapture methods-useful or misleading? Int J Epidemiol 2001;30:12-4. doi:10.1093/ije/30.1.12 PMID:11171841
17. Karanth KU. Estimating tiger Panthera tigris populations from camera-trap data using capture- recapture models. Biol Conserv 1995;71:333-8. doi:10.1016/0006-3207(94)00057-W
18. Aaron DJ, Chang Y-F, Markovic N, LaPorte RE. Estimating the lesbian population: a capture-recapture approach. J Epidemiol Community Health 2003;57:207-9. doi:10.1136/jech.57.3.207 PMID:12594197
19. Eisele TP, Lindblade KA, Rosen DH, Odhiambo F, Vulule JM, Slutsker L. Evaluating the completeness of demographic surveillance of children less than five years old in western Kenya: a capture-recapture approach. Am J Trop Med Hyg 2003;69:92-7. PMID:12932105
Submitted: 27 October 2010
Revised version received: 16 July 2011
Accepted: 19 July 2011
Published online: 23 August 2011
* Correspondence to Patama Vapattanawong (e-mail: prpvp@mahidol.ac.th).