3,970 research outputs found
Stochastic Algorithm For Parameter Estimation For Dense Deformable Template Mixture Model
Estimating probabilistic deformable template models is a new approach in the
fields of computer vision and probabilistic atlases in computational anatomy. A
first coherent statistical framework modelling the variability as a hidden
random variable has been given by Allassonni\`ere, Amit and Trouv\'e in [1] in
simple and mixture of deformable template models. A consistent stochastic
algorithm has been introduced in [2] to face the problem encountered in [1] for
the convergence of the estimation algorithm for the one component model in the
presence of noise. We propose here to go on in this direction of using some
"SAEM-like" algorithm to approximate the MAP estimator in the general Bayesian
setting of mixture of deformable template model. We also prove the convergence
of this algorithm toward a critical point of the penalised likelihood of the
observations and illustrate this with handwritten digit images
DPO - Denoising, Deconvolving, and Decomposing Photon Observations
The analysis of astronomical images is a non-trivial task. The D3PO algorithm
addresses the inference problem of denoising, deconvolving, and decomposing
photon observations. Its primary goal is the simultaneous but individual
reconstruction of the diffuse and point-like photon flux given a single photon
count image, where the fluxes are superimposed. In order to discriminate
between these morphologically different signal components, a probabilistic
algorithm is derived in the language of information field theory based on a
hierarchical Bayesian parameter model. The signal inference exploits prior
information on the spatial correlation structure of the diffuse component and
the brightness distribution of the spatially uncorrelated point-like sources. A
maximum a posteriori solution and a solution minimizing the Gibbs free energy
of the inference problem using variational Bayesian methods are discussed.
Since the derivation of the solution is not dependent on the underlying
position space, the implementation of the D3PO algorithm uses the NIFTY package
to ensure applicability to various spatial grids and at any resolution. The
fidelity of the algorithm is validated by the analysis of simulated data,
including a realistic high energy photon count image showing a 32 x 32 arcmin^2
observation with a spatial resolution of 0.1 arcmin. In all tests the D3PO
algorithm successfully denoised, deconvolved, and decomposed the data into a
diffuse and a point-like signal estimate for the respective photon flux
components.Comment: 22 pages, 8 figures, 2 tables, accepted by Astronomy & Astrophysics;
refereed version, 1 figure added, results unchanged, software available at
http://www.mpa-garching.mpg.de/ift/d3po
Neural Connectivity with Hidden Gaussian Graphical State-Model
The noninvasive procedures for neural connectivity are under questioning.
Theoretical models sustain that the electromagnetic field registered at
external sensors is elicited by currents at neural space. Nevertheless, what we
observe at the sensor space is a superposition of projected fields, from the
whole gray-matter. This is the reason for a major pitfall of noninvasive
Electrophysiology methods: distorted reconstruction of neural activity and its
connectivity or leakage. It has been proven that current methods produce
incorrect connectomes. Somewhat related to the incorrect connectivity
modelling, they disregard either Systems Theory and Bayesian Information
Theory. We introduce a new formalism that attains for it, Hidden Gaussian
Graphical State-Model (HIGGS). A neural Gaussian Graphical Model (GGM) hidden
by the observation equation of Magneto-encephalographic (MEEG) signals. HIGGS
is equivalent to a frequency domain Linear State Space Model (LSSM) but with
sparse connectivity prior. The mathematical contribution here is the theory for
high-dimensional and frequency-domain HIGGS solvers. We demonstrate that HIGGS
can attenuate the leakage effect in the most critical case: the distortion EEG
signal due to head volume conduction heterogeneities. Its application in EEG is
illustrated with retrieved connectivity patterns from human Steady State Visual
Evoked Potentials (SSVEP). We provide for the first time confirmatory evidence
for noninvasive procedures of neural connectivity: concurrent EEG and
Electrocorticography (ECoG) recordings on monkey. Open source packages are
freely available online, to reproduce the results presented in this paper and
to analyze external MEEG databases
Hierarchical Cosmic Shear Power Spectrum Inference
We develop a Bayesian hierarchical modelling approach for cosmic shear power
spectrum inference, jointly sampling from the posterior distribution of the
cosmic shear field and its (tomographic) power spectra. Inference of the shear
power spectrum is a powerful intermediate product for a cosmic shear analysis,
since it requires very few model assumptions and can be used to perform
inference on a wide range of cosmological models \emph{a posteriori} without
loss of information. We show that joint posterior for the shear map and power
spectrum can be sampled effectively by Gibbs sampling, iteratively drawing
samples from the map and power spectrum, each conditional on the other. This
approach neatly circumvents difficulties associated with complicated survey
geometry and masks that plague frequentist power spectrum estimators, since the
power spectrum inference provides prior information about the field in masked
regions at every sampling step. We demonstrate this approach for inference of
tomographic shear -mode, -mode and -cross power spectra from a
simulated galaxy shear catalogue with a number of important features; galaxies
distributed on the sky and in redshift with photometric redshift uncertainties,
realistic random ellipticity noise for every galaxy and a complicated survey
mask. The obtained posterior distributions for the tomographic power spectrum
coefficients recover the underlying simulated power spectra for both - and
-modes.Comment: 16 pages, 8 figures, accepted by MNRA
- …