Statistical inference for the optimal approximating model

In the setting of high-dimensional linear models with Gaussian noise, we investigate the possibility of confidence statements connected to model selection. Although there exist numerous procedures for adaptive (point) estimation, the construction of adaptive confidence regions is severely limited (cf. Li in Ann Stat 17:1001–1008, 1989). The present paper sheds new light on this gap. We develop exact and adaptive confidence regions for the best approximating model in terms of risk. One of our constructions is based on a multiscale procedure and a particular coupling argument. Utilizing exponential inequalities for noncentral χ2-distributions, we show that the risk and quadratic loss of all models within our confidence region are uniformly bounded by the minimal risk times a factor close to one.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Ecosystem service trade-offs and stakeholder claims - a meso-scale study from a Madagascar biodiversity hotspot

The north-eastern escarpment of Madagascar has been labelled a global biodiversity hotspot due to its extremely high rates of endemic species which are heavily threatened by accelerated deforestation rates and landscape change. The traditional practice of shifting cultivation or "tavy" used by the majority of land users in this area to produce subsistence rice is commonly blamed for these threats. A wide range of stakeholders ranging from conservation to development agencies, and from the private to the public sector has therefore been involved in trying to find solutions to protect the remaining forest fragments and to increase agricultural production. Consequently, provisioning, regulating and socio-cultural services of this forest-mosaic landscape are fundamentally altered leading to trade-offs between them and consequently new winners and losers amongst the stakeholders at different scales. However, despite a growing amount of evidence from case studies analysing local changes, the regional dynamics of the landscape and their contribution to such trade-offs remain poorely understood. This study therefore aims at using generalised landscape units as a base for the assessment of multi-level stakeholder claims on ecosystem services to inform negotiation, planning and decision making at a meso-scale. The presented study applies a mixed-method approach combining remote sensing, GIS and socio-economic methods to reveal current landscape dynamics, their change over time and the corresponding ecosystem service trade-offs induced by diverse stakeholder claims on the regional level. In a first step a new regional land cover classification for three points in time (1995, 2005 and 2011) was conducted including agricultural classes characteristic for shifting cultivation systems. Secondly, a novel GIS approach, termed “landscape mosaics approach” originally developed to assess dynamics of shifting cultivation landscapes in Laos was applied. Through this approach generalised landscape mosaics were generated allowing for a better understanding of changes in land use intensities instead of land cover. As a next step we will try to use these landscape units as proxies to map provisioning and regulating ecosystem services throughout the region. Through the overlay with other regional background data such as accessibility and population density and information from a region-wide stakeholder analysis, multiscale trade-offs between different services will be highlighted. The trade-offs observed on the regional scale will then be validated through a socio-economic ground-truthing within selected sites at the local scale. We propose that such meso-scale knowledge is required by all stakeholders involved in decision making towards sustainable development of north-eastern Madagascar.

MorphoGraphX: A platform for quantifying morphogenesis in 4D

Morphogenesis emerges from complex multiscale interactions between genetic and mechanical processes. To understand these processes, the evolution of cell shape, proliferation and gene expression must be quantified. This quantification is usually performed either in full 3D, which is computationally expensive and technically challenging, or on 2D planar projections, which introduces geometrical artifacts on highly curved organs. Here we present MorphoGraphX (www.MorphoGraphX.org), a software that bridges this gap by working directly with curved surface images extracted from 3D data. In addition to traditional 3D image analysis, we have developed algorithms to operate on curved surfaces, such as cell segmentation, lineage tracking and fluorescence signal quantification. The software’s modular design makes it easy to include existing libraries, or to implement new algorithms. Cell geometries extracted with MorphoGraphX can be exported and used as templates for simulation models, providing a powerful platform to investigate the interactions between shape, genes and growth.DOI: http://dx.doi.org/10.7554/eLife.05864.001Author keywordsResearch organism

Modeling the formation and aging of secondary organic aerosols during CalNex 2010

