Overfitting Bayesian mixture models with an unknown number of components

This paper proposes solutions to three issues pertaining to the estimation of finite mixture models with an unknown number of components: the non-identifiability induced by overfitting the number of components, the mixing limitations of standard Markov Chain Monte Carlo (MCMC) sampling techniques, and the related label switching problem. An overfitting approach is used to estimate the number of components in a finite mixture model via a Zmix algorithm. Zmix provides a bridge between multidimensional samplers and test based estimation methods, whereby priors are chosen to encourage extra groups to have weights approaching zero. MCMC sampling is made possible by the implementation of prior parallel tempering, an extension of parallel tempering. Zmix can accurately estimate the number of components, posterior parameter estimates and allocation probabilities given a sufficiently large sample size. The results will reflect uncertainty in the final model and will report the range of possible candidate models and their respective estimated probabilities from a single run. Label switching is resolved with a computationally light-weight method, Zswitch, developed for overfitted mixtures by exploiting the intuitiveness of allocation-based relabelling algorithms and the precision of label-invariant loss functions. Four simulation studies are included to illustrate Zmix and Zswitch, as well as three case studies from the literature. All methods are available as part of the R package Zmix, which can currently be applied to univariate Gaussian mixture models.

Accelerating Pseudo-Marginal MCMC using Gaussian Processes

Pseudo-marginal methods such as the grouped independence Metropolis-Hastings (GIMH) and Markov chain within Metropolis (MCWM) algorithms have been introduced in the literature as an approach to perform Bayesian inference in latent variable models. These methods replace intractable likelihood calculations with unbiased estimates within Markov chain Monte Carlo algorithms. The GIMH method has the posterior of interest as its limiting distribution, but suffers from poor mixing if it is too computationally intensive to obtain high-precision likelihood estimates. The MCWM algorithm has better mixing properties, but less theoretical support. In this paper we propose to use Gaussian processes (GP) to accelerate the GIMH method, whilst using a short pilot run of MCWM to train the GP. Our new method, GP-GIMH, is illustrated on simulated data from a stochastic volatility and a gene network model.

The impact of inter-fraction set-up errors on the probability of pulmonary and cardiac complication in left-sided breast cancer patients

Purpose This study evaluated the impact of patient set-up errors on the probability of pulmonary and cardiac complications in the irradiation of left-sided breast cancer. Methods and Materials Using the CMS XiO Version 4.6 (CMS Inc., St Louis, MO) radiotherapy planning system's NTCP algorithm and the Lyman -Kutcher-Burman (LKB) model, we calculated the DVH indices for the ipsilateral lung and heart and the resultant normal tissue complication probabilities (NTCP) for radiation-induced pneumonitis and excess cardiac mortality in 12 left-sided breast cancer patients. Results Isocenter shifts in the posterior direction had the greatest effect on the lung V20, heart V25, mean and maximum doses to the lung and the heart. Dose volume histograms (DVH) results show that the ipsilateral lung V20 tolerance was exceeded in 58% of the patients after 1cm posterior shifts. Similarly, the heart V25 tolerance was exceeded after 1cm antero-posterior and left-right isocentric shifts in 70% of the patients. The baseline NTCPs for radiation-induced pneumonitis ranged from 0.73% - 3.4% with a mean value of 1.7%. The maximum reported NTCP for radiation-induced pneumonitis was 5.8% (mean 2.6%) after 1cm posterior isocentric shift. The NTCP for excess cardiac mortality were 0 % in 100% of the patients (n=12) before and after setup error simulations. Conclusions Set-up errors in left sided breast cancer patients have a statistically significant impact on the Lung NTCPs and DVH indices. However, with a central lung distance of 3cm or less (CLD <3cm), and a maximum heart distance of 1.5cm or less (MHD<1.5cm), the treatment plans could tolerate set-up errors of up to 1cm without any change in the NTCP to the heart.

CBCT image guidance in head and neck irradiation: The impact of daily and weekly imaging protocols

