Sequential Monte Carlo methods for Bayesian elliptic inverse problems

In this article, we consider a Bayesian inverse problem associated to elliptic partial differential equations in two and three dimensions. This class of inverse problems is important in applications such as hydrology, but the complexity of the link function between unknown field and measurements can make it difficult to draw inference from the associated posterior. We prove that for this inverse problem a basic sequential Monte Carlo (SMC) method has a Monte Carlo rate of convergence with constants which are independent of the dimension of the discretization of the problem; indeed convergence of the SMC method is established in a function space setting. We also develop an enhancement of the SMC methods for inverse problems which were introduced in Kantas et al. (SIAM/ASA J Uncertain Quantif 2:464–489, 2014); the enhancement is designed to deal with the additional complexity of this elliptic inverse problem. The efficacy of the methodology and its desirable theoretical properties, are demonstrated for numerical examples in both two and three dimensions

Endovenous laser ablation therapy in children: applications and outcomes

BACKGROUND: Endovenous laser ablation is well recognized as the first-line treatment for superficial venous reflux with varicose veins in adults. It is not widely reported and is not an established practice in pediatric patients. OBJECTIVE: To illustrate a variety of pediatric venous conditions in which endovenous laser ablation can be utilized and to demonstrate its feasibility and safety in children. MATERIALS AND METHODS: We conducted a retrospective review of endovenous laser ablation procedures performed between January 2007 and July 2014 at two large pediatric institutions. RESULTS: We included 35 patients (17 males) who underwent endovenous laser ablation to 43 veins. Median age at first treatment was 14 years (range: 3-18 years). Median weight was 56 kg (range: 19-97 kg). Underlying diagnoses were common venous malformation (15), Klippel-Trenaunay syndrome (8), superficial venous reflux with varicose veins (5), verrucous hemangioma-related phlebectasia (4), venous varix (2) and arteriovenous fistula (1). The most common aim of treatment was to facilitate sclerotherapy. Thirty-four patients had treatment in the lower limbs and one patient in an upper limb. Ten of the veins treated with endovenous laser ablation had an additional procedure performed to close the vein. Complications attributable to endovenous laser ablation occurred in two patients (6%). One patient experienced post-procedural pain and one patient developed a temporary sensory nerve injury. Median clinical follow-up was 13 months (range: 28 days-5.7 years). The aim of the treatment was achieved in 29 of the 35 (83%) patients. CONCLUSION: Endovenous laser ablation is technically feasible and safe in children. It can be used in the management of a range of pediatric venous diseases with good outcomes

What is the correct cost functional for variational data assimilation?

Variational approaches to data assimilation, and weakly constrained four dimensional variation (WC-4DVar) in particular, are important in the geosciences but also in other communities (often under different names). The cost functions and the resulting optimal trajectories may have a probabilistic interpretation, for instance by linking data assimilation with maximum aposteriori (MAP) estimation. This is possible in particular if the unknown trajectory is modelled as the solution of a stochastic differential equation (SDE), as is increasingly the case in weather forecasting and climate modelling. In this situation, the MAP estimator (or “most probable path” of the SDE) is obtained by minimising the Onsager–Machlup functional. Although this fact is well known, there seems to be some confusion in the literature, with the energy (or “least squares”) functional sometimes been claimed to yield the most probable path. The first aim of this paper is to address this confusion and show that the energy functional does not, in general, provide the most probable path. The second aim is to discuss the implications in practice. Although the mentioned results pertain to stochastic models in continuous time, they do have consequences in practice where SDE’s are approximated by discrete time schemes. It turns out that using an approximation to the SDE and calculating its most probable path does not necessarily yield a good approximation to the most probable path of the SDE proper. This suggest that even in discrete time, a version of the Onsager–Machlup functional should be used, rather than the energy functional, at least if the solution is to be interpreted as a MAP estimator

Mapping species distributions: A comparison of skilled naturalist and lay citizen science recording