Four different literature parameterizations for the formation and evolution of urban secondary organic aerosol (SOA) frequently used in 3-D models are evaluated using a 0-D box model representing the Los Angeles metropolitan region during the California Research at the Nexus of Air Quality and Climate Change (CalNex) 2010 campaign. We constrain the model predictions with measurements from several platforms and compare predictions with particle- and gas-phase observations from the CalNex Pasadena ground site. That site provides a unique opportunity to study aerosol formation close to anthropogenic emission sources with limited recirculation. The model SOA that formed only from the oxidation of VOCs (V-SOA) is insufficient to explain the observed SOA concentrations, even when using SOA parameterizations with multi-generation oxidation that produce much higher yields than have been observed in chamber experiments, or when increasing yields to their upper limit estimates accounting for recently reported losses of vapors to chamber walls. The Community Multiscale Air Quality (WRF-CMAQ) model (version 5.0.1) provides excellent predictions of secondary inorganic particle species but underestimates the observed SOA mass by a factor of 25 when an older VOC-only parameterization is used, which is consistent with many previous model–measurement comparisons for pre-2007 anthropogenic SOA modules in urban areas. Including SOA from primary semi-volatile and intermediate-volatility organic compounds (P-S/IVOCs) following the parameterizations of Robinson et al. (2007), Grieshop et al. (2009), or Pye and Seinfeld (2010) improves model–measurement agreement for mass concentration. The results from the three parameterizations show large differences (e.g., a factor of 3 in SOA mass) and are not well constrained, underscoring the current uncertainties in this area. Our results strongly suggest that other precursors besides VOCs, such as P-S/IVOCs, are needed to explain the observed SOA concentrations in Pasadena. All the recent parameterizations overpredict urban SOA formation at long photochemical ages (3 days) compared to observations from multiple sites, which can lead to problems in regional and especially global modeling. However, reducing IVOC emissions by one-half in the model to better match recent IVOC measurements improves SOA predictions at these long photochemical ages. Among the explicitly modeled VOCs, the precursor compounds that contribute the greatest SOA mass are methylbenzenes. Measured polycyclic aromatic hydrocarbons (naphthalenes) contribute 0.7% of the modeled SOA mass. The amounts of SOA mass from diesel vehicles, gasoline vehicles, and cooking emissions are estimated to be 16–27, 35–61, and 19–35 %, respectively, depending on the parameterization used, which is consistent with the observed fossil fraction of urban SOA, 71(+-3) %. The relative contribution of each source is uncertain by almost a factor of 2 depending on the parameterization used. In-basin biogenic VOCs are predicted to contribute only a few percent to SOA. A regional SOA background of approximately 2.1 μgm-3 is also present due to the long-distance transport of highly aged OA, likely with a substantial contribution from regional biogenic SOA. The percentage of SOA from diesel vehicle emissions is the same, within the estimated uncertainty, as reported in previous work that analyzed the weekly cycles in OA concentrations (Bahreini et al., 2012; Hayes et al., 2013). However, the modeling work presented here suggests a strong anthropogenic source of modern carbon in SOA, due to cooking emissions, which was not accounted for in those previous studies and which is higher on weekends. Lastly, this work adapts a simple two-parameter model to predict SOA concentration and O/C from urban emissions. This model successfully predicts SOA concentration, and the optimal parameter combination is very similar to that found for Mexico City. This approach provides a computationally inexpensive method for predicting urban SOA in global and climate models. We estimate pollution SOA to account for 26 Tg yr-1 of SOA globally, or 17% of global SOA, one third of which is likely to be non-fossil.

Optimal Confidence Bands for Shape-Restricted Curves

