The Node-wise Pseudo-marginal Method

by   Denishrouf Thesingarajah, et al.

Motivated by problems from neuroimaging in which existing approaches make use of "mass univariate" analysis which neglects spatial structure entirely, but the full joint modelling of all quantities of interest is computationally infeasible, a novel method for incorporating spatial dependence within a (potentially large) family of model-selection problems is presented. Spatial dependence is encoded via a Markov random field model for which a variant of the pseudo-marginal Markov chain Monte Carlo algorithm is developed and extended by a further augmentation of the underlying state space. This approach allows the exploitation of existing unbiased marginal likelihood estimators used in settings in which spatial independence is normally assumed thereby facilitating the incorporation of spatial dependence using non-spatial estimates with minimal additional development effort. The proposed algorithm can be realistically used for analysis of subsets of large image moderately sized data sets such as 2D slices of whole 3D dynamic PET brain images or other regions of interest. Principled approximations of the proposed method, together with simple extensions based on the augmented spaces, are investigated and shown to provide similar results to the full pseudo-marginal method. Such approximations and extensions allow the improved performance obtained by incorporating spatial dependence to be obtained at negligible additional cost. An application to measured PET image data shows notable improvements in revealing underlying spatial structure when compared to current methods that assume spatial independence.



page 18

page 20

page 26

page 30


Markov chain Monte Carlo importance samplers for Bayesian models with intractable likelihoods

We consider the efficient use of an approximation within Markov chain Mo...

A subsampling approach for Bayesian model selection

It is common practice to use Laplace approximations to compute marginal ...

A Bayesian semi-parametric hybrid model for spatial extremes with unknown dependence structure

The max-stable process is an asymptotically justified model for spatial ...

Large Sample Asymptotics of the Pseudo-Marginal Method

The pseudo-marginal algorithm is a variant of the Metropolis-Hastings al...

bssm: Bayesian Inference of Non-linear and Non-Gaussian State Space Models in R

We present an R package bssm for Bayesian non-linear/non-Gaussian state ...

Structure Learning of Contextual Markov Networks using Marginal Pseudo-likelihood

Markov networks are popular models for discrete multivariate systems whe...

Modelling Extremes of Spatial Aggregates of Precipitation using Conditional Methods

Inference on the extremal behaviour of spatial aggregates of precipitati...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.