Nonlinear Mixed-Effects Model Development and Simulation Using nlmixr and Related R Open-Source Packages.
Ontology highlight
ABSTRACT:
SUBMITTER: Fidler M
PROVIDER: S-EPMC6765694 | biostudies-literature | 2019 Sep
REPOSITORIES: biostudies-literature
Publications
Nonlinear Mixed-Effects Model Development and Simulation Using nlmixr and Related R Open-Source Packages.
Fidler Matthew M Wilkins Justin J JJ Hooijmaijers Richard R Post Teun M TM Schoemaker Rik R Trame Mirjam N MN Xiong Yuan Y Wang Wenping W
CPT: pharmacometrics & systems pharmacology 20190716 9
nlmixr is a free and open-source R package for fitting nonlinear pharmacokinetic (PK), pharmacodynamic (PD), joint PK-PD, and quantitative systems pharmacology mixed-effects models. Currently, nlmixr is capable of fitting both traditional compartmental PK models as well as more complex models implemented using ordinary differential equations. We believe that, over time, it will become a capable, credible alternative to commercial software tools, such as NONMEM, Monolix, and Phoenix NLME. ...[more]
Similar Datasets
Project description:A nonlinear mixed-effects approach is developed for disease progression models that incorporate variation in age in a Bayesian framework. We further generalize the probability model for sensitivity to depend on age at diagnosis, time spent in the preclinical state and sojourn time. The developed models are then applied to the Johns Hopkins Lung Project data and the Health Insurance Plan for Greater New York data using Bayesian Markov chain Monte Carlo and are compared with the estimation method that does not consider random-effects from age. Using the developed models, we obtain not only age-specific individual-level distributions, but also population-level distributions of sensitivity, sojourn time and transition probability.
| S-EPMC4718583 | biostudies-literature
Project description:Rift Valley fever (RVF) is one of the major viral zoonoses in Africa, affecting humans and several domestic animal species. The epidemics in eastern Africa occur in a 5-15 year cycle coinciding with abnormally high rainfall generally associated to the warm phase of the El Niño event. However, recently, evidence has been gathered of inter-epidemic transmission. An open-source, easily applicable, accessible and modifiable model was built to simulate the transmission dynamics of RVF. The model was calibrated using data collected in the Kilombero Valley in Tanzania with people and cattle as host species and Ædes mcintoshi, Æ. ægypti and two Culex species as vectors. Simulations were run over a period of 27 years using standard parameter values derived from two previous studies in this region. Our model predicts low-level transmission of RVF, which is in line with epidemiological studies in this area. Emphasis in our simulation was put on both the dynamics and composition of vector populations in three ecological zones, in order to elucidate the respective roles played by different vector species: the model output did indicate the necessity of Culex involvement and also indicated that vertical transmission in Ædes mcintoshi may be underestimated. This model, being built with open-source software and with an easy-to-use interface, can be adapted by researchers and control program managers to their specific needs by plugging in new parameters relevant to their situation and locality.
| S-EPMC6326482 | biostudies-literature
Project description:Bioturbation is one of the most widespread forms of ecological engineering and has significant implications for the structure and functioning of ecosystems, yet our understanding of the processes involved in biotic mixing remains incomplete. One reason is that, despite their value and utility, most mathematical models currently applied to bioturbation data tend to neglect aspects of the natural complexity of bioturbation in favour of mathematical simplicity. At the same time, the abstract nature of these approaches limits the application of such models to a limited range of users. Here, we contend that a movement towards process-based modelling can improve both the representation of the mechanistic basis of bioturbation and the intuitiveness of modelling approaches. In support of this initiative, we present an open source modelling framework that explicitly simulates particle displacement and a worked example to facilitate application and further development. The framework combines the advantages of rule-based lattice models with the application of parameterisable probability density functions to generate mixing on the lattice. Model parameters can be fitted by experimental data and describe particle displacement at the spatial and temporal scales at which bioturbation data is routinely collected. By using the same model structure across species, but generating species-specific parameters, a generic understanding of species-specific bioturbation behaviour can be achieved. An application to a case study and comparison with a commonly used model attest the predictive power of the approach.
| S-EPMC3230619 | biostudies-literature
Project description:We present a Bayesian nonlinear mixed-effects location scale model (NL-MELSM). The NL-MELSM allows for fitting nonlinear functions to the location, or individual means, and the scale, or within-person variance. Specifically, in the context of learning, this model allows the within-person variance to follow a nonlinear trajectory, where it can be determined whether variability reduces during learning. It incorporates a sub-model that can predict nonlinear parameters for both the location and scale. This specification estimates random effects for all nonlinear location and scale parameters that are drawn from a common multivariate distribution. This allows estimation of covariances among the random effects, within and across the location and the scale. These covariances offer new insights into the interplay between individual mean structures and intra-individual variability in nonlinear parameters. We take a fully Bayesian approach, not only for ease of estimation but also for inference because it provides the necessary and consistent information for use in psychological applications, such as model selection and hypothesis testing. To illustrate the model, we use data from 333 individuals, consisting of three age groups, who participated in five learning trials that assessed verbal memory. In an exploratory context, we demonstrate that fitting a nonlinear function to the within-person variance, and allowing for individual variation therein, improves predictive accuracy compared to customary modeling techniques (e.g., assuming constant variance). We conclude by discussing the usefulness, limitations, and future directions of the NL-MELSM.
| S-EPMC6800615 | biostudies-literature
Project description:The free and open-source package nlmixr implements pharmacometric nonlinear mixed effects model parameter estimation in R. It provides a uniform language to define pharmacometric models using ordinary differential equations. Performances of the stochastic approximation expectation-maximization (SAEM) and first order-conditional estimation with interaction (FOCEI) algorithms in nlmixr were compared with those found in the industry standards, Monolix and NONMEM, using the following two scenarios: a simple model fit to 500 sparsely sampled data sets and a range of more complex compartmental models with linear and nonlinear clearance fit to data sets with rich sampling. Estimation results obtained from nlmixr for FOCEI and SAEM matched the corresponding output from NONMEM/FOCEI and Monolix/SAEM closely both in terms of parameter estimates and associated standard errors. These results indicate that nlmixr may provide a viable alternative to existing tools for pharmacometric parameter estimation.
| S-EPMC6930853 | biostudies-literature
Project description:Numerical body models of children are used for designing medical devices, including but not limited to optical imaging, ultrasound, CT, EEG/MEG, and MRI. These models are used in many clinical and neuroscience research applications, such as radiation safety dosimetric studies and source localization. Although several such adult models have been reported, there are few reports of full-body pediatric models, and those described have several limitations. Some, for example, are either morphed from older children or do not have detailed segmentations. Here, we introduce a 29-month-old male whole-body native numerical model, "MARTIN", that includes 28 head and 86 body tissue compartments, segmented directly from the high spatial resolution MRI and CT images. An advanced auto-segmentation tool was used for the deep-brain structures, whereas 3D Slicer was used to segment the non-brain structures and to refine the segmentation for all of the tissue compartments. Our MARTIN model was developed and validated using three separate approaches, through an iterative process, as follows. First, the calculated volumes, weights, and dimensions of selected structures were adjusted and confirmed to be within 6% of the literature values for the 2-3-year-old age-range. Second, all structural segmentations were adjusted and confirmed by two experienced, sub-specialty certified neuro-radiologists, also through an interactive process. Third, an additional validation was performed with a Bloch simulator to create synthetic MR image from our MARTIN model and compare the image contrast of the resulting synthetic image with that of the original MRI data; this resulted in a "structural resemblance" index of 0.97. Finally, we used our model to perform pilot MRI safety simulations of an Active Implantable Medical Device (AIMD) using a commercially available software platform (Sim4Life), incorporating the latest International Standards Organization guidelines. This model will be made available on the Athinoula A. Martinos Center for Biomedical Imaging website.
| S-EPMC7806143 | biostudies-literature
Project description:We introduce June, an open-source framework for the detailed simulation of epidemics on the basis of social interactions in a virtual population constructed from geographically granular census data, reflecting age, sex, ethnicity and socio-economic indicators. Interactions between individuals are modelled in groups of various sizes and properties, such as households, schools and workplaces, and other social activities using social mixing matrices. June provides a suite of flexible parametrizations that describe infectious diseases, how they are transmitted and affect contaminated individuals. In this paper, we apply June to the specific case of modelling the spread of COVID-19 in England. We discuss the quality of initial model outputs which reproduce reported hospital admission and mortality statistics at national and regional levels as well as by age strata.
| S-EPMC8261230 | biostudies-literature
Project description:Virtual representations of the knee joint can provide clinicians, scientists, and engineers the tools to explore mechanical functions of the knee and its tissue structures in health and disease. Modeling and simulation approaches such as finite element analysis also provide the possibility to understand the influence of surgical procedures and implants on joint stresses and tissue deformations. A large number of knee joint models are described in the biomechanics literature. However, freely accessible, customizable, and easy-to-use models are scarce. Availability of such models can accelerate clinical translation of simulations, where labor-intensive reproduction of model development steps can be avoided. Interested parties can immediately utilize readily available models for scientific discovery and clinical care. Motivated by this gap, this study aims to describe an open source and freely available finite element representation of the tibiofemoral joint, namely Open Knee, which includes the detailed anatomical representation of the joint's major tissue structures and their nonlinear mechanical properties and interactions. Three use cases illustrate customization potential of the model, its predictive capacity, and its scientific and clinical utility: prediction of joint movements during passive flexion, examining the role of meniscectomy on contact mechanics and joint movements, and understanding anterior cruciate ligament mechanics. A summary of scientific and clinically directed studies conducted by other investigators are also provided. The utilization of this open source model by groups other than its developers emphasizes the premise of model sharing as an accelerator of simulation-based medicine. Finally, the imminent need to develop next-generation knee models is noted. These are anticipated to incorporate individualized anatomy and tissue properties supported by specimen-specific joint mechanics data for evaluation, all acquired in vitro from varying age groups and pathological states.
| S-EPMC4876308 | biostudies-literature
Project description:Precise modeling of disease progression in neurodegenerative disorders may enable early intervention before clinical manifestation of a disease, which is crucial since early intervention at the premanifest stage is expected to be more effective. Neuroimaging biomarkers are indicative of the underlying disease pathology and may be used to predict future disease occurrence at the premanifest stage. As observed in many pivotal studies, longitudinal measurements of clinical outcomes, such as motor or cognitive symptoms, often present nonlinear sigmoid shapes over time, where the inflection points of the trajectories mark a meaningful time in disease progression. Therefore, to identify neuroimaging biomarkers predicting disease progression, we propose a nonlinear mixed effects model based on a sigmoid function to predict longitudinal clinical outcomes, and associate a linear combination of neuroimaging biomarkers with subject-specific inflection points. Based on an expectation-maximization (EM) algorithm, we propose a method that can fit a nonlinear model with many potentially correlated biomarkers for random inflection points while achieving computational stability. Variable selection is introduced in the algorithm in order to identify important biomarkers of disease progression and to reduce prediction variability. We apply the proposed method to the data from the Predictors of Huntington's Disease study to select brain subcortical regional volumes predictive of the inflection points of the motor and cognitive function trajectories. Our results reveal that brain atrophy in the striatum and expansion of the ventricular system are highly predictive of the inflection points. Furthermore, these inflection points may precede clinically defined disease onset by as early as a decade and thus may be useful biomarkers as early signs of Huntington's Disease onset.
| S-EPMC5790569 | biostudies-literature
Project description:Owing to the increase in freely available software and data for cheminformatics and structural bioinformatics, research for computer-aided drug design (CADD) is more and more built on modular, reproducible, and easy-to-share pipelines. While documentation for such tools is available, there are only a few freely accessible examples that teach the underlying concepts focused on CADD, especially addressing users new to the field. Here, we present TeachOpenCADD, a teaching platform developed by students for students, using open source compound and protein data as well as basic and CADD-related Python packages. We provide interactive Jupyter notebooks for central CADD topics, integrating theoretical background and practical code. TeachOpenCADD is freely available on GitHub: https://github.com/volkamerlab/TeachOpenCADD .
| S-EPMC6454689 | biostudies-literature