Sediment concentration prediction and statistical evaluation for annual load estimation

We consider the development of statistical models for prediction of constituent concentration of riverine pollutants, which is a key step in load estimation from frequent flow rate data and less frequently collected concentration data. We consider how to capture the impacts of past flow patterns via the average discounted flow (ADF) which discounts the past flux based on the time lapsed - more recent fluxes are given more weight. However, the effectiveness of ADF depends critically on the choice of the discount factor which reflects the unknown environmental cumulating process of the concentration compounds. We propose to choose the discount factor by maximizing the adjusted R-2 values or the Nash-Sutcliffe model efficiency coefficient. The R2 values are also adjusted to take account of the number of parameters in the model fit. The resulting optimal discount factor can be interpreted as a measure of constituent exhaustion rate during flood events. To evaluate the performance of the proposed regression estimators, we examine two different sampling scenarios by resampling fortnightly and opportunistically from two real daily datasets, which come from two United States Geological Survey (USGS) gaging stations located in Des Plaines River and Illinois River basin. The generalized rating-curve approach produces biased estimates of the total sediment loads by -30% to 83%, whereas the new approaches produce relatively much lower biases, ranging from -24% to 35%. This substantial improvement in the estimates of the total load is due to the fact that predictability of concentration is greatly improved by the additional predictors.

Delineating environmental envelopes to improve mapping of species distributions, via a hurdle model with CART &/or MaxEnT

Species distribution modelling (SDM) typically analyses species’ presence together with some form of absence information. Ideally absences comprise observations or are inferred from comprehensive sampling. When such information is not available, then pseudo-absences are often generated from the background locations within the study region of interest containing the presences, or else absence is implied through the comparison of presences to the whole study region, e.g. as is the case in Maximum Entropy (MaxEnt) or Poisson point process modelling. However, the choice of which absence information to include can be both challenging and highly influential on SDM predictions (e.g. Oksanen and Minchin, 2002). In practice, the use of pseudo- or implied absences often leads to an imbalance where absences far outnumber presences. This leaves analysis highly susceptible to ‘naughty-noughts’: absences that occur beyond the envelope of the species, which can exert strong influence on the model and its predictions (Austin and Meyers, 1996). Also known as ‘excess zeros’, naughty noughts can be estimated via an overall proportion in simple hurdle or mixture models (Martin et al., 2005). However, absences, especially those that occur beyond the species envelope, can often be more diverse than presences. Here we consider an extension to excess zero models. The two-staged approach first exploits the compartmentalisation provided by classification trees (CTs) (as in O’Leary, 2008) to identify multiple sources of naughty noughts and simultaneously delineate several species envelopes. Then SDMs can be fit separately within each envelope, and for this stage, we examine both CTs (as in Falk et al., 2014) and the popular MaxEnt (Elith et al., 2006). We introduce a wider range of model performance measures to improve treatment of naughty noughts in SDM. We retain an overall measure of model performance, the area under the curve (AUC) of the Receiver-Operating Curve (ROC), but focus on its constituent measures of false negative rate (FNR) and false positive rate (FPR), and how these relate to the threshold in the predicted probability of presence that delimits predicted presence from absence. We also propose error rates more relevant to users of predictions: false omission rate (FOR), the chance that a predicted absence corresponds to (and hence wastes) an observed presence, and the false discovery rate (FDR), reflecting those predicted (or potential) presences that correspond to absence. A high FDR may be desirable since it could help target future search efforts, whereas zero or low FOR is desirable since it indicates none of the (often valuable) presences have been ignored in the SDM. For illustration, we chose Bradypus variegatus, a species that has previously been published as an exemplar species for MaxEnt, proposed by Phillips et al. (2006). We used CTs to increasingly refine the species envelope, starting with the whole study region (E0), eliminating more and more potential naughty noughts (E1–E3). When combined with an SDM fit within the species envelope, the best CT SDM had similar AUC and FPR to the best MaxEnt SDM, but otherwise performed better. The FNR and FOR were greatly reduced, suggesting that CTs handle absences better. Interestingly, MaxEnt predictions showed low discriminatory performance, with the most common predicted probability of presence being in the same range (0.00-0.20) for both true absences and presences. In summary, this example shows that SDMs can be improved by introducing an initial hurdle to identify naughty noughts and partition the envelope before applying SDMs. This improvement was barely detectable via AUC and FPR yet visible in FOR, FNR, and the comparison of predicted probability of presence distribution for pres/absence.

Genetic studies of body mass index yield new insights for obesity biology

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in up to 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 × 10−8), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for ~2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous system in obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

A novel approach to assessing road-curve crash severity

