A pseudo-marginal sequential Monte Carlo algorithm for random effects models in Bayesian sequential design

A computationally efficient sequential Monte Carlo algorithm is proposed for the sequential design of experiments for the collection of block data described by mixed effects models. The difficulty in applying a sequential Monte Carlo algorithm in such settings is the need to evaluate the observed data likelihood, which is typically intractable for all but linear Gaussian models. To overcome this difficulty, we propose to unbiasedly estimate the likelihood, and perform inference and make decisions based on an exact-approximate algorithm. Two estimates are proposed: using Quasi Monte Carlo methods and using the Laplace approximation with importance sampling. Both of these approaches can be computationally expensive, so we propose exploiting parallel computational architectures to ensure designs can be derived in a timely manner. We also extend our approach to allow for model uncertainty. This research is motivated by important pharmacological studies related to the treatment of critically ill patients.

Understanding public response to a congestion charge : a random-effects ordered logit approach

Public acceptance is consistently listed as having an enormous impact on the implementation and success of a congestion charge scheme. This paper investigates public acceptance of such a scheme in Australia. Surveys were conducted in Brisbane and Melbourne, the two fastest growing Australian cities. Using an ordered logit modeling approach, the survey data including stated preferences were analyzed to pinpoint the important factors influencing people’s attitudes to a congestion charge and, in turn, to their transport mode choices. To accommodate the nature of, and to account for the resulting heterogeneity of the panel data, random effects were considered in the models. As expected, this study found that the amount of the congestion charge and the financial benefits of implementing it have a significant influence on respondents’ support for the charge and on the likelihood of their taking a bus to city areas. However, respondents’ current primary transport mode for travelling to the city areas has a more pronounced impact. Meanwhile, respondents’ perceptions of the congestion charge’s role in protecting the environment by reducing vehicle emissions, and of the extent to which the charge would mean that they travelled less frequently to the city for shopping or entertainment, also have a significant impact on their level of support for its implementation. We also found and explained notable differences across two cities. Finally, findings from this study have been fully discussed in relation to the literature.

Can measures of limb loading and dynamic stability during the squat maneuver provide an index of early functional recovery following unilateral total hip arthroplasty?

Objective: To investigate limb loading and dynamic stability during squatting in the early functional recovery of total hip arthroplasty (THA) patients. Design: Cohort study Setting: Inpatient rehabilitation clinic. Participants: A random sample of 61 THA patients (34♂/27♀; 62±9 yrs, 77±14 kg, 174±9 cm) was assessed twice, 13.2±3.8 days (PRE) and 26.6±3.3 days post-surgery (POST), and compared with a healthy reference group (REF) (22♂/16♀; 47±12yrs; 78±20kg; 175±10cm). Interventions: THA patients received two weeks of standard in-patient rehabilitation. Main Outcome Measure(s): Inter-limb vertical force distribution and dynamic stability during the squat maneuver, as defined by the root mean square (RMS) of the center of pressure in antero-posterior and medio-lateral directions, of operated (OP) and non-operated (NON)limbs. Self-reported function was assessed via FFb-H-OA 2.0 questionnaire. Results: At PRE, unloading of the OP limb was 15.8% greater (P<.001, d=1.070) and antero-posterior and medio-lateral center of pressure RMS were 30-34% higher in THA than REF P<.05). Unloading was reduced by 12.8% towards a more equal distribution from PRE to POST (P<.001, d=0.874). Although medio-lateral stability improved between PRE and POST (OP: 14.8%, P=.024, d=0.397; NON: 13.1%, P=.015, d=0.321), antero-posterior stability was not significantly different. Self-reported physical function improved by 15.8% (P<.001, d=0.965). Conclusion(s): THA patients unload the OP limb and are dynamically more unstable during squatting in the early rehabilitation phase following total hip replacement than healthy adults. Although loading symmetry and medio-lateral stability improved to the level of healthy adults with rehabilitation, antero-posterior stability remained impaired. Measures of dynamic stability and load symmetry during squatting provide quantitative information that can be used to clinically monitor early functional recovery from THA.

Thrust allocation with linear constrained quadratic cost function

