Performance in population models for count data, part II: a new SAEM algorithm.

Radojka Savic 1, * Marc Lavielle 2
* Auteur correspondant
2 SELECT - Model selection in statistical learning
Inria Saclay - Ile de France, LMO - Laboratoire de Mathématiques d'Orsay, CNRS - Centre National de la Recherche Scientifique : UMR
Abstract : Analysis of count data from clinical trials using mixed effect analysis has recently become widely used. However, algorithms available for the parameter estimation, including LAPLACE and Gaussian quadrature (GQ), are associated with certain limitations, including bias in parameter estimates and the long analysis runtime. The stochastic approximation expectation maximization (SAEM) algorithm has proven to be a very efficient and powerful tool in the analysis of continuous data. The aim of this study was to implement and investigate the performance of a new SAEM algorithm for application to count data. A new SAEM algorithm was implemented in MATLAB for estimation of both, parameters and the Fisher information matrix. Stochastic Monte Carlo simulations followed by re-estimation were performed according to scenarios used in previous studies (part I) to investigate properties of alternative algorithms (Plan et al., 2008, Abstr 1372 [ http://wwwpage-meetingorg/?abstract=1372 ]). A single scenario was used to explore six probability distribution models. For parameter estimation, the relative bias was less than 0.92% and 4.13% for fixed and random effects, for all models studied including ones accounting for over- or under-dispersion. Empirical and estimated relative standard errors were similar, with distance between them being <1.7% for all explored scenarios. The longest CPU time was 95 s for parameter estimation and 56 s for SE estimation. The SAEM algorithm was extended for analysis of count data. It provides accurate estimates of both, parameters and standard errors. The estimation is significantly faster compared to LAPLACE and GQ. The algorithm is implemented in Monolix 3.1, (beta-version available in July 2009).
Type de document :
Article dans une revue
Journal of Pharmacokinetics and Pharmacodynamics, Springer Verlag, 2009, 36 (4), pp.367-79. 〈10.1007/s10928-009-9127-7〉
Liste complète des métadonnées

Littérature citée [6 références]  Voir  Masquer  Télécharger

http://www.hal.inserm.fr/inserm-00470358
Contributeur : Thu Thuy Nguyen <>
Soumis le : mardi 6 avril 2010 - 13:02:34
Dernière modification le : jeudi 11 janvier 2018 - 06:22:14
Document(s) archivé(s) le : mardi 14 septembre 2010 - 16:58:01

Fichiers

Savic_JPKPD_2009.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

Collections

Citation

Radojka Savic, Marc Lavielle. Performance in population models for count data, part II: a new SAEM algorithm.. Journal of Pharmacokinetics and Pharmacodynamics, Springer Verlag, 2009, 36 (4), pp.367-79. 〈10.1007/s10928-009-9127-7〉. 〈inserm-00470358〉

Partager

Métriques

Consultations de la notice

205

Téléchargements de fichiers

728