Validation of a container ship model for parametric rolling

The objective of this work is to formulate a nonlinear, coupled model of a container ship during parametric roll resonance, and to validate the model using experimental data.

A new approach to spatial data interpolation using higher-order statistics

Interpolation techniques for spatial data have been applied frequently in various fields of geosciences. Although most conventional interpolation methods assume that it is sufficient to use first- and second-order statistics to characterize random fields, researchers have now realized that these methods cannot always provide reliable interpolation results, since geological and environmental phenomena tend to be very complex, presenting non-Gaussian distribution and/or non-linear inter-variable relationship. This paper proposes a new approach to the interpolation of spatial data, which can be applied with great flexibility. Suitable cross-variable higher-order spatial statistics are developed to measure the spatial relationship between the random variable at an unsampled location and those in its neighbourhood. Given the computed cross-variable higher-order spatial statistics, the conditional probability density function (CPDF) is approximated via polynomial expansions, which is then utilized to determine the interpolated value at the unsampled location as an expectation. In addition, the uncertainty associated with the interpolation is quantified by constructing prediction intervals of interpolated values. The proposed method is applied to a mineral deposit dataset, and the results demonstrate that it outperforms kriging methods in uncertainty quantification. The introduction of the cross-variable higher-order spatial statistics noticeably improves the quality of the interpolation since it enriches the information that can be extracted from the observed data, and this benefit is substantial when working with data that are sparse or have non-trivial dependence structures.

Exact parametric confidence intervals for Bland-Altman limits of agreement

Purpose The previous literature on Bland-Altman analysis only describes approximate methods for calculating confidence intervals for 95% Limits of Agreement (LoAs). This paper describes exact methods for calculating such confidence intervals, based on the assumption that differences in measurement pairs are normally distributed. Methods Two basic situations are considered for calculating LoA confidence intervals: the first where LoAs are considered individually (i.e. using one-sided tolerance factors for a normal distribution); and the second, where LoAs are considered as a pair (i.e. using two-sided tolerance factors for a normal distribution). Equations underlying the calculation of exact confidence limits are briefly outlined. Results To assist in determining confidence intervals for LoAs (considered individually and as a pair) tables of coefficients have been included for degrees of freedom between 1 and 1000. Numerical examples, showing the use of the tables for calculating confidence limits for Bland-Altman LoAs, have been provided. Conclusions Exact confidence intervals for LoAs can differ considerably from Bland and Altman’s approximate method, especially for sample sizes that are not large. There are better, more precise methods for calculating confidence intervals for LoAs than Bland and Altman’s approximate method, although even an approximate calculation of confidence intervals for LoAs is likely to be better than none at all. Reporting confidence limits for LoAs considered as a pair is appropriate for most situations, however there may be circumstances where it is appropriate to report confidence limits for LoAs considered individually.

Poverty and non-communicable diseases in South Africa

Background: High levels of wealth inequality with improved health statistics in South Africa (SA) provide an important opportunity to investigate non-communicable diseases (NCDs) among the poor. Aims: This paper uses two distinct national data sets to contrast patterns of mortality in rich and poor areas and explore the associations between poverty, risk factors, health care and selected NCDs diseases in South African adults. Methods: Causes of premature mortality in 1996 experienced in the poorest magisterial districts are compared with those in the richest, using average household wealth to classify districts. Logistic and multinomial regression are used to investigate the association of a household asset index and selected chronic conditions, related risk factors and healthcare indicators using data from the 1998 South African Demographic and Health Survey. Results: NCDs accounted for 39% and 33% of premature mortality in rich and poor districts respectively. The household survey data showed that the risk factors hypertension and obesity increased with increasing wealth, while most of the lifestyle factors, such as light smoking, domestic exposure to ``smoky'' fuels and alcohol dependence were associated with poverty. Treatment status for hypertension and asthma was worse for poor people than for rich people. Conclusions: The study suggests that NCDs and lifestyle-related risk factors are prevalent among the poor in SA and treatment for chronic diseases is lacking for poor people. The observed increase in hypertension and obesity with wealth suggests that unless comprehensive health promotion strategies are implemented, there will be an unmanageable chronic disease epidemic with future socioeconomic development in SA.

Demographics and discharge outcomes of dysvascular and non-vascular lower limb amputees at a subacute rehabilitation unit: A 7-year series

