Email updates

Keep up to date with the latest news and content from BMC Pregnancy and Childbirth and BioMed Central.

Open Access Research article

Risk adjustment models for interhospital comparison of CS rates using Robson’s ten group classification system and other socio-demographic and clinical variables

Paola Colais1*, Maria P Fantini2, Danilo Fusco1, Elisa Carretta2, Elisa Stivanello2, Jacopo Lenzi2, Giulia Pieri2 and Carlo A Perucci3

Author Affiliations

1 Department of Epidemiology, Regional Health Service, Lazio Region, Italy

2 Department of Medicine and Public Health, University of Bologna, Bologna, Italy

3 National Agency of Regional Health Services, Italy

For all author emails, please log on.

BMC Pregnancy and Childbirth 2012, 12:54  doi:10.1186/1471-2393-12-54

The electronic version of this article is the complete one and can be found online at:

Received:1 April 2011
Accepted:23 May 2012
Published:21 June 2012

© 2012 Colais et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



Caesarean section (CS) rate is a quality of health care indicator frequently used at national and international level. The aim of this study was to assess whether adjustment for Robson’s Ten Group Classification System (TGCS), and clinical and socio-demographic variables of the mother and the fetus is necessary for inter-hospital comparisons of CS rates.


The study population includes 64,423 deliveries in Emilia-Romagna between January 1, 2003 and December 31, 2004, classified according to theTGCS. Poisson regression was used to estimate crude and adjusted hospital relative risks of CS compared to a reference category. Analyses were carried out in the overall population and separately according to the Robson groups (groups I, II, III, IV and V–X combined). Adjusted relative risks (RR) of CS were estimated using two risk-adjustment models; the first (M1) including the TGCS group as the only adjustment factor; the second (M2) including in addition demographic and clinical confounders identified using a stepwise selection procedure. Percentage variations between crude and adjusted RRs by hospital were calculated to evaluate the confounding effect of covariates.


The percentage variations from crude to adjusted RR proved to be similar in M1 and M2 model. However, stratified analyses by Robson’s classification groups showed that residual confounding for clinical and demographic variables was present in groups I (nulliparous, single, cephalic, ≥37 weeks, spontaneous labour) and III (multiparous, excluding previous CS, single, cephalic, ≥37 weeks, spontaneous labour) and IV (multiparous, excluding previous CS, single, cephalic, ≥37 weeks, induced or CS before labour) and to a minor extent in groups II (nulliparous, single, cephalic, ≥37 weeks, induced or CS before labour) and IV (multiparous, excluding previous CS, single, cephalic, ≥37 weeks, induced or CS before labour).


The TGCS classification is useful for inter-hospital comparison of CS section rates, but residual confounding is present in the TGCS strata.


Caesarean section (CS) rate is one of the most frequently used indicators of health care quality at the national and international level for clinical governance and outcome research. Hospitals and health-care systems are often compared on the basis of this indicator with the implicit assumption that lower rates reflect more appropriate practice, although the rate that defines optimum quality of care is undefined and seems to depend on the characteristics of the populations under study [1]. The World Health Organization has indicated that a CS rate greater than 10–15% is not justified in any region of the world. Rates are higher in developed countries, Latin America and the Caribbean, and lower in other developing countries [2-6].

In 2005, the Italian CS rate was the highest in Europe (38.5% vs. an average European rate of 23.7%) and one of the highest in the world [6]. In Italy, national CS rates have increased from 32% in 2001 to 38.5% in 2005. This increase was found both for primary CS and repeated CS. Primary caesarean deliveries contribute 2/3 to the overall CS rate, although the contribution of repeated CS is higher in regions with high overall CS rates [7]. Primary caesarean deliveries are an important target for reduction, because they lead to an increased risk for repeated caesarean delivery [8-10]. Therefore, some authors suggested to focus on primary CS for inter-hospital comparison and quality improvement [11], and others, based on evidence suggesting that non-vertex and multiple births may have better outcomes with cesarean deliveries [12], omitted these categories from the calculation of CS rates and focused on nulliparous term cephalic singleton (NTCS) deliveries.

