Detecting Unspecified Structure in Low-Count Images
File(s)lira_testing_rev.pdf (1.17 MB) pdf (1).pdf (1.23 MB)
Accepted version
Published version
OA Location
Author(s)
Stein, NM
van Dyk, DA
Kashyap, VL
Siemiginowska, A
Type
Journal Article
Abstract
Unexpected structure in images of astronomical sources often presents itself upon
visual inspection of the image, but such apparent structure may either correspond to
true features in the source or be due to noise in the data. This paper presents a
method for testing whether inferred structure in an image with Poisson noise represents a
significant departure from a baseline (null) model of the image. To infer image structure,
we conduct a Bayesian analysis of a full model that uses a multiscale component to
allow flexible departures from the posited null model. As a test statistic, we use a
tail probability of the posterior distribution under the full model. This choice of test
statistic allows us to estimate a computationally efficient upper bound on a p-value
that enables us to draw strong conclusions even when there are limited computational
resources that can be devoted to simulations under the null model. We demonstrate
the statistical performance of our method on simulated images. Applying our method
to an X-ray image of the quasar 0730+257, we find significant evidence against the null
model of a single point source and uniform background, lending support to the claim of
an X-ray jet.
visual inspection of the image, but such apparent structure may either correspond to
true features in the source or be due to noise in the data. This paper presents a
method for testing whether inferred structure in an image with Poisson noise represents a
significant departure from a baseline (null) model of the image. To infer image structure,
we conduct a Bayesian analysis of a full model that uses a multiscale component to
allow flexible departures from the posited null model. As a test statistic, we use a
tail probability of the posterior distribution under the full model. This choice of test
statistic allows us to estimate a computationally efficient upper bound on a p-value
that enables us to draw strong conclusions even when there are limited computational
resources that can be devoted to simulations under the null model. We demonstrate
the statistical performance of our method on simulated images. Applying our method
to an X-ray image of the quasar 0730+257, we find significant evidence against the null
model of a single point source and uniform background, lending support to the claim of
an X-ray jet.
Date Issued
2015-11-01
Date Acceptance
2015-08-31
Citation
Astrophysical Journal, 2015, 813 (1)
ISSN
1538-4357
Publisher
American Astronomical Society
Journal / Book Title
Astrophysical Journal
Volume
813
Issue
1
Copyright Statement
© 2015. The American Astronomical Society. All rights reserved.
Sponsor
The Royal Society
Commission of the European Communities
Science and Technology Facilities Council (STFC)
National Science Foundation (US)
Grant Number
WM110023
FP7-PEOPLE-2012-CIG-321865
ST/K001051/1
DMS 15-13484
Subjects
Science & Technology
Physical Sciences
Astronomy & Astrophysics
galaxies: jets
methods: data analysis
methods: statistical
quasars: individual (0730+257)
techniques: image processing
X-rays: general
HIERARCHICAL-MODELS
RECONSTRUCTION
RESTORATION
STATISTICS
CHECKING
Publication Status
Published
Article Number
66