Robust linear mixed models with normal/independent distributions and Bayesian MCMC implementation

Linear mixed effects models have been widely used in analysis of data where responses are clustered around some random effects, so it is not reasonable to assume independence between observations in the same cluster. In most biological applications, it is assumed that the distributions of the random effects and of the residuals are Gaussian. This makes inferences vulnerable to the presence of outliers. Here, linear mixed effects models with normal/independent residual distributions for robust inferences are described. Specific distributions examined include univariate and multivariate versions of the Student-t, the slash and the contaminated normal. A Bayesian framework is adopted and Markov chain Monte Carlo is used to carry out the posterior analysis. The procedures are illustrated using birth weight data on rats in a texicological experiment. Results from the Gaussian and robust models are contrasted, and it is shown how the implementation can be used for outlier detection. The thick-tailed distributions provide an appealing robust alternative to the Gaussian process in linear mixed models, and they are easily implemented using data augmentation and MCMC techniques.

Approximations of the Fisher information for the construction of efficient experimental designs in nonlinear mixed effects models

Magdeburg, Univ., Fak. für Mathematik, Diss., 2012

The utility and application of mixed-effects models in second language research

Second language acquisition researchers often face particular challenges when attempting to generalize study findings to the wider learner population. For example, language learners constitute a heterogeneous group, and it is not always clear how a study’s findings may generalize to other individuals who may differ in terms of language background and proficiency, among many other factors. In this paper, we provide an overview of how mixed-effects models can be used to help overcome these and other issues in the field of second language acquisition. We provide an overview of the benefits of mixed-effects models and a practical example of how mixed-effects analyses can be conducted. Mixed-effects models provide second language researchers with a powerful statistical tool in the analysis of a variety of different types of data.

An overview of mixed-effects statistical models for second language researchers

As in any field of scientific inquiry, advancements in the field of second language acquisition (SLA) rely in part on the interpretation and generalizability of study findings using quantitative data analysis and inferential statistics. While statistical techniques such as ANOVA and t-tests are widely used in second language research, this review article provides a review of a class of newer statistical models that have not yet been widely adopted in the field, but have garnered interest in other fields of language research. The class of statistical models called mixed-effects models are introduced, and the potential benefits of these models for the second language researcher are discussed. A simple example of mixed-effects data analysis using the statistical software package R (R Development Core Team, 2011) is provided as an introduction to the use of these statistical techniques, and to exemplify how such analyses can be reported in research articles. It is concluded that mixed-effects models provide the second language researcher with a powerful tool for the analysis of a variety of types of second language acquisition data.

Bayesian longitudinal data analysis with mixed models and thick-tailed distributions using MCMC

Linear mixed effects models are frequently used to analyse longitudinal data, due to their flexibility in modelling the covariance structure between and within observations. Further, it is easy to deal with unbalanced data, either with respect to the number of observations per subject or per time period, and with varying time intervals between observations. In most applications of mixed models to biological sciences, a normal distribution is assumed both for the random effects and for the residuals. This, however, makes inferences vulnerable to the presence of outliers. Here, linear mixed models employing thick-tailed distributions for robust inferences in longitudinal data analysis are described. Specific distributions discussed include the Student-t, the slash and the contaminated normal. A Bayesian framework is adopted, and the Gibbs sampler and the Metropolis-Hastings algorithms are used to carry out the posterior analyses. An example with data on orthodontic distance growth in children is discussed to illustrate the methodology. Analyses based on either the Student-t distribution or on the usual Gaussian assumption are contrasted. The thick-tailed distributions provide an appealing robust alternative to the Gaussian process for modelling distributions of the random effects and of residuals in linear mixed models, and the MCMC implementation allows the computations to be performed in a flexible manner.

Prevalence of mixed neuropathologies in autopsied older adults and associations with clinical disease progression

Thesis (Ph.D.)--University of Washington, 2016-06

Effects of introduced moose (Alces alces) on vegetation composition, nutrient dynamics, and decomposition rates in boreal forest ecosystems in Newfoundland, Canada

