JDR JDR Most Cited Articles
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
QUICK SEARCH:   [advanced]


     

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chuang, S.-K.
Right arrow Articles by Cai, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chuang, S.-K.
Right arrow Articles by Cai, T.
J Dent Res 85(12):1147-1151, 2006
© 2006 International and American Associations for Dental Research

RESEARCH REPORT
Clinical

Predicting Clustered Dental Implant Survival Using Frailty Methods

S.-K. Chuang1,*, and T. Cai2

1 Department of Oral and Maxillofacial Surgery, Massachusetts General Hospital and Harvard School of Dental Medicine, 55 Fruit Street, Warren 1201, Boston, MA 02114, USA; and
2 Department of Biostatistics, Harvard School of Public Health, 677 Huntington Avenue, Boston, MA 02115, USA

* corresponding author, PO Box 67376, Chestnut Hill Station, Chestnut Hill, MA 02467, USA, schuang{at}hsph.harvard.edu


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The purpose of this study was to predict future implant survival using information on risk factors and on the survival status of an individual’s existing implant(s). We considered a retrospective cohort study with 677 individuals having 2349 implants placed. We proposed to predict the survival probabilities using the Cox proportional hazards frailty model, with three important risk factors: smoking status, timing of placement, and implant staging. For a non-smoking individual with 2 implants placed, an immediate implant and in one stage, the marginal probability that 1 implant would survive 12 months was 85.8% (95%CI: 77%, 91.7%), and the predicted joint probability of surviving for 12 months was 75.1% (95%CI: 62.1%, 84.7%). If 1 implant was placed earlier and had survived for 12 months, then the second implant had an 87.5% (95%CI: 80.3%, 92.4%) chance of surviving 12 months. Such conditional and joint predictions can assist in clinical decision-making for individuals.

KEY WORDS: clustered data • survival predictions • frailty • correlated survival analysis • proportional hazards model • dental implants


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
In clinical dental research, each participant often contributes multiple, potentially clustered observations to a database. For example, in dental implants, multiple implants placed within the same individual will be clustered or correlated. It is common practice for those conducting dental research to assume independent outcome observations in survival analyses (Higuchi et al., 1995; Wheeler, 1996; Buser et al., 1997; Brocard et al., 2000; Testori et al., 2001). These approaches tend to underestimate the variability of the estimated covariate effects or survival probabilities, because the possible within-cluster correlations are overlooked. To account for correlation, many authors have recommended randomly selecting 1 implant per person, and the use of only the subset for analysis (Lekholm et al., 1999; Herrmann et al., 1999, 2005; Gomez-Roman et al., 2001; Weibrich et al., 2001). However, this approach does not fully utilize all the available information, and thus may be inefficient. This gives rise to the need for the application of novel survival methods that can efficiently analyze clustered dental survival data.

There have been considerable statistical methods reported for the analysis of clustered survival data (Klein, 1992; Lee et al., 1992; Lin, 1994; Hougard, 1995, 2000; Parner, 1998; Spiekerman and Lin, 1998; Cai et al., 2000, 2002; Fine et al., 2003; Yin and Cai, 2005). Among these, two types of models have been proposed: frailty models and marginal models. Frailty models explicitly specify the within-cluster correlation structure, which allows for joint inference about the survival times within the same cluster. Marginal models leave the correlation structure completely unspecified, but adjust for the correlation by using robust sandwich-type estimators for the variance. Applications of these important methodologies to dental research have not yet been adequately explored.

We have previously applied the marginal approach (Chuang et al., 2002b) to predict covariate specific implant failure. The marginal approach, while effective, can predict only marginal survival probabilities of a single implant. In this research, we proposed to apply the frailty approach to predict implant survival accounting for within-individual correlation. In contrast to the marginal approach, the frailty approach allowed us to make joint predictions for several implants in the same individual. The primary purpose of this research was to predict the joint survival probability of multiple implants within the same individual, given their covariate information, and to predict the conditional survival probability of a future implant, given the survival status of an existing implant. The secondary aim was to compare and contrast the marginal and frailty approaches by comparing the point and interval estimates of the survival probabilities obtained based on the 2 methods.


   MATERIALS & METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