Objective To examine personal and social demographics, and rehabilitation discharge outcomes of dysvascular and non-vascular lower limb amputees. Methods In total, 425 lower limb amputation inpatient rehabilitation admissions (335 individuals) from 2005 to 2011 were examined. Admission and discharge descriptive statistics (frequency, percentages) were calculated and compared by aetiology. Results Participants were male (74%), aged 65 years (s.d. 14), born in Australia (72%), had predominantly dysvascular aetiology (80%) and a median length of stay 48 days (interquartile range (IQR): 25–76). Following amputation, 56% received prostheses for mobility, 21% (n = 89) changed residence and 28% (n = 116) required community services. Dysvascular amputees were older (mean 67 years, s.d. 12 vs 54 years, s.d. 16; P < 0.001) and recorded lower functional independence measure – motor scores at admission (z = 3.61, P < 0.001) and discharge (z = 4.52, P < 0.001). More nonvascular amputees worked before amputation (43% vs 11%; P < 0.001), were prescribed a prosthesis by discharge (73% vs 52%; P < 0.001) and had a shorter length of stay (7 days, 95% confidence interval: –3 to 17), although this was not statistically significant. Conclusions Differences exist in social and demographic outcomes between dysvascular and non-vascular lower limb amputees.

Young and unaffected by road policing strategies: Using deterrence theory to explain provisional drivers' (non)compliance

Newly licenced drivers are disproportionately represented in traffic injuries and crash statistics. Despite the implementation of countermeasures designed to improve safety, such as graduated driver licencing (GDL) schemes, many young drivers do not comply with road rules. This study used a reconceptualised deterrence theory framework to investigate young drivers’ perceptions of the enforcement of road rules in general and those more specifically related to GDL. A total of 236 drivers aged 17–24 completed a questionnaire assessing their perceptions of various deterrence mechanisms (personal and vicarious) and their compliance with both GDL-specific and general road rules. Hierarchical multiple regressions conducted to explore noncompliant behaviour revealed that, contrary to theoretical expectations, neither personal nor vicarious punishment experiences affected compliance in the expected direction. Instead, the most influential factors contributing to noncompliance were licence type (P2) and, counterintuitively, having previously been exposed to enforcement. Parental enforcement was also significant in the prediction of transient rule violations, but not fixed rule violations or overall noncompliance. Findings are discussed in light of several possibilities, including an increase in violations due to more time spent on the road, an ‘emboldening effect’ noted in prior studies and possible conceptual constraints regarding the deterrence variables examined in this study.

Stochastic analysis of a radial-inflow turbine in the presence of parametric uncertainties

This paper presents an uncertainty quantification study of the performance analysis of the high pressure ratio single stage radial-inflow turbine used in the Sundstrand Power Systems T-100 Multi-purpose Small Power Unit. A deterministic 3D volume-averaged Computational Fluid Dynamics (CFD) solver is coupled with a non-statistical generalized Polynomial Chaos (gPC) representation based on a pseudo-spectral projection method. One of the advantages of this approach is that it does not require any modification of the CFD code for the propagation of random disturbances in the aerodynamic and geometric fields. The stochastic results highlight the importance of the blade thickness and trailing edge tip radius on the total-to-static efficiency of the turbine compared to the angular velocity and trailing edge tip length. From a theoretical point of view, the use of the gPC representation on an arbitrary grid also allows the investigation of the sensitivity of the blade thickness profiles on the turbine efficiency. The gPC approach is also applied to coupled random parameters. The results show that the most influential coupled random variables are trailing edge tip radius coupled with the angular velocity.

Motor response suppression and the prepotent tendency to respond: A parametric fMRI study

In the present study we utilised functional magnetic resonance imaging (fMRI) to examine cerebral activation during performance of a classic motor task in which response suppression load was parametrically varied. Linear increases in activity were observed in a distributed network of regions across both cerebral hemispheres, although with more extensive involvement of the right prefrontal cortex. Activated regions included prefrontal, parietal and occipitotemporal cortices. Decreasing activation was similarly observed in a distributed network of regions. These response forms are discussed in terms of an increasing requirement for visual cue discrimination and suppression/selection of motor responses, and a decreasing probability of the occurrence of non-target stimuli and attenuation of a prepotent tendency to respond. The results support recent proposals for a dominant role for the right-hemisphere in performance of motor response suppression tasks that emphasise the importance of the right prefrontal cortex.

Bayesian parametric bootstrap for models with intractable likelihoods

In this paper it is demonstrated how the Bayesian parametric bootstrap can be adapted to models with intractable likelihoods. The approach is most appealing when the semi-automatic approximate Bayesian computation (ABC) summary statistics are selected. After a pilot run of ABC, the likelihood-free parametric bootstrap approach requires very few model simulations to produce an approximate posterior, which can be a useful approximation in its own right. An alternative is to use this approximation as a proposal distribution in ABC algorithms to make them more efficient. In this paper, the parametric bootstrap approximation is used to form the initial importance distribution for the sequential Monte Carlo and the ABC importance and rejection sampling algorithms. The new approach is illustrated through a simulation study of the univariate g-and- k quantile distribution, and is used to infer parameter values of a stochastic model describing expanding melanoma cell colonies.

Whole-genome screening in ankylosing spondylitis: Evidence of non-MHC genetic-susceptibility loci