Purpose This study evaluated the impact of a daily and weekly image-guided radiotherapy protocols in reducing setup errors and setting of appropriate margins in head and neck cancer patients. Materials and methods Interfraction and systematic shifts for the hypothetical day 1–3 plus weekly imaging were extrapolated from daily imaging data from 31 patients (964 cone beam computed tomography (CBCT) scans). In addition, residual setup errors were calculated by taking the average shifts in each direction for each patient based on the first three shifts and were presumed to represent systematic setup error. The clinical target volume (CTV) to planning target volume (PTV) margins were calculated using van Herk formula and analysed for each protocol. Results The mean interfraction shifts for daily imaging were 0·8, 0·3 and 0·5 mm in the S-I (superior-inferior), L-R (left-right) and A-P (anterior-posterior) direction, respectively. On the other hand the mean shifts for day 1–3 plus weekly imaging were 0·9, 1·8 and 0·5 mm in the S-I, L-R and A-P direction, respectively. The mean day 1–3 residual shifts were 1·5, 2·1 and 0·7 mm in the S-I, L-R and A-P direction, respectively. No significant difference was found in the mean setup error for the daily and hypothetical day 1–3 plus weekly protocol. However, the calculated CTV to PTV margins for the daily interfraction imaging data were 1·6, 3·8 and 1·4 mm in the S-I, L-R and A-P directions, respectively. Hypothetical day 1–3 plus weekly resulted in CTV–PTV margins of 5, 4·2 and 5 mm in the S-I, L-R and A-P direction. Conclusions The results of this study show that a daily CBCT protocol reduces setup errors and allows setup margin reduction in head and neck radiotherapy compared to a weekly imaging protocol.

Atmospheric visibility and PM10 as indicators of new particle formation in an urban environment

It is well-known that new particle formation (NPF) in the atmosphere is inhibited by pre-existing particles in the air that act as condensation sinks to decrease the concentration and, thus, the supersaturation of precursor gases. In this study, we investigate the effects of two parameters - atmospheric visibility, expressed as the particle back-scatter coefficient (BSP), and PM10 particulate mass concentration, on the occurrences of NPF events in an urban environment where the majority of precursor gases originate from motor vehicle and industrial sources. This is the first attempt to derive direct relationships between each of these two parameters and the occurrence of NPF. NPF events were identified from data obtained with a neutral cluster and air ion spectrometer over 245 days within a calendar year. Bayesian logistic regression was used to determine the probability of observing NPF as functions of BSP and PM10. We show that the BSP at 08 h on a given day is a reliable indicator of an NPF event later that day. The posterior median probability of observing an NPF event was greater than 0.5 (95%) when the BSP at 08 h was less than 6.8 Mm-1.

Hypoxic corneal changes following 8 hours of scleral contact lens wear

- Purpose To examine the change in corneal thickness and posterior curvature following 8 hours of miniscleral contact lens wear. - Methods Scheimpflug imaging (Pentacam HR, Oculus) was captured before, and immediately following, 8 hours of miniscleral contact lens wear for 15 young (mean age 22 ± 3 years), healthy participants with normal corneae. Natural diurnal variations were considered by measuring baseline corneal changes obtained on a separate control day without contact lens wear. - Results Over the central 6 mm of the cornea, a small, but highly statistically significant amount of edema was observed following 8 hours of miniscleral lens wear, after accounting for normal diurnal fluctuations (mean ± standard deviation percentage swelling 1.70 ± 0.98%, p < 0.0001). Posterior corneal topography remained stable following lens wear (-0.01 ± 0.07 mm steepening over the central 6 mm, p = 0.60). The magnitude of posterior corneal topographical changes following lens wear did not correlate with the extent of lens-related corneal edema (r = -0.16, p = 0.57). Similarly, the initial central corneal vault (maximum post-lens tear layer depth) was not associated with corneal swelling following lens removal (r = 0.27, p = 0.33). - Conclusions While a small amount of corneal swelling was induced following 8 hours of miniscleral lens wear (on average <2%), modern high Dk miniscleral contact lenses that vault the cornea do not induce clinically significant corneal edema or hypoxic related posterior corneal curvature changes during short-term wear. Longer-term studies of compromised eyes (e.g. corneal ectasia) are still required to inform the optimum lens and fitting characteristics for safe scleral lens wear to minimize corneal hypoxia.

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.