To assess the ability of traditional biological recording schemes and lay citizen science approaches to gather data on species distributions and changes therein, we examined bumblebee records from the UK’s national repository (National Biodiversity Network) and from BeeWatch. The two recording approaches revealed similar relative abundances of bumblebee species but different geographical distributions. For the widespread common carder (Bombus pascuorum), traditional recording scheme data were patchy, both spatially and temporally, reflecting active record centre rather than species distribution. Lay citizen science records displayed more extensive geographic coverage, reflecting human population density, thus offering better opportunities to account for recording effort. For the rapidly spreading tree bumblebee (Bombus hypnorum), both recording approaches revealed similar distributions due to a dedicated mapping project which overcame the patchy nature of naturalist records. We recommend, where possible, complementing skilled naturalist recording with lay citizen science programmes to obtain a nation-wide capability, and stress the need for timely uploading of data to the national repository

A demonstration of 'broken' visual space

It has long been assumed that there is a distorted mapping between real and ‘perceived’ space, based on demonstrations of systematic errors in judgements of slant, curvature, direction and separation. Here, we have applied a direct test to the notion of a coherent visual space. In an immersive virtual environment, participants judged the relative distance of two squares displayed in separate intervals. On some trials, the virtual scene expanded by a factor of four between intervals although, in line with recent results, participants did not report any noticeable change in the scene. We found that there was no consistent depth ordering of objects that can explain the distance matches participants made in this environment (e.g. A > B > D yet also A < C < D) and hence no single one-to-one mapping between participants’ perceived space and any real 3D environment. Instead, factors that affect pairwise comparisons of distances dictate participants’ performance. These data contradict, more directly than previous experiments, the idea that the visual system builds and uses a coherent 3D internal representation of a scene

Parameter identification problems in the modelling of cell motility

We present a novel parameter identification algorithm for the estimation of parameters in models of cell motility using imaging data of migrating cells. Two alternative formulations of the objective functional that measures the difference between the computed and observed data are proposed and the parameter identification problem is formulated as a minimisation problem of nonlinear least squares type. A Levenberg–Marquardt based optimisation method is applied to the solution of the minimisation problem and the details of the implementation are discussed. A number of numerical experiments are presented which illustrate the robustness of the algorithm to parameter identification in the presence of large deformations and noisy data and parameter identification in three dimensional models of cell motility. An application to experimental data is also presented in which we seek to identify parameters in a model for the monopolar growth of fission yeast cells using experimental imaging data. Our numerical tests allow us to compare the method with the two different formulations of the objective functional and we conclude that the results with both objective functionals seem to agree

Characterizations of how species mediate ecosystem properties require more comprehensive functional effect descriptors

The importance of individual species in mediating ecosystem process and functioning is generally accepted, but categorical descriptors that summarize species-specific contributions to ecosystems tend to reference a limited number of biological traits and underestimate the importance of how organisms interact with their environment. Here, we show how three functionally contrasting sediment-dwelling marine invertebrates affect fluid and particle transport - important processes in mediating nutrient cycling - and use high-resolution reconstructions of burrow geometry to determine the extent and nature of biogenic modification. We find that individual functional effect descriptors fall short of being able to adequately characterize how species mediate the stocks and flows of important ecosystem properties and that, in contrary to common practice and understanding, they are not substitutable with one another because they emphasize different aspects of species activity and behavior. When information derived from these metrics is combined with knowledge of how species behave and modify their environment, however, detailed mechanistic information emerges that increases the likelihood that a species functional standing will be appropriately summarized. Our study provides evidence that more comprehensive functional effect descriptors are required if they are to be of value to those tasked with projecting how altered biodiversity will influence future ecosystems

SCAMP:standardised, concentrated, additional macronutrients, parenteral nutrition in very preterm infants: a phase IV randomised, controlled exploratory study of macronutrient intake, growth and other aspects of neonatal care