Ankylosing spondylitis (AS) is a common inflammatory arthritis predominantly affecting the axial skeleton. Susceptibility to the disease is thought to be oligogenic. To identify the genes involved, we have performed a genomewide scan in 185 families containing 255 affected sibling pairs. Two-point and multipoint nonparametric linkage analysis was performed. Regions were identified showing "suggestive" or stronger linkage with the disease on chromosomes 1p, 2q, 6p, 9q, 10q, 16q, and 19q. The MHC locus was identified as encoding the greatest component of susceptibility, with an overall LOD score of 15.6. The strongest non-MHC linkage lies on chromosome 16q (overall LOD score 4.7). These results strongly support the presence of non-MHC genetic-susceptibility factors in AS and point to their likely locations.

Development of near infrared analysis of faeces to estimate non-grass proportions in diets selected by cattle grazing tropical pastures

Grass (monocots) and non-grass (dicots) proportions in ruminant diets are important nutritionally because the non-grasses are usually higher in nutritive value, particularly protein, than the grasses, especially in tropical pastures. For ruminants grazing tropical pastures where the grasses are C-4 species and most non-grasses are C-3 species, the ratio of C-13/C-12 in diet and faeces, measured as delta C-13 parts per thousand, is proportional to dietary non-grass%. This paper describes the development of a faecal near infrared (NIR) spectroscopy calibration equation for predicting faecal delta C-13 from which dietary grass and non-grass proportions can be calculated. Calibration development used cattle faeces derived from diets containing only C-3 non-grass and C-4 grass components, and a series of expansion and validation steps was employed to develop robustness and predictive reliability. The final calibration equation contained 1637 samples and faecal delta C-13 range (parts per thousand) of [12.27]-[27.65]. Calibration statistics were: standard error of calibration (SEC) of 0.78, standard error of cross-validation (SECV) of 0.80, standard deviation (SD) of reference values of 3.11 and R-2 of 0.94. Validation statistics for the final calibration equation applied to 60 samples were: standard error of prediction (SEP) of 0.87, bias of -0.15, R-2 of 0.92 and RPD of 3.16. The calibration equation was also tested on faeces from diets containing C-4 non-grass species or temperate C-3 grass species. Faecal delta C-13 predictions indicated that the spectral basis of the calibration was not related to C-13/C-12 ratios per se but to consistent differences between grasses and non-grasses in chemical composition and that the differences were modified by photosynthetic pathway. Thus, although the calibration equation could not be used to make valid faecal delta C-13 predictions when the diet contained either C-3 grass or C-4 non-grass, it could be used to make useful estimates of dietary non-grass proportions. It could also be ut :sed to make useful estimates of non-grass in mixed C-3 grass/non-grass diets by applying a modified formula to calculate non-grass from predicted faecal delta C-13. The development of a robust faecal-NIR calibration equation for estimating non-grass proportions in the diets of grazing cattle demonstrated a novel and useful application of NIR spectroscopy in agriculture.

Dynamic Bayesian network inferencing for non-homogeneous complex systems

Dynamic Bayesian Networks (DBNs) provide a versatile platform for predicting and analysing the behaviour of complex systems. As such, they are well suited to the prediction of complex ecosystem population trajectories under anthropogenic disturbances such as the dredging of marine seagrass ecosystems. However, DBNs assume a homogeneous Markov chain whereas a key characteristics of complex ecosystems is the presence of feedback loops, path dependencies and regime changes whereby the behaviour of the system can vary based on past states. This paper develops a method based on the small world structure of complex systems networks to modularise a non-homogeneous DBN and enable the computation of posterior marginal probabilities given evidence in forwards inference. It also provides an approach for an approximate solution for backwards inference as convergence is not guaranteed for a path dependent system. When applied to the seagrass dredging problem, the incorporation of path dependency can implement conditional absorption and allows release from the zero state in line with environmental and ecological observations. As dredging has a marked global impact on seagrass and other marine ecosystems of high environmental and economic value, using such a complex systems model to develop practical ways to meet the needs of conservation and industry through enhancing resistance and/or recovery is of paramount importance.

Predicting Additive and Non-additive Genetic Effects from Trials Where Traits Are Affected by Interplot Competition

There are two key types of selection in a plant breeding program, namely selection of hybrids for potential commercial use and the selection of parents for use in future breeding. Oakey et al. (in Theoretical and Applied Genetics 113, 809-819, 2006) showed how both of these aims could be achieved using pedigree information in a mixed model analysis in order to partition genetic effects into additive and non-additive effects. Their approach was developed for field trial data subject to spatial variation. In this paper we extend the approach for data from trials subject to interplot competition. We show how the approach may be used to obtain predictions of pure stand additive and non-additive effects. We develop the methodology in the context of a single field trial using an example from an Australian sorghum breeding program.