In 2001, a new classification for CS known as the “ten group” (TGCS) or Robson classification was proposed. This classification system categorizes women into 10 mutually exclusive groups, considering the following a priori criteria: parity, the previous obstetric record of the woman, the course of labour, including pre-labour CS, and gestational age [11]. Several studies used the Robson classification system to compare CS rates within specific subsets of an obstetric population to overcome many of the historic controversies that have arisen when comparing overall caesarean rates among different populations [13-17]. A recent systematic review [18] supported the use of TGCS classification over other classifications based on characteristics of parturients for auditing purposes and comparison of CS rates across different settings.

The TGCS classifies CS according to the characteristics of each woman and her pregnancy. However, caesarean delivery has many other indications such as fetal distress, dystocia, placenta previa, HIV, and other conditions of the mother and foetus [19]. The failure to account for such patient-specific risk factors may lead to biased interhospital comparisons [20,21]. The possible confounding effect is caused by the heterogeneous distribution of CS risk factors across hospitals, that is not taken into account in the a priori Robson classification.

A recent study addressing interhospital comparison of CS rates in women with a primary CS and those with NTCS deliveries emphasized a differential need to adjust for clinical variables of the mother and the foetus and for socio-demographic characteristics of the mother [22]. Adjustment proved to be warranted when the indicator of interest was primary CS rate but not when the indicator was NTCS CS rate. The aim of the present study is to define a risk-adjustment model for interhospital comparison using TGCS classification and clinical and socio-demographic characteristics of the mother and the fetus not included in the TGCS classification that are indications for CS.


Study population

Deliveries in Emilia-Romagna Region (Italy) from January 1, 2003 through December 31, 2004 were extracted from Hospital Discharge Abstracts of mothers and their babies (SDO), and Birth Certificates (CedaP). Record linkage was performed between the SDO and CedaP databases.

The SDO data includes demographics (ID number, sex, date and place of birth, place of residence), discharge ID, admission and discharge dates, discharge diagnoses and procedures (International Classification of Diseases, 9th Revision, Clinical Modification ICD-IX-CM), ward(s) of hospitalization, date(s) of in-hospital transfer, and the regional code of the admitting facility.

Birth Certificates include demographic data of the mother, information on presentation and multiple pregnancy (singleton cephalic, singleton breech, trasverse or oblique lie, etc.), parity (nulliparous, multiparous), the course of labour and delivery (spontaneous labour, induced labour or CS before labour) and gestational age (defined as the number of completed weeks at the time of birth).

To identify the delivery by Hospital Discharge Abstracts of mothers and their babies (SDO), we used three different sources of information: DRGs 370–375, procedure codes ICD-9-CM 72.*–74.* and diagnosis codes ICD-9-CM 640.*–676.*, V27.

DRG 370, procedure codes ICD-9-CM 74.0, 74.1, 74.2, 74.4, 74, 99 and diagnosis codes ICD-9-CM 669.7 were used to identify caesarean deliveries.

Mothers under the age of 11 or over the age of 50 years; mothers who were discharged from a hospital without an operating room; and infants with a birth weight under 550 g or over 6000 g were excluded.

CS rates were calculated as the ratio of caesarean deliveries to total deliveries.

Deliveries were retrospectively classified according to Robson’s Ten Group Classification System (TGCS) using information in the databases

The following socio-demographic variables were considered as potential risk factors for caesarean sections: maternal age (classified as <20, 20–34, or ≥35 years), citizenship (Italian, from developing countries, from developed countries other than Italy) and educational level (≤ 5, 6–8, 9–13, or ≥14 years). Maternal and neonatal clinical factors that constitute indication for CS were extracted using primary and secondary discharge diagnoses of Hospital Discharge Abstracts (see Additional file 1: Appendix A for the ICD-9-CM codes).

Additional file 1. Appendix A. A) ICD 9-CM codes to identify clinical variables from mothers’ discharge records. B) ICD 9-CM codes to identify variables from neonatal discharge records. C) ) ICD 9-CM codes to identify variables from both neonatal and maternal discharge records.

Format: DOC Size: 47KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

