Abstract
Background
Although the general statistical advice is to keep continuous exposure variables as continuous in statistical analyses, categorisation is still a common approach in medical research. In a recent paper from the Hyperglycaemia and Adverse Pregnancy Outcome (HAPO) Study, categorisation of body mass index (BMI) was used when analysing the effect of BMI on adverse pregnancy outcomes. The lowest category, labelled "underweight", was used as the reference category.
Methods
The present paper gives a summary of reasons for categorisation and methodological drawbacks of this approach. We also discuss the choice of reference category and alternative analyses. We exemplify our arguments by a reanalysis of results from the HAPO paper.
Results
Categorisation of continuous exposure data results in loss of power and other methodological challenges. An unfortunate choice of reference category can give additional lack of precision and obscure the interpretation of risk estimates. A highlighted odds ratio (OR) in the HAPO study is the OR for birth weight >90^{th }percentile for women in the highest compared to the lowest BMI category ("obese class III" versus "underweight"). This estimate was OR = 4.55 and OR = 3.52, with two different multiple logistic regression models. When using the "normal weight" category as the reference, our corresponding estimates were OR = 2.03 and OR = 1.62, respectively. Moreover, our choice of reference category also gave narrower confidence intervals.
Summary
Due to several methodological drawbacks, categorisation should be avoided. Modern statistical analyses should be used to analyse continuous exposure data, and to explore nonlinear relations. If continuous data are categorised, special attention must be given to the choice of reference category.
Background
Although the general statistical advice is to keep continuous exposure variables as continuous in statistical analyses [1,2], categorisation of continuous variables is still a common approach in medical research [3,4]. In the present paper we revisit this issue, and we have chosen to exemplify our arguments through a discussion and reanalysis of data from the recently published paper "Hyperglycaemia and Adverse Pregnancy Outcome (HAPO) Study: associations with maternal body mass index." [5].
The main objective of the international, multicentre HAPO study was to clarify the risk of adverse outcomes associated with various degrees of glucose intolerance during pregnancy, including degrees less than overt diabetes [6]. With its impressive sample size, doubleblinding of glucose levels and thoroughly developed standardized routines for data collection, this is already a wellcited and highly valued study in the fields of obstetrics and endocrinology. In total, 25 505 women were enrolled in the study from 2000 to 2006. Glucose tolerance was measured by a 75g 2h oral glucose tolerance test (OGTT) in a large, heterogeneous, multinational, ethnically diverse cohort of women at approximately 28 weeks of gestation. The participating women, caregivers and HAPO staff were blinded to glucose tolerance values, except when predefined thresholds were met. Height and weight were also measured by a standardized procedure at the OGTT visit. Both glucose tolerance and body mass index (BMI, kg/m^{2}) were measured on continuous scales. In the study protocol it is argued for a possible continuous relation between glucose and pregnancy outcomes [6]. However, in the statistical analyses, BMI was categorised.
Our aim is to discuss arguments for categorisation of exposure variables and the corresponding methodological challenges of doing so. We also discuss issues that need attention if categorisation is the strategy of choice. Furthermore, we show how different choices of reference category can influence the results presented in the HAPO paper. Finally, we present alternative strategies to that of categorisation of continuous variables.
Regression analysis and linearity assumptions
The primary analytic tool in the HAPO paper [5] was logistic regression analysis. In standard logistic regression analyses, linearity in the logit is assumed [7]. Similarly, in linear regression, a linear relation is assumed between the exposure variable and the response variable. These assumptions give easily interpretable results for continuous exposure variables. However, a linearity check should always be carried out. If the observed data fulfil the linearity assumptions, the general advice is to keep continuous variables as continuous in such analyses [1,2]. A wide range of linearity diagnostics is available for different regression models [7]. Although we will use logistic regression analyses from the HAPO paper as our example, our arguments apply to other regression models as well. The authors of the HAPO paper have discarded a linear model with BMI as a continuous exposure variable because of a statistically significant squared term of BMI. This is not in itself a valid argument for rejection of a linear model, unless additional linearity diagnostics have also been used. In any case, presence of nonlinearity must be reflected in regression models. Several approaches are possible. The most commonly used approach, also taken in the HAPO paper, is to categorise the data [3,4].
The cost of categorisation of continuous data
Categorisation of continuous exposure data is intriguing because it seems to represent a simple solution to a methodological challenge. It is a common belief that categorical exposure data are both more robust, i.e. not easily affected by small departures from model assumptions, easier to interpret and more feasible to present than their continuous counterparts. Categorisation also mirrors clinical practice. Research results are crucial in the development of clinical medical practice where decisions typical are categorical in nature. The similarities between clinical decisionmaking and a categorical data analysis and presentation approach can make reading and understanding scientific papers with categorical data easier for clinicians. However, categorisation comes with some costs. Firstly, there is a loss of power due to reduced variability in the data. Several papers have addressed this issue and simulation studies have been used to quantify the loss [1,2]. In a study as large as the HAPO study, there is enough information to retain sufficient power even after the categorising of data. Indeed, the sample size calculations for the HAPO study were based on categorised data [6]. Paradoxically, one consequence of large sample sizes is that even small deviations from linearity can be detected, which again leads to categorisation of the originally continuous variables, resulting in loss of power and the need for larger sample sizes. Secondly, categorisation conceals information about the details of the nonlinear relation. When categorising, one assumes that the relation between the independent and the dependent variable is constant within intervals. This implies that any change in effect within an interval will be lost, also biologically plausible ones. Less obvious, but also important consequences of categorisation are those related to residual confounding [1,2], increased risk of false positive results [1,2] and issues about misclassification [8], which all may give unintentionally biased results. It has been shown that less of the confounding will be removed if one adjusts for a categorised variable, instead of the original continuous variable [1,2]. Further, choice of cutoff values can significantly impact effect measures, e.g. odd ratios, and can thereby lead to artificial associations between variables [1,2]. Finally, most exposure variables are prone to measurement errors. When categorising a continuous variable, nondifferential errors can lead to differential misclassification [8].
Categorisation in the HAPO paper
Assuming that the choice of categorising data is made despite of the shortcomings of such a strategy, cutoff values and reference category must be chosen with care. It is common to use internationally and clinically accepted categories, if such exist. The HAPO paper [5] focuses on BMI as the exposure variable for several pregnancy outcomes. The international classification of BMI (kg/m^{2}) in adults suggested by WHO [9] is "underweight" (<18.5), "normal weight" (18.524.9), "overweight" (25.029.9) and "obese class IIII" (30.034.9, 35.039.9 and ≥40, respectively). In the HAPO study, BMI was measured only at gestational weeks 2432 (at the OGTT). Due to the expected weight gain during the first 2432 weeks of the pregnancy the WHO categories did not readily apply. Instead, "comparable category limits for BMI" at 28 weeks of gestation were estimated by a regression analysis of BMI measured at the OGTT on recalled prepregnancy BMI and gestational age at the OGTT: <22.6, 22.628.4, 28.532.9, 33.037.4, 37.441.9 and ≥42.0 kg/m^{2}.
This usage of regression modelling deserves a remark. A regression model consists of two main elements; a deterministic term, linking the exposure to the outcome, and an error term which incorporates the residual variance of the observations, which is not expressed by the link function. This implies that one cannot merely transform the midpregnancy BMI to a predicted prepregnancy BMI, based on the regression estimates, as part of the variance has been discarded when the prepregnant BMI was estimated. Additionally, it is not easy to deduct how this could influence misclassification issues.
In statistical analysis of categorical data, the reference category is the category against which the others are compared [7]. Both the biological interpretation of the estimated association and the number of observations should be considered when choosing the reference category. In the HAPO paper, the estimated "normal range" BMI category at gestational week 2432 was 22.628.4 kg/m^{2}. With its almost 12 000 (51%) women it was also the largest group, making it the preferred reference category both by means of a natural reference and by considerations of statistical power. However, the lowest BMI group, labelled "underweight", with 2989 (13%) women was chosen as the reference. This resulted in broader confidence intervals, and somewhat misleading effect estimates. Clinically, the underweight group is susceptible to having additional health problems, like eating disorders or chronic diseases influencing body weight. Thus, there are several reasons why using the "underweight" group as the reference could influence the estimated associations in the HAPO study.
The HAPO paper estimated associations between BMI and 11 different outcomes. The results were presented as frequencies and adjusted effect estimates (odds ratios, ORs, and 95% confidence intervals, CIs) from two different multiple logistic regression models. The crude ORs were not given, but can easily be calculated from the frequencies in the tables. We have focused on the macrosomia (high birth weight) results, and Table 1 shows results from logistic regression analyses of the effect of BMI on macrosomia, with two different choices of reference category. We first considered the crude results with reference categories BMI <22.6 kg/m^{2 }or BMI 22.628.4 kg/m^{2}. Although the same increasing trend is seen, the two analyses give different impressions to the reader. By using "underweight" as the reference group, the increasing trend is striking, with a nearly fivefold increased risk for macrosomia in the highest BMI group ("obese class III"). When "normal range" is used as the reference group, the similarities of the three highest BMI groups are clearly seen, in addition to the statistically significantly lower OR for the "underweight" group. The lack of increase in ORs for the three highest BMI groups is even more convincing for the adjusted ORs. In addition, the table shows the standard errors of the betas in the logistic regressions and the widths of the unadjusted 95% CIs, illustrating the loss in statistical power by using a small group as the reference, compared to using the larger normal range group as the reference. It should be noted that using SE for the betas instead of CI for the ORs as measures of precision of the estimates show much less exaggeration of lack of precision when changing reference category. However, most scientific papers that present results from logistic regression focus on ORs rather than betas. We have therefore chosen to present both, also the widths of the CIs. By changing the reference category, the main conclusions in the paper would not be altered, but the effect sizes would be relative to pregnant women with a normal range BMI, as well as being more precisely estimated.
Table 1. The effect of BMI on macrosomia: Different reference categories influence effect estimates and precision.
Alternatives to categorisation
Application of more complex modern statistical techniques is necessary if we want to keep variables continuous in situations where nonlinearity is evident. Suitable statistical methods include piecewise linear regression [10] and spline based techniques [11], for example generalised additive models (GAMs). By applying such techniques, one will gain statistical power, avoid unintended methodological challenges based on classification problems, as well as allowing the data to provide more explicit details about the nonlinear relations. This could include apparent changepoints or thresholds values where a change in the outcome variable can be observed, piecewise linear relations, "plateau" effects and so on. The existence and locations of changepoints can further be formally tested for and estimated [10]. In the HAPO paper, a plateau effect for the higher BMI categories is discussed. Both clinical practice and physiological theories would benefit from a statistical test for the existence of a plateau effect, and an estimate of the BMI value where the risk levelled off. Application of more complex analytical tools will to some extent complicate the technical analysis and the task of applying and communicating the research results. However, this should not prevent a more optimal analysis.
Summary
Applying and communicating results from statistical analyses is a challenging aspect of the scientific process. The results should be presented in a simple, yet balanced manner, and should be based on adequate statistical methods. Researchers should be aware of the tradeoff between simplicity in presentation and unintended methodological challenges arising from simplifications. Particularly, if one chooses to categorise continuous exposure variables, special attention must be given both to the definition of categories and the choice of reference category. If we apply modern statistical analysis techniques, we might extract more detailed knowledge and get a better utilisation of our data.
Abbreviations
HAPO Study: Hyperglycaemia and Adverse Pregnancy Outcome Study; BMI: Body Mass Index defined as weight/height^{2}, where weight is measured in kilograms (kg) and height is measured in meters (m); OGTT: Oral glucose tolerance test; GAMs: Generalised additive models.
Competing interests
There were no competing interests, and no personal relationships, academic competition, or intellectual commitments that might bias the work or interfere with objective judgment.
Authors' contributions
KFF did the statistical analyses. All authors contributed to the discussions about the topic, revision the manuscript and to the final approval of the manuscript.
Acknowledgements
The authors KFF and EQ were funded by the National Resource Centre for Women's Health. The Centre funded a PhD position in biostatistics for KFF and a Post Doc position in gestational endocrinology for EQ.
References