Curves are a common feature of road infrastructure; however crashes on road curves are associated with increased risk of injury and fatality to vehicle occupants. Countermeasures require the identification of contributing factors. However, current approaches to identifying contributors use traditional statistical methods and have not used self-reported narrative claim to identify factors related to the driver, vehicle and environment in a systemic way. Text mining of 3434 road-curve crash claim records filed between 1 January 2003 and 31 December 2005 at a major insurer in Queensland, Australia, was undertaken to identify risk levels and contributing factors. Rough set analysis was used on insurance claim narratives to identify significant contributing factors to crashes and their associated severity. New contributing factors unique to curve crashes were identified (e.g., tree, phone, over-steer) in addition to those previously identified via traditional statistical analysis of Police and licensing authority records. Text mining is a novel methodology to improve knowledge related to risk and contributing factors to road-curve crash severity. Future road-curve crash countermeasures should more fully consider the interrelationships between environment, the road, the driver and the vehicle, and education campaigns in particular could highlight the increased risk of crash on road-curves.

Gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

Background Adolescent Idiopathic Scoliosis is the most common type of spinal deformity, and whilst the risk of progression appears to be biomechanically mediated (larger deformities are more likely to progress), the detailed biomechanical mechanisms driving progression are not well understood. Gravitational forces in the upright position are the primary sustained loads experienced by the spine. In scoliosis they are asymmetrical, generating moments about the spinal joints which may promote asymmetrical growth and deformity progression. Using 3D imaging modalities to estimate segmental torso masses allows the gravitational loading on the scoliotic spine to be determined. The resulting distribution of joint moments aids understanding of the mechanics of scoliosis progression. Methods Existing low-dose CT scans were used to estimate torso segment masses and joint moments for 20 female scoliosis patients. Intervertebral joint moments at each vertebral level were found by summing the moments of each of the torso segment masses above the required joint. Results The patients’ mean age was 15.3 years (SD 2.3; range 11.9 – 22.3 years); mean thoracic major Cobb angle 52° (SD 5.9°; range 42°-63°) and mean weight 57.5 kg (SD 11.5 kg; range 41 – 84.7 kg). Joint moments of up to 7 Nm were estimated at the apical level. No significant correlation was found between the patients’ major Cobb angles and apical joint moments. Conclusions Patients with larger Cobb angles do not necessarily have higher joint moments, and curve shape is an important determinant of joint moment distribution. These findings may help to explain the variations in progression between individual patients. This study suggests that substantial corrective forces are required of either internal instrumentation or orthoses to effectively counter the gravity-induced moments acting to deform the spinal joints of idiopathic scoliosis patients.

Dark-energy constraints and correlations with systematics from CFHTLS weak lensing, SNLS supernovae Ia and WMAP5

Aims We combine measurements of weak gravitational lensing from the CFHTLS-Wide survey, supernovae Ia from CFHT SNLS and CMB anisotropies from WMAP5 to obtain joint constraints on cosmological parameters, in particular, the dark-energy equation-of-state parameter w. We assess the influence of systematics in the data on the results and look for possible correlations with cosmological parameters. Methods We implemented an MCMC algorithm to sample the parameter space of a flat CDM model with a dark-energy component of constant w. Systematics in the data are parametrised and included in the analysis. We determine the influence of photometric calibration of SNIa data on cosmological results by calculating the response of the distance modulus to photometric zero-point variations. The weak lensing data set is tested for anomalous field-to-field variations and a systematic shape measurement bias for high-redshift galaxies. Results Ignoring photometric uncertainties for SNLS biases cosmological parameters by at most 20% of the statistical errors, using supernovae alone; the parameter uncertainties are underestimated by 10%. The weak-lensing field-to-field variance between 1 deg2-MegaCam pointings is 5-15% higher than predicted from N-body simulations. We find no bias in the lensing signal at high redshift, within the framework of a simple model, and marginalising over cosmological parameters. Assuming a systematic underestimation of the lensing signal, the normalisation increases by up to 8%. Combining all three probes we obtain -0.10 < 1 + w < 0.06 at 68% confidence ( -0.18 < 1 + w < 0.12 at 95%), including systematic errors. Our results are therefore consistent with the cosmological constant . Systematics in the data increase the error bars by up to 35%; the best-fit values change by less than 0.15.

Oral medication delivery in impaired swallowing: Thickening liquid medications for safe swallowing alters dissolution characteristics