The study was carried out in compliance with the Italian law on privacy (Art. 20–21, DL 196/2003) and the regulations of the Regional Health Authority of Emilia-Romagna on data management. Access to the data was approved by the hospital trust administration. Data were anonymized at the regional statistical office where a unique identifier, the same for all databases, was assigned to each patient. This identifier does not allow to trace the patient’s identity and other sensitive data. When anonymized administrative data are used to inform health care planning activities, the study is exempt from notification to the Ethics Committee and no specific written consent is needed to use patient information stored in the hospital databases.

Statistical analysis

Analyses were initially carried out on the entire population. We then analyzed the I–IV Robson groups separately, and the V–X Robson groups taken together. The last six Robson groups accounted for about 20% of all deliveries and could not be analyzed separately because of the small number of deliveries in these groups by hospital. Crude and adjusted relative risks of CS for each hospital were calculated using as the reference category hospitals with the lowest CS rates. These hospitals were identified by means of a recursive procedure developed by P.Re.Val.E. Project [23]. Adjusted RR of CS (caesarean section risk for patients admitted to a specific hospital vs. caesarean section risk for patients admitted to the reference category) were obtained by using modified Poisson regression models based on the Huber sandwich estimate, that improves efficiency in mean-variance relationship.

Specifically, two risk adjustment models were set up. The first model (M1) was built using TGCS as the only potential confounding factor. The second model (M2) included, in addition to TGCS, a number of potential confounders (demographic and clinical variables related to the mother and foetus) selected according to available scientific evidence. These include age, citizenship, severe co-morbid illness of the mother, diabetes, hypertension, HIV, lung diseases, ante-partum haemorrhage/abruption placentae/placenta praevia, eclampsia/pre-eclampsia, foetal-pelvic disproportion/excessive development of the infant, polyhydramnios, oligohydramnios, isoimmunisation, premature rupture of the membranes, abortion threats/assisted fecundation, congenital malformation, problem of the amnios, post-maturity and macrosomia, intrauterine growth retardation (see Additional file 1 for the ICD-9-CM codes). A stepwise selection procedure (significance level for entry of 0.10 and 0.05 for stay), was used to remove variables unrelated to CS.

Stratum-specific models were defined for groups I to IV and V-X that included only clinical and demographic variables selected using a stepwise selection procedure.

Adjusted relative risks (RRs) and percentage variations between crude and adjusted RRs by hospital were then calculated to evaluate the amount of confounding. The presence of confounding was defined as a percentage variation ≥10% between crude and adjusted RRs [24,25]. The same hospitals selected as reference group in the overall population were used as the reference group in the stratified analyses. The significance level for the RR was set at 5% (p < 0.05). All analyses were performed using SAS Version 8.02.


A total of 64,423 deliveries in Emilia-Romagna occurred between January 1, 2003, and December 31, 2004. The overall crude CS rate was 30.4%, and the CS rate in the reference group was 23.1%. Figure 1 shows the TGCS distribution in the study population and the proportion of CS in each group.

thumbnailFigure 1. Percentage distribution of deliveries by TGCS over the total number of deliveries. Light gray: non-CS; gray: CS.

The first Robson group was the most frequent (28.0% of deliveries), while group VIII was the less frequent (1.1%). The first four groups constitute approximately the 80% of the entire population. Group V had the highest CS rate (93.6%), while group III showed the lowest recourse to CS, with a rate of 6.3%. Table 1 reports the number of caesarean deliveries, crude and adjusted caesarean section RRs, their statistical significance, and the percentage variation by hospital. The RR percentage variations by hospital estimated using the M1 model were very similar with those estimated using M2 model. The M1 adjusted RRs led to percentage variations greater than 10% in 16 out of 25 hospitals, while the M2 adjusted RRs led to percentage variations greater than 10% in 15 out of 25 hospitals.

Table 1. Crude and adjusted c-section RRs of hospitals

Tables 2, 3, 4, 5 provide the number of caesarean deliveries, crude and adjusted caesarean section RRs, their statistical significance, and the percentage variation by hospital in groups I, II, III, IV. In group I (Table 2), 10 hospitals had a percentage variation greater than 10%, with a very high reduction in adjusted compared to the crude RR for hospital M (36.9%). In group II no percentage variation ≥10% was observed. In group III (Table 4), 12 hospitals had a percentage variation in relative risk higher than 10 %, and again hospital M proved to have the largest value (36.6%). In group IV 7 hospitals exhibited variations ≥10 % (Table 5). Lastly, Table 6 reports the number of caesarean deliveries, crude and adjusted caesarean section RRs, their statistical significance, and the percentage variation by hospital for Robson groups V–X. We found that only 2 hospitals (L, T) had a percentage variation greater than 10%.

Table 2. Crude and adjusted c-section RRs of hospitals in group I

Table 3. Crude and adjusted c-section RRs of hospitals in group II

Table 4. Crude and adjusted c-section RRs of hospitals in group III

Table 5. Crude and adjusted c-section RRs of hospitals in group IV

Table 6. Crude and adjusted c-section RRs of hospitals in groups V–X


The aim of the present study was to define a risk-adjustment model for inter-hospital comparison, using TGCS classification and variables not included in the TGCS classification that are indications for CS. Our results indicate that TGCS classification should be used to control for the hospital case mix in terms of parity, presentation, gestational age and multiple pregnancy. However, some residual variability, in the overall population and in the I and III groups, is accounted for by clinical and socio-demographic confounders.

Several studies used different risk adjustment [26-28] techniques to compare CS rates across hospitals. The stratified analyses proposed by Robson with “a priori” criteria of classification seem to overcome the problem of risk adjustment. However, our results suggest that the TGCS classification is not sufficient to remove case mix differences present in the first four TGCS groups, especially groups I and III. Brennan et al. [29], recently demonstrated a wide variation of caesarean section rates in women in spontaneous cephalic term labour (I and III TGCS groups) among 9 international “third level” hospitals and suggested the need to verify the role of other potential confounding factors when comparative evaluation is carried out. Our results incorporate in the analyses some clinical and socio-demographic CS risk factors related to mother and foetus, not included in the TGCS, but do not consider the organizational and process variables such as midwifery care, use of oxytocin to correct dystocia, intrapartum foetal blood sampling mentioned by Brennan et al. [29].

One of the strength of the present study is the opportunity to use two current administrative databases with a very good record linkage (higher then 95%) and to take advantage of data collected from two different sources. In this study, caesarean section occurrences were evaluated using discharge record data. Accuracy, completeness, and quality of records may differ from hospital to hospital, however the CS rates and proportion of all patients in Robson’s groups (data not shown) are similar to those reported in other studies [30,31]. The potential for inconsistencies in coding discharge records may challenge the accuracy of the assessment of outcome and of risk factors. Missing information on important risk factors and errors in coding may in fact lead to subsequent errors in adjustment and this represent a limit of the study. Finally, part of the limitations of administrative data may be due to the basic tension which exists between using the same data for reimbursement and for measuring quality. “When the use is reimbursement, there is a tendency to perform coding quickly and to maximize the coding of complications and co-morbidities. When the use is to assess quality, however, it is important for coders to have a complete record and to restrict diagnosis coding to conditions that affect patient care [32].” For instance, hypertension and diabetes may intervene in the algorithm used to determine the case mix of an admission and thus be rewarding in financial terms, whereas this may not be the case for labour induction and history of a previous caesarean. Nevertheless, administrative databases are widely available at the national and regional level, and are currently utilized to compare outcomes, including CS, of inpatient care in Italy [33]. Risk adjustment models should be time and population specific, and have recently proved to be useful for monitoring caesarean section rates and for interhospital comparison [26]. Methods used to develop models based on administrative information can potentially be generalized to other populations.

The TGCS is a good tool for clinical audit practices as it enables professionals to compare their CS practice with homogeneous a priori risk populations. Since low-risk deliveries both in nulliparous and multiparous women are an important target for reduction because in these categories the large majority of inappropriate CS can be found, our analytical method may be useful to partial out the effect of clinical and demographic variables in the TGCS groups. Furthermore, it is important to focus on the first four TGCS groups that in our country represent, given the low fertility rate (1.28 for years 2000–2005) [33] and current obstetric practice in relation to the management of deliveries, about two thirds of all deliveries [34] because in these groups it is most likely to find inappropriate CS.

Reducing the number of unnecessary CS in low-risk women is also a good strategy to indirectly reduce the CS in women with previous CS.