Bayesian estimation of the dynamics of pandemic (H1N1) 2009 influenza transmission in Queensland: A space–time SIR-based model

Background A pandemic strain of influenza A spread rapidly around the world in 2009, now referred to as pandemic (H1N1) 2009. This study aimed to examine the spatiotemporal variation in the transmission rate of pandemic (H1N1) 2009 associated with changes in local socio-environmental conditions from May 7–December 31, 2009, at a postal area level in Queensland, Australia. Method We used the data on laboratory-confirmed H1N1 cases to examine the spatiotemporal dynamics of transmission using a flexible Bayesian, space–time, Susceptible-Infected-Recovered (SIR) modelling approach. The model incorporated parameters describing spatiotemporal variation in H1N1 infection and local socio-environmental factors. Results The weekly transmission rate of pandemic (H1N1) 2009 was negatively associated with the weekly area-mean maximum temperature at a lag of 1 week (LMXT) (posterior mean: −0.341; 95% credible interval (CI): −0.370–−0.311) and the socio-economic index for area (SEIFA) (posterior mean: −0.003; 95% CI: −0.004–−0.001), and was positively associated with the product of LMXT and the weekly area-mean vapour pressure at a lag of 1 week (LVAP) (posterior mean: 0.008; 95% CI: 0.007–0.009). There was substantial spatiotemporal variation in transmission rate of pandemic (H1N1) 2009 across Queensland over the epidemic period. High random effects of estimated transmission rates were apparent in remote areas and some postal areas with higher proportion of indigenous populations and smaller overall populations. Conclusions Local SEIFA and local atmospheric conditions were associated with the transmission rate of pandemic (H1N1) 2009. The more populated regions displayed consistent and synchronized epidemics with low average transmission rates. The less populated regions had high average transmission rates with more variations during the H1N1 epidemic period.

Assessing the combined effect of asbestos exposure and smoking on lung cancer: a Bayesian approach

We review the literature on the combined association between lung cancer and two environmental exposures, asbestos exposure and smoking, and explore a Bayesian approach to assess evidence of interaction between the exposures. The meta-analysis combines separate indices of additive and multiplicative relationships and multivariate relative risk estimates. By making inferences on posterior probabilities we can explore both the form and strength of interaction. This analysis may be more informative than providing evidence to support one relation over another on the basis of statistical significance. Overall, we find evidence for a more than additive and less than multiplicative relation.

Morphological changes in the conjunctiva, episclera and sclera following short-term miniscleral contact lens wear in rigid lens neophytes

PURPOSE To quantify the influence of short-term wear of miniscleral contact lenses on the morphology of the corneo-scleral limbus, the conjunctiva, episclera and sclera. METHODS OCT images of the anterior eye were captured before, immediately following 3h of wear and then 3h after removal of a miniscleral contact lens for 10 young (27±5 years) healthy participants (neophyte rigid lens wearers). The region of analysis encompassed 1mm anterior, to 3.5mm posterior to the scleral spur. Natural diurnal variations in thickness were measured on a separate day and compensated for in subsequent analyses. RESULTS Following 3h of lens wear, statistically significant tissue thinning was observed across all quadrants, with a mean decrease in thickness of -24.1±3.6μm (p<0.001), which diminished, but did not return to baseline 3h after lens removal (-16.9±1.9μm, p<0.001). The largest tissue compression was observed in the superior quadrant (-49.9±8.5μm, p<0.01) and in the annular zone 1.5mm from the scleral spur (-48.2±5.7μm), corresponding to the approximate edge of the lens landing zone. Compression of the conjunctiva/episclera accounted for about 70% of the changes. CONCLUSIONS Optimal fitting miniscleral contact lenses worn for three hours resulted in significant tissue compression in young healthy eyes, with the greatest thinning observed superiorly, potentially due to the additional force of the eyelid, with a partial recovery of compression 3h after lens removal. Most of the morphological changes occur in the conjunctiva/episclera layers.

