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

The role of ongoing dendritic oscillations in single-neuron dynamics

The dendritic tree contributes significantly to the elementary computations a neuron performs while converting its synaptic inputs into action potential output. Traditionally, these computations have been characterized as temporally local, near-instantaneous mappings from the current input of the cell to its current output, brought about by somatic summation of dendritic contributions that are generated in spatially localized functional compartments. However, recent evidence about the presence of oscillations in dendrites suggests a qualitatively different mode of operation: the instantaneous phase of such oscillations can depend on a long history of inputs, and under appropriate conditions, even dendritic oscillators that are remote may interact through synchronization. Here, we develop a mathematical framework to analyze the interactions of local dendritic oscillations, and the way these interactions influence single cell computations. Combining weakly coupled oscillator methods with cable theoretic arguments, we derive phase-locking states for multiple oscillating dendritic compartments. We characterize how the phase-locking properties depend on key parameters of the oscillating dendrite: the electrotonic properties of the (active) dendritic segment, and the intrinsic properties of the dendritic oscillators. As a direct consequence, we show how input to the dendrites can modulate phase-locking behavior and hence global dendritic coherence. In turn, dendritic coherence is able to gate the integration and propagation of synaptic signals to the soma, ultimately leading to an effective control of somatic spike generation. Our results suggest that dendritic oscillations enable the dendritic tree to operate on more global temporal and spatial scales than previously thought

Cellular expression, trafficking, and function of two isoforms of human ULBP5/RAET1G

Background: The activating immunoreceptor NKG2D is expressed on Natural Killer (NK) cells and subsets of T cells. NKG2D contributes to anti-tumour and anti-viral immune responses in vitro and in vivo. The ligands for NKG2D in humans are diverse proteins of the MIC and ULBP/RAET families that are upregulated on the surface of virally infected cells and tumours. Two splicing variants of ULBP5/RAET1G have been cloned previously, but not extensively characterised. Methodology/Principal Findings: We pursue a number of approaches to characterise the expression, trafficking, and function of the two isoforms of ULBP5/RAET1G. We show that both transcripts are frequently expressed in cell lines derived from epithelial cancers, and in primary breast cancers. The full-length transcript, RAET1G1, is predicted to encode a molecule with transmembrane and cytoplasmic domains that are unique amongst NKG2D ligands. Using specific anti-RAET1G1 antiserum to stain tissue microarrays we show that RAET1G1 expression is highly restricted in normal tissues. RAET1G1 was expressed at a low level in normal gastrointestinal epithelial cells in a similar pattern to MICA. Both RAET1G1 and MICA showed increased expression in the gut of patients with celiac disease. In contrast to healthy tissues the RAET1G1 antiserum stained a wide variety or different primary tumour sections. Both endogenously expressed and transfected RAET1G1 was mainly found inside the cell, with a minority of the protein reaching the cell surface. Conversely the truncated splicing variant of RAET1G2 was shown to encode a soluble molecule that could be secreted from cells. Secreted RAET1G2 was shown to downregulate NKG2D receptor expression on NK cells and hence may represent a novel tumour immune evasion strategy. Conclusions/Significance: We demonstrate that the expression patterns of ULBP5RAET1G are very similar to the well-characterised NKG2D ligand, MICA. However the two isoforms of ULBP5/RAET1G have very different cellular localisations that are likely to reflect unique functionality

"It's making contacts" : notions of social capital and implications for widening access to medical education

Acknowledgements Our thanks to the Medical Schools Council (MSC) of the UK for funding Study A; REACH Scotland for funding Study B; and Queen Mary University of London, and to the medical school applicants and students who gave their time to be interviewed. Our thanks also to Dr Sean Zhou and Dr Sally Curtis, and Manjul Medhi, for their help with data collection for studies A and B respectively. Our thanks also to Dr Lara Varpio, Uniformed Services University of the USA, for her advice and guidance on collating data sets and her comments on the draft manuscript.Peer reviewedPublisher PD

Intervening to improve health indicators among Australian farm families