<p>Abstract</p> <p>Background</p> <p>Infants born <29 weeks gestation are at high risk of neurocognitive disability. Early postnatal growth failure, particularly head growth, is an important and potentially reversible risk factor for impaired neurodevelopmental outcome. Inadequate nutrition is a major factor in this postnatal growth failure, optimal protein and calorie (macronutrient) intakes are rarely achieved, especially in the first week. Infants <29 weeks are dependent on parenteral nutrition for the bulk of their nutrient needs for the first 2-3 weeks of life to allow gut adaptation to milk digestion. The prescription, formulation and administration of neonatal parenteral nutrition is critical to achieving optimal protein and calorie intake but has received little scientific evaluation. Current neonatal parenteral nutrition regimens often rely on individualised prescription to manage the labile, unpredictable biochemical and metabolic control characteristic of the early neonatal period. Individualised prescription frequently fails to translate into optimal macronutrient delivery. We have previously shown that a standardised, concentrated neonatal parenteral nutrition regimen can optimise macronutrient intake.</p> <p>Methods</p> <p>We propose a single centre, randomised controlled exploratory trial of two standardised, concentrated neonatal parenteral nutrition regimens comparing a standard macronutrient content (maximum protein 2.8 g/kg/day; lipid 2.8 g/kg/day, dextrose 10%) with a higher macronutrient content (maximum protein 3.8 g/kg/day; lipid 3.8 g/kg/day, dextrose 12%) over the first 28 days of life. 150 infants 24-28 completed weeks gestation and birthweight <1200 g will be recruited. The primary outcome will be head growth velocity in the first 28 days of life. Secondary outcomes will include a) auxological data between birth and 36 weeks corrected gestational age b) actual macronutrient intake in first 28 days c) biomarkers of biochemical and metabolic tolerance d) infection biomarkers and other intravascular line complications e) incidence of major complications of prematurity including mortality f) neurodevelopmental outcome at 2 years corrected gestational age</p> <p>Trial registration</p> <p>Current controlled trials: <a href="http://www.controlled-trials.com/ISRCTN76597892">ISRCTN76597892</a>; EudraCT Number: 2008-008899-14</p

Concurrent use of prescription drugs and herbal medicinal products in older adults: A systematic review

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 International License (http://creativecommons.org/licenses/by-nc/4.0/), which permits any noncommercial use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.The use of herbal medicinal products (HMPs) is common among older adults. However, little is known about concurrent use with prescription drugs as well as the potential interactions associated with such combinations. Objective Identify and evaluate the literature on concurrent prescription and HMPs use among older adults to assess prevalence, patterns, potential interactions and factors associated with this use. Methods Systematic searches in MEDLINE, PsycINFO, EMBASE, CINAHL, AMED, Web of Science and Cochrane from inception to May 2017 for studies reporting concurrent use of prescription medicines with HMPs in adults (≥65 years). Quality was assessed using the Joanna Briggs Institute checklists. The Evidence for Policy and Practice Information and Co-ordinating Centre (EPPI-Centre) three stage approach to mixed method research was used to synthesise data. Results Twenty-two studies were included. A definition of HMPs or what was considered HMP was frequently missing. Prevalence of concurrent use by older adults varied widely between 5.3% and 88.3%. Prescription medicines most combined with HMPs were antihypertensive drugs, beta blockers, diuretics, antihyperlipidemic agents, anticoagulants, analgesics, antihistamines, antidiabetics, antidepressants and statins. The HMPs most frequently used were: ginkgo, garlic, ginseng, St John’s wort, Echinacea, saw palmetto, evening primrose oil and ginger. Potential risks of bleeding due to use of ginkgo, garlic or ginseng with aspirin or warfarin was the most reported herb-drug interaction. Some data suggests being female, a lower household income and less than high school education were associated with concurrent use. Conclusion Prevalence of concurrent prescription drugs and HMPs use among older adults is substantial and potential interactions have been reported. Knowledge of the extent and manner in which older adults combine prescription drugs will aid healthcare professionals can appropriately identify and manage patients at risk.Peer reviewedFinal Published versio

Neurogenesis Drives Stimulus Decorrelation in a Model of the Olfactory Bulb

The reshaping and decorrelation of similar activity patterns by neuronal networks can enhance their discriminability, storage, and retrieval. How can such networks learn to decorrelate new complex patterns, as they arise in the olfactory system? Using a computational network model for the dominant neural populations of the olfactory bulb we show that fundamental aspects of the adult neurogenesis observed in the olfactory bulb -- the persistent addition of new inhibitory granule cells to the network, their activity-dependent survival, and the reciprocal character of their synapses with the principal mitral cells -- are sufficient to restructure the network and to alter its encoding of odor stimuli adaptively so as to reduce the correlations between the bulbar representations of similar stimuli. The decorrelation is quite robust with respect to various types of perturbations of the reciprocity. The model parsimoniously captures the experimentally observed role of neurogenesis in perceptual learning and the enhanced response of young granule cells to novel stimuli. Moreover, it makes specific predictions for the type of odor enrichment that should be effective in enhancing the ability of animals to discriminate similar odor mixtures