Analyzing the emergence times of permanent teeth: an example of modeling the covariance matrix with interval-censored data.

*(English)*Zbl 07257142Summary: Based on a data set obtained in a large dental longitudinal study, conducted in Flanders (Belgium), the joint emergence distribution of seven teeth was modeled as a function of gender and caries experience on primary teeth. Besides establishing the marginal dependence of emergence on the covariates, there was also interest in examining the impact of the covariates on the association among emergence times. This allows the establishment of the preferred rankings of emergence and their dependence on covariates. To this end, the covariance matrix was modeled as a function of covariates. Modeling the covariance matrix in this way needs to ensure the positive definiteness of the covariance matrix and it is preferable that the regression parameters of the model are interpretable. The modified Cholesky decomposition of the covariance matrix, as suggested by Pourahmadi, splits up the covariance matrix into two parts where the parameters can be interpreted, given a natural ranking of the responses. This approach was used here taking into account that the emergence times are interval-censored. Hence, we opted for a Bayesian implementation of the data augmentation algorithm.

##### MSC:

62-XX | Statistics |

##### Keywords:

bivariate survival; covariance matrix; data augmentation; emergence times; interval-censored; multivariate normality
Full Text:
DOI

**OpenURL**

##### References:

[1] | Adler P (1963) Effect of some environmental factors on sequence of permanent tooth eruption . Journal of Dental Research 42, 605-616 . |

[2] | Barnard J , McCulloch R and Meng X (2000) Modelling covariance matrices in terms of standard deviations and correlations, with application to shrinkage . Statistica Sinica 10, 1281-1311 . · Zbl 0980.62045 |

[3] | Bogaerts K and Lesaffre E (2006) Estimating measures of association for bivariate interval censored data in a flexible way, available from the author. |

[4] | Bogaerts K , Leroy R , Declerck D and Lesaffre E (2002) Modelling tooth emergence data based on multivariate interval censored-data . Statistics in Medicine 21, 3375-3387 . |

[5] | Browne W , Draper D , Goldstein H and Rasbash, J (2002) Bayesian and likelihood methods for fitting multilevel models with complex level-1 variation . Computational Statistics and Data Analysis 39, 203-225 . · Zbl 1132.62312 |

[6] | Carlin BP and Louis TA (2000) Bayes and empirical bayes methods for data analysis. Chapman & Hall/CRC . · Zbl 1017.62005 |

[7] | Carvalho J , Ekstrand K and Thylstrup A (1989) Dental plaque and caries on occlusal surfaces of first permanent molars in relation to stage of eruption . Journal of Dental Research 68, 773-779 . |

[8] | Cecere S , Jara A and Lesaffre E (2006) Analyzing the emergence times of permanent teeth: an example of Bayesian modeling of a covariance matrix with interval-censored data. Technical report, IAP-statistics Technical Report Series. Report TR0608. |

[9] | Chen M-H and Shao Q-M (1999). Monte Carlo estimation of Bayesian Credible and HPD intervals . Journal of Computational and Graphical Statistics 8, 69-92 . |

[10] | Daniels M and Pourahmadi M (2002) Bayesian analysis of covariance matrices and dynamic models for longitudinal data . Biometrika 89, 553-566 . · Zbl 1036.62019 |

[11] | Demirjian A (1978) Dentition. In Falkner F and Tanner JM , editors, Human Growth 2. Baillère Tindahl , 413-444. |

[12] | Eskeli R , Laine-Avala M , Hausen H and Pahkala R (1999) Standards for permanent tooth emergence in Finnish children . The Angle Orthodontist, 69, 529-533 . |

[13] | Genz A , Bretz F and Hothorn T (2005) mvtnorm: Multivariate Normal and T Distribution. R package version 0.7-2. |

[14] | Hickernell F and Owen A (2005) Control variates for Quasi-Monte Carlo . Statistical Science 20, 1-31 . · Zbl 1100.65006 |

[15] | Kochhar R and Richardson A (1998) The chronology and sequence of eruption of human permanent teeth in Northern Ireland . International Journal of Paediatric Dentistry 8(4), 243-252 . |

[16] | Leroy R , Bogaerts K , Lesaffre E and Declerck D (2003a) The emergence of permanent teeth in Flemish children . Community Dentistry and Oral Epidemiology 31, 30-39 . |

[17] | Leroy R , Bogaerts K , Lesaffre E and Declerck D (2003b) Impact of caries experience in the deciduous molars on the emergence of the successors . European Journal of Oral Science 111, 106-110 . |

[18] | Mungonzibwa E , Kuijpers-Jagtman A , Laine-Avala M and van’t Hof M (2002) Emergence of permanent teeth in Tanzanian children . Community Dentistry and Oral Epidemiology 30, 455-462 . |

[19] | Parner E (2001) Non-parametric maximum likelihood estimation for bivariate interval censored data. Technical report, Department of Biostatistics, University of Aarhus . |

[20] | Parner E , Heidmann J , Vaeth M and Poulsen S (2001) A longitudinal study of time trends in the eruption of permanent teeth in Danish children . Archives of Oral Biology 46, 425-431 . |

[21] | Parner E , Heidmann J , Kjaer I , Vaeth M and Poulsen S (2002) Biological interpretation of the correlation of emergence times of permanent teeth . Journal of Dental Research 81(7), 451-454 . |

[22] | Pinheiro J and Bates D (1996) Unconstrained parametrizations for variance-covariance matrices . Statistics and Computing 6, 289-296 . |

[23] | Plummer M , Best N , Cowles K and Vines K (2005) Output analysis and diagnostics for MCMC. The Coda Package Version 0.9-5. |

[24] | Pourahmadi M (1999) Joint mean-covariance models with applications to longitudinal data: unconstrained parametrisation . Biometrika 86, 677-690 . · Zbl 0949.62066 |

[25] | Quinn K and Martin A (2002) Scythe Statistical Library Version 0.3. |

[26] | Ripley B (1987) Stochastic simulation. Wiley . · Zbl 0613.65006 |

[27] | Savara B and Steen J (1973) Timing and sequence of eruption of permanent teeth in a longitudinal sample of children from Oregon . Journal of American Dental Association 97, 209-214 . |

[28] | Shannon CE and Weaver W (1963) Mathematical theory of communication. University of Illinois Press . · Zbl 0126.35701 |

[29] | Smith B (2005) Bayesian Output Analysis Program (BOA) for MCMC, 2005 Version 1.1.5-2. |

[30] | Song J , Barnhart H and Lyles R (2004) A GEE approach for estimating correlation coefficients involving left-censored variables . Journal of Data Science 2, 245-257 . |

[31] | Tanner MA and Wong WH (1987) The calculation of posterior density by data augmentation . Journal of the American Statistical Association 82, 528-540 . · Zbl 0619.62029 |

[32] | Turnbull B (1976) The empirical distribution function with arbitrarily grouped, censored and truncated data . Journal of The Royal Statistical Society 38, 290-295 . · Zbl 0343.62033 |

[33] | Vanobbergen J , Martens L , Lesaffre E and Declerck D (2000) The Signal-Tandmobiel® project, a longitudinal intervention health promotion study in Flanders (Belgium): baseline and first year results . European Journal of Paediatric Dentistry 2, 87-96 . |

[34] | Virtanen J , Bloigu R and Larmas M (1994) Timing of eruption of permanent teeth: standard Finnish patient documents . Community Dentistry and Oral Epidemiology, 22, 286-288 . |

[35] | Wuertz D (2005) fOptions: Financial Software Collection -fOptions. R package version 220.10063. Retrieved from http://www.rmetrics.org. |

[36] | Yan J and Fine J (2004) Estimating equations for association measures . Statistics in Medicine 23, 859-874 . |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.