Dose-response analysis can be carried out using multi-purpose commercial statistical software, but except for a few special cases the analysis easily becomes cumbersome as relevant, non-standard output requires manual programming. The extension package drc for the statistical environment R provides a flexible and versatile infrastructure for dose-response analyses in general. The present version of the package, reflecting extensions and modifications over the last decade, provides a user-friendly interface to specify the model assumptions about the dose-response relationship and comes with a number of extractors for summarizing fitted models and carrying out inference on derived parameters. The aim of the present paper is to provide an overview of state-of-the-art dose-response analysis, both in terms of general concepts that have evolved and matured over the years and by means of concrete examples.

/pdf/dose-response-analysis-using-r-1rjy8qez6o.pdf

Dose-Response Analysis Using R

Elucidating the factors that impinge on the stability of bacterial communities in the vagina may help in predicting the risk of diseases that affect women’s health. Here, we describe the temporal dynamics of the composition of vaginal bacterial communities in 32 reproductive-age women over a 16-week period. The analysis revealed the dynamics of five major classes of bacterial communities and showed that some communities change markedly over short time periods, whereas others are relatively stable. Modeling community stability using new quantitative measures indicates that deviation from stability correlates with time in the menstrual cycle, bacterial community composition, and sexual activity. The women studied are healthy; thus, it appears that neither variation in community composition per se nor higher levels of observed diversity (co-dominance) are necessarily indicative of dysbiosis.

Temporal Dynamics of the Human Vaginal Microbiota

Kwashiorkor, an enigmatic form of severe acute malnutrition, is the consequence of inadequate nutrient intake plus additional environmental insults. To investigate the role of the gut microbiome, we studied 317 Malawian twin pairs during the first 3 years of life. During this time, half of the twin pairs remained well nourished, whereas 43% became discordant, and 7% manifested concordance for acute malnutrition. Both children in twin pairs discordant for kwashiorkor were treated with a peanut-based, ready-to-use therapeutic food (RUTF). Time-series metagenomic studies revealed that RUTF produced a transient maturation of metabolic functions in kwashiorkor gut microbiomes that regressed when administration of RUTF was stopped. Previously frozen fecal communities from several discordant pairs were each transplanted into gnotobiotic mice. The combination of Malawian diet and kwashiorkor microbiome produced marked weight loss in recipient mice, accompanied by perturbations in amino acid, carbohydrate, and intermediary metabolism that were only transiently ameliorated with RUTF. These findings implicate the gut microbiome as a causal factor in kwashiorkor.

http://www.humeco.m.u-tokyo.ac.jp/individuals/umezaki/PDF_files/Meeting/Malawi_Gutmicrobiome_Kwashiorkor_Science2013.pdf

Gut microbiomes of Malawian twin pairs discordant for kwashiorkor

In mountainous regions, climate warming is expected to shift species' ranges to higher altitudes. Evidence for such shifts is still mostly from revisitations of historical sites. We present recent (2001 to 2008) changes in vascular plant species richness observed in a standardized monitoring network across Europe's major mountain ranges. Species have moved upslope on average. However, these shifts had opposite effects on the summit floras' species richness in boreal-temperate mountain regions (+3.9 species on average) and Mediterranean mountain regions (-1.4 species), probably because recent climatic trends have decreased the availability of water in the European south. Because Mediterranean mountains are particularly rich in endemic species, a continuation of these trends might shrink the European mountain flora, despite an average increase in summit species richness across the region.

https://digital.csic.es/bitstream/10261/50439/1/Villar_Benito-Science2012.pdf

Recent Plant Diversity Changes on Europe’s Mountain Summits

Determining the form of key predator-prey relationships is critical for understanding marine ecosystem dynamics. Using a comprehensive global database, we quantified the effect of fluctuations in food abundance on seabird breeding success. We identified a threshold in prey (fish and krill, termed “forage fish”) abundance below which seabirds experience consistently reduced and more variable productivity. This response was common to all seven ecosystems and 14 bird species examined within the Atlantic, Pacific, and Southern Oceans. The threshold approximated one-third of the maximum prey biomass observed in long-term studies. This provides an indicator of the minimal forage fish biomass needed to sustain seabird productivity over the long term.

/pdf/global-seabird-response-to-forage-fish-depletion-one-third-w1c0wdyr8g.pdf

Global Seabird Response to Forage Fish Depletion-One-Third for the Birds

Douglas M. Bates Department of Statistics University of Wisconsin Madison Jose C. Pinheiro Bell Laboratories Lucent Technologies 1 Recent developments in computational methods for maximum likelihood (ML) or restricted maximum likelihood (REML) estimation of parameters in general linear mixed-effects models have made the analysis of data in typical agricultural settings much easier. With software such as SAS PROC MIXED we are able to handle da~ from random-effects one-way classifications, from blocked designs including incomplete blocked designs, from hierarchical designs such as splitplot designs, and other types of data that may be described as repeated measures or longitudinal data or growth-curve data. It is especially helpful that the new computational methods do not depend on balance in the data so we are able to deal more easily with observational studies or with randomly missing data in a designed experiment. We describe some of the new computational approaches and how they are implemented in the nlme3.0 library for the S-PLUS language. One of the most powerful features of this language is the graphics capabilities, especially the trellis graphics facilities developed by Bill Cleveland and his coworkers at Bell Labs. Although most participants in this conference may be more familiar with SAS, and most of the models described here can be fit with PROC MIXED or the NLiNMIX macro or new P ROC N LM IXED in SAS version 7, some exposure to the combination of graphical display and model-fitting approaches from S-PLUS may be informative. 1 Annual Conference on Applied Statistics in Agriculture Kansas State University New Prairie Press http://newprairiepress.org/agstatconference/1998/proceedings/2 2 Kansas State University We show how data exploration with trellis graphics, followed by fitting and comparing mixedeffects models, followed by graphical assessment of the fitted model can be used in a variety of situations. On some occasions, such as modeling growth curves, a linear trend or polynomial trend or other types of linear statistical models for the within-subject time dependence are just not going to do an adequate job of representing the data. In those cases, a nonlinear model is more appropriate. We show how the concept of a random coefficient model can be extended to nonlinear models so as to fit nonlinear mixed-effects models.

R-core

Papers

Linear and nonlinear mixed-effects models

Linear and Nonlinear Mixed Effects Models [R package nlme version 3.1-149]

Linear and Nonlinear Mixed Effects Models [R package nlme version 3.1-152]

Linear and Nonlinear Mixed Effects Models [R package nlme version 3.1-151]