Globally, consumers affect ecosystem processes including nutrient dynamics. Herbivores have been known to slow nutrient flow in boreal forest ecosystems. I examined the effects of introduced moose on disturbed forests of Newfoundland, Canada by conducting a field experiment during August - November 2014 in 20 paired moose exclosure-control plots. I tested whether moose browsing directly and indirectly affected forests by measuring plant species composition, litter quality and quantity, soil quality, and decomposition rates in areas moose exclosure-control plots. I analyzed moose effects using linear mixed effects models and found evidence indicating that moose reduce plant height and litter biomass affecting the availability of carbon, nitrogen, and phosphorus. However, plant diversity, soil quality, and litter decomposition did not differ between moose exclosures and controls. Moose in Newfoundland directly influence plant regeneration and litter biomass while indirect effects on soil ecosystems may be limited by time, disturbance, and climate.

Effect of the exposure to maternal smoking during pregnancy and childhood on the body mass index until adolescence

OBJECTIVE Investigate the effect of exposure to smoking during pregnancy and early childhood on changes in the body mass index (BMI) from birth to adolescence.METHODS A population-based cohort of children (0-5 years old) from Cuiabá, Midwest Brazil, was assessed in 1999-2000 (n = 2,405). Between 2009 and 2011, the cohort was re-evaluated. Information about birth weight was obtained from medical records, and exposure to smoking during pregnancy and childhood was assessed at the first interview. Linear mixed effects models were used to estimate the association between exposure to maternal smoking during pregnancy and preschool age, and the body mass index of children at birth, childhood and adolescence.RESULTS Only 11.3% of the mothers reported smoking during pregnancy, but most of them (78.2%) also smoked during early childhood. Among mothers who smoked only during pregnancy (n = 59), 97.7% had smoked only in the first trimester. The changes in body mass index at birth and in childhood were similar for children exposed and those not exposed to maternal smoking. However, from childhood to adolescence the rate of change in the body mass index was higher among those exposed only during pregnancy than among those who were not exposed.CONCLUSIONS Exposure to smoking only during pregnancy, especially in the first trimester, seems to affect changes in the body mass index until adolescence, supporting guidelines that recommend women of childbearing age to stop smoking.

Blood Pressure Variation Throughout Pregnancy According to Early Gestational BMI: A Brazilian Cohort

Background: The maternal cardiovascular system undergoes progressive adaptations throughout pregnancy, causing blood pressure fluctuations. However, no consensus has been established on its normal variation in uncomplicated pregnancies. Objective: To describe the variation in systolic blood pressure (SBP) and diastolic blood pressure (DBP) levels during pregnancy according to early pregnancy body mass index (BMI). Methods: SBP and DBP were measured during the first, second and third trimesters and at 30-45 days postpartum in a prospective cohort of 189 women aged 20-40 years. BMI (kg/m2) was measured up to the 13th gestational week and classified as normal-weight (<25.0) or excessive weight (≥25.0). Longitudinal linear mixed-effects models were used for statistical analysis. Results: A decrease in SBP and DBP was observed from the first to the second trimester (βSBP=-0.394; 95%CI: -0.600- -0.188 and βDBP=-0.617; 95%CI: -0.780- -0.454), as was an increase in SBP and DBP up to 30-45 postpartum days (βSBP=0.010; 95%CI: 0.006-0.014 and βDBP=0.015; 95%CI: 0.012-0.018). Women with excessive weight at early pregnancy showed higher mean SBP in all gestational trimesters, and higher mean DBP in the first and third trimesters. Excessive early pregnancy BMI was positively associated with prospective changes in SBP (βSBP=7.055; 95%CI: 4.499-9.610) and in DBP (βDBP=3.201; 95%CI: 1.136-5.266). Conclusion: SBP and DBP decreased from the first to the second trimester and then increased up to the postpartum period. Women with excessive early pregnancy BMI had higher SBP and DBP than their normal-weight counterparts throughout pregnancy, but not in the postpartum period.