Acetaminophen (paracetamol) is available in a wide range of oral formulations designed to meet the needs of the population across the age-spectrum, but for people with impaired swallowing, i.e. dysphagia, both solid and liquid medications can be difficult to swallow without modification. The effect of a commercial polysaccharide thickener, designed to be added to fluids to promote safe swallowing by dysphagic patients, on rheology and acetaminophen dissolution was tested using crushed immediate-release tablets in water, effervescent tablets in water, elixir and suspension. The inclusion of the thickener, comprised of xanthan gum and maltodextrin, had a considerable impact on dissolution; acetaminophen release from modified medications reached 12-50% in 30 minutes, which did not reflect the pharmacopeia specification for immediate release preparations. Flow curves reflect the high zero-shear viscosity and the apparent yield stress of the thickened products. The weak gel nature, in combination with high G’ values compared to G” (viscoelasticity) and high apparent yield stress, impact drug release. The restriction on drug release from these formulations is not influenced by the theoretical state of the drug (dissolved or dispersed), and the approach typically used in clinical practice (mixing crushed tablets into pre-prepared thickened fluid) cannot be improved by altering the order of incorporation or mixing method.

A curve shortening flow rule for closed embedded plane curves with a prescribed rate of change in enclosed area

Motivated by a problem from fluid mechanics, we consider a generalization of the standard curve shortening flow problem for a closed embedded plane curve such that the area enclosed by the curve is forced to decrease at a prescribed rate. Using formal asymptotic and numerical techniques, we derive possible extinction shapes as the curve contracts to a point, dependent on the rate of decreasing area; we find there is a wider class of extinction shapes than for standard curve shortening, for which initially simple closed curves are always asymptotically circular. We also provide numerical evidence that self-intersection is possible for non-convex initial conditions, distinguishing between pinch-off and coalescence of the curve interior.

Nitrification (DMPP) and urease (NBPT) inhibitors had no effect on pasture yield, nitrous oxide emissions nor nitrate leaching under irrigation in a hot-dry climate

Modern dairy farming in Australia relies on substantial inputs of fertiliser nitrogen (N) to underpin economic production. However, N lost from dairy systems represents an opportunity cost and can pose a number of environmental risks. Nitrogen cycle inhibitors can be co-applied with N fertilisers to slow the conversion of urea to NH4+ to reduce losses via volatilisation, and slow the conversion of NH4+ to NO3- to minimize leaching of NO3- and gaseous losses via nitrification and denitrification. In a field campaign in a high input ryegrass-kikuyu pasture system we compared the soil N pools, losses and pasture production between a) urea coated with the nitrification inhibitor (3,4-dimethyl pyrazole phosphate - DMPP) b) urea coated with the urease inhibitor (N-(n-butyl) thiophosphoric triamide - NBPT) and c) standard urea. There was no treatment effect (P>0.05) on soil mineral N, pasture yield, N2O flux nor leaching of NO3- cf. standard urea. We hypothesise that at our site, because gaseous losses were highly episodic (rainfall was erratic and displayed no seasonal rainfall nor soil wetting pattern) that there was a lack of coincidence of N application and conditions conducive to gaseous losses, thus the effectiveness of the inhibitor products was minimal and did not result in an increase in pasture yield. There remains a paucity of knowledge on N cycle inhibitors in relation to their effective use in field system to increase N use efficiency. Further research is required to define under what field conditions inhibitor products are effective in order to be able to provide accurate advice to managers of nitrogen in production systems.

Appendiceal sonography in children: A retrospective analysis as a platform for potential implementation of diagnostic categories

- Background Sonography is an important diagnostic tool in children with suspected appendicitis. Reported accuracy and appendiceal visualisation rates vary significantly, as does the management of equivocal ultrasound findings. The aim of this study was to audit appendiceal sonography at a tertiary children's hospital, and provide baseline data for a future prospective study. - Summary of work Records of children who underwent ultrasound studies for possible appendicitis between January 2008 and December 2010 were reviewed. Variables included patient demographics, sonographic appendix characteristics, and secondary signs. Descriptive statistics and analysis using ANOVA, Mann-Whitney U test, and ROC curves were performed. Mater Human Research Ethic Committee approval was granted. - Summary of results There were 457 eligible children. Using a dichotomous diagnostic model (including equivocal results), sensitivity was 89.6%, specificity 91.6%, and diagnostic yield of 40.7%. ROC curve analysis of a 6mm diameter cut-off was 0.88 AUC (95% CI 0.80 to 0.95). - Discussion and conclusions Sonography is an accurate test for acute appendicitis in children, with a high sensitivity and negative predictive value. A diameter of 6mm as an absolute cut-off in a binary model can lead to false findings. Results were compared with available literature. Recent publications propose categorising diameter1 and integrating secondary signs2 to improve accuracy and provide more meaningful results to clinicians. This study will be a benchmark for future studies with multiple diagnostic categorisation.