A SAEM algorithm for fused Lasso penalized nonlinear mixed effect models: application to group comparison in pharmacokinetics. (English) Zbl 1468.62152
Summary: Nonlinear mixed effect models are classical tools to analyze nonlinear longitudinal data in many fields such as population pharmacokinetics. Groups of observations are usually compared by introducing the group affiliations as binary covariates with a reference group that is stated among the groups. This approach is relatively limited as it allows only the comparison of the reference group to the others. The proposed method compares groups using a penalized likelihood approach. Groups are described by the same structural model but with parameters that are group specific. The likelihood is penalized with a fused lasso penalty that induces sparsity in the differences between groups for both fixed effects and variances of random effects. A penalized Stochastic Approximation EM algorithm is proposed that is coupled to Alternating Direction Method Multipliers to solve the maximization step. An extensive simulation study illustrates the performance of this algorithm when comparing more than two groups. Then the approach is applied to real data from two pharmacokinetic drug-drug interaction trials.
MSC:
| 62-08 | Computational methods for problems pertaining to statistics |
| 62P10 | Applications of statistics to biology and medical sciences; meta analysis |
| 62J07 | Ridge regression; shrinkage estimators (Lasso) |
Keywords:
nonlinear mixed effect model; SAEM algorithm; fused Lasso; group comparison; pharmacokineticsReferences:
| [1] | Arribas-Gil, A.; Bertin, K.; Meza, C.; Rivoirard, V., Lasso-type estimators for semiparametric nonlinear mixed-effects models estimation, Stat. Comput., 24, 3, 443-460, (2014) · Zbl 1325.62070 |
| [2] | Atchade, Y.F., Fort, G., Moulines, E., 2014. On stochastic proximal gradient algorithms. ArXiv Preprint arXiv:1402.2365. · Zbl 1433.90199 |
| [3] | Bach, F.; Jenatton, R.; Mairal, J.; Obozinski, G., Convex optimization with sparsity-inducing norms, (Optimization for Machine Learning, (2011)), 19-53 |
| [4] | Bell, E. T., Exponential numbers, Amer. Math. Monthly, 411-419, (1934) · JFM 60.0116.01 |
| [5] | Bertrand, J.; De Iorio, M.; Balding, D. J., Integrating dynamic mixed-effect modelling and penalized regression to explore genetic association with pharmacokinetics, Pharmacogenetics Genom., 25, 5, 231, (2015) |
| [6] | Bien, J.; Tibshirani, R. J., Sparse estimation of a covariance matrix, Biometrika, 98, 4, 807-820, (2011) · Zbl 1228.62063 |
| [7] | Blech, S.; Ebner, T.; Ludwig-Schwellinger, E.; Stangier, J.; Roth, W., The metabolism and disposition of the oral direct thrombin inhibitor, dabigatran, in humans, Drug Metab. Dispos., 36, 2, 386-399, (2008) |
| [8] | Bondell, H. D.; Krishna, A.; Ghosh, S. K., Joint variable selection for fixed and random effects in linear mixed-effects models, Biometrics, 66, 4, 1069-1077, (2010) · Zbl 1233.62134 |
| [9] | Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J., Distributed optimization and statistical learning via the alternating direction method of multipliers, Found. Trends Mach. Learn., 3, 1, 1-122, (2011) · Zbl 1229.90122 |
| [10] | Colby, E., Bair, E., 2013. Cross-validation for nonlinear mixed effects models, Cross-Validation for Nonlinear Mixed Effects Models. |
| [11] | Danaher, P.; Wang, P.; Witten, D. M., The joint graphical lasso for inverse covariance estimation across multiple classes, J. R. Stat. Soc. Ser. B, (2013) |
| [12] | Delattre, M.; Lavielle, M.; Poursat, M.-A., A note on BIC in mixed-effects models, Electron. J. Stat., 8, 456-475, (2014) · Zbl 1348.62186 |
| [13] | Delavenne, X.; Ollier, E.; Basset, T.; Bertoletti, L.; Accassat, S.; Garcin, A.; Laporte, S.; Zufferey, P.; Mismetti, P., A semi-mechanistic absorption model to evaluate drug-drug interaction with dabigatran: application with clarithromycin, Br. J. Clin. Pharmacol., 76, 1, 107-113, (2013) |
| [14] | Delyon, B.; Lavielle, M.; Moulines, E., Convergence of a stochastic approximation version of the EM algorithm, Ann. Statist., 94-128, (1999) · Zbl 0932.62094 |
| [15] | Dubois, A.; Lavielle, M.; Gsteiger, S.; Pigeolet, E.; Mentré, F., Model-based analyses of bioequivalence crossover trials using the stochastic approximation expectation maximisation algorithm, Stat. Med., 30, 21, 2582-2600, (2011) |
| [16] | Efron, B.; Hastie, T.; Johnstone, I.; Tibshirani, R., Least angle regression, Ann. Statist., 32, 2, 407-499, (2004) · Zbl 1091.62054 |
| [17] | Höfling, H., Binder, H., Schumacher, M., 2010. A coordinate-wise optimization algorithm for the fused lasso. ArXiv Preprint arXiv:1011.6409. |
| [18] | Kuhn, E.; Lavielle, M., Maximum likelihood estimation in nonlinear mixed effects models, Comput. Statist. Data Anal., 49, 4, 1020-1038, (2005) · Zbl 1429.62279 |
| [19] | Meinshausen, N., Relaxed lasso, Comput. Statist. Data Anal., 52, 1, 374-393, (2007) · Zbl 1452.62522 |
| [20] | Oelker, M.-R.; Gertheiss, J.; Tutz, G., Regularization and model selection with categorical predictors and effect modifiers in generalized linear models, Stat. Model., 14, 2, 157-177, (2014) · Zbl 07257900 |
| [21] | Ollier, E.; Hodin, S.; Basset, T.; Accassat, S.; Bertoletti, L.; Mismetti, P.; Delavenne, X., In vitro and in vivo evaluation of drug-drug interaction between dabigatran and proton pump inhibitors, Fundam. Clin. Pharmacol., (2015) |
| [22] | Ollier, E., Viallon, V., 2015. Regression modeling on stratified data: automatic and covariate-specific selection of the reference stratum with simple l_1-norm penalties. ArXiv Preprint arXiv:1508.05476. |
| [23] | Panhard, X.; Samson, A., Extension of the SAEM algorithm for nonlinear mixed models with 2 levels of random effects, Biostatistics, 10, 1, 121-135, (2009) · Zbl 1437.62573 |
| [24] | Rohart, F.; San Cristobal, M.; Laurent, B., Selection of fixed effects in high dimensional linear mixed models using a multicycle ECM algorithm, Comput. Statist. Data Anal., (2014) · Zbl 1506.62156 |
| [25] | Samson, A.; Lavielle, M.; Mentré, F., The SAEM algorithm for group comparison tests in longitudinal data analysis based on non-linear mixed-effects model, Stat. Med., 26, 27, 4860-4875, (2007) |
| [26] | Schelldorfer, J.; Buhlmann, P.; De Geer, S., Estimation for high-dimensional linear mixed-effects models using 1-penalization, Scand. J. Statist., 38, 197-214, (2011) · Zbl 1246.62161 |
| [27] | Tibshirani, R.; Saunders, M.; Rosset, S.; Zhu, J.; Knight, K., Sparsity and smoothness via the fused lasso, J. R. Stat. Soc. Ser. B Stat. Methodol., 67, 1, 91-108, (2005) · Zbl 1060.62049 |
| [28] | Viallon, V.; Lambert-Lacroix, S.; Hoefling, H.; Picard, F., On the robustness of the generalized fused lasso to prior specifications, Stat. Comput., 1-17, (2014) |
| [29] | Wang, H., Coordinate descent algorithm for covariance graphical lasso, Stat. Comput., 1-9, (2013) |
| [30] | Weiss, M.; Sermsappasuk, P.; Siegmund, W., Modeling the kinetics of digoxin absorption: enhancement by p-glycoprotein inhibition, J. Clin. Pharmacol., 52, 3, 381-387, (2012) |
| [31] | Witten, D. M.; Tibshirani, R., Covariance-regularized regression and classification for high dimensional problems, J. R. Stat. Soc. Ser. B Stat. Methodol., 71, 3, 615-636, (2009) · Zbl 1250.62033 |
| [32] | Zou, H., The adaptive lasso and its oracle properties, J. Amer. Statist. Assoc., 101, 476, 1418-1429, (2006) · Zbl 1171.62326 |
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. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
