Understanding Uncertainties in Non-Linear Population Trajectories: A Bayesian Semi-Parametric Hierarchical Approach to Large-Scale Surveys of Coral Cover

Recently, attempts to improve decision making in species management have focussed on uncertainties associated with modelling temporal fluctuations in populations. Reducing model uncertainty is challenging; while larger samples improve estimation of species trajectories and reduce statistical errors, they typically amplify variability in observed trajectories. In particular, traditional modelling approaches aimed at estimating population trajectories usually do not account well for nonlinearities and uncertainties associated with multi-scale observations characteristic of large spatio-temporal surveys. We present a Bayesian semi-parametric hierarchical model for simultaneously quantifying uncertainties associated with model structure and parameters, and scale-specific variability over time. We estimate uncertainty across a four-tiered spatial hierarchy of coral cover from the Great Barrier Reef. Coral variability is well described; however, our results show that, in the absence of additional model specifications, conclusions regarding coral trajectories become highly uncertain when considering multiple reefs, suggesting that management should focus more at the scale of individual reefs. The approach presented facilitates the description and estimation of population trajectories and associated uncertainties when variability cannot be attributed to specific causes and origins. We argue that our model can unlock value contained in large-scale datasets, provide guidance for understanding sources of uncertainty, and support better informed decision making

A novel population balance model for the dilute acid hydrolysis of hemicellulose

Acid hydrolysis is a popular pretreatment for removing hemicellulose from lignocelluloses in order to produce a digestible substrate for enzymatic saccharification. In this work, a novel model for the dilute acid hydrolysis of hemicellulose within sugarcane bagasse is presented and calibrated against experimental oligomer profiles. The efficacy of mathematical models as hydrolysis yield predictors and as vehicles for investigating the mechanisms of acid hydrolysis is also examined. Experimental xylose, oligomer (degree of polymerisation 2 to 6) and furfural yield profiles were obtained for bagasse under dilute acid hydrolysis conditions at temperatures ranging from 110C to 170C. Population balance kinetics, diffusion and porosity evolution were incorporated into a mathematical model of the acid hydrolysis of sugarcane bagasse. This model was able to produce a good fit to experimental xylose yield data with only three unknown kinetic parameters ka, kb and kd. However, fitting this same model to an expanded data set of oligomeric and furfural yield profiles did not successfully reproduce the experimental results. It was found that a ``hard-to-hydrolyse'' parameter, $\alpha$, was required in the model to ensure reproducibility of the experimental oligomer profiles at 110C, 125C and 140C. The parameters obtained through the fitting exercises at lower temperatures were able to be used to predict the oligomer profiles at 155C and 170C with promising results. The interpretation of kinetic parameters obtained by fitting a model to only a single set of data may be ambiguous. Although these parameters may correctly reproduce the data, they may not be indicative of the actual rate parameters, unless some care has been taken to ensure that the model describes the true mechanisms of acid hydrolysis. It is possible to challenge the robustness of the model by expanding the experimental data set and hence limiting the parameter space for the fitting parameters. The novel combination of ``hard-to-hydrolyse'' and population balance dynamics in the model presented here appears to stand up to such rigorous fitting constraints.

Colonoscopy requirements of population screening for colorectal cancer in New Zealand

Aim: To estimate the colonoscopy burden of introducing population screening for colorectal cancer in New Zealand. Methods: Screening for colorectal cancer using biennial immunochemical faecal occult blood tests offered to people aged 50-74 years of age was modelled using population estimates from Statistics New Zealand for 2011-2031. Modelling to determine colonoscopy requirements was based on participation and test positivity rates from published results of screening programmes. Estimates of the number of procedures required for ongoing adenoma surveillance were calculated using screening literature results of adenoma yield, and New Zealand Guidelines for Adenoma Surveillance. Sensitivity analysis was undertaken on key parameters. Results: For a test positivity of 6.4%, biennial screening using immunochemical faecal occult blood testing with a 60% participation rate, would require 18,000 colonoscopies nationally, increasing to 28,000 by 2031. The majority of procedures are direct referrals from a positive FOBT, with surveillance colonoscopy numbers building over time. Conclusion: Colonoscopy requirements for immunochemical faecal occult blood based population screening for colorectal cancer are high. Significant expansion of services is required and careful management of surveillance procedures to ensure timely delivery of initial colonoscopies whilst maintaining symptomatic services. A model re-run informed by data from the screening pilot will allow improved estimates for the New Zealand setting.

Estimation of cosmological parameters using adaptive importance sampling

We present a Bayesian sampling algorithm called adaptive importance sampling or population Monte Carlo (PMC), whose computational workload is easily parallelizable and thus has the potential to considerably reduce the wall-clock time required for sampling, along with providing other benefits. To assess the performance of the approach for cosmological problems, we use simulated and actual data consisting of CMB anisotropies, supernovae of type Ia, and weak cosmological lensing, and provide a comparison of results to those obtained using state-of-the-art Markov chain Monte Carlo (MCMC). For both types of data sets, we find comparable parameter estimates for PMC and MCMC, with the advantage of a significantly lower wall-clock time for PMC. In the case of WMAP5 data, for example, the wall-clock time scale reduces from days for MCMC to hours using PMC on a cluster of processors. Other benefits of the PMC approach, along with potential difficulties in using the approach, are analyzed and discussed.

Bayesian model comparison in cosmology with Population Monte Carlo

We use Bayesian model selection techniques to test extensions of the standard flat LambdaCDM paradigm. Dark-energy and curvature scenarios, and primordial perturbation models are considered. To that end, we calculate the Bayesian evidence in favour of each model using Population Monte Carlo (PMC), a new adaptive sampling technique which was recently applied in a cosmological context. The Bayesian evidence is immediately available from the PMC sample used for parameter estimation without further computational effort, and it comes with an associated error evaluation. Besides, it provides an unbiased estimator of the evidence after any fixed number of iterations and it is naturally parallelizable, in contrast with MCMC and nested sampling methods. By comparison with analytical predictions for simulated data, we show that our results obtained with PMC are reliable and robust. The variability in the evidence evaluation and the stability for various cases are estimated both from simulations and from data. For the cases we consider, the log-evidence is calculated with a precision of better than 0.08. Using a combined set of recent CMB, SNIa and BAO data, we find inconclusive evidence between flat LambdaCDM and simple dark-energy models. A curved Universe is moderately to strongly disfavoured with respect to a flat cosmology. Using physically well-motivated priors within the slow-roll approximation of inflation, we find a weak preference for a running spectral index. A Harrison-Zel'dovich spectrum is weakly disfavoured. With the current data, tensor modes are not detected; the large prior volume on the tensor-to-scalar ratio r results in moderate evidence in favour of r=0.

The tired foot. The effects of physical activity, energy expenditure and sleep on the diabetic foot population

Background The most common pathway to development of diabetes foot ulcers is repetitive daily activity stress on the plantar surface of the neuropathic foot. Studies suggest an association between different diabetic foot complications and physical activity. However, to the best of the authors knowledge the steps/day and sleep patterns of people with diabetic foot ulcers has yet to be investigated. This observational study aims to investigate the physical activity and sleep patterns of three groups of adults with type 2 diabetes and different foot complications Methods Participants with type 2 diabetes were recruited into three groups: 1. those with no reported foot complications (DNIL), 2. those with diagnosis of neuropathy (DPN) and 3. those with a neuropathic ulcer (DFU). Exclusion criteria included peripheral arterial disease and mobility aid use. Participants wore a SenseWear Pro 3 Armband continuously for 7 days and completed an Epworth Sleepiness Scale. The Armband is a validated automated measure of activity (walking steps, average Metabolic Equivalent Task (MET), physical activity (>3 METs) duration), energy expenditure(kJ) (total and physical activity (>3 METs)) and sleep (duration). Data on age, sex, BMI, diabetes duration and HbA1c were also collected. Results Sixty-Six (14 DNIL, 22 DPN and 30 DFU's participants were recruited; 71% males, mean age 61(±12) years, diabetes duration 13(±9) years, HbA1c 8.3(±2.8), BMI 32.6(±5.9), average METs 1.2(0.2). Significant differences were reported in mean(SD) steps/day (5,859(±2,381) in DNIL; 5,007(±3,349) in DPN and 3,271(±2,417) in DFU's and daily energy expenditure (10,868(±1,307)kJ in DNIL; 11,060(±1,916)kJ in DPN and 13,006(± 3,559) in DFU's(p <0.05). No significant differences were reported for average METs, physical activity duration or energy expenditure, sleep time or Epworth score (p>0.1). Conclusions Preliminary findings suggest people with diabetes are sedentary. Results indicate that patients with a diabetic foot ulcer work significantly less than those with neuropathy or nil complications and use significantly more energy to do so. Sleep Parameters showed no differences. Recruitment is still on going.

Reducing the global burden of depression : population level analysis of intervention cost-effectiveness in 14 world regions

International evidence on the cost and effects of interventions for reducing the global burden of depression remain scarce. Aims: To estimate the population-level cost-effectiveness of evidence-based depression interventions and their contribution towards reducing current burden. Method: Primary-care-based depression interventions were modelled at the level of whole populations in 14 epidemiological subregions of the world. Total population-level costs (in international dollars or I$) and effectiveness (disability adjusted life years (DALYs) averted) were combined to form average and incremental cost-effectiveness ratios. Results: Evaluated interventions have the potential to reduce the current burden of depression by 10–30%. Pharmacotherapy with older antidepressant drugs, with or without proactive collaborative care, are currently more cost-effective strategies than those using newer antidepressants, particularly in lower-income subregions. Conclusions: Even in resource-poor regions, each DALYaverted by efficient depression treatments in primary care costs less than 1 year of average per capita income, making such interventions a cost-effective use of health resources. However, current levels of burden can only be reduced significantlyif there is a substantialincrease substantial increase intreatment coverage.