Bayesian data analysis
The goal of Bayesian [1] [2] or integrated data analysis (IDA) is to combine the information from a set of diagnostics providing complementary information in order to recover the best possible reconstruction of the actual state of the system subjected to measurement. [3] [4] [5] [6] [7] [8] Like Function parametrization (FP), this technique requires having a forward model to predict the measurement readings for any given state of the physical system; however
- instead of computing an estimate of the inverse of the forward model (as with FP), IDA finds the best model state corresponding to a specific measurement by a maximization procedure (maximization of the likelihood);
- the handling of error propagation is more sophisticated within IDA, allowing non-Gaussian error distributions and absolutely general parameter interdependencies; and
- additionally, it provides a systematic way to include prior knowledge into the analysis.
Bayes' Theorem
The method is based on Bayes' Theorem, expressed as follows:
Here, α are a set of model parameters, and P is the probability distribution of these parameters, given the experimental data d, their errors σ, and additional information I.
Bayes' Theorem expresses this probability distribution as a product of the likelihood L of obtaining the cited experimental data given some values of the model parameters α as well as σ and I, and the prior distribution π that expresses the knowledge concerning the model parameters preceding any measurement.
The likelihood L is computed using a forward model of the experiment, returning the value of simulated measurements while assuming a given physical state of the experimental system. It should be noted that this forward model (from system parameters to measurements) is often much easier to compute than the reverse mapping (from measurements to system parameters), as the latter is often the inverse of a projection, which is therefore typically ill-determined.
The normalization of the equation serves to maintain its character of a probability distribution, although it is not important for the determination of the best values of the parameters and their errors. The optimum reconstruction is determined by maximizing the posterior P, varying the parameters α.
See also
References
- ↑ D.S. Sivia, Data Analysis: A Bayesian Tutorial, Oxford University Press, USA (1996) ISBN 0198518897
- ↑ P. Gregory, Bayesian Logical Data Analysis for the Physical Sciences, Cambridge University Press, Cambridge (2005) ISBN 052184150X
- ↑ R. Fischer, C. Wendland, A. Dinklage, et al, Thomson scattering analysis with the Bayesian probability theory, Plasma Phys. Control. Fusion 44 (2002) 1501
- ↑ R. Fischer, A. Dinklage, and E. Pasch, Bayesian modelling of fusion diagnostics, Plasma Phys. Control. Fusion 45 (2003) 1095-1111
- ↑ R. Fischer, A. Dinklage, Integrated data analysis of fusion diagnostics by means of the Bayesian probability theory, Rev. Sci. Instrum. 75 (2004) 4237
- ↑ A. Dinklage, R. Fischer, and J. Svensson, Topics and Methods for Data Validation by Means of Bayesian Probability Theory, Fusion Sci. Technol. 46 (2004) 355
- ↑ J. Svensson, A. Werner, Large Scale Bayesian Data Analysis for Nuclear Fusion Experiments, IEEE International Symposium on Intelligent Signal Processing (2007) 1
- ↑ R. Fischer, C.J. Fuchs, B. Kurzan, et al., Integrated Data Analysis of Profile Diagnostics at ASDEX Upgrade, Fusion Sci. Technol. 58 (2010) 675