In this paper, a method of thrust allocation based on a linearly constrained quadratic cost function capable of handling rotating azimuths is presented. The problem formulation accounts for magnitude and rate constraints on both thruster forces and azimuth angles. The advantage of this formulation is that the solution can be found with a finite number of iterations for each time step. Experiments with a model ship are used to validate the thrust allocation system.

Visual function and the ageing visual system

The population is ageing. Globally, the number of older adults (aged 60 years or over) is expected to more than double, from 841 million people in 2013 to more than 2 billion in 2050.1 In light of the increasing size of the older adult population, there is a pressing need to better identify the nature of, and mechanisms underlying, age-related vision impairment and the functional impact it has on the performance of everyday activities in older adults. The content of this feature issue reflects the diversity of research currently being undertaken on the topic of the ageing visual system and the important visual challenges that this presents for our ageing patient population. The scope is broad and includes topics relating to three main related themes: 1) The treatment of age-related ocular disorders and diseases and their consequences, including presbyopia and AMD; 2) The impact of age-related visual changes on everyday activities in older people, including mobility, driving and falls, and; 3) The interaction of age-related visual impairments and other age-related impairments including hearing and cognitive changes.

Theoretical modelling of momentum transfer function of bi-disperse porous media

The aim of this paper is to obtain the momentum transfer coefficient between the two phases, denoted by f and p, occupying a bi-disperse porous medium by mapping the available experimental data to the theoretical model proposed by Nield and Kuznetsov. Data pertinent to plate-fin heat exchangers, as bi-disperse porous media, were used. The measured pressure drops for such heat exchangers are then used to give the overall permeability which is linked to the porosity and permeability of each phase as well as the interfacial momentum transfer coefficient between the two phases. Accordingly, numerical values are obtained for the momentum transfer coefficient for three different fin spacing values considered in the heat exchanger experiments.

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.

Factors influencing robustness and effectiveness of conditional random fields in active learning frameworks

Active learning approaches reduce the annotation cost required by traditional supervised approaches to reach the same effectiveness by actively selecting informative instances during the learning phase. However, effectiveness and robustness of the learnt models are influenced by a number of factors. In this paper we investigate the factors that affect the effectiveness, more specifically in terms of stability and robustness, of active learning models built using conditional random fields (CRFs) for information extraction applications. Stability, defined as a small variation of performance when small variation of the training data or a small variation of the parameters occur, is a major issue for machine learning models, but even more so in the active learning framework which aims to minimise the amount of training data required. The factors we investigate are a) the choice of incremental vs. standard active learning, b) the feature set used as a representation of the text (i.e., morphological features, syntactic features, or semantic features) and c) Gaussian prior variance as one of the important CRFs parameters. Our empirical findings show that incremental learning and the Gaussian prior variance lead to more stable and robust models across iterations. Our study also demonstrates that orthographical, morphological and contextual features as a group of basic features play an important role in learning effective models across all iterations.

Comparison of ocular modulation transfer function determined by a ray-tracing aberrometer and a double-pass system in early cataract patients

BACKGROUND: The evaluation of retinal image quality in cataract eyes has gained importance and the clinical modulation transfer functions (MTF) can obtained by aberrometer and double pass (DP) system. This study aimed to compare MTF derived from a ray tracing aberrometer and a DP system in early cataractous and normal eyes. METHODS: There were 128 subjects with 61 control eyes and 67 eyes with early cataract defined according to the Lens Opacities Classification System III. A laser ray-tracing wavefront aberrometer (iTrace) and a double pass (DP) system (OQAS) assessed ocular MTF for 6.0 mm pupil diameters following dilation. Areas under the MTF (AUMTF) and their correlations were analyzed. Stepwise multiple regression analysis assessed factors affecting the differences between iTrace- and OQAS-derived AUMTF for the early cataract group. RESULTS: For both early cataract and control groups, iTrace-derived MTFs were higher than OQAS-derived MTFs across a range of spatial frequencies (P < 0.01). No significant difference between the two groups occurred for iTrace-derived AUMTF, but the early cataract group had significantly smaller OQAS-derived AUMTF than did the control group (P < 0.01). AUMTF determined from both the techniques demonstrated significant correlations with nuclear opacities, higher-order aberrations (HOAs), visual acuity, and contrast sensitivity functions, while the OQAS-derived AUMTF also demonstrated significant correlations with age and cortical opacity grade. The factors significantly affecting the difference between iTrace and OQAS AUMTF were root-mean-squared HOAs (standardized beta coefficient = -0.63, P < 0.01) and age (standardized beta coefficient = 0.26, P < 0.01). CONCLUSIONS: MTFs determined from a iTrace and a DP system (OQAS) differ significantly in early cataractous and normal subjects. Correlations with visual performance were higher for the DP system. OQAS-derived MTF may be useful as an indicator of visual performance in early cataract eyes.

