|
|
||||||||
RESEARCH REPORT |
1 Biostatistics Unit, Academic Unit of Epidemiology and Health Services Research, Medical School, University of Leeds, 24 Hyde Terrace, Leeds, LS2 9LN, UK;
2 Dental Training School Malaysia, No 3 Sepoy Lines Road, 10450 Penang, Malaysia;
3 Department of Periodontology &
4 Department of Paediatric Dentistry, Eastman Dental Institute, University College London, 256 Grays Inn Road, London, UK; and
5 Department of Oral and Maxillofacial Medicine and Pathology, Guys Kings and St Thomas Schools of Medicine, Dentistry and Biomedical Sciences, Caldecot Rd, London, UK;
*corresponding author, m.s.gilthorpe{at}leeds.ac.uk
| ABSTRACT |
|---|
|
|
|---|
KEY WORDS: multilevel modeling hierarchical linear modeling random coefficients lifetime cumulative attachment loss pocket probing depth
| INTRODUCTION |
|---|
|
|
|---|
Aggregation of information to the subject level has led to the impression that progression of periodontal diseases is a continuous phenomenon, occurring at a relatively constant pace, only varying between population groups (Löe et al., 1978). However, frequent examinations of individual sites, along with stringent criteria for what constitutes disease rather than measurement error, led Socransky et al. (1984) to postulate that periodontal disease progressed through bursts of activity, accompanied by remission, repair, or further activity. This gave rise to two competing theories: the "linear" (continuous-rate) model, where, overall, sites slowly and progressively lose attachment; and the (random) "burst" model, where multiple sites show breakdown within a short period, with periods of remission that might last for months, years, or decades.
It has been argued that measurement error, and the possibility of different disease patterns at different sites, could produce erroneous evidence of burst progression (Ralls and Cohen, 1986; Cohen and Ralls, 1988). Yang et al. (1991) advocated regression methods of analysis, though Allen and Hausmann (1995), using simulated measurements, found model-fitting assessed by the least-squares criterion unreliable. Hujoel and Leroux (1998) found that lack-of-fit methods were unable to resolve the dilemma and that burst sizes of 3-5 mm would be necessary to be reliably distinguished from linear progressionan unlikely clinical scenario. For the "burst" model, statistical methods have utilized the distribution of the sizes of bursts over time as a normal random variable, or as a uniform random variable (Yang et al., 1993). The fit of these models to longitudinal data was assessed (Sterne et al., 1990, 1992; Jeffcoat and Reddy, 1991; Machtei et al., 1993), but findings were not conclusive.
Several studies have therefore attempted to model periodontal disease progression, though none has satisfactorily overcome the complexity of the data structure. This study proposes the use of multilevel modeling (MLM) to exploit periodontal data hierarchy, delivering a comprehensive model that describes the underlying progression of periodontal disease. Lifetime cumulative attachment loss (LCAL) (Eaton et al., 2001; Griffiths et al., 2001) and pocket probing depth (PPD) are modeled over time, and these models are further developed to consider potential periodontal disease risk factors and confounders.
| MATERIALS & METHODS |
|---|
|
|
|---|
Model Considerations
LCAL and PPD, measured in millimeters, were treated as positive continuous outcomes, since their underlying distributions were sufficiently close to Normal for standard assumptions to be adopted. Two "trajectory" models were sought to describe the values of LCAL and PPD, respectively. A trajectory in this instance was a hypothetical "path" of outcomes across all measurement occasions during the study. Any parameterization of change would suffice, provided that it facilitated the description of non-linear transitions over time and provided an evaluation of how change (linear progression/recession) related to acceleration/deceleration (non-linear progression/recession). With only three measurement occasions, trajectories were described by a quadratic expression of time through two time-dependent covariates: linear time (t) and quadratic time (t2), where t was the duration between initial and subsequent measurements, centered on the second occasion for improved model estimation (Gilthorpe and Cunningham, 2000).
Periodontal "Risk Factors" and Confounders
Initially, multilevel models were developed without any consideration of factors that might potentially affect either outcome. Clinical factors recorded at the time of each inspection reflected oral hygiene and oral health status, and were considered as potential "risk factors" and/or confounders. Calculus, plaque, and bleeding were coded as absent (0) or present (1) for each site. Aggregate variables were constructed, such as the proportion of sites with plaque present within each subject. The proportions of sites with LCAL > 1 mm and PPD > 2 mm were also derived. Subject-based proportions were centered. Subjects were coded as smokers (1) or non-smokers (0) and drinkers/non-drinkers of alcoholic beverages. Tooth type was coded such that central incisors were the reference against which all other teeth were compared.
Some covariates were correlated, thereby giving rise to collinearity, potentially yielding unreliable model coefficients and elevated standard errors (Slinker and Glantz, 1985). Each covariate was modeled separately before being combined in one model. Contrasts between coefficients of these models yielded insight into how robust model estimates were in the presence of multicollinearity. Non-significant covariates were not removed from any model, since variable subset selection is not straightforward (Miller, 1990), and retaining these covariates only affected other parts of the model through collinearity, and this was addressed directly.
Multilevel Models
The basic form of a multilevel model has been described previously within a series of methodological articles (Gilthorpe et al., 2000a,b,c, 2001, 2002; Gilthorpe and Cunningham, 2000; Lewsey et al., 2000, 2001). In brief, a multilevel model is a regression model for data that form a hierarchy, with variation at each level determined separately. Coefficients for explanatory variables may exhibit random variation about their mean estimate (at any level), leading to the inclusion of additional variance and covariance terms. Each covariance provides a measure of correlation between two variance terms. These models, known as random coefficient models, were evaluated for each periodontal outcome, yielding an overall mean trajectory across all sites, teeth, and subjects. All possible changes are summarized in the Fig.
Data were analyzed within the statistical package MLwiN (Rasbash et al., 2000), and model assumptions were confirmed by analysis of residuals (Goldstein, 1995). Model fit and improvement were assessed by examination of the -2 Log Likelihood statistic (McCullagh and Nelder, 1989). For nested models (i.e., where covariates are only added or subtracted, but not both), changes in this parameter followed a chi-squared distribution, with the number of degrees of freedom equal to the number of covariates added or subtracted.
| RESULTS |
|---|
|
|
|---|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
While cyclical change might be described as regression to the mean (RTM) (Blomqvist, 1987; Egelberg, 1989), results at the tooth and subject levels were not an artefact of measurement error. At the site level, the consistently negative correlation between the intercept and linear-time indicates that sites with initially higher-than-average scores experienced reduced linear change, such that some sites improved while others deteriorated at a much slower rate. This could be due to RTM. However, for higher levels of the hierarchy, RTM yields much less influence due to the averaging of site-level errors across teeth and subjects. While minor discrepancies in the correlations estimated might occur due to measurement error, the inference that change is oscillatory depends only upon the correlations sign.
The proposed theory needs verification for other cohorts, especially for more severely diseased subjects. It is necessary to validate the principle that disease patterns are cyclical, i.e., verifying the sign of time-dependent covariate correlations.
The Role of Covariates
Determination of potential risk factors and confounders is important for our understanding of disease prevention and treatment. After accounting for these factors, there remained residual variation, indicating that not all potentially important factors had been included. Trajectories for both outcomes were attenuated to reveal less improvement once risk factors were included. The effect of covariates was similar, though changes in the time-dependent coefficients differed between models, since overall changes in each outcome initially differed.
Not all risk factors were informative. Alcohol consumption had no impact on either outcome, and tobacco smoking demonstrated only a modest protective effect for PPD. Along with the usual difficulties with self-reported variables, these factors were recorded only at baseline, so there was no knowledge of changing behaviors. The impacts of smoking and alcohol therefore remain ambiguous.
Multicollinearity was observed through inconsistent subject-based coefficients between models that considered covariates in isolation and in combination (Table 3
). These covariates were aggregated variables, derived from site-based measures. Collinearity was greatest for the two most-correlated aggregate measures: LCAL > 1 mm and PPD > 2 mm. There were sizeable correlations between aggregates and their site-based measures. However, it is an advantage of MLM that measures related through aggregation can be combined in the same model, as illustrated by site-based covariates yielding robust coefficient estimates when subject-based aggregate coefficients were included (Table 3
). It is erroneous to place any meaningful interpretation on the subject-based covariate coefficientsan important observation given that subject-level (aggregate) covariates are often used in periodontal studies where less appropriate statistical methods are adopted.
Sites that experienced fluctuating LCAL/PPD (non-constant sites) were expected to experience greater disease on average, because many sites at baseline had little or no disease; for these sites, there was only one direction to change. Disease was progressively worse from anterior to posterior sites, while supragingival calculus moderated the progression of LCAL/PPD, mainly because supragingival calculus accumulated more in the lower anterior teeth.
The Role of Multilevel Modeling
The most powerful feature of MLM is the facility to investigate the underlying complexity of hierarchical systems, simultaneously modeling fixed effects and complex variation. Methods that accommodate hierarchy but fail to model variation explicitly (e.g., Generalized Estimating Equations) could not explore periodontal data in this manner. MLM provided a novel insight into the dynamic hierarchical system of periodontal breakdown and its progression, illustrating complex patterns of disease variation that had not been described by other methodologies. Many concepts attributed to the "linear" and "burst" models of periodontal disease progression are satisfied when examined in a multilevel context, since aspects of both theories have validity at different levels of the hierarchical system.
|
| ACKNOWLEDGMENTS |
|---|
Received July 5, 2001; Last revision October 30, 2002; Accepted November 7, 2002
| REFERENCES |
|---|
|
|
|---|
Blomqvist N (1987). On the bias caused by regression toward the mean in studying the relation between change and initial value. J Clin Periodontol 14:3437.[ISI][Medline]
Cohen ME, Ralls SA (1988). Distributions of periodontal attachment levels. Mathematical models and implications. J Periodontol 59:254258.[Medline]
Eaton KA, Duffy S, Griffiths GS, Gilthorpe MS, Johnson NW (2001). The influence of partial and full-mouth recordings on estimates of prevalence and extent of lifetime cumulative attachment loss: a study in a population of young male military recruits. J Periodontol 72:140145.[Medline]
Egelberg J (1989). The impact of regression towards the mean on probing changes in studies on the effect of periodontal therapy. J Clin Periodontol 16:120123.[Medline]
Gilthorpe MS, Cunningham SJ (2000). The application of multilevel, multivariate modelling to orthodontic research data. Community Dent Health 17:236242.[ISI][Medline]
Gilthorpe MS, Maddick IH, Petrie A (2000a). Introduction to Bayesian modelling in dental research. Community Dent Health 17:218221.[ISI][Medline]
Gilthorpe MS, Maddick IH, Petrie A (2000b). Introduction to multilevel modelling in dental research. Community Dent Health 17:222226.[ISI][Medline]
Gilthorpe MS, Griffiths GS, Maddick IH, Zamzuri AT (2000c). The application of multilevel modelling to periodontal research data. Community Dent Health 17:227235.[ISI][Medline]
Gilthorpe MS, Griffiths GS, Maddick IH, Zamzuri AT (2001). An application of multilevel modelling to longitudinal periodontal research data. Community Dent Health 18:7986.[ISI][Medline]
Gilthorpe MS, Mayhew MT, Bulman JS (2002). Multilevel survival analysis of amalgam restorations amongst RAF personnel. Community Dent Health 19:311.[Medline]
Goldstein H (1995). Multilevel statistical models. London: Edward Arnold.
Griffiths GS, Sterne JAC, Wilton JM, Eaton KA, Johnson NW (1992). Associations between volume and flow rate of gingival crevicular fluid and clinical assessments of gingival inflammation in a population of British male adolescents. J Clin Periodontol 19:464470.[Medline]
Griffiths GS, Duffy S, Eaton KA, Gilthorpe MS, Johnson NW (2001). Prevalence and extent of lifetime cumulative attachment loss (LCAL) at different thresholds and associations with clinical variables: changes in a population of young male military recruits over 3 years. J Periodontol 28:961969.
Hujoel PP, Leroux BG (1998). Evaluating the burst hypothesis at a site-specific level using the lack-of-fit test. J Periodontol 69:357362.[Medline]
Jeffcoat MK, Reddy MS (1991). Progression of probing attachment loss in adult periodontitis. J Periodont Res 62:185189.
Lewsey JD, Gilthorpe MS, Bulman JS, Bedi R (2000). Is modelling dental caries a normal thing to do? Community Dent Health 17:212217.[Medline]
Lewsey JD, Gilthorpe MS, Gulabivala K (2001). An introduction to meta-analysis within the framework of multilevel modelling using the probability of success of root canal treatment as an illustration. Community Dent Health 18:131137.[Medline]
Löe H, Anerud A, Boysen H, Smith M (1978). The natural history of periodontal disease in man. The rate of periodontal disease destruction before 40 years of age. J Periodontol 49:607620.[ISI][Medline]
Machtei EE, Norderyd J, Koch G, Dunford R, Grossi S, Genco RJ (1993). The rate of periodontal attachment loss in subjects with established periodontitis. J Periodontol 64:713718.[Medline]
McCullagh P, Nelder JA (1989). Generalized linear models. 2nd ed. London: Chapman & Hall.
Miller AJ (1990). Subset selection in regression. London: Chapman & Hall.
Ralls SA, Cohen ME (1986). Problems in identifying "bursts" of periodontal attachment loss. J Periodontol 57:746752.[Medline]
Rasbash J, Browne W, Cameron B, Charlton C (2000). MLwiN version 1.10.0006. London: Multilevel Models Project, Institute of Education.
Slinker BY, Glantz SA (1985). Multiple regression for physiological data analysis: the problem of multicollinearity. Am J Physiol 249:R1R12.
Socransky SS, Haffajee AD, Goodson JM, Lindhe J (1984). New concepts of destructive periodontal disease. J Clin Periodontol 11:2132.[ISI][Medline]
Sterne JA, Curtis MA, Gillett IR, Griffiths GS, Maiden MF, Wilton JM, et al. (1990). Statistical models for data from periodontal research. J Clin Periodontol 17:129137.[Medline]
Sterne JAC, Kingman A, Löe H (1992). Assessing the nature of periodontal disease progressionan application of covariance structure estimation. Appl Statist 41:539552.
Yang MC, Marks RG, Clark WB, Magnusson I (1991). Evaluation of statistical methods for monitoring periodontal disease. Stat Med 10:10891097.[Medline]
Yang MC, Namgung YY, Marks RG, Magnusson I, Clark WB (1993). Change detection on longitudinal data in periodontal research. J Periodontal Res 28:152160.[Medline]
This article has been cited by other articles:
![]() |
W. M. Thomson, R. Poulton, J. M. Broadbent, T. E. Moffitt, A. Caspi, J. D. Beck, D. Welch, and R. J. Hancox Cannabis Smoking and Periodontal Disease Among Young Adults JAMA, February 6, 2008; 299(5): 525 - 531. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Blance, Y.-K. Tu, and M. S Gilthorpe A multilevel modelling solution to mathematical coupling Statistical Methods in Medical Research, December 1, 2005; 14(6): 553 - 565. [Abstract] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| IADR Journals | Advances in Dental Research ® |
| Journal of Dental Research ® | Critical Reviews (1990-2004) |