Altman DG, Royston P: Statistics notes  The cost of dichotomising continuous variables.

van Walraven C, Hart RG: Leave 'em alone  Why continuous variables should be analyzed as such.

Pocock SJ, Collier TJ, Dandreo KJ, De Stavola BL, Goldman MB, Kalish LA, et al.: Issues in the reporting of epidemiological studies: a survey of recent practice.

Turner E, Dobson J, Pocock S: Categorisation of continuous risk factors in epidemiological publications: a survey of current practice.

HAPO Study Cooperative Research Group: Hyperglycaemia and Adverse Pregnancy Outcome (HAPO) Study: associations with maternal body mass index.

Contreras M, Sacks DA, Watson W, Dooley SL, Foderaro M, Niznik C, Polk DH, Bjaloncik J, Catalano PM, Fox S, Gullion L, Johnson C, Lindsay CA, Makovos H, Saker F, Carpenter MW, Somers MH, Amankwah KS, Chan PC, Kapur B, Kenshole A, Lawrence G, Matheson K, Mayes L, Owen H, Basdeo P, Cave C, Fenty G, Gibson E, Hennis A, Rotchell YE, Spooner C, Thomas HAR, Fox J, Hadden DR, Halliday H, Irwin J, Kearney O, McCance DR, Smye M, Traub AI, Cruickshank JK, Dry J, Holt AC, Khan F, Lambert C, Maresh M, Prichard F, van Haeften TW, van de Hengel AMR, Visser GHA, Zwart A, Chaovarindr U, Chotigeat U, Deerochanawong C, Panyasiri I, Sanguanpong P, Amichay D, Golan A, Marks K, Mazor M, Ronen J, Wiznitzer A, Chen R, Faktorovich A, Harel D, Hoter N, Pardo J, Bowling F, Cowley D, Liley HG, McIntyre HD, Morrison B, Tudehope D, Kong SL, Li CY, Ng KF, Ng PC, Rogers MS, Beverly J, Edgar M, Giles W, Gill A, Lowe J, Verma J, Koh A, Tan E, Teo CC, Rajadurai V, Wee HY, Yeo GSH, Coustan D, Haydon B, Alexander A, Hadden DR, AttiasRaved O, Hod M, Oats JJN, Parry AF, Frank AS, Lowe LP, Metzger BE, Thomas A, Case T, Cholod P, Dyer AR, Engelman L, Xiao M, Burgess CI, Lappin TRJ, Nesbitt GS, Trimble ER, Coustan D, Dyer AR, Hadden DR, Hod M, Metzger BE, Lowe LP, Oats JJN, Persson B, Trimble ER, Cutter GR, Gabbe SG, Hare JW, Wagenknecht LE: The Hyperglycemia and Adverse Pregnancy Outcome (HAPO) Study.
International Journal of Gynecology & Obstetrics 2002, 78:6977.

Hosmer DW, Lemeshow S: Applied logistic regression. 2nd edition. John Wiley & Sons. U.S.A; 2000.

Flegal KM, Keyl PM, Nieto FJ: Differential Misclassification Arising from Nondifferential Errors in Exposure Measurement.

World Health Organization: Obesity: Preventing and managing the global epidemic  Introduction. Geneva: World Health Organization; 2000.

Muggeo VMR: Estimating regression models with unknown breakpoints.

Gilboa SM, Correa A, Alverson CJ: Use of spline regression in an analysis of maternal prepregnancy body mass index and adverse birth outcomes: Does it tell us more than we already know?
Prepublication history
The prepublication history for this paper can be accessed here: