Greenhouse gas emissions from sub-tropical agricultural soils after addition of organic by-products

As the cost of mineral fertilisers increases globally, organic soil amendments (OAs) from agricultural sources are increasingly being used as substitutes for nitrogen. However, the impact of OAs on the production of greenhouse gases (CO2 and N2O) is not well understood. A 60-day laboratory incubation experiment was conducted to investigate the impacts of applying OAs (equivalent to 296 kg N ha−1 on average) on N2O and CO2 emissions and soil properties of clay and sandy loam soils from sugar cane production. The experiment included 6 treatments, one being an un-amended (UN) control with addition of five OAs being raw mill mud (MM), composted mill mud (CM), high N compost (HC), rice husk biochar (RB), and raw mill mud plus rice husk biochar (MB). These OAs were incubated at 60, 75 and 90% water-filled pore space (WFPS) at 25°C with urea (equivalent to 200 kg N ha−1) added to the soils thirty days after the incubation commenced. Results showed WFPS did not influence CO2 emissions over the 60 days but the magnitude of emissions as a proportion of C applied was RB < CM < MB < HC linear models being developed to predict CO2 and N2O emissions as a function of the dry matter and C/N ratio of the OAs, WFPS, and the soil CEC. Application of RB reduced N2O emissions by as much as 42-64% depending on WFPS. The reductions in both CO2 and N2O emissions after application of RB were due to a reduced bioavailability of C and not immobilisation of N. These findings show that the effect of OAs on soil GHG emissions can vary substantially depending on their chemical properties. OAs with a high availability of labile C and N can lead to elevated emissions of CO2 and N2O, while rice husk biochar showed potential in reducing overall soil GHG emissions.

Contributions to Bayesian experimental design

This thesis progresses Bayesian experimental design by developing novel methodologies and extensions to existing algorithms. Through these advancements, this thesis provides solutions to several important and complex experimental design problems, many of which have applications in biology and medicine. This thesis consists of a series of published and submitted papers. In the first paper, we provide a comprehensive literature review on Bayesian design. In the second paper, we discuss methods which may be used to solve design problems in which one is interested in finding a large number of (near) optimal design points. The third paper presents methods for finding fully Bayesian experimental designs for nonlinear mixed effects models, and the fourth paper investigates methods to rapidly approximate the posterior distribution for use in Bayesian utility functions.

Development of a log-quadratic model to describe microbial inactivation, illustrated by thermal inactivation of Clostridium botulinum

In the commercial food industry, demonstration of microbiological safety and thermal process equivalence often involves a mathematical framework that assumes log-linear inactivation kinetics and invokes concepts of decimal reduction time (DT), z values, and accumulated lethality. However, many microbes, particularly spores, exhibit inactivation kinetics that are not log linear. This has led to alternative modeling approaches, such as the biphasic and Weibull models, that relax strong log-linear assumptions. Using a statistical framework, we developed a novel log-quadratic model, which approximates the biphasic and Weibull models and provides additional physiological interpretability. As a statistical linear model, the log-quadratic model is relatively simple to fit and straightforwardly provides confidence intervals for its fitted values. It allows a DT-like value to be derived, even from data that exhibit obvious "tailing." We also showed how existing models of non-log-linear microbial inactivation, such as the Weibull model, can fit into a statistical linear model framework that dramatically simplifies their solution. We applied the log-quadratic model to thermal inactivation data for the spore-forming bacterium Clostridium botulinum and evaluated its merits compared with those of popular previously described approaches. The log-quadratic model was used as the basis of a secondary model that can capture the dependence of microbial inactivation kinetics on temperature. This model, in turn, was linked to models of spore inactivation of Sapru et al. and Rodriguez et al. that posit different physiological states for spores within a population. We believe that the log-quadratic model provides a useful framework in which to test vitalistic and mechanistic hypotheses of inactivation by thermal and other processes. Copyright © 2009, American Society for Microbiology. All Rights Reserved.

The contribution of three components of nutrition knowledge to socio-economic differences in food purchasing choices

Objectives: To assess socio-economic differences in three components of nutrition knowledge, i.e. knowledge of (i) the relationship between diet and disease, (ii) the nutrient content of foods and (iii) dietary guideline recommendations; furthermore, to determine if socio-economic differences in nutrition knowledge contribute to inequalities in food purchasing choices. Design: The cross-sectional study considered household food purchasing,nutrition knowledge, socio-economic and demographic information. Household food purchasing choices were summarised by three indices, based on self-reported purchasing of sixteen groceries, nineteen fruits and twenty-one vegetables. Socio-economic position (SEP) was measured by household income and education. Associations between SEP, nutrition knowledge and food purchasing were examined using general linear models adjusted for age, gender, household type and household size. Setting: Brisbane, Australia in 2000. Subjects: Main household food shoppers (n 1003, response rate 66?4 %), located in fifty small areas (Census Collectors Districts). Results: Shoppers in households of low SEP made food purchasing choices that were less consistent with dietary guideline recommendations: they were more likely to purchase grocery foods comparatively higher in salt, sugar and fat, and lower in fibre, and they purchased a narrower range of fruits and vegetables. Those of higher SEP had greater nutrition knowledge and this factor attenuated most associations between SEP and food purchasing choices. Among nutrition knowledge factors, knowledge of the relationship between diet and disease made the greatest and most consistent contribution to explaining socio-economic differences in food purchasing. Conclusions: Addressing inequalities in nutrition knowledge is likely to reduce socio-economic differences in compliance with dietary guidelines. Improving knowledge of the relationship between diet and disease appears to be a particularly relevant focus for health promotion aimed to reduce socio-economic differences in diet and related health inequalities.

Can perceptual indices estimate physiological strain across a range of environments and metabolic workloads when wearing explosive ordnance disposal and chemical protective clothing?

Objective Explosive ordnance disposal (EOD) often requires technicians to wear multiple protective garments in challenging environmental conditions. The accumulative effect of increased metabolic cost coupled with decreased heat dissipation associated with these garments predisposes technicians to high levels of physiological strain. It has been proposed that a perceptual strain index (PeSI) using subjective ratings of thermal sensation and perceived exertion as surrogate measures of core body temperature and heart rate, may provide an accurate estimation of physiological strain. Therefore, this study aimed to determine if the PeSI could estimate the physiological strain index (PSI) across a range of metabolic workloads and environments while wearing heavy EOD and chemical protective clothing. Methods Eleven healthy males wore an EOD and chemical protective ensemble while walking on a treadmill at 2.5, 4 and 5.5 km·h− 1 at 1% grade in environmental conditions equivalent to wet bulb globe temperature (WBGT) 21, 30 and 37 °C. WBGT conditions were randomly presented and a maximum of three randomised treadmill walking trials were completed in a single testing day. Trials were ceased at a maximum of 60-min or until the attainment of termination criteria. A Pearson's correlation coefficient, mixed linear model, absolute agreement and receiver operating characteristic (ROC) curves were used to determine the relationship between the PeSI and PSI. Results A significant moderate relationship between the PeSI and the PSI was observed [r = 0.77; p < 0.001; mean difference = 0.8 ± 1.1 a.u. (modified 95% limits of agreement − 1.3 to 3.0)]. The ROC curves indicated that the PeSI had a good predictive power when used with two, single-threshold cut-offs to differentiate between low and high levels of physiological strain (area under curve: PSI three cut-off = 0.936 and seven cut-off = 0.841). Conclusions These findings support the use of the PeSI for monitoring physiological strain while wearing EOD and chemical protective clothing. However, future research is needed to confirm the validity of the PeSI for active EOD technicians operating in the field.

Mandatory naptimes in child care and children's nighttime sleep

Objectives To examine the relationship between mandatory naptimes in child care and children's nighttime sleep duration, both concurrently and 12 months later once in school. Methods A sample of 168 children (50-72 months; 55% males) attending licensed child care centers were observed across their morning and throughout their scheduled naptime. Mandatory naptime was determined as the period in which children were not permitted any alternative activity except lying on their bed. Teachers reported each child's napping in child care. Nighttime and total sleep duration was reported by parents at 2 time points, in child care and in the second semester of their first school year. General linear models were used to examine group differences in sleep duration between children experiencing 0 to 60 minutes and >60 minutes of mandatory naptime, adjusting for key confounders. Path analysis was conducted to test a mediation model in which mandatory naptime is associated with nighttime sleep duration through increased napping in child care. Results Children who experienced >60 minutes of mandatory naptime in child care had significantly less nighttime sleep than those with 0 to 60 minutes of mandatory naptime. This difference persisted at 12-month follow-up, once children were in school. Napping in child care mediated the relationship between mandatory naptime and duration of nighttime sleep. Conclusions Exposure to mandatory naptimes of >60 minutes in child care is associated with decreased duration of nighttime sleep that endures beyond child care attendance. Given the large number of children who attend child care, sleep practices within these settings present an important focus for child health.

A robust optimization approach for index tracking problem

Index tracking is an investment approach where the primary objective is to keep portfolio return as close as possible to a target index without purchasing all index components. The main purpose is to minimize the tracking error between the returns of the selected portfolio and a benchmark. In this paper, quadratic as well as linear models are presented for minimizing the tracking error. The uncertainty is considered in the input data using a tractable robust framework that controls the level of conservatism while maintaining linearity. The linearity of the proposed robust optimization models allows a simple implementation of an ordinary optimization software package to find the optimal robust solution. The proposed model of this paper employs Morgan Stanley Capital International Index as the target index and the results are reported for six national indices including Japan, the USA, the UK, Germany, Switzerland and France. The performance of the proposed models is evaluated using several financial criteria e.g. information ratio, market ratio, Sharpe ratio and Treynor ratio. The preliminary results demonstrate that the proposed model lowers the amount of tracking error while raising values of portfolio performance measures.

Developments of the total entropy utility function for the dual purpose of model discrimination and parameter estimation in Bayesian design

The total entropy utility function is considered for the dual purpose of Bayesian design for model discrimination and parameter estimation. A sequential design setting is proposed where it is shown how to efficiently estimate the total entropy utility for a wide variety of data types. Utility estimation relies on forming particle approximations to a number of intractable integrals which is afforded by the use of the sequential Monte Carlo algorithm for Bayesian inference. A number of motivating examples are considered for demonstrating the performance of total entropy in comparison to utilities for model discrimination and parameter estimation. The results suggest that the total entropy utility selects designs which are efficient under both experimental goals with little compromise in achieving either goal. As such, the total entropy utility is advocated as a general utility for Bayesian design in the presence of model uncertainty.

Generalizing the use of geographical weights in biodiversity modelling

Aim Determining how ecological processes vary across space is a major focus in ecology. Current methods that investigate such effects remain constrained by important limiting assumptions. Here we provide an extension to geographically weighted regression in which local regression and spatial weighting are used in combination. This method can be used to investigate non-stationarity and spatial-scale effects using any regression technique that can accommodate uneven weighting of observations, including machine learning. Innovation We extend the use of spatial weights to generalized linear models and boosted regression trees by using simulated data for which the results are known, and compare these local approaches with existing alternatives such as geographically weighted regression (GWR). The spatial weighting procedure (1) explained up to 80% deviance in simulated species richness, (2) optimized the normal distribution of model residuals when applied to generalized linear models versus GWR, and (3) detected nonlinear relationships and interactions between response variables and their predictors when applied to boosted regression trees. Predictor ranking changed with spatial scale, highlighting the scales at which different species–environment relationships need to be considered. Main conclusions GWR is useful for investigating spatially varying species–environment relationships. However, the use of local weights implemented in alternative modelling techniques can help detect nonlinear relationships and high-order interactions that were previously unassessed. Therefore, this method not only informs us how location and scale influence our perception of patterns and processes, it also offers a way to deal with different ecological interpretations that can emerge as different areas of spatial influence are considered during model fitting.

A Model to Estimate the Population Contributing to the Wastewater Using Samples Collected on Census Day

An important uncertainty when estimating per capita consumption of, for example, illicit drugs by means of wastewater analysis (sometimes referred to as “sewage epidemiology”) relates to the size and variability of the de facto population in the catchment of interest. In the absence of a day-specific direct population count any indirect surrogate model to estimate population size lacks a standard to assess associated uncertainties. Therefore, the objective of this study was to collect wastewater samples at a unique opportunity, that is, on a census day, as a basis for a model to estimate the number of people contributing to a given wastewater sample. Mass loads for a wide range of pharmaceuticals and personal care products were quantified in influents of ten sewage treatment plants (STP) serving populations ranging from approximately 3500 to 500 000 people. Separate linear models for population size were estimated with the mass loads of the different chemical as the explanatory variable: 14 chemicals showed good, linear relationships, with highest correlations for acesulfame and gabapentin. De facto population was then estimated through Bayesian inference, by updating the population size provided by STP staff (prior knowledge) with measured chemical mass loads. Cross validation showed that large populations can be estimated fairly accurately with a few chemical mass loads quantified from 24-h composite samples. In contrast, the prior knowledge for small population sizes cannot be improved substantially despite the information of multiple chemical mass loads. In the future, observations other than chemical mass loads may improve this deficit, since Bayesian inference allows including any kind of information relating to population size.

Postnatal depressive symptoms amongst women in Central Vietnam: A cross-sectional study investigating prevalence and associations with social, cultural and infant factors

Background This study investigated the prevalence and socio-cultural correlates of postnatal mood disturbance amongst women 18–45 years old in Central Vietnam. Son preference and traditional confinement practices were explored as well as factors such as poverty, parity, family and intimate partner relationships and infant health. Methods A cross-sectional study was conducted in twelve randomly selected Commune Health Centres from urban and rural districts of Thua Thien Hue Province, Vietnam. Mother-infant dyads one to six months postpartum were invited to participate. Questionnaires from 431 mothers (urban n = 216; rural n = 215) assessed demographic and family characteristics, traditional confinement practices, son preference, infant health and social capital. The Edinburgh Postnatal Depression Scale (EPDS) and WHO5 Wellbeing Index indicated depressive symptoms and emotional wellbeing. Data were analysed using general linear models. Results Using an EPDS cut-off of 12/13, 18.1 % (n = 78, 95 % CI 14.6 - 22.1) of women had depressive symptoms (20.4 % urban; 15.8 % rural). Contrary to predictions, infant gender and traditional confinement were unrelated to depressive symptoms. Poverty, food insecurity, being frightened of family members, and intimate partner violence increased both depressive symptoms and lowered wellbeing. The first model accounted for 30.2 % of the variance in EPDS score and found being frightened of one’s husband, husband’s unemployment, breastfeeding difficulties, infant diarrhoea, and cognitive social capital were associated with higher EPDS scores. The second model had accounted for 22 % of the variance in WHO5 score. Living in Hue city, low education, poor maternal competence and a negative family response to the baby lowered maternal wellbeing. Conclusions Traditional confinement practices and son preference were not linked to depressive symptoms among mothers, but were correlates of family relationships and wellbeing. Poverty, food insecurity, violence, infant ill health, and discordant intimate and family relationships were linked with depressive symptoms in Central Vietnam.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P combined = 2.76 × 10 -17 and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species. © 2013 Nature America, Inc. All rights reserved.

Intra-cluster correlation structure in longitudinal data analysis: Selection criteria and misspecification tests

Selection criteria and misspecification tests for the intra-cluster correlation structure (ICS) in longitudinal data analysis are considered. In particular, the asymptotical distribution of the correlation information criterion (CIC) is derived and a new method for selecting a working ICS is proposed by standardizing the selection criterion as the p-value. The CIC test is found to be powerful in detecting misspecification of the working ICS structures, while with respect to the working ICS selection, the standardized CIC test is also shown to have satisfactory performance. Some simulation studies and applications to two real longitudinal datasets are made to illustrate how these criteria and tests might be useful.

Efficient estimation for rank-based regression with clustered data

Rank-based inference is widely used because of its robustness. This article provides optimal rank-based estimating functions in analysis of clustered data with random cluster effects. The extensive simulation studies carried out to evaluate the performance of the proposed method demonstrate that it is robust to outliers and is highly efficient given the existence of strong cluster correlations. The performance of the proposed method is satisfactory even when the correlation structure is misspecified, or when heteroscedasticity in variance is present. Finally, a real dataset is analyzed for illustration.

Rank regression for accelerated failure time model with clustered and censored data

For clustered survival data, the traditional Gehan-type estimator is asymptotically equivalent to using only the between-cluster ranks, and the within-cluster ranks are ignored. The contribution of this paper is two fold: - (i) incorporating within-cluster ranks in censored data analysis, and; - (ii) applying the induced smoothing of Brown and Wang (2005, Biometrika) for computational convenience. Asymptotic properties of the resulting estimating functions are given. We also carry out numerical studies to assess the performance of the proposed approach and conclude that the proposed approach can lead to much improved estimators when strong clustering effects exist. A dataset from a litter-matched tumorigenesis experiment is used for illustration.