Mixtures of stochastic differential equations with random effects: Application to data clustering

We consider N independent stochastic processes (X-i(t), t is an element of [0, T-i]), i = 1,..., N, defined by a stochastic differential equation with drift term depending on a random variable phi(i). The distribution of the random effect phi(i) is a Gaussian mixture distribution, depending on unknown parameters which are to be estimated from the continuous observation of the processes X-i. The likelihood of the observation is explicit. When the number of components is known, we prove the consistency of the exact maximum likelihood estimators and use the EM algorithm to compute it. When the number of components is unknown, BIC (Bayesian Information Criterion) is applied to select it. To assign each individual to a class, we define a classification rule based on estimated posterior probabilities. A simulation study illustrates our estimation and classification method on various models. A real data analysis is performed on growth curves with convincing results. (C) 2015 Elsevier B.V. All rights reserved.