Objective Bayesian analysis for the normal compositional model

The issue of objective prior specification for the parameters in the normal compositional model is considered within the context of statistical analysis of linearly mixed structures in image processing. In particular, the Jeffreys prior for the vector of fractional abundances in case of a known covariance matrix is derived. If an additional unknown variance parameter is present, the Jeffreys prior and the reference prior are computed and it is proven that the resulting posterior distributions are proper. Markov chain Monte Carlo strategies are proposed to efficiently sample from the posterior distributions and the priors are compared on the grounds of the frequentist properties of the resulting Bayesian inferences. The default Bayesian analysis is illustrated by a dataset taken from fluorescence spectroscopy.

[1]  D. Stein,et al.  Application of the normal compositional model to the analysis of hyperspectral imagery , 2003, IEEE Workshop on Advances in Techniques for Analysis of Remotely Sensed Data, 2003.

[2]  H. Jeffreys,et al.  Theory of probability , 1896 .

[3]  Yincai Tang,et al.  Objective Bayesian analysis of accelerated competing failure models under Type-I censoring , 2011, Comput. Stat. Data Anal..

[4]  Rui Paulo Default priors for Gaussian processes , 2005 .

[5]  Michael T. Eismann,et al.  Stochastic Mixture Modeling , 2006 .

[6]  P. Diggle,et al.  Model‐based geostatistics , 2007 .

[7]  Roger N. Clark,et al.  The US Geological Survey, digital spectral reflectance library: version 1: 0.2 to 3.0 microns , 1993 .

[8]  Chein-I Chang,et al.  Semi-Supervised Linear Spectral Unmixing Using a Hierarchical Bayesian Model for Hyperspectral Imagery , 2008, IEEE Transactions on Signal Processing.

[9]  John F. Mustard,et al.  Spectral unmixing , 2002, IEEE Signal Process. Mag..

[10]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[11]  E. L. Lehmann,et al.  Theory of point estimation , 1950 .

[12]  A. Fiacco A Finite Algorithm for Finding the Projection of a Point onto the Canonical Simplex of R " , 2009 .

[13]  James O. Berger,et al.  Objective Bayesian Analysis for the Multivariate Normal Model , 2006 .

[14]  L. Wasserman,et al.  The Selection of Prior Distributions by Formal Rules , 1996 .

[15]  Jean-Yves Tourneret,et al.  Bayesian Estimation of Linear Mixtures Using the Normal Compositional Model. Application to Hyperspectral Imagery , 2010, IEEE Transactions on Image Processing.

[16]  J. Bernardo Reference Posterior Distributions for Bayesian Inference , 1979 .

[17]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[18]  Alfonso Fernández-Manso,et al.  Spectral unmixing , 2012 .

[19]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[20]  J. Berger,et al.  Objective Bayesian Analysis of Spatially Correlated Data , 2001 .

[21]  José M. Bioucas-Dias,et al.  Vertex component analysis: a fast algorithm to unmix hyperspectral data , 2005, IEEE Transactions on Geoscience and Remote Sensing.

[22]  Chein-I. Chang Hyperspectral Data Exploitation: Theory and Applications , 2007 .

[23]  Hannes Kazianka,et al.  A Bayesian approach to estimating linear mixtures with unknown covariance structure , 2011 .

[24]  M. Aitkin Posterior Bayes Factors , 1991 .

[25]  E. Ziegel,et al.  Geostatistics Wollongong '96 , 1997 .