Let Y be a stochastic process on [0,1] satisfying dY(t)=n 1/2 f(t)dt+dW(t) , where n≥1 is a given scale parameter (`sample size'), W is standard Brownian motion and f is an unknown function. Utilizing suitable multiscale tests, we construct confidence bands for f with guaranteed given coverage probability, assuming that f is isotonic or convex. These confidence bands are computationally feasible and shown to be asymptotically sharp optimal in an appropriate sense.

Confidence Bands for Isotonic Median Curves using Sign-Tests

Suppose that one observes independent random variables (X1, Y1), (X2, Y2), …, (Xn, Yn) in R2 with unknown distributions, except that Median(Yi | Xi = M(x) for some unknown isotonic function M. We describe an explicit algorithm for the computation of confidence bands for the median function M whose running time is of order O(n2). The bands rely on multiscale sign tests and are shown to have desirable asymptotic properties.

Structure–mechanics relationships in mineralized tendons

In this paper, we review the hierarchical structure and the resulting elastic properties of mineralized tendons as obtained by various multiscale experimental and computational methods spanning from nano- to macroscale. The mechanical properties of mineralized collagen fibres are important to understand the mechanics of hard tissues constituted by complex arrangements of these fibres, like in human lamellar bone. The uniaxial mineralized collagen fibre array naturally occurring in avian tendons is a well studied model tissue for investigating various stages of tissue mineralization and the corresponding elastic properties. Some avian tendons mineralize with maturation, which results in a graded structure containing two zones of distinct morphology, circumferential and interstitial. These zones exhibit different amounts of mineral, collagen, pores and a different mineral distribution between collagen fibrillar and extrafibrillar space that lead to distinct elastic properties. Mineralized tendon cells have two phenotypes: elongated tenocytes placed between fibres in the circumferential zone and cuboidal cells with lower aspect ratios in the interstitial zone. Interestingly some regions of avian tendons seem to be predestined to mineralization, which is exhibited as specific collagen cross-linking patterns as well as distribution of minor tendon constituents (like proteoglycans) and loss of collagen crimp. Results of investigations in naturally mineralizing avian tendons may be useful in understanding the pathological mineralization occurring in some human tendons.

Matching Social and Ecological Systems in Complex Ocean Fisheries

This paper considers ocean fisheries as complex adaptive systems and addresses the question of how human institutions might be best matched to their structure and function. Ocean ecosystems operate at multiple scales, but the management of fisheries tends to be aimed at a single species considered at a single broad scale. The paper argues that this mismatch of ecological and management scale makes it difficult to address the fine-scale aspects of ocean ecosystems, and leads to fishing rights and strategies that tend to erode the underlying structure of populations and the system itself. A successful transition to ecosystem-based management will require institutions better able to economize on the acquisition of feedback about the impact of human activities. This is likely to be achieved by multiscale institutions whose organization mirrors the spatial organization of the ecosystem and whose communications occur through a polycentric network. Better feedback will allow the exploration of fine-scale science and the employment of fine-scale fishing restraints, better adapted to the behavior of fish and habitat. The scale and scope of individual fishing rights also needs to be congruent with the spatial structure of the ecosystem. Place-based rights can be expected to create a longer private planning horizon as well as stronger incentives for the private and public acquisition of system relevant knowledge.

Contribution of activity and skin temperature to the risk of developing pressure ulcers in nursing facility residents

Purpose. The central concepts in pressure ulcer risk are exposure to external pressure caused by inactivity and tissue tolerance to pressure, a factor closely related to blood flow. Inactivity measures are effective in predicting pressure ulcer risk. The purpose of the study is to evaluate whether a physiological measure of skin blood flow improves pressure ulcer risk prediction. Skin temperature regularity and self-similarity, as proxy measures of blood flow, and not previously described, may be undefined pressure ulcer risk factors. The specific aims were to determine whether a sample of nursing facility residents at high risk of pressure ulcers classified using the Braden Scale for Pressure Sore Risk© differ from a sample of low risk residents according to (1) exposure to external pressure as measured by resident activity, (2) tissue tolerance to external pressure as measured by skin temperature, and (3) skin temperature fluctuations and recovery in response to a commonly occurring stressor, bathing and additionally whether (4) scores on the Braden Scale mobility subscale score are related to entropy and the spectral exponent. ^ Methods. A two group observational time series design was used to describe activity and skin temperature regularity and self-similarity, calculating entropy and the spectral exponent using detrended fluctuation analysis respectively. Twenty nursing facility residents wore activity and skin temperature monitors for one week. One bathing episode was observed as a commonly occurring stressor for skin temperature.^ Results. Skin temperature multiscale entropy (MSE), F(1, 17) = 5.55, p = .031, the skin temperature spectral exponent, F(1, 17) = 6.19, p = .023, and the activity mean MSE, F(1, 18) = 4.52, p = .048 differentiated the risk groups. The change in skin temperature entropy during bathing was significant, t(16) = 2.55, p = .021, (95% CI, .04-.40). Multiscale entropy for skin temperature was lowest in those who developed pressure ulcers, F(1, 18) = 35.14, p < .001.^ Conclusions. This study supports the tissue tolerance component of the Braden and Bergstrom conceptual framework and shows differences in skin temperature multiscale entropy between pressure ulcer risk categories, pressure ulcer outcome, and during a commonly occurring stressor. ^

Spectroscopic and functional investigations of the AMPA subtype of ionotropic glutamate receptors

Ion channels play a crucial role in the functioning of different systems of the body because of their ability to bridge the cell membrane and allow ions to pass in and out of the cell. Ionotropic glutamate receptors are one class of these important proteins and have been shown to be critical in propagating synaptic transmission in the central nervous system and in other diverse functions throughout the body. Because of their wide-ranging effects, this family of receptors is an important target for structure-function investigations to understand their mechanism of action. ^ α-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors are one subtype of glutamate receptors and have been shown to be the primary receptors involved in rapid excitatory signaling in the central nervous system. Agonist binding to the extracellular ligand binding domain of these receptors causes various conformational changes that culminate in formation of the ion channel. Previous structural investigations have provided important information about their mechanism of action, including uncovering a relationship between the degree of cleft closure in the binding domain and activation of the receptor. However, what question remains unanswered is how specific interactions between the agonist and the protein interplay with cleft closure to mediate receptor activation. ^ To investigate this question, I applied a multiscale approach to investigate the effects of agonist binding on various levels. Vibrational spectroscopy was utilized to investigate molecular-level interactions in the binding pocket, and fluorescence resonance energy transfer (FRET) was employed to measure cleft closure in the isolated ligand binding domain. The results of these studies in the isolated binding domain were then correlated to activation of the full receptor. These investigations showed a relationship between the strength of the interaction at the α-amine group of the agonist and extent of receptor activation, where a stronger interaction correlated to a larger activation, which was upheld even when the extent of cleft closure did not correlate to activation. These results show that this interaction at the α-amine group is critical in mediating the allosteric mechanism of activation and provide a bit more insight into how agonist binding is coupled to channel gating in AMPA receptors. ^