Semi-supervised linear spectral unmixing using a hierarchical Bayesian model for hyperspectral imagery

This paper proposes a hierarchical Bayesian model that can be used for semi-supervised hyperspectral image unmixing. The model assumes that the pixel reflectances result from linear combinations of pure component spectra contaminated by an additive Gaussian noise. The abundance parameters appearing in this model satisfy positivity and additivity constraints. These constraints are naturally expressed in a Bayesian context by using appropriate abundance prior distributions. The posterior distributions of the unknown model parameters are then derived. A Gibbs sampler allows one to draw samples distributed according to the posteriors of interest and to estimate the unknown abundances. An extension of the algorithm is finally studied for mixtures with unknown numbers of spectral components belonging to a know library. The performance of the different unmixing strategies is evaluated via simulations conducted on synthetic and real data

Development and characterisation of a novel three-dimensional inter-kingdom wound biofilm model

Chronic diabetic foot ulcers are frequently colonised and infected by polymicrobial biofilms that ultimately prevent healing. This study aimed to create a novel in vitro inter-kingdom wound biofilm model on complex hydrogel-based cellulose substrata to test commonly used topical wound treatments. Inter-kingdom triadic biofilms composed of Candida albicans, Pseudomonas aeruginosa, and Staphylococcus aureus were shown to be quantitatively greater in this model compared to a simple substratum when assessed by conventional culture, metabolic dye and live dead qPCR. These biofilms were both structurally complex and compositionally dynamic in response to topical therapy, so when treated with either chlorhexidine or povidone iodine, principal component analysis revealed that the 3-D cellulose model was minimally impacted compared to the simple substratum model. This study highlights the importance of biofilm substratum and inclusion of relevant polymicrobial and inter-kingdom components, as these impact penetration and efficacy of topical antiseptics

WSRT Faraday tomography of the Galactic ISM at \lambda \sim 0.86 m

We investigate the distribution and properties of Faraday rotating and synchrotron emitting regions in the Galactic ISM in the direction of the Galactic anti-centre. We apply Faraday tomography to a radio polarization dataset that we obtained with the WSRT. We developed a new method to calculate a linear fit to periodic data, which we use to determine rotation measures from our polarization angle data. From simulations of a Faraday screen + noise we could determine how compatible the data are with Faraday screens. An unexpectedly large fraction of 14% of the lines-of-sight in our dataset show an unresolved main component in the Faraday depth spectrum. For lines-of-sight with a single unresolved component we demonstrate that a Faraday screen in front of a synchrotron emitting region that contains a turbulent magnetic field component can explain the data.Comment: 5 pages, 5 figures. Accepted for publication as a Letter to the Editor in A&

Unsupervised Bayesian linear unmixing of gene expression microarrays

Background: This paper introduces a new constrained model and the corresponding algorithm, called unsupervised Bayesian linear unmixing (uBLU), to identify biological signatures from high dimensional assays like gene expression microarrays. The basis for uBLU is a Bayesian model for the data samples which are represented as an additive mixture of random positive gene signatures, called factors, with random positive mixing coefficients, called factor scores, that specify the relative contribution of each signature to a specific sample. The particularity of the proposed method is that uBLU constrains the factor loadings to be non-negative and the factor scores to be probability distributions over the factors. Furthermore, it also provides estimates of the number of factors. A Gibbs sampling strategy is adopted here to generate random samples according to the posterior distribution of the factors, factor scores, and number of factors. These samples are then used to estimate all the unknown parameters. Results: Firstly, the proposed uBLU method is applied to several simulated datasets with known ground truth and compared with previous factor decomposition methods, such as principal component analysis (PCA), non negative matrix factorization (NMF), Bayesian factor regression modeling (BFRM), and the gradient-based algorithm for general matrix factorization (GB-GMF). Secondly, we illustrate the application of uBLU on a real time-evolving gene expression dataset from a recent viral challenge study in which individuals have been inoculated with influenza A/H3N2/Wisconsin. We show that the uBLU method significantly outperforms the other methods on the simulated and real data sets considered here. Conclusions: The results obtained on synthetic and real data illustrate the accuracy of the proposed uBLU method when compared to other factor decomposition methods from the literature (PCA, NMF, BFRM, and GB-GMF). The uBLU method identifies an inflammatory component closely associated with clinical symptom scores collected during the study. Using a constrained model allows recovery of all the inflammatory genes in a single factor