Cold water immersion enhances recovery of submaximal muscle function after resistance exercise

We investigated the effect of cold water immersion (CWI) on the recovery of muscle function and physiological responses following high-intensity resistance exercise. Using a randomized, cross-over design, 10 physically active men performed high-intensity resistance exercise, followed by one of two recovery interventions: 10 min of cold water immersion at 10°C, or 10 min active recovery (low-intensity cycling). After the recovery interventions, maximal muscle function was assessed after 2 h and 4 h by measuring jump height and isometric squat strength. Submaximal muscle function was assessed after 6 h by measuring the average load lifted during six sets of 10 squats at 80% 1RM. Intramuscular temperature (1 cm) was also recorded, and venous blood samples were analyzed for markers of metabolism, vasoconstriction and muscle damage. CWI did not enhance recovery of maximal muscle function. However, during the final three sets of the submaximal muscle function test, the participants lifted a greater load (p<0.05; 38%; Cohen’s d 1.3) following CWI compared with active recovery. During CWI, muscle temperature decreased 6°C below post-exercise values, and remained below pre-exercise values for another 35 min. Venous blood O2 saturation decreased below pre-exercise values for 1.5 h after CWI. Serum endothelin-1 concentration did not change after CWI, whereas it decreased after active recovery. Plasma myoglobin concentration was lower, whereas plasma interleukin-6 concentration was higher after CWI compared with active recovery. These results suggest that cold water immersion after resistance exercise allow athletes to complete more work during subsequent training sessions, which could enhance long-term training adaptations.

Effect of tocopherol on atherosclerosis, vascular function and inflammation in Apolipoprotein E knockout mice with subtotal nephrectomy

Background/Aims: Inflammation and endothelial dysfunction contribute to cardiovascular disease, prevalent in chronic kidney disease (CKD). Antioxidant supplements such as tocopherols may reduce inflammation and atherosclerosis. This study aimed to investigate the effect of tocopherol supplementation on vascular function, aortic plaque formation, and inflammation in apolipoprotein E−/− mice with 5/6 nephrectomy as a model of combined cardiovascular and kidney disease. Methods: Nephrectomized mice were assigned to a normal chow diet group (normal chow), a group receiving 1000 mg/kg diet of α-tocopherol supplementation or a group receiving 1000 mg/kg diet mixed-tocopherol (60% γ-tocopherol). Results: Following 12 weeks, in vitro aortic endothelial-independent relaxation was enhanced with both α-tocopherol and mixed-tocopherol (P < 0.05), while mixed-tocopherol enhanced aortic contraction at noradrenaline concentrations of 3 × 10−7 M to 3 × 10−5 M (P < 0.05), when compared to normal chow. Supplementation with α- and mixed-tocopherol reduced systemic concentrations of IL-6 (P < 0.001 and P < 0.001, respectively) and IL-10 (P < 0.05 and P < 0.001, respectively), while α-tocopherol also reduced MCP-1 (P < 0.05) and tumor necrosis factor (TNF)-α (P < 0.05). Aortic sinus plaque area was significantly reduced with α-tocopherol supplementation when compared to normal chow (P < 0.01). Conclusion: Tocopherol supplementation favorably influenced vascular function and cytokine profile, while it was also effective in reducing atherosclerosis in the apolipoprotein E−/− mouse with CKD.

Pre-processing for approximate Bayesian computation in image analysis