RESULTS
 DISCUSSION
 REFERENCES
 
Methods
The research methods have been described in detail in other studies (Chuang et al., 2002a; Vehemente et al., 2002). We conducted a retrospective cohort study design composed of individuals who had at least 1 Bicon® implant placed between May 20, 1992, and July 6, 2000, at the Implant Dentistry Center based at Faulkner Hospital (IDC-FH), Boston, MA, USA. We reviewed the participants’ records to collect information regarding exposures associated with implant failures grouped into the following main categories: demographic, health status, anatomic, implant- and abutment-specific, anticipated restoration, peri-operative chemotherapy, reconstructive, and operator. The major outcome variable of interest was implant failure, which was defined as the removal of the implant for any reason (Dental Implant Clinical Research Group, 1997). We estimated survival time by computing the difference in time (mos) between implant placement and implant explantation, or the date of the last follow-up visit, for individuals whose implants had not been removed. The use of human participants in this study followed a protocol approved by the human research committee.

Statistical Issues
In this report, our interest focused on survival prediction for multiple implants under the Cox proportional hazards frailty model (Hougaard, 1995, 2000). Specifically, we assumed that, conditional on some unobserved frailty random variable, the time to dental implant failure would follow a Cox proportional hazards model. Multiple dental implants from each individual share a common frailty random variable, which accounts for the within-individual correlation. Under this model, we described procedures for predicting the conditional and joint survival probabilities for dental implants with given covariate information. In particular, the non-parametric maximum likelihood estimator for the unknown parameters in the model (Parner, 1998) was used to construct plug-in estimators for the survival probabilities. The bootstrap method (Efron and Tibshirani, 1993) was used to obtain the interval estimates. The set of covariates selected to be used in this study—i.e., current tobacco use, timing of implant placement, and implant staging—was previously identified by means of the Cox frailty regression model, adjusted for clustered observations (Chuang et al., 2005). Databases were stored in SAS-PC Version 8.2 (SAS Institute, Inc., 1999) files. The statistical program for clustered survival predictions was written by the co-author (Cai) in the S-plus environment (Version 6.0, 2000). [We have made the S-plus computer codes available on the Internet at the following addresses: http://biosun1.harvard.edu/~tcai/FUN_CoxFrailtyPred.q and http://biosun1.harvard.edu/~tcai/CoxFrailtyPred.q.]

Statistical Notation
Let Tik denote the failure time of the kth implant for the ith individual, k = 1,2,3,....,Ki; i = 1,2,3,...,n, where n is the sample size. Let Zik be the corresponding covariate vector for Tik and Cik be the censoring time. For Tik, one can observe only Xik = min(Tik, Cik) and {delta}ik = 1(Tik = Xik). Conditional on the covariate vector Zi = (Zi1,...,ZiKi)', the censoring vector, Ci = (Ci1,...,CiKi)', is assumed to be independent of the failure time vector, Ti = (Ti1,...,TiKi)'. We assume that, conditional on Zi and an unobserved frailty wi, Tik follows a Cox proportional hazards model:


Formula 1(1)

where {lambda}ik (t) is the hazard function for Tik, {lambda}0 (t) is an unspecified baseline hazard function, and ß is the true regression coefficient vector. The frailty wi accounts for the within-individual correlation, due to some unobserved common covariate information. The unobserved wi’s are assumed to be independent and identically distributed, with unit mean and some unknown variance {sigma}2. Different individuals could have different values of frailty, and the variability in the wi’s reflects the heterogeneity of risks between individuals. For computational convenience, the frailty distribution is often taken to be a Gamma, i.e., wi~{sigma}2gamma({sigma}–2), and thus we assume a gamma-frailty throughout.

Various inference procedures have been proposed for the covariate effect ß, the baseline cumulative hazard function, {Lambda}0(t) = {int}0t{lambda}0(u)du, and the variance component {sigma}2 (Klein, 1992; Nielsen et al., 1992; Parner, 1998). One popular estimator is the non-parametric maximum likelihood estimator (NPMLE) studied by Parner (1998). The NPMLE can be implemented through an EM-algorithm (Therneau and Grambsch, 2000), which alternates between two steps: (1) the M-step, in which ß and {Lambda}0 are updated as in normal Cox regression by treating the current estimates of {w1, ..., wn} as fixed values and offsets; and (2) the E-step, in which {w1,..., wn} are computed as the expected value, given the current values of ß and {Lambda}0 and the data. These estimates depend on {sigma}2 and are denoted by Formula 1 ({sigma}2) and Formula 10(|{sigma}2). {sigma}2 can then be estimated by maximizing the profile log-likelihood (Parner, 1998). [See Hougaard (2000) and Therneau and Grambsch (2000) for more details on the estimation procedures for ß, {Lambda}0, and the variance component {sigma}2.]

To predict the survival probabilities for the implants based on the Cox gamma frailty model, we noted that the survival probability for an implant with covariate Z can be expressed as:


Formula 2(2)

Note that S(t|Z) here is the population average survival probability, averaged across all implants with covariate Z and all possible frailty levels w. A plug-in estimator for S(t|Z) may be obtained as Formula 2(t|Z) = {1 + Formula 22Formula 20(t)eFormula 2'Z}–{sigma}–2, where Formula 22, Formula 2, and Formula 20. are the respective estimators for {sigma}2, ß, and {Lambda}0. The main strength of the frailty model is in its ability to make simultaneous joint inference about the survival probabilities of multiple implants from the same individual. For example, for K implants of an individual with covariate levels (Z11,...,Z1K), the joint survival probability, S(t1,...,tK) = P(T11 ≥ t1,..., T1K ≥ tK | Z11,..., Z1K), can be estimated by


Formula 3(3)

This also allowed us to predict the conditional survival probability for a future implant, given the information on the survival of existing implants and their covariate levels. This can be achieved based on the definition of the conditional probability:


Formula 4(4)

To construct confidence intervals for these predicted survival probabilities, we needed to account for the variabilities in all the estimated parameters, and thus required the estimation of the covariance of Formula 4 = (Formula 42, Formula 4, Formula 40 (t1),..., Formula 40 (tK))'. Due to the difficulty in obtaining an analytical variance-covariance estimator for Formula 4, we approximated the limiting distribution of Formula 4 by the bootstrap method (Efron and Tibshirani, 1993; Monaco et al., 2005), which is widely used in practice when analytical variance of an estimator is difficult to obtain. To account for within-individual correlations, we treated multiple observations from the same individual as a unit, and thus the original dataset was treated as n datapoints. We obtained each bootstrap sample by sampling a size of n datapoints with replacement from the original data; we used a total of 1000 bootstrap samples. For the mth bootstrap sample, we obtained a bootstrap estimator for {theta}, denoted by Formula 4(m). Plugging Formula 4(m) into Eq. (3), we obtained a corresponding bootstrap estimator for the survival probability Formula 4(m)(t1,...tK). Then the distribution of Formula 4, Formula 4(t1,...tK), and the estimated conditional probabilities could be approximated by the empirical distribution of their respective bootstrap replicates. For example, the variance of Formula 4(t1,...tK) can be estimated by


Formula 5(5)

A 95% confidence interval for S(t1,...,tK) can be obtained based on the 2.5% and 97.5% quantiles of {Formula 5(1)(t1,...tK), ...,Formula 51000(t1,...tK)}.


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
DISCUSSION
 REFERENCES
 
The study group had 677 participants, who had 2349 dental implants placed. On average, 3.5 implants (range, 1–22) were placed per person. There were 137 implants that failed. There were 57 smokers (10.3%, N = 554). Out of 2349 implants, 339 (14.4%) were placed in one stage, and 243 (10.3%) were placed immediately after tooth extraction. The details of the descriptive statistics have been described (Chuang et al., 2002a).

We previously identified 3 key variables associated with implant survival: current tobacco use (yes or no), timing of implant placement (delayed vs. immediate), and implant staging (one- or two-stage) (Chuang et al., 2002a). Based on a clustered marginal approach (Chuang et al., 2002b), we found that, among those who had a delayed procedure, did not smoke, and underwent a two-stage implant procedure (best-case scenario), the predicted one- and five-year survival rates were 97.2% and 93.4%, respectively. For those who had an immediate implant placed, smoked, and underwent a one-stage procedure (worst-case scenario), the predicted survival rates at 1 and 5 yrs were 58.5% and 27.6%, respectively. We compared the predicted survival probabilities based on the marginal approach (Chuang et al., 2002b) and the current frailty prediction method for non-smokers with immediate implant placement in one stage (Table 1Go). The predicted 12-month (1 yr) survival was 83.7%, based on the marginal method, and 85.8% with the frailty method. Estimates from these 2 methods were not drastically different, suggesting that both models may fit the data well.


View this table:
[in this window]
[in a new window]
  		Table 1. Predicted S(t), the Survival Probability of T ≥ t (mos), for an Immediate Implant Placed in One Stage in a Non-smoker, along with 95% Confidence Intervals, Based on Both the Marginal Method and the Frailty Method
 
In this study, we were particularly interested in making joint inference about multiple implants in the same individual. First, we considered the subgroup of non-smoking persons having 2 immediate implants placed in 1 stage. The predicted survival probability (95% confidence interval) was 75.1% (62.1%, 84.7%) for both implants surviving 12 mos (1 yr), and was 55.4% (36.2%, 71.6%) for both implants surviving 60 mos (5 yrs) (Table 2Go). In an individual with an existing implant, one may be interested in predicting the survival for a future implant, given the survival status of the current implant. For example, for a non-smoker with an existing immediate one-stage implant that has survived 1 yr, the predicted probability for a future implant with the same covariate surviving one year is 87.5%, with a 95% confidence interval (80.3%, 92.4%) (Table 3Go). We next considered those who were non-smokers having the first implant placed in a delayed fashion in 1 stage, and having the second implant placed in an immediate fashion in 1 stage. The predicted one- and five- survival rates for such 2 implants were 79.6% (70.3%, 86.5%) and 61.6% (44.0%, 74.4%), respectively (Table 4Go). The results of the above analyses were specific and unique only to the frailty approach prediction method. This cannot be performed with the marginal approach method.


View this table:
[in this window]
[in a new window]
  		Table 2. Predicted S(t1,t2), the Joint Survival Probability of T1 ≥ t1 (mos) and T2 ≥ t2 (mos) for Two Immediate Implants Placed in One Stage in a Non-smoker
 

View this table:
[in this window]
[in a new window]
  		Table 3. Predicted Conditional Probability of T2 ≥ t2 (mos) Given T1 ≥ t1 (mos) for Two Immediate Implants Placed in Stage One in a Non-smoker
 

View this table:
[in this window]
[in a new window]
  		Table 4. Predicted Joint Probabilities of T1 and T2 Both Exceeding t (mos) Given the Different Exposures for the Two Implants
 

   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
REFERENCES
 
Multivariate failure times with clustered or dependent data observations are commonly encountered in dental research. To make valid and efficient statistical inference in such settings, one needs to account for within-cluster correlations. Useful tools for analyzing clustered survival data include the marginal and frailty survival models. The marginal approach leaves the correlation structure unspecified, and adjusts for the correlation by using a sandwich-type variance estimator (Lee et al., 1992). The frailty approach specifies the correlation structure through the distribution of the unobserved shared frailty random variable w. One main advantage of the frailty model is in its ability to make joint prediction of the survival of multiple implants from the same individual. It also allows us to predict the survival of a future implant based on the survival status of current implants and covariate information. The purposes of this study were two-fold: to produce joint predictions of the survival of multiple implants, based on valid and efficient survival methods, accounting for within-cluster correlation; and to compare and contrast how the marginal approach and the frailty approach differ in survival predictions, given important exposure information. The estimators used in this paper for the unknown parameters in the frailty model were based on the non-parametric maximum likelihood method (Parner, 1998), and the implementation of this method has been discussed in detail (Therneau and Grambsch, 2000).

To illustrate the strength of the frailty model, we investigated within-cluster implant prediction using the ’other implant(s) information’ from the same individual. Our first model (marginal method) consisted of results from our previous study (Chuang et al., 2002b) (Table 1Go). In our second model, we utilized the frailty approach using implants from the same individual for survival predictions (frailty method). When comparing the two valid analytic strategies, we found that the one- and five-year survival point estimates did not differ drastically, but the variance estimates did differ. The 95% confidence intervals for the marginal method were wider than those of the frailty method. This can be partly attributed to the fact that the frailty method assumes a stronger model by specifying the correlation structure. In the setting of clustered observations, we believe that both methods are valid and efficient, as long as the assumed models fit the data well. The ultimate choice of method depends on the research question and the goals for the particular study.

In summary, clustered survival observations are frequently encountered in many different areas of person-oriented dental research. In particular, we are interested in making joint predictions for the survival of multiple dental implants in the same individual, given various clinical parameters. The Cox frailty model used in this study produced valid and efficient estimates, with the strength of allowing for joint prediction. Using available information about clinical parameters, along with the survival status of existing implants, we were able to make more precise predictions about the survival of future implants. Such predictions can assist in surgical treatment-planning and thus lead to better clinical decision-making and patient care.


   ACKNOWLEDGMENTS
 
This research was supported by NIH/NIDCR R03 DE016919 (SKC, TC), and by NIH grant R37 AI024643 (TC). We thank the Harvard School of Public Health’s Department of Biostatistics for access to statistical computation. The authors also recognize the clinicians and staff of the Implant Dentistry Centre at the Faulkner Hospital, Boston, MA, for their cooperation in this study, and for free and unfettered access to patient records.

Received November 29, 2005; Last revision June 13, 2006; Accepted September 25, 2006


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Brocard D, Barthet P, Baysse E, Duffort JF, Eller P, Justumus P, et al. (2000). A multicenter report on 1,022 consecutively placed ITI implants: a 7-year longitudinal study. Int J Oral Maxillofac Implants 15:691–700.[ISI][Medline]

Buser D, Mericske-Stern R, Bernard JP, Behneke A, Behneke N, Hirt HP, et al. (1997). Long-term evaluation of non-submerged ITI implants. Part 1: 8-year life table analysis of a prospective multi-center study with 2359 implants. Clin Oral Implants Res 8:161–172.[ISI][Medline]

Cai T, Wei LJ, Wilcox M (2000). Semiparametric regression analysis for clustered failure time data. Biometrika 87:867–878.[Abstract/Free Full Text]

Cai T, Cheng SC, Wei LJ (2002). Semiparametric mixed-effects models for clustered failure time data. J Am Stat Assoc 97:514–522.

Chuang SK, Wei LJ, Douglass CW, Dodson TB. (2002a). Risk factors for dental implant failure: a strategy for the analysis of clustered failure-time observations. J Dent Res 81:572–577.[Abstract/Free Full Text]

Chuang SK, Tian L, Wei LJ, Dodson TB (2002b). Predicting dental implant survival by use of the marginal approach of the semi-parametric survival methods for clustered observations. J Dent Res 81:851–855.[Abstract/Free Full Text]

Chuang SK, Cai T, Douglass CW, Wei LJ, Dodson TB (2005). Frailty approach for the analysis of clustered failure time observations in dental research. J Dent Res 84:54–58.[Abstract/Free Full Text]

Dental Implant Clinical Research Group (1997). Clinical factors of importance to successful implant therapy. Introduction. J Oral Maxillofac Surg 55(12 Suppl 5):7–11.[ISI][Medline]

Efron B, Tibshirani RJ (1993). An introduction to the bootstrap. New York: Chapman & Hall.

Fine JP, Glidden DV, Lee KE (2003). A simple estimator for a shared frailty regression model. J R Stat Soc B 65:317–329.

Gomez-Roman G, Kruppenbacher M, Weber H, Schulte W (2001). Immediate postextraction implant placement with root-analog stepped implants: surgical procedure and statistical outcome after 6 years. Int J Oral Maxillofac Implants 16:503–513.[ISI][Medline]

Herrmann I, Lekholm U, Holm S, Karlsson S (1999). Impact of implant interdependency when evaluating success rates: a statistical analysis of multicenter results. Int J Prosthodont 12:160–166.[ISI][Medline]

Herrmann I, Lekholm U, Holm S, Kultje C (2005). Evaluation of patient and implant characteristics as potential prognostic factors for oral implant failures. Int J Oral Maxillofac Implants 20:220–230.[ISI][Medline]

Higuchi KW, Folmer T, Kultje C (1995). Implant survival rates in partially edentulous patients: a 3-year prospective multicenter study. J Oral Maxillofac Surg 53:264–268.[ISI][Medline]

Hougaard P (1995). Frailty models for survival data. Lifetime Data Anal 1:255–273.[Medline]

Hougaard P (2000). Analysis of multivariate survival data. New York: Springer-Verlag.

Klein JP (1992). Semiparametric estimation of random effects using the Cox model based on the EM algorithm. Biometrics 48:795–806.[ISI][Medline]

Lee EW, Wei LJ, Amato DA (1992). Cox-type regression analysis for large numbers of small groups of correlated failure time observations. In: Survival analysis: state of the art. Klein JP, Goel PK, editors. Dordrecht, Germany: Kluwer Academic, pp. 237–247.

Lekholm U, Gunne J, Henry P, Higuchi K, Linden U, Bergström C, et al. (1999). Survival of the Brånemark implant in partially edentulous jaws: a 10-year prospective multicenter study. Int J Oral Maxillofac Implants 14:639–645.

Lin DY (1994). Cox regression analysis of multivariate failure time data: the marginal approach. Stat Med 13:2233–2247.[ISI][Medline]

Monaco J, Cai J, Grizzle J (2005). Bootstrap analysis of multivariate failure time data. Stat Med 24:3387–3400.[ISI][Medline]

Nielsen G, Gill R, Andersen P, Sorensen T (1992). A counting process approach to maximum likelihood estimation in frailty models. Scand J Stat 19:25–43.

Parner E (1998). Asymptotic theory for the correlated gamma-frailty model. Ann Statist 26:183–214.

SAS Institute Inc. (1999). SAS System. Version 8.2. Cary, NC: SAS Institute Inc.

S-Plus (2000). StatSci, A division of MathSoft, Inc., Seattle, Washington.

Spiekerman CF, Lin DY (1998). Marginal regression models for multivariate failure time data. J Am Stat Assoc 93:1164–1175.

Testori T, Wiseman L, Woolfe S, Porter SS (2001). A prospective multicenter clinical study of the Osseotite implant: four-year interim report. Int J Oral Maxillofac Implants 16:193–200.[ISI][Medline]

Therneau TM, Grambsch PM (2000). Modeling survival data. Extending the Cox model. New York: Springer-Verlag.

Vehemente VA, Chuang SK, Daher S, Muftu A, Dodson TB (2002). Risk factors affecting dental implant survival. J Oral Implantol 28:74–81.[Medline]

Weibrich G, Buch RS, Wegener J, Wagner W (2001). Five-year prospective follow-up report of the Astra Tech standard dental implant in clinical treatment. Int J Oral Maxillofac Implants 16:557–562.[ISI][Medline]

Wheeler SL (1996). Eight-year clinical retrospective study of titanium plasma-sprayed and hydroxyapatite-coated cylinder implants. Int J Oral Maxillofac Implants 11:340–350.[Medline]

Yin G, Cai J (2005). Quantile regression models with multivariate failure time data. Biometrics 61:151–161.[ISI][Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chuang, S.-K.
Right arrow Articles by Cai, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chuang, S.-K.
Right arrow Articles by Cai, T.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
IADR Journals Advances in Dental Research ®
Journal of Dental Research ® Critical Reviews (1990-2004)