Evidence from molecules and morphology expands Podonomopsis Brundin (Diptera: Chironomidae: Podonominae) to include ‘genus Chile'

The informal taxon ‘genus Chile’ of Brundin, based solely on pupal exuviae of a podonomine Chironomidae, has remained inadequately known for half a century. New collections reveal life associations, and provide molecular data to hypothesise a precise phylogenetic placement in the austral Podonominae. A densely sampled molecular phylogeny based on two nuclear and one mitochondrial DNA markers shows ‘genus Chile’ to be the sister group to Podonomopsis Brundin, 1966. Within Podonomopsis a clade of South American species is sister to all Australian species. We discuss how to rank such a sister group taxon and treat ‘genus Chile’ as a new subgenus Araucanopsis, subg. nov. with the new species, Podonomopsis (Araucanopsis) avelasse, sp. nov. from Chile and Argentina as genotype of the monotypic subgenus. We describe P. (A.) avelasse in all stages and provide an expanded diagnosis and description of Podonomopsis to include Araucanopsis. A dated biogeographic hypothesis (chronogram) infers the most recent common ancestor (tmcra) of expanded Podonomopsis at 95 million years ago (Mya) (68–122 Mya 95% highest posterior density), ‘core’ Podonomopsis at 83 Mya (58–108) and Australian Podonomopsis at 65 Mya (44–87). All dates are before the South America–Australia geological separation through Antarctica, supporting previous conclusions that the taxon distribution is ‘Gondwanan’ in origin. Podonomopsis, even as expanded here, remains unknown from New Zealand or elsewhere on extant Zealandia.

The neuroscience of inspirational leadership: the importance of collective-oriented language and shared group membership

Effective leaders are believed to inspire followers by providing inclusive visions of the future that followers can identify with. In the present study, we examined the neural mechanisms underlying this process, testing key hypotheses derived from transformational and social identity approaches to leadership. While undergoing functional MRI, supporters from the two major Australian political parties (Liberal vs. Labor) were presented with inspirational collective-oriented and noninspirational personal-oriented statements made by in-group and out-group leaders. Imaging data revealed that inspirational (rather than noninspirational) statements from in-group leaders were associated with increased activation in the bilateral rostral inferior parietal lobule, pars opercularis, and posterior midcingulate cortex: brain areas that are typically implicated in controlling semantic information processing. In contrast, for out-group leaders, greater activation in these areas was associated with noninspirational statements. In addition, noninspirational statements by in-group (but not out-group) leaders resulted in increased activation in the medial prefrontal cortex, an area typically associated with reasoning about a person’s mental state. These results show that followers processed identical statements qualitatively differently as a function of leaders’ group membership, thus demonstrating that shared identity acts as an amplifier for inspirational leadership communication.

Tissue thickness calculation in ocular optical coherence tomography

Thickness measurements derived from optical coherence tomography (OCT) images of the eye are a fundamental clinical and research metric, since they provide valuable information regarding the eye’s anatomical and physiological characteristics, and can assist in the diagnosis and monitoring of numerous ocular conditions. Despite the importance of these measurements, limited attention has been given to the methods used to estimate thickness in OCT images of the eye. Most current studies employing OCT use an axial thickness metric, but there is evidence that axial thickness measures may be biased by tilt and curvature of the image. In this paper, standard axial thickness calculations are compared with a variety of alternative metrics for estimating tissue thickness. These methods were tested on a data set of wide-field chorio-retinal OCT scans (field of view (FOV) 60° x 25°) to examine their performance across a wide region of interest and to demonstrate the potential effect of curvature of the posterior segment of the eye on the thickness estimates. Similarly, the effect of image tilt was systematically examined with the same range of proposed metrics. The results demonstrate that image tilt and curvature of the posterior segment can affect axial tissue thickness calculations, while alternative metrics, which are not biased by these effects, should be considered. This study demonstrates the need to consider alternative methods to calculate tissue thickness in order to avoid measurement error due to image tilt and curvature.