Most of the existing algorithms for approximate Bayesian computation (ABC) assume that it is feasible to simulate pseudo-data from the model at each iteration. However, the computational cost of these simulations can be prohibitive for high dimensional data. An important example is the Potts model, which is commonly used in image analysis. Images encountered in real world applications can have millions of pixels, therefore scalability is a major concern. We apply ABC with a synthetic likelihood to the hidden Potts model with additive Gaussian noise. Using a pre-processing step, we fit a binding function to model the relationship between the model parameters and the synthetic likelihood parameters. Our numerical experiments demonstrate that the precomputed binding function dramatically improves the scalability of ABC, reducing the average runtime required for model fitting from 71 hours to only 7 minutes. We also illustrate the method by estimating the smoothing parameter for remotely sensed satellite imagery. Without precomputation, Bayesian inference is impractical for datasets of that scale.

Acoustic identification of 12 species of echolocating bat by discriminant function analysis and artificial neural networks

We recorded echolocation calls from 14 sympatric species of bat in Britain. Once digitised, one temporal and four spectral features were measured from each call. The frequency-time course of each call was approximated by fitting eight mathematical functions, and the goodness of fit, represented by the mean-squared error, was calculated. Measurements were taken using an automated process that extracted a single call from background noise and measured all variables without intervention. Two species of Rhinolophus were easily identified from call duration and spectral measurements. For the remaining 12 species, discriminant function analysis and multilayer back-propagation perceptrons were used to classify calls to species level. Analyses were carried out with and without the inclusion of curve-fitting data to evaluate its usefulness in distinguishing among species. Discriminant function analysis achieved an overall correct classification rate of 79% with curve-fitting data included, while an artificial neural network achieved 87%. The removal of curve-fitting data improved the performance of the discriminant function analysis by 2 %, while the performance of a perceptron decreased by 2 %. However, an increase in correct identification rates when curve-fitting information was included was not found for all species. The use of a hierarchical classification system, whereby calls were first classified to genus level and then to species level, had little effect on correct classification rates by discriminant function analysis but did improve rates achieved by perceptrons. This is the first published study to use artificial neural networks to classify the echolocation calls of bats to species level. Our findings are discussed in terms of recent advances in recording and analysis technologies, and are related to factors causing convergence and divergence of echolocation call design in bats.

A hybrid quality function deployment and cybernetic analytic network process model for project manager selection

Identifying appropriate decision criteria and making optimal decisions in a structured way is a complex process. This paper presents an approach for doing this in the form of a hybrid Quality Function Deployment (QFD) and Cybernetic Analytic Network Process (CANP) model for project manager selection. This involves the use of QFD to translate the owner's project management expectations into selection criteria and the CANP to weight the expectations and selection criteria. The supermatrix approach then prioritises the candidates with respect to the overall decision-making goal. A case study is used to demonstrate the use of the model in selecting a renovation project manager. This involves the development of 18 selection criteria in response to the owner's three main expectations of time, cost and quality.

Building indifferentiable compression functions from the PGV compression functions

Preneel, Govaerts and Vandewalle (PGV) analysed the security of single-block-length block cipher based compression functions assuming that the underlying block cipher has no weaknesses. They showed that 12 out of 64 possible compression functions are collision and (second) preimage resistant. Black, Rogaway and Shrimpton formally proved this result in the ideal cipher model. However, in the indifferentiability security framework introduced by Maurer, Renner and Holenstein, all these 12 schemes are easily differentiable from a fixed input-length random oracle (FIL-RO) even when their underlying block cipher is ideal. We address the problem of building indifferentiable compression functions from the PGV compression functions. We consider a general form of 64 PGV compression functions and replace the linear feed-forward operation in this generic PGV compression function with an ideal block cipher independent of the one used in the generic PGV construction. This modified construction is called a generic modified PGV (MPGV). We analyse indifferentiability of the generic MPGV construction in the ideal cipher model and show that 12 out of 64 MPGV compression functions in this framework are indifferentiable from a FIL-RO. To our knowledge, this is the first result showing that two independent block ciphers are sufficient to design indifferentiable single-block-length compression functions.