Clinical pharmacokinetics of mycophenolic acid and its metabolites in solid organ transplant recipient

Mycophenolate mofetil (MMF), an ester prodrug of the immunosuppressant mycophenolic acid (MPA), is widely used for maintenance immunosuppressive therapy and prevention of renal allograft rejection in renal transplant recipients.MPA inhibits inosine monophosphate dehydrogenase (IMPDH), an enzyme involved in the “de novo” synthesis of purine nucleotides, thus suppressing both T-cell and B-cell proliferation. MPA shows a complex pharmacokinetics with considerable interand intra- patient by between- and within patient variabilities associated to MPA exposure. Several factors may contribute to it. The pharmacokinetic modeling according to the population pharmacokinetic approach with the non-linear mixed effects models has shown to be a powerful tool to describe the relationships between MMF doses and the MPA exposures and also to identify potential predictive patients’ demographic and clinical characteristics for dose tailoring during the post-transplant immunosuppresive treatment.

Clinical pharmacokinetics of mycophenolic acid and its metabolites in solid organ transplant recipient

Mycophenolate mofetil (MMF), an ester prodrug of the immunosuppressant mycophenolic acid (MPA), is widely used for maintenance immunosuppressive therapy and prevention of renal allograft rejection in renal transplant recipients.MPA inhibits inosine monophosphate dehydrogenase (IMPDH), an enzyme involved in the “de novo” synthesis of purine nucleotides, thus suppressing both T-cell and B-cell proliferation. MPA shows a complex pharmacokinetics with considerable interand intra- patient by between- and within patient variabilities associated to MPA exposure. Several factors may contribute to it. The pharmacokinetic modeling according to the population pharmacokinetic approach with the non-linear mixed effects models has shown to be a powerful tool to describe the relationships between MMF doses and the MPA exposures and also to identify potential predictive patients’ demographic and clinical characteristics for dose tailoring during the post-transplant immunosuppresive treatment.

Clinical pharmacokinetics of mycophenolic acid and its metabolites in solid organ transplant recipient

Mycophenolate mofetil (MMF), an ester prodrug of the immunosuppressant mycophenolic acid (MPA), is widely used for maintenance immunosuppressive therapy and prevention of renal allograft rejection in renal transplant recipients.MPA inhibits inosine monophosphate dehydrogenase (IMPDH), an enzyme involved in the “de novo” synthesis of purine nucleotides, thus suppressing both T-cell and B-cell proliferation. MPA shows a complex pharmacokinetics with considerable interand intra- patient by between- and within patient variabilities associated to MPA exposure. Several factors may contribute to it. The pharmacokinetic modeling according to the population pharmacokinetic approach with the non-linear mixed effects models has shown to be a powerful tool to describe the relationships between MMF doses and the MPA exposures and also to identify potential predictive patients’ demographic and clinical characteristics for dose tailoring during the post-transplant immunosuppresive treatment.

Clinical pharmacokinetics of mycophenolic acid and its metabolites in solid organ transplant recipient

Mycophenolate mofetil (MMF), an ester prodrug of the immunosuppressant mycophenolic acid (MPA), is widely used for maintenance immunosuppressive therapy and prevention of renal allograft rejection in renal transplant recipients.MPA inhibits inosine monophosphate dehydrogenase (IMPDH), an enzyme involved in the “de novo” synthesis of purine nucleotides, thus suppressing both T-cell and B-cell proliferation. MPA shows a complex pharmacokinetics with considerable interand intra- patient by between- and within patient variabilities associated to MPA exposure. Several factors may contribute to it. The pharmacokinetic modeling according to the population pharmacokinetic approach with the non-linear mixed effects models has shown to be a powerful tool to describe the relationships between MMF doses and the MPA exposures and also to identify potential predictive patients’ demographic and clinical characteristics for dose tailoring during the post-transplant immunosuppresive treatment.