Bayesian estimation of small effects in exercise and sports science

The aim of this paper is to provide a Bayesian formulation of the so-called magnitude-based inference approach to quantifying and interpreting effects, and in a case study example provide accurate probabilistic statements that correspond to the intended magnitude-based inferences. The model is described in the context of a published small-scale athlete study which employed a magnitude-based inference approach to compare the effect of two altitude training regimens (live high-train low (LHTL), and intermittent hypoxic exposure (IHE)) on running performance and blood measurements of elite triathletes. The posterior distributions, and corresponding point and interval estimates, for the parameters and associated effects and comparisons of interest, were estimated using Markov chain Monte Carlo simulations. The Bayesian analysis was shown to provide more direct probabilistic comparisons of treatments and able to identify small effects of interest. The approach avoided asymptotic assumptions and overcame issues such as multiple testing. Bayesian analysis of unscaled effects showed a probability of 0.96 that LHTL yields a substantially greater increase in hemoglobin mass than IHE, a 0.93 probability of a substantially greater improvement in running economy and a greater than 0.96 probability that both IHE and LHTL yield a substantially greater improvement in maximum blood lactate concentration compared to a Placebo. The conclusions are consistent with those obtained using a ‘magnitude-based inference’ approach that has been promoted in the field. The paper demonstrates that a fully Bayesian analysis is a simple and effective way of analysing small effects, providing a rich set of results that are straightforward to interpret in terms of probabilistic statements.

Avian diversification patterns across the K-Pg boundary: Influence of calibrations, datasets, and model misspecification

Birds represent the most diverse extant tetrapod clade, with ca. 10,000 extant species, and the timing of the crown avian radiation remains hotly debated. The fossil record supports a primarily Cenozoic radiation of crown birds, whereas molecular divergence dating analyses generally imply that this radiation was well underway during the Cretaceous. Furthermore, substantial differences have been noted between published divergence estimates. These have been variously attributed to clock model, calibration regime, and gene type. One underappreciated phenomenon is that disparity between fossil ages and molecular dates tends to be proportionally greater for shallower nodes in the avian Tree of Life. Here, we explore potential drivers of disparity in avian divergence dates through a set of analyses applying various calibration strategies and coding methods to a mitochondrial genome dataset and an 18-gene nuclear dataset, both sampled across 72 taxa. Our analyses support the occurrence of two deep divergences (i.e., the Palaeognathae/Neognathae split and the Galloanserae/Neoaves split) well within the Cretaceous, followed by a rapid radiation of Neoaves near the K-Pg boundary. However, 95% highest posterior density intervals for most basal divergences in Neoaves cross the boundary, and we emphasize that, barring unreasonably strict prior distributions, distinguishing between a rapid Early Paleocene radiation and a Late Cretaceous radiation may be beyond the resolving power of currently favored divergence dating methods. In contrast to recent observations for placental mammals, constraining all divergences within Neoaves to occur in the Cenozoic does not result in unreasonably high inferred substitution rates. Comparisons of nuclear DNA (nDNA) versus mitochondrial DNA (mtDNA) datasets and NT- versus RY-coded mitochondrial data reveal patterns of disparity that are consistent with substitution model misspecifications that result in tree compression/tree extension artifacts, which may explain some discordance between previous divergence estimates based on different sequence types. Comparisons of fully calibrated and nominally calibrated trees support a correlation between body mass and apparent dating error. Overall, our results are consistent with (but do not require) a Paleogene radiation for most major clades of crown birds.