In conclusion, our results indicate that parity and type of labour should be taken into account in risk adjustment models for interhospital comparison. Moreover, Robson’s classification proved to be useful to compare caesarean rates among hospitals even though the presence of residual confounding related to clinical and socio-demographic variables within strata may lead to potential bias, especially in low-risk nulliparous and multiparous women with spontaneous labour (groups I and III). Only after eliminating confounders in comparative evaluation of hospital performance we may be confident that we are considering unnecessary variability and inappropriate procedures. Unnecessary variability must be the target for health care quality improvement activities.

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

CP: study conception and design, data analysis, interpretation of data. FMP: study conception and design, interpretation of data. FD: study design, data analysis, interpretation of data. CE: data analysis. SE: data analysis. LJ: statistical analysis. PG: interpretation of data. PCA: study conception and design, interpretation of data. All authors read and approved the final manuscript.


  1. Canadian Institute for Health Information: Giving Birth in Canada: Providers of Maternity and Infart Care. Ottawa. Canada: CIHI; 2004. OpenURL

  2. MacDorman M, Menacker F, Declercq E: Cesarean birth in the United States: epidemiology, trends, and outcomes.

    Clin Perinatol 2008, 35:293-307. PubMed Abstract | Publisher Full Text OpenURL

  3. Belizan J, Althabe F, Cafferata ML: Health consequences of the increasing caesarean section rates.

    Epidemiology 2007, 18:485-486. PubMed Abstract | Publisher Full Text OpenURL

  4. Betran AP, Merialdi M, Lauer A, et al.: Rates of caesarean section: analysis of global, regional and national estimates.

    Paediatr Perinat Epidemiol 2007, 21:98-113. PubMed Abstract | Publisher Full Text OpenURL

  5. Belizan JM, Althabe F, Barros S, et al.: Rates and implications of caesarean section in Latin America: ecological study.

    BMJ 1999, 519:1397-1402. OpenURL

  6. World Health Organization: Appropriate technology for birth.

    Lancet 1985, 2:436-437. PubMed Abstract | Publisher Full Text OpenURL

  7. Osservatorio Nazionale sulla Salute nelle Regioni Italiane: Rapporto Osservasalute 2008. In Stato di salute e qualità dell’assistenza nelle regioni italiane. Milano: Press; 2008. OpenURL

  8. Royal College of Obstetricians and Gynaecologists: Birth after previous caesarean birth. London (UK): RCOG Press; 2007.

    Green-top guideline; n°45


  9. Robson MS: Can we reduce the caesarean section rate?

    Best Pract Res Clin Obstet Gynaecol 2001, 15:179-194. PubMed Abstract | Publisher Full Text OpenURL

  10. Meikle SF, Steiner CA, Zhang J, Lawrence WL: A national estimate of the elective primary caesarean delivery rate.

    Obstet Gynecol 2005, 105(4):751-756. PubMed Abstract | Publisher Full Text OpenURL

  11. Main EK, Bloomfield L, Hunt G: Sutter Health, First Pregnancy and delivery Clinical Initiative Committee: Development of large-scale obstetric quality-improvement program that focused on the nulliparous patient at term.

    Am J Obstet Gynecol 2004, 190:1747-56. PubMed Abstract | Publisher Full Text OpenURL

  12. Kontopoulos EV, Ananth CV, Smulian JC, Vintzileos AM: The impact of route of delivery and presentation on twin neonatal and infant mortality: a population-based study in the USA, 1995–97.

    J Matern Fetal Neonatal Med 2004, 15(4):219-24. PubMed Abstract | Publisher Full Text OpenURL

  13. Tamburrano MT, Iaccarino C, Castagnaro C: Un primo tentativo di analisi dei parti cesarei attraverso le classi di Robson: limiti e potenzialità dei CeDAP.

    Rapporto Osservasalute 2007, 1:244-250. OpenURL

  14. Regione Emilia-Romagna: La nascita in Emilia-Romagna. 5° Rapporto sui dati del Certificato di Assistenza al Parto (CedAP) - Anno 2007. Bologna: Regione Emilia-Romagna; 2008:77-107. OpenURL

  15. Fischer A, LaCoursiere DY, Barnard P, et al.: Differences Between Hospitals in Cesarean Rates for Term Primigravidas With Cephalic Presentation.

    Obstet Gynecol 2005, 105(4):816-821. PubMed Abstract | Publisher Full Text OpenURL

  16. Robson M: Classification of cesarean sections.

    Fetal Matern Med Rev 2001, 12:23-9. OpenURL

  17. Betrán AP, Gulmezoglu AM, Robson M, et al.: WHO global survey on maternal and perinatal health in Latin America: classifying caesarean sections.

    Reprod Health 2009, 6:18. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Torloni MR, Betran AP, Souza JP, et al.: Classification for Cesarean Section: A Systemac Review.

    PLoS One 2011, 6(1):e14566. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Gregory KD, Korst LM, Lawrence DP: Variation in elective primary caesarean delivery by patient and hospital factors.

    Am J Obstet Gynecol 2001, 184:1521-1531. PubMed Abstract | Publisher Full Text OpenURL

  20. Royal College of Obstetricians and Gynaecologists: National Collaborating Centre for Women’s and Children’s Health Commissioned by the National Institute for Clinical Excellence. Cesarean Section, clinical guideline. London: RCOG Press; 2004. PubMed Abstract OpenURL

  21. Shearer El: Caesarean section: medical benefits and costs.

    Soc Sci Med 1993, 37(10):1223-31.


    PubMed Abstract | Publisher Full Text OpenURL

  22. Stivanello E, Rucci P, Carretta E, et al.: Risk Adjustment for Inter-Hospital Comparison of Caesarean Delivery Rates in Low-risk Deliveries.

    PLoS One 2011, 6(11):e28060.


    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Lazio Regional Health Service (RHS):

    Outcome Evaluation Program 2006–2008. 2011.

    Available at: webcite. Accessed March 16, 2011

    PubMed Abstract | Publisher Full Text OpenURL

  24. Rothman KJ: Epidemiology. An Introduction. Cary, NC: Oxford University Press; 2002. OpenURL

  25. Grayson DA: Confounding confounding.

    Am J Epidemiol 1987, 126:546-53. PubMed Abstract | Publisher Full Text OpenURL

  26. Iezzoni LI: Risk adjustment for measuring healthcare outcomes. 2nd edition. Chicago, IL: Health Administration Press; 1997. OpenURL

  27. Arcà M, Fusco D, Barone AP, Perucci CA: Risk adjustment and outcome research. Part I.

    J Cardiovasc Med 2006, 7(9):682-90. Publisher Full Text OpenURL

  28. Fantini MP, Stivanello E, Frammartino B, et al.: Risk adjustment for inter-hospital comparison of primary cesarean section rates: need, validity and parsimony.

    BMC Health Serv Res 2006, 6:100. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Brennan DJ, Robson MS, Murphy M, et al.: Comparative analysis of international cesarean delivery rates using 10-group classification identifies significant variation in spontaneous labor.

    Am J Obstet Gynecol 2009, 201(3):308.e1-8. Publisher Full Text OpenURL

  30. McCarthy FP, Rigg L, Cady L, Cullinane F: A new way of looking at caesarean section births.

    Aust N Z J Obstet Gynaecol 2007, 47:316-320. PubMed Abstract | Publisher Full Text OpenURL

  31. Scarella A, Chamy V, Sepúlveda M, et al.: Medical audit using the Ten Group Classification System and its impact on the cesarean section rate.

    Eur J Obstet Gynecol Reprod Biol 2011, 154(2):136-40. PubMed Abstract | Publisher Full Text OpenURL

  32. AHRQ:

    AHRQuality Indicators e-Newsletter: AHRQ-QI-TIPS: ICD-9-CM Coding Issues. 2005.

    [http://www.qualityindica ter.htm webcite


  33. Lazio Regional Health Service (RHS):

    Outcome Evaluation Program 2006–2008. 2011.

    Available at: webcite. Accessed March 16, 2011

    PubMed Abstract | Publisher Full Text OpenURL

  34. United Nations Department of Economic and Social Affairs, Population Division: “United Nations World Population Prospects: 2006 revision, Table A.15”. New York: UN; 2007.

    Retrieved 7 December 2009


Pre-publication history

The pre-publication history for this paper can be accessed here: