Modè<strong>le</strong> nonlinéaire <strong>à</strong> processus latent 741022 Biometrics, December 2006has previously been discussed by Mo<strong>le</strong>nberghs and Verbeke(2001).5.5 Multivariate Model versus Univariate ModelsFor each test, the univariate model detected an associationbetween apoE genotype and cognition with a larger p-value(p-value from the likelihood ratio test for apoE × t 2 parameter:p = 0.0043 for MMSE, p = 0.018 for IST, p = 0.020for BVRT, p = 0.054 for the DSSTW) than for the multivariatemodel (p = 0.0018). By using a multivariate modelcompared to four univariate models, we had a gain of powerin assessing the association between apoE genotype and cognition.Moreover, note that interpretation of the associationwith the latent process and with each psychometric test isdifferent.The gain in efficiency can also be evaluated by comparingthe AIC from the multivariate model and the AIC computedby pooling the likelihoods from the four univariate modelswith the total number of parameters in these four models.In our case, even if we added in the multivariate model theconstraint that apoE had a common effect on the four tests,the AIC from the multivariate model was markedly better(39503.1 vs. 40203.8 for the four univariate models).6. DiscussionWe proposed a nonlinear model for multivariate longitudinalnon-Gaussian quantitative outcomes when the outcomes areindirect measures of a common underlying continuous-timeprocess. Such data are very frequent in psychometrics, but themethodology has many other potential areas of application,as, for instance, the study of the course of chronic illnessesevaluated by several biological markers.In this work, psychometric tests are analyzed by consideringtheir sum scores as quantitative variab<strong>le</strong>s. This is themost frequent way to consider psychometric tests in gerontology:the summary scores are used to evaluate cognitive<strong>le</strong>vel and risk of dementia (Hall et al., 2001; Sliwinski etal., 2003; Amieva et al., 2005). From a neuropsychologicalperspective, the alternative approach, which consists of analyzingitem responses using SEM or item response model(Skrondal and Rabe-Hesketh, 2004, Chapter 3), could be usefulif the objective was to understand the underlying componentsof the tests. This methodology is interesting when thetests consist of a limited number of items evaluating differentcognitive domains (such as the MMSE) or exhibiting different<strong>le</strong>vels of difficulty (such as the BVRT). On the contrary,this methodology would be difficult to apply to the IST score,which is the number of words cited by the subjects (exceptfor considering the four subscores for each semantic category)and the DSSTW score, which is the count of symbols correctlyassigned to a sequence of numbers.Given that the summary scores are quantitative discretevariab<strong>le</strong>s, we could either consider them as continuous variab<strong>le</strong>sor as ordinal variab<strong>le</strong>s. However, threshold models forordinal data require estimation of one threshold for each possib<strong>le</strong>value of the scores, which would be very chal<strong>le</strong>nging formultivariate modeling of scores with so many different values.In our application including four tests, this would haveimplied estimation of more than 150 additional parameters.Thus, we decided to analyze the scores as continuous variab<strong>le</strong>sand to use nonlinear transformations depending on a limitednumber of parameters as link functions between the Gaussianlatent process and the observed outcomes. Various linkfunctions have been considered, but we found that the betaCDF was f<strong>le</strong>xib<strong>le</strong> enough with only two parameters. Withmany fewer parameters than threshold models, the estimatedcurves provide interesting information on the relationship betweenevolution of the latent cognitive <strong>le</strong>vel and evolution ofthe observed scores. Moreover, goodness-of-fit analyses showthat the beta transforms of the four test scores fitted wella Gaussian distribution. Neverthe<strong>le</strong>ss, if necessary for otherapplications, it would be easy to include a different familyof continuous transformation for each test. This model couldalso be extended to allow a mixture of continuous, binary, andordinal outcomes with few categories as in Dunson (2000),Dunson (2003), or Rabe-Hesketh et al. (2004).Another asset of this model is the way of accounting forcovariate effects. Using fixed contrasts, we were ab<strong>le</strong> to distinguishbetween the association with the latent process and thedifferential association with the various psychometric tests.Being ab<strong>le</strong> to compare covariate effects over the tests is alsoan advantage of multivariate modeling. With a large numberof tests it may also be possib<strong>le</strong> and advantageous to have theeffects of the covariates on the tests be random rather thanfixed.We assumed the data to be missing at random and thusignorab<strong>le</strong> using a maximum likelihood approach. Even if wehave shown that missing data were associated with randomeffects, this does not preclude missing data from being ignorab<strong>le</strong>.Indeed, if missing data at the last three tests depend onthe observed MMSE score or on the observed evolution of theMMSE, it induces a dependency on the random effects, butthe missing data are ignorab<strong>le</strong> because missing values maybe predicted using observed data. This is an advantage ofmultivariate modeling, in that by using more observed informationit is more robust to missing data. Neverthe<strong>le</strong>ss, as itis not excluded that missing data are informative, it could beuseful to jointly model time to dropout using a shared randomeffect model as in Roy and Lin (2002). However, thiswould increase comp<strong>le</strong>xity of the estimation process and estimateswould depend on uncheckab<strong>le</strong> parametric assumptions.Another useful extension would be to jointly model dementia,defining dementia diagnosis as the time at which the latentprocess first reaches an estimated threshold (Hashemi,Jacqmin-Gadda, and Commenges, 2003).Acknow<strong>le</strong>dgementsThe authors thank Jean-François Dartigues, Luc Letenneur,and the PAQUID research team for providing the data.ReferencesAmieva, H., Jacqmin-Gadda, H., Orgogozo, J., LeCarret,N., Helmer, C., Letenneur, L., Barberger-Gateau, P.,Fabrigou<strong>le</strong>, C., and Dartigues, J. (2005). The 9 year cognitivedecline before dementia of the Alzheimer type: Aprospective population-based study. Brain 128, 1093–1101.Arminger, G. and Muthén, B. (1998). A Bayesian approach tononlinear latent variab<strong>le</strong> models using the Gibbs samp<strong>le</strong>r
Modè<strong>le</strong> nonlinéaire <strong>à</strong> processus latent 75Nonlinear Model with Latent Process for Cognitive Evolution 1023and the Metropolis-Hastings algorithm. Psychometrika63, 271–300.Dunson, D. (2000). Bayesian latent variab<strong>le</strong> models for clusteredmixed outcomes. Journal of the Royal StatisticalSociety, Series B 62, 355–366.Dunson, D. (2003). Dynamic latent trait models for multidimensionallongitudinal data. Journal of the AmericanStatistical Association 98, 555–563.Fabrigou<strong>le</strong>, C., Rouch, I., Taberly, A., Letenneur, L.,Commenges, D., Mazaux, J.-M., Orgogozo, J.-M., andDartigues, J.-F. (1998). Cognitive process in preclinicalphase of dementia. Brain 121, 135–141.Farrer, L., Cupp<strong>le</strong>s, L., Haines, J., Hyman, B., Kukull, W.,Mayeux, R., Myers, R., Pericak-Vance, M., Risch, N.,and Van Duijn, C. (1997). Effects of age, sex and ethnicityon the association between apolipoprotein E genotypeand Alzheimer disease. A meta-analysis. APOE andAlzheimer Disease Meta Analysis Consortium. Journal ofthe American Medical Association 278, 1349–1356.Gray, S. M. and Brookmeyer, R. (1998). Estimating a treatmenteffect from multidimensional longitudinal data.Biometrics 54, 976–988.Hall, C., Ying, J., Kuo, L., Sliwinski, M., Buschke, H., Katz,M., and Lipton, R. (2001). Estimation of bivariate measurementshaving different change points, with applicationto cognitive ageing. Statistics in Medicine 20, 3695–3714.Harvey, D., Beckett, L., and Mungas, D. (2003). Multivariatemodeling of two associated cognitive outcomes in a longitudinalstudy. Journal of Alzheimer’s Disease 5, 357–365.Hashemi, R., Jacqmin-Gadda, H., and Commenges, D. (2003).A latent process model for joint modeling of events andmarker. Lifetime Data Analysis 9, 331–343.Jacqmin-Gadda, H., Fabrigou<strong>le</strong>, C., Commenges, D., andDartigues, J.-F. (1997). A 5-year longitudinal study ofthe Mini Mental State Examination in normal aging.American Journal of Epidemiology 145, 498–506.Jöreskog, K. and Yang, F. (1996). Nonlinear structural equationmodels: The Kenny-Judd model with interaction effects.In Advanced Structural Equation Modeling: Issuesand Techniques, G. A. Marcoulides and R. E. Schumacker(eds), 57–88. Mahwah, New Jersey: Lawrence ErlbaumAssociates.Kendall, M. and Stuart, A. (1977). The Advanced Theory ofStatistics, Volume 1. New York: MacMillan Publishing.Lee, S. and Song, X. (2004). Bayesian model comparison ofnonlinear structural equation models with missing continuousand ordinal categorical data. British Journal ofMathematical and Statistical Psychology 57, 131–150.Letenneur, L., Commenges, D., Dartigues, J.-F., andBarberger-Gateau, P. (1994). Incidence of dementia andAlzheimer’s disease in the elderly community residents ofsouth-western France. International Journal of Epidemiology23, 577–590.Longford, N. and Muthén, B. (1992). Factor analysis for clusteredobservations. Psychometrika 57, 581–597.Marquardt, D. (1963). An algorithm for <strong>le</strong>ast-squares estimationof nonlinear parameters. SIAM Journal on AppliedMathematics 11, 431–441.Mo<strong>le</strong>nberghs, G. and Verbeke, G. (2001). A review on linearmixed models for longitudinal data, possibly subject todropout. Statistical Modelling 1, 235–269.Muthén, B. (2002). Beyond SEM: General latent variab<strong>le</strong>modeling. Behaviormetrika 29, 81–117.Rabe-Hesketh, S., Skrondal, A., and Pick<strong>le</strong>s, A. (2004). Generalizedmulti<strong>le</strong>vel structural equation modelling. Psychometrika69, 167–190.Roy, J. and Lin, X. (2000). Latent variab<strong>le</strong> models for longitudinaldata with multip<strong>le</strong> continuous outcomes. Biometrics56, 1047–1054.Roy, J. and Lin, X. (2002). Analysis of multivariate longitudinaloutcomes with nonignorab<strong>le</strong> dropouts and missingcovariates: Changes in methadone treatment practices.Journal of the American Statistical Association 97, 40–52.Salthouse, T., Hancock, H., Meinz, E., and Hambrick, D.(1996). Interrelations of age, visual acuity, and cognitivefunctioning. Journal of Gerontology: PsychologicalSciences 51B, 317–330.Sammel, M. and Ryan, L. (1996). Latent variab<strong>le</strong> models withfixed effects. Biometrics 52, 650–663.Sánchez, B. N., Budtz-Jorgensen, E., Ryan, L. M., and Hu,H. (2005). Structural equation models: A review withapplications to environmental epidemiology. Journal ofthe American Statistical Association 100, 1443–1455.Skrondal, A. and Rabe-Hesketh, S. (2004). Generalized LatentVariab<strong>le</strong> Modelling: Multi<strong>le</strong>vel, Longitudinal and StructuralEquation Models. Boca Raton, Florida: Chapman& Hall/CRC.Sliwinski, M., Hofer, S., and Hall, C. (2003). Correlated andcoup<strong>le</strong>d cognitive change in older adults with or withoutpreclinical dementia. Psychology and Aging 18, 672–683.Song, X. and Lee, S. (2004). Bayesian analysis of two-<strong>le</strong>velnonlinear structural equation models with continuousand polytomous data. British Journal of Mathematicaland Statistical Psychology 57, 29–52.Wall, M. and Amemiya, Y. (2000). Estimation for polynomialstructural equation models. Journal of the American StatisticalAssociation 95, 929–940.Winnock, M., Letenneur, L., Jacqmin-Gadda, H.,Dallongevil<strong>le</strong>, J., Amouyel, P., and Dartigues, J.(2002). Longitudinal analysis of the effect of apolipoproteinE ɛ4 and education on cognitive performancein elderly subjects: The PAQUID study. Journal ofNeurology, Neurosurgery and Psychiatry 72, 794–797.Yalcin, I. and Amemiya, Y. (2001). Nonlinear factor analysisas a statistical method. Statistical Science 16, 275–294.Received January 2005. Revised December 2005.Accepted January 2006.AppendixDetails on the TransformationsFor subject i and test k, we give the expressions of the functionand its Jacobian for the linear transformation, the combinationof a linear and a logit transformation, and the WeibullCDF.The linear transformation is defined for y ∈ R,η 1k ∈ R, andη 2k ∈ R ∗ as
- Page 1 and 2:
Université Victor Segalen Bordeaux
- Page 3 and 4:
3RemerciementsA Monsieur Jean-Louis
- Page 5 and 6:
5Un immense merci à tous ceux que
- Page 7 and 8:
7A Delphine et ses mille et une his
- Page 9 and 10:
TABLE DES MATIÈRES 92.3.2 Estimati
- Page 11 and 12:
TABLE DES MATIÈRES 116.3.1 Adéqua
- Page 13 and 14:
Introduction 13ans atteintes d’un
- Page 15 and 16:
Introduction 15Mais, pour l’insta
- Page 17 and 18:
Introduction 171.2 Problèmes méth
- Page 19 and 20:
Introduction 191.2.3 Association en
- Page 21 and 22:
Introduction 21QUID suggère que le
- Page 23 and 24: Chapitre 2Etat des connaissancesCe
- Page 25 and 26: Etat des connaissances 25Plusieurs
- Page 27 and 28: Etat des connaissances 27Prise en c
- Page 29 and 30: Etat des connaissances 29Extensions
- Page 31 and 32: Etat des connaissances 31données s
- Page 33 and 34: Etat des connaissances 33variables
- Page 35 and 36: Etat des connaissances 35(2000) ne
- Page 37 and 38: Etat des connaissances 372.3 Modél
- Page 39 and 40: Etat des connaissances 39chapitre 3
- Page 41 and 42: Etat des connaissances 41de mélang
- Page 43 and 44: Etat des connaissances 43Hawkins et
- Page 45 and 46: Etat des connaissances 45tique de d
- Page 47 and 48: Etat des connaissances 47L’estima
- Page 49 and 50: Etat des connaissances 49normales s
- Page 51 and 52: Etat des connaissances 512.4 Modél
- Page 53 and 54: Etat des connaissances 53en tant qu
- Page 55 and 56: Etat des connaissances 55interactio
- Page 57 and 58: Etat des connaissances 57décrite p
- Page 59 and 60: Etat des connaissances 59entre l’
- Page 61 and 62: Etat des connaissances 612.4.3 Cas
- Page 63 and 64: Etat des connaissances 63l’inform
- Page 65 and 66: Modèle nonlinéaire à processus l
- Page 67 and 68: Modèle nonlinéaire à processus l
- Page 69 and 70: Modèle nonlinéaire à processus l
- Page 71 and 72: Modèle nonlinéaire à processus l
- Page 73: Modèle nonlinéaire à processus l
- Page 77 and 78: Modèle nonlinéaire à processus l
- Page 79 and 80: Modèle nonlinéaire à processus l
- Page 81 and 82: Modèle nonlinéaire à processus l
- Page 83 and 84: Modèle nonlinéaire à processus l
- Page 85 and 86: Modèle nonlinéaire à processus l
- Page 87 and 88: Modèle nonlinéaire à processus l
- Page 89 and 90: Modèle nonlinéaire à processus l
- Page 91 and 92: Modèle nonlinéaire à processus l
- Page 93 and 94: Chapitre 4Modèle nonlinéaire à c
- Page 95 and 96: Modèle nonlinéaire à classes lat
- Page 97 and 98: Modèle nonlinéaire à classes lat
- Page 99 and 100: Modèle nonlinéaire à classes lat
- Page 101 and 102: Modèle nonlinéaire à classes lat
- Page 103 and 104: Modèle nonlinéaire à classes lat
- Page 105 and 106: Modèle nonlinéaire à classes lat
- Page 107 and 108: Modèle nonlinéaire à classes lat
- Page 109 and 110: Modèle nonlinéaire à classes lat
- Page 111 and 112: Modèle nonlinéaire à classes lat
- Page 113 and 114: Modèle nonlinéaire à classes lat
- Page 115 and 116: Modèle nonlinéaire à classes lat
- Page 117 and 118: Modèle nonlinéaire à classes lat
- Page 119 and 120: Modèle nonlinéaire à classes lat
- Page 121 and 122: Modèle nonlinéaire à classes lat
- Page 123 and 124: Modèle nonlinéaire à classes lat
- Page 125 and 126:
Modèle nonlinéaire à classes lat
- Page 127 and 128:
Modèle nonlinéaire à classes lat
- Page 129 and 130:
Modèle nonlinéaire à classes lat
- Page 131 and 132:
Modèle nonlinéaire à classes lat
- Page 133 and 134:
Modèle nonlinéaire à classes lat
- Page 135 and 136:
Modèle nonlinéaire à classes lat
- Page 137 and 138:
Modèle nonlinéaire à classes lat
- Page 139 and 140:
Modèle nonlinéaire à classes lat
- Page 141 and 142:
Modèle nonlinéaire à classes lat
- Page 143 and 144:
Modèle nonlinéaire à classes lat
- Page 145 and 146:
Modèle nonlinéaire à classes lat
- Page 147 and 148:
Modèle nonlinéaire à classes lat
- Page 149 and 150:
Chapitre 6Discussion et perspective
- Page 151 and 152:
Discussion et perspectives 151la pr
- Page 153 and 154:
Discussion et perspectives 153Varia
- Page 155 and 156:
Discussion et perspectives 1556.2 M
- Page 157 and 158:
Discussion et perspectives 157à cl
- Page 160 and 161:
Discussion et perspectives 160effet
- Page 162 and 163:
Discussion et perspectives 162Dans
- Page 164 and 165:
Chapitre 7BibliographieAmieva, H.,
- Page 166 and 167:
Bibliographie 1661, S19-25.Brown, E
- Page 168 and 169:
Bibliographie 168Sons, New-York.Fol
- Page 170 and 171:
Bibliographie 170Hogan, J. W. et La
- Page 172 and 173:
Bibliographie 172Lin, H., McCulloch
- Page 174 and 175:
Bibliographie 174Park, CA.Muthén,
- Page 176 and 177:
Bibliographie 176Schlattmann, P. (2
- Page 178 and 179:
Bibliographie 178in the random-effe
- Page 180 and 181:
Chapitre 8Annexes8.1 Liste des publ
- Page 182 and 183:
Annexes 182Pau (France)Communicatio
- Page 184 and 185:
Annexes 184166 C. Proust, H. Jacqmi
- Page 186 and 187:
Annexes 186168 C. Proust, H. Jacqmi
- Page 188 and 189:
Annexes 188170 C. Proust, H. Jacqmi
- Page 190 and 191:
Annexes 190172 C. Proust, H. Jacqmi
- Page 192 and 193:
AnnexesAmerican Journal of Epidemio
- Page 194 and 195:
Annexes 194Psychometric Tests’ Se
- Page 196 and 197:
Annexes 196Psychometric Tests’ Se
- Page 198 and 199:
Annexes 198Psychometric Tests’ Se
- Page 200 and 201:
Annexes 200Abstract :When investiga
- Page 202 and 203:
Annexes 202educated subjects. This
- Page 204 and 205:
Annexes 204(iii) the recognition fo
- Page 206 and 207:
Annexes 206Explanatory variablesIn
- Page 208 and 209:
Annexes 208IST15 (median=28, IQR=24
- Page 210 and 211:
Annexes 210and on the mean evolutio
- Page 212 and 213:
Annexes 212These findings should be
- Page 214 and 215:
Annexes 214Appendix : model specifi
- Page 216 and 217:
Annexes 216References[1] Amieva H,
- Page 218 and 219:
Annexes 218[19] Letenneur L, Commen
- Page 220 and 221:
Annexes 220Figure 1 : (A) Predicted
- Page 222 and 223:
Annexes 222Table 1: demographic and
- Page 224:
Annexes 224MMSE 0.0037 * 0.0013,0.0