The Sustainable Farm Families project (http://www.sustainablefarmfamilies.org.au/) was a 3-year demonstration and education project designed to influence farmer behavior with respect to family health and well-being among cropping and grazing farmers in Victoria, New South Wales, and South Australia, Australia. The project was conducted by the Western District Health Service, Hamilton, Australia, in partnership with farmers; Farm Management 500 (peer discussion group); the Victorian Farmers Federation; Royal Melbourne Institute of Technology; and Land Connect. During the 3 years of the project, 128 farmers&mdash;men (70) and women (58)&mdash;were enrolled. The project utilized a combination of small group workshops, individualized health action plans, and health education opportunities to encourage farm safety and health behavior changes and to elicit sustained improvements in the following health indicators: body mass index (BMI), total cholesterol, fasting blood glucose, and blood pressure. Mean changes in these health indicators were analyzed using repeated measures analysis of variance (ANOVA) and McNemar\u27s test compared the proportion of individuals with elevated indicators. Among participants with elevated values at baseline, the following average reductions were observed: BMI 0.44 kg/m2 (p = .0034), total cholesterol 48.7 mg/dl (p &lt; .0001), blood glucose 10.1 mg/dl (p = .0016), systolic blood pressure 12.5 mm Hg (p &lt; .0001), and diastolic blood pressure 5.0 mm Hg (p = .0007). The proportion of participants with elevated total cholesterol at baseline decreased after 24 months (p &lt; .001). Such findings suggest that proactive intervention by farmer associations, rural health services, and government agencies may be an effective vehicle for promoting voluntary farm safety and health behavior change while empowering farm families to achieve measurable reductions in important health risk factors. <br /

A class of extremising sphere-valued maps with inherent maximal tori symmetries in SO(n)

In this paper we consider an energy functional depending on the norm of the gradient and seek to extremise it over an admissible class of Sobolev maps defined on an annulus and taking values on the unit sphere whilst satisfying suitable boundary conditions. We establish the existence of an infinite family of solutions with certain symmetries to the associated nonlinear Euler-Lagrange system in even dimensions and discuss the stability of such extremisers by way of examining the positivity of the second variation of the energy at these solutions

Using ordinal logistic regression to evaluate the performance of laser-Doppler predictions of burn-healing time

Background Laser-Doppler imaging (LDI) of cutaneous blood flow is beginning to be used by burn surgeons to predict the healing time of burn wounds; predicted healing time is used to determine wound treatment as either dressings or surgery. In this paper, we do a statistical analysis of the performance of the technique. Methods We used data from a study carried out by five burn centers: LDI was done once between days 2 to 5 post burn, and healing was assessed at both 14 days and 21 days post burn. Random-effects ordinal logistic regression and other models such as the continuation ratio model were used to model healing-time as a function of the LDI data, and of demographic and wound history variables. Statistical methods were also used to study the false-color palette, which enables the laser-Doppler imager to be used by clinicians as a decision-support tool. Results Overall performance is that diagnoses are over 90% correct. Related questions addressed were what was the best blood flow summary statistic and whether, given the blood flow measurements, demographic and observational variables had any additional predictive power (age, sex, race, % total body surface area burned (%TBSA), site and cause of burn, day of LDI scan, burn center). It was found that mean laser-Doppler flux over a wound area was the best statistic, and that, given the same mean flux, women recover slightly more slowly than men. Further, the likely degradation in predictive performance on moving to a patient group with larger %TBSA than those in the data sample was studied, and shown to be small. Conclusion Modeling healing time is a complex statistical problem, with random effects due to multiple burn areas per individual, and censoring caused by patients missing hospital visits and undergoing surgery. This analysis applies state-of-the art statistical methods such as the bootstrap and permutation tests to a medical problem of topical interest. New medical findings are that age and %TBSA are not important predictors of healing time when the LDI results are known, whereas gender does influence recovery time, even when blood flow is controlled for. The conclusion regarding the palette is that an optimum three-color palette can be chosen 'automatically', but the optimum choice of a 5-color palette cannot be made solely by optimizing the percentage of correct diagnoses

Search for Second-Generation Scalar Leptoquarks in ppˉ\bm{p \bar{p}} Collisions at s\sqrt{s}=1.96 TeV

Results on a search for pair production of second generation scalar leptoquark in $p \bar{p}$ collisions at $\sqrt{s}$=1.96 TeV are reported. The data analyzed were collected by the CDF detector during the 2002-2003 Tevatron Run II and correspond to an integrated luminosity of 198 pb$^{-1}$. Leptoquarks (LQ) are sought through their decay into (charged) leptons and quarks, with final state signatures represented by two muons and jets and one muon, large transverse missing energy and jets. We observe no evidence for $LQ$ production and derive 95% C.L. upper limits on the $LQ$ production cross sections as well as lower limits on their mass as a function of $\beta$, where $\beta$ is the branching fraction for $LQ \to \mu q$.Comment: 9 pages (3 author list) 5 figure