Calibration concordance for astronomical instruments via multiplicative shrinkage
File(s)18-jasa-calib.pdf (1.25 MB)
Working paper
Author(s)
Type
Working Paper
Abstract
Calibration data are often obtained by observing several well-understood
objects simultaneously with multiple instruments, such as satellites for
measuring astronomical sources. Analyzing such data and obtaining proper
concordance among the instruments is challenging when the physical source
models are not well understood, when there are uncertainties in "known"
physical quantities, or when data quality varies in ways that cannot be fully
quantified. Furthermore, the number of model parameters increases with both the
number of instruments and the number of sources. Thus, concordance of the
instruments requires careful modeling of the mean signals, the intrinsic source
differences, and measurement errors. In this paper, we propose a log-Normal
hierarchical model and a more general log-t model that respect the
multiplicative nature of the mean signals via a half-variance adjustment, yet
permit imperfections in the mean modeling to be absorbed by residual variances.
We present analytical solutions in the form of power shrinkage in special cases
and develop reliable MCMC algorithms for general cases. We apply our method to
several data sets obtained with a variety of X-ray telescopes such as Chandra.
We demonstrate that our method provides helpful and practical guidance for
astrophysicists when adjusting for disagreements among instruments.
objects simultaneously with multiple instruments, such as satellites for
measuring astronomical sources. Analyzing such data and obtaining proper
concordance among the instruments is challenging when the physical source
models are not well understood, when there are uncertainties in "known"
physical quantities, or when data quality varies in ways that cannot be fully
quantified. Furthermore, the number of model parameters increases with both the
number of instruments and the number of sources. Thus, concordance of the
instruments requires careful modeling of the mean signals, the intrinsic source
differences, and measurement errors. In this paper, we propose a log-Normal
hierarchical model and a more general log-t model that respect the
multiplicative nature of the mean signals via a half-variance adjustment, yet
permit imperfections in the mean modeling to be absorbed by residual variances.
We present analytical solutions in the form of power shrinkage in special cases
and develop reliable MCMC algorithms for general cases. We apply our method to
several data sets obtained with a variety of X-ray telescopes such as Chandra.
We demonstrate that our method provides helpful and practical guidance for
astrophysicists when adjusting for disagreements among instruments.
Date Issued
2018-09-26
Citation
2018
Publisher
arXiv
Copyright Statement
© 2018 The Author(s).
Identifier
http://arxiv.org/abs/1711.09429v3
Subjects
stat.AP
stat.AP
Publication Status
Published