## The Proton Radius from Bayesian Inference

Krzysztof M. Graczyk ∗ and Cezary Juszczak

Institute for Theoretical Physics, University of Wroc/suppress law, pl. M. Borna 9, 50-204, Wroc/suppress law, Poland

The methods of Bayesian statistics are used to extract the value of the proton radius from the elastic ep scattering data in a model independent way. To achieve that goal a large number of parametrizations (equivalent to neural network schemes) are considered and ranked by their conditional probability P (parametrization | data) instead of using the minimal error criterion. As a result the most probable proton radii values ( r p E = 0 899 . ± 0 003 fm, . r p M = 0 879 . ± 0 007 fm) are obtained . and systematic error due to freedom in the choice of parametrization is estimated. Correcting the data for the two photon exchange effect leads to smaller difference between the extracted values of r p E and r p M . The results disagree with recent muonic atom measurements.

PACS numbers: 13.40.Gp, 25.30.Bf, 14.20.Dh, 84.35.+i

Keywords: proton radius, proton form-factors, two-photon exchange correction

The problem of the proton radius is being discussed in the particle and atomic physics for many years but recently it has received even more attention due to the new results coming form the Lamb shift measurements of the muonic hydrogen atom ( µp ) [1]. It is a very accurate estimate of the proton radius but it is inconsistent with the CODATA value (compilation of the measurements in hydrogen atom and analysis of the electron-proton scattering data) [2]. This disagreement, called 'the proton radius puzzle', has not been explained yet. There are many proposals to resolve the problem including considering different types of interaction for ep and µp [3] or 'weakening the assumptions of perturbative formulation of quantum electrodynamics within the proton' [4]. In this paper we introduce a novel approach to the extraction of the proton radius from the ep scattering data, which allows one to control the model dependence of the result within a Bayesian objective algorithm.

The electromagnetic (E-M) structure of the proton is encoded in the electric G E and magnetic G M form factors, which are also the crucial input for the atomic (hydrogen-like) calculations [5]. In the low Q 2 limit in the Breit frame they are related with the distributions of the electric charge and magnetic momentum inside the nucleon [6]. In particular the nucleon radius is expressed by the slope of the form factors at vanishing Q 2 :

$$r _ { E, M } ^ { p } \equiv \left ( - \frac { 6 } { G _ { E, M } ( 0 ) } \, \frac { d G _ { E, M } ( Q ^ { 2 } ) } { d Q ^ { 2 } } \Big | _ { Q ^ { 2 } = 0 } ^ { \frac { 1 } { 2 } } \cdot \quad ( 1 ) \quad \text{direct} \\ \text{The value of $r_{E,M}^{p}$ can be obtained from the scattering} \\ \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{or} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \,\text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is}\, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{is} \, \text{\ }$$

∣ The value of r p E,M can be obtained from the scattering data (mainly elastic ep ). Recent results include [7-11] and a more complete list can be found in Refs. [3, 12, 13]. These analyses have been performed in the spirit of the frequentistic statistics.

The values of the proton radius (related to the electric charge distribution) obtained by different authors during

∗ Electronic address: krzysztof.graczyk@ift.uni.wroc.pl

FIG. 1: (Color online) The scheme of the neural network (with H = 4 hidden units) used in the analysis, for reference see Eq. 6. Every line (connection) corresponds to one parameter w i . Gray circles denote the sigmoid functions, open circles are constants.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0150v3-figure-4beb615a-6e22-4b72-ba2f-019b2c557a3a.jpg)

the last fifty years, range from about 0.8 fm to about 0.9 fm [14]. Indeed, as it was demonstrated in [8, 9] the results depend the choice of the form factor parametrization but it is obvious that the collection of datasets used in the analysis also matters. Even the right choice of the number of parameters in a given class of parametrizations can be a challenge. More parameters in the fit can always reduce the χ 2 error but at some point the fit tends to reproduce statical fluctuations of the experimental measurements. Such a model overfits the data i.e. describes existing points with unrealistic precision but has no predictive power - introducing new data points drastically increases the χ 2 error and spoils the quality of the fit. This is connected with the so called bias variance tradeoff. Some attempt to resolve this problem has been made in Ref. [9] where conformal mapping was used in order to exploit the analytic properties of the form factors.

Our philosophy is different. In contrast to the frequentistic statistics methods used in the above papers, we use the Bayesian statistics which is well suited to model comparison. Within this framework, adapted for neural networks, we select the most optimal model (according to the data) and extract the value of the proton radius.

The extraction of the nucleon form factors in [15] was our first use of the Bayesian neural network framework. In the next paper [16] the approach was significantly im-

FIG. 2: (Color online) Logarithm of evidence for the best neural networks in each scheme. Solid/dashed lines correspond to the analyses of the data corrected/not corrected by the TPE effect. Points mark the best models for each analysis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0150v3-figure-98e29052-a871-4af3-8749-a86d4f22abf1.jpg)

proved to obtain the two-photon exchange (TPE) correction to unpolarized elastic ep cross section. However, because of some physical assumptions these analyses were not dedicated to the lowQ 2 region. In the present paper we overcome this difficulty because: (i) the TPE contribution (elastic and inelastic described by box diagrams with nucleon and ∆(1232) resonance as intermediate hadronic states) calculated in [17] is subtracted from the data and only E-M form factors are extracted from the cross section and polarization transfer data; (ii) the form factors are re-parameterized in order to automatically fulfil G E (0) = G M (0) /µ p = 1 ( µ p is the proton magnetic moment); (iii) our software has been significantly improved in terms of accuracy and efficiency.

Here we briefly outline the Bayesian framework, for more details see [15-17]. The main idea is to consider a large class of models. By a model we mean the function N used to fit the data D and two conditional probabilities: P ( { w i }|N ) - prior which accommodates the initial assumption about the { w i } , and the likelihood P ( D|{ w i } N , ); { w i } is the set of parameters of the function N . Then the models are ranked by the conditional probability P ( N|D ), which can be replaced by the evidence P ( D|N ) if the data is the same and no model is preferred at the beginning of the analysis.

For any given model we denote with { w i } MP the configuration of the parameters which maximizes the posterior

$$P ( \{ w _ { i } \} | \mathcal { D }, \mathcal { N } ) = \frac { P ( \mathcal { D } | \{ w _ { i } \}, \mathcal { N } ) P ( \{ w _ { i } \} | \mathcal { N } ) } { P ( \mathcal { D } | \mathcal { N } ) }. \quad ( 2 ) \quad \text{by o} \\ \sigma G _ { \varepsilon } ^ { 2 }.$$

Finding { w i } MP in one of the steps of the analysis.

Usually the evidence is peaked at { w i } MP and it simplifies to the likelihood at the maximum multiplied by the Occam factor which makes too complex models less

FIG. 3: (Color online) The form factor G /G E D ( G D = 1 (1 + / Q / . 2 0 71GeV ) ). 2 2 Red, blue and magenta lines correspond to the best models in the analysis with data limited by Q 2 cutoff = 1, 3 and 10 GeV 2 respectively. The areas colored with magenta, blue and red denote 1 σ uncertainty due to the change in the fit parameters (calculated within the Hessian approximation). The gray lines show models which are best for neural networks with definite number of hidden units, but are not globally best.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0150v3-figure-83c780a2-b528-4523-a0b3-9df794b3a803.jpg)

likely [18]

$$\begin{array} { c } \Phi _ { \text{$\mathcal{ } } } & \cdot & \cdot \\ \text{$P(\mathcal{ }D|\mathcal{N})$} \approx P(\mathcal{ }D|\{w_{i}\}_{M P},\mathcal{N})$} \underbrace { ( 2 \pi ) ^ { \frac { W } { 2 } } ( \det A ) ^ { - \frac { 1 } { 2 } } } _ { O c c a m \, f a c t o r }, \\ \intertext { \text{te} } & \text{where $W$ is the number of parameters $\{w_{i}$} \} \text{ and } A _ { i j } = \\ & - \nabla \dots \nabla \dots \, \ln P ( \i u u \, \downarrow | \mathcal{ }D \, \mathcal{ }M \, | } \end{array}$$

∣ For the function N we take the feed-forward neural networks with one hidden layer of units (Fig. 1), which correspond to linear combinations of sigmoid functions. This class of functions can be used to approximate any continuous function with arbitrary precision [19].

where W is the number of parameters { w i } and A ij = - ∇ w i ∇ w j ln P ( { w i }|D N , ) ∣ { w i } = { w i } MP .

The prior function for the neural network parameters { w i } has the standard form P { ( w i }| α, N ) ∼ e -αE w where E w = 1 2 ∑ W i =1 w 2 i . The parameter α is the regularizer which determines the width of the initial Gaussian distribution for the parameters. In principle it should be treated as another parameter of the model and its optimal value α MP is also obtained during the analysis. The likelihood is defined by the χ 2 distribution

$$\mathbb { N } _ { e _ { - } } \quad \mathcal { P } \left ( \mathcal { D } | \left \{ w _ { i } \right \}, \alpha, \mathcal { N } \right ) \sim e ^ { - \chi _ { c r. } ^ { 2 } \left ( \mathcal { D }, \{ w _ { i } \} \right ) - \chi _ { P T } ^ { 2 } \left ( \mathcal { D }, \{ w _ { i } \} \right ) }. \quad ( 4 )$$

The χ 2 cr. includes the unpolarized elastic ep scattering cross section data, which are assumed to be parametrized by one-photon exchange formula σ R = Q G 2 2 M / 4 M 2 p + εG 2 E ( M p is the proton mass, ε denotes photon polarizability). Each independent set of cross section measurements is characterized by its systematic normalization uncertainty so we introduce a separate normalization parameter for every data set. Values of these parameters are

obtained by an iterative algorithm described in [16]. The χ 2 PT includes the form factor ratio µ G p E /G M data from the polarization transfer measurements. We consider the same selection of data as in [16] with one additional data set from [20] and make the prior assumption that all data sets are equally relevant. The MAMI data [7] are not included in the analysis because of the difficulties described in [21].

The optimal configuration { α, w i } MP is found by an iterative procedure based on the fact that the maximum of the posterior corresponds to the minimum over { w i } of the error function S ( D { , w i } ) ≡ χ 2 cr. ( D { , w i } ) + χ 2 PT ( D { , w i } ) + α MP E w and also ∂ ∂α P D| ( α, N ) ∣ α = α MP = 0, where

$$\begin{array} { c } \chi _ { c r. } ( \mathcal { V }, \underline { \mathcal { W } } _ { i J } ) \rightarrow \chi _ { P T } ( \mathcal { V }, \underline { \mathcal { W } } _ { i J } ) \rightarrow \alpha _ { M P } \mathcal { L } _ { w } \quad \text{and} \quad \text{as} \quad \cdots \rightarrow \\ \frac { \partial } { \partial \alpha } \mathcal { P } \left ( \mathcal { D } \right | \alpha, \mathcal { N } \right ) \Big | _ { \alpha = \alpha _ { M P } } = 0, \text{ where} & & \text{$T$} \right | \\ & & \text{our} \\ \mathcal { P } \left ( \mathcal { D } \right | \alpha, \mathcal { N } \right ) = \int d w \mathcal { P } \left ( \mathcal { D } \right | \left \{ w _ { i } \right \}, \alpha, \mathcal { N } \right ) \mathcal { P } \left ( \left \{ w _ { i } \right \} \right | \alpha, \mathcal { N } \right ). \quad \text{num} \\ \text{for} \quad & \text{dep} \\ \text{Eventually the algorithm of evidence for the model } \mathcal { M } & \text{$\Delta$} \quad \cdots \right. \end{array}$$

Eventually the logarithm of evidence for the model N (in the Hessian approximation) reads:

$$\ln P ( \mathcal { D } | \mathcal { N } ) \ \approx \ - S ( \mathcal { D }, \{ w _ { i } \} _ { M P } ) + \frac { W } { 2 } \ln \alpha _ { M P } - \frac { 1 } { 2 } \ln \det A \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\mathcal { N } \right ) \ & \approx \ - S ( \mathcal { D }, \{ w _ { i } \} _ { M P } ) + \frac { W } { 2 } \ln \alpha _ { M P } - \frac { 1 } { 2 } \ln \det A \quad \text{cons} \\ & \quad - \frac { 1 } { 2 } \ln \sum _ { i = 1 } ^ { W } \frac { 1 } { 2 } \frac { \lambda _ { i } } { \lambda _ { i } + \alpha _ { M P } } \quad ( 5 ) \quad \text{izati} \\ A. \, \cdot \ = \ \nabla \quad \nabla \quad \mathcal { S } ( \mathcal { D } \ \mathcal { I } w. \, \mathfrak { U } \ \text{and} \ \lambda. \, \mathfrak { s } \ \text{ are the} \quad \text{This}$$

where A kj = ∇ w k ∇ w j S ( D { , w i } ) and λ k -s are the eigenvalues of the matrix ∇ w i ∇ w j ( χ cross ( D { , w i } ) + χ PT ( D { , w i } )).

As the form factors G E , G M describe the properties of the same object, they must be related, so we parametrize them with a single feed forward neural network with one input ( Q 2 ) and two outputs, o E and o M (see Fig. 1),

$$G _ { E } = ( 1 - Q ^ { 2 } o _ { E } ) G _ { D }, \ \ G _ { M } = \mu _ { p } ( 1 - Q ^ { 2 } o _ { M } ) G _ { D }, \ ( 6 ) \quad \ \ln ^ { \mu }$$

where G D ( Q 2 ) = ( 1 + Q /M 2 2 V ) -2 and M 2 V = 0 71 GeV . . 2 Then the value of the proton radius is given by

$$\left ( r _ { E, M } ^ { p } \right ) ^ { 2 } = \left ( r _ { d i p o l e } ^ { p } \right ) ^ { 2 } + 6 o _ { E, M } ( Q ^ { 2 } = 0 ), \quad ( 7 ) \quad \text{also} \\ \text{where} \, ( r ^ { p } \quad \rangle ^ { 2 } = 1 2 / M ^ { 2 } \colon \text{We see that the network out} \quad \text{$\tan$}$$

where ( r p dipole ) 2 = 12 /M 2 V . We see that the network outputs at Q 2 = 0 directly contribute to the deviation of the proton radius from the dipole form factor value.

In Refs. [14, 22] it was shown that modifying the cross section data by the Coulomb distortion (CD) is important in the extraction of the proton radius. On the other hand, the CD is a part of the TPE effect 1 , which has to be subtracted from the data to obtain the form factors which agree with the polarization transfer measurements

1 The relevance of the TPE effect and higher order Born corrections in the proton radius extraction is also discussed in Refs. [23] and [24, 25] respectively.

TABLE I: The values of the proton radii obtained with 1 σ uncertainty due to variation in the parameter space. H is the number of hidden units of the best model.

|   Q 2 cutoff (GeV | r p M (fm)        | r p E (fm)        |   H |
|-------------------|-------------------|-------------------|-----|
|                 1 | 0 . 879 ± 0 . 007 | 0 . 899 ± 0 . 003 |   1 |
|                 3 | 0 . 883 ± 0 . 007 | 0 . 899 ± 0 . 003 |   4 |
|                10 | 0 . 953 ± 0 . 065 | 0 . 897 ± 0 . 005 |   6 |

[26]. Similarly as in [27] we correct the cross section data by the TPE contribution (calculated in [17] in similar way as in [28, 29]).

The set of all possible parametrizations is infinite. In our approach the models are grouped according to the number of units in the hidden layer. For each group we found the { w , α i } which maximize the evidence. The dependence of the maximal evidence on the number of hidden units is plotted in Fig. 2. For every analysis the evidence reaches a peak (which defines the best model) and then starts to fall. Because of this fact we do not consider the networks with more then 40 hidden units. The total number of obtained neural network parameterizations exceeded half million.

It needs to be underlined that the models which maximize the evidence are not those which minimize the χ 2 . This is true in each class of models and also globally. The χ 2 min decreases with the number of parameters possibly saturating at some point. But the models with the lowest χ 2 tend to overfit the data, which is one of the reasons why the χ 2 min is not a suitable criterion for ranking the models. The other problems of χ 2 -based criteria are discussed in [3]. Objective mathematical methods for model comparison are provided by the Bayesian statistics.

In Fig. 3 it is shown how the form factor plots depend on the choice of the parametrization. The gray lines correspond to the best models in each group, while the globally best for each analysis are shown with color. Certainly the choice of the form factor parametrization has also strong impact on the obtained proton radius value (Fig. 4). Hence it is crucial to have criterion for finding the model which is the most favorable by the measurements. On the other hand each model can be true with the probability P ( N|D ) (proportional to the evidence) so it is possible to calculate the expected value of the proton radius and the systematic uncertainty due to the choice of the model. In practice the expected value is very close to the most probable value (see Tables I and II).

In Fig. 4 we show the proton radii values obtained from the analyses where data points with Q 2 above Q 2 cutoff = 1, 3 and 10 GeV 2 are rejected (for the data uncorrected by the TPE we considered only Q 2 cutoff = 1 , 3 GeV ). 2 The choice of the value for the cutoff has small impact on the extraction of the proton radii. For the result of our analysis we take those from the lowest cut, namely r p E = 0 899 . ± 0 003 fm and . r p M = 0 879 . ± 0 007 fm. .

TABLE II: The expected value (according to evidence probability distribution) of the proton radii with systematic uncertainty due to the choice of the parametrization (given by the variance).

|   Q 2 cutoff (GeV | r p M (fm)          | r p E (fm)          |
|-------------------|---------------------|---------------------|
|                 1 | 0 . 8792 ± 0 . 0006 | 0 . 8989 ± 0 . 0001 |
|                 3 | 0 . 8828 ± 0 . 0063 | 0 . 8988 ± 0 . 0003 |
|                10 | 0 . 9205 ± 0 . 0606 | 0 . 8968 ± 0 . 0029 |

On the other hand the r p E 's obtained based on the minimal error criterion depend on the cutoff choice. These results are smaller than those indicated by the evidence and have larger uncertainties, see Fig. 4. We leave for future studies within the Bayesian framework the quantitative investigation of the dependence of the results on the data set selection and Q 2 cutoff value. Eventually it can be seen in Fig. 4 that using the data uncorrected by TPE effect leads to larger uncertainties of r p E and r p M and larger difference r p E -r p M (0 02 fm with TPE and 0 09 fm . . without).

We would like to emphasize that according to our knowledge the present work is the first extraction of the proton radius based on the Bayesian methods. The extracted value of the proton radius r p E agrees with previous results [9, 14, 27] but disagrees with the muonic atom measurements [1]. The subtracting of the TPE correction from the cross section data plays an important role in the analysis.

The calculations have been carried out in Wroclaw Centre for Networking and Supercomputing ( http://www.wcss.wroc.pl ), grant No. 268.

- [1] R. Pohl, A. Antognini, F. Nez, F. D. Amaro, F. Biraben, et al. , Nature 466 , 213 (2010).
- [2] P. J. Mohr, B. N. Taylor, and D. B. Newell, Rev.Mod.Phys. 84 , 1527 (2012), arXiv:1203.5425 [physics.atom-ph] .
- [3] E. Kraus, K. Mesick, A. White, R. Gilman, and S. Strauch, (2014), arXiv:1405.4735 [nucl-ex] .
- [4] K. Pachucki and K. A. Meissner, (2014), arXiv:1405.6582 [hep-ph] .
- [5] K. Pachucki, Phys.Rev. A53 , 2092 (1996).
- [6] F. Ernst, R. Sachs, and K. Wali, Phys.Rev. 119 , 1105 (1960).
- [7] J. Bernauer et al. (A1 Collaboration), Phys.Rev.Lett. 105 , 242001 (2010), arXiv:1007.5076 [nucl-ex] .
- [8] J. Bernauer et al. (A1 Collaboration), (2013), 10.1103/PhysRevC.90.015206, arXiv:1307.6227 [nucl-ex] .
- [9] R. J. Hill and G. Paz, Phys.Rev. D82 , 113005 (2010), arXiv:1008.4619 [hep-ph] .
- [10] I. Lorenz, H.-W. Hammer, and U.-G. Meissner, Eur.Phys.J. A48 , 151 (2012), arXiv:1205.6628 [hep-ph] .
- [11] I. Lorenz and U. G. Meissner, (2014), arXiv:1406.2962 [hep-ph] .
- [12] Z. Epstein, G. Paz, and J. Roy, (2014), arXiv:1407.5683 [hep-ph] .
- [13] S. G. Karshenboim, (2014), arXiv:1405.6515 [hep-ph] .
- [14] I. Sick, Phys.Lett. B576 , 62 (2003), arXiv:nucl-ex/0310008 [nucl-ex] .
- [15] K. M. Graczyk, P. Plonski, and R. Sulej, JHEP 1009 , 053 (2010), arXiv:1006.0342 [hep-ph]

FIG. 4: (Color online) The proton radii corresponding to the models which are the best within particular data selection and neural network scheme. The results for the data corrected by the TPE are marked with black /triangle , /square , and © for Q 2 cutoff = 1, 3 and 10 GeV 2 respectively. The results for the uncorrected data are marked with blue ∇ and ♦ for Q 2 cutoff = 1 and 3 GeV 2 respectively. The filled points mark the best result for each case. The three leftmost red points are the best results according to the minimum of the error function. For clarity their x-coordinate has been changed. Their ln(evidence) values are -718, -1070, and -1311. The picture is based on the analysis of 39 · 5 = 195 models.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0150v3-figure-b7d1c667-893e-4211-9e48-3a77cd107391.jpg)

.

- [16] K. M. Graczyk, Phys.Rev. C84 , 034314 (2011), arXiv:1106.1204 [hep-ph] .
- [17] K. M. Graczyk, Phys.Rev. C88 , 065205 (2013), arXiv:1306.5991 [hep-ph] .
- [18] D. MacKay, Bayesian Methods for Adaptive Models , Ph.D. thesis, California Institute of Technology (1991).
- [19] G. Cybenko, Math Control, Signal 2 , 303 (1989).

- [20] X. Zhan, K. Allada, D. Armstrong, J. Arrington, W. Bertozzi, et al. , Phys.Lett. B705 , 59 (2011), arXiv:1102.0318 [nucl-ex] .
- [21] J. Arrington, Phys.Rev.Lett. 107 , 119101 (2011), arXiv:1108.3058 [nucl-ex] .
- [22] R. Rosenfelder, Phys.Lett. B479 , 381 (2000), arXiv:nucl-th/9912031 [nucl-th] .
- [23] M. Gorchtein, (2014), arXiv:1406.1612 [nucl-th] .
- [24] J. Arrington, J.Phys. G40 , 115003 (2013), arXiv:1210.2677 [nucl-ex] .
- [25] J. Arrington and I. Sick, Phys.Rev. C70 , 028203 (2004), arXiv:nucl-ex/0406014 [nucl-ex] .
- [26] J. Arrington, W. Melnitchouk, and
- J. Tjon, Phys.Rev. C76 , 035205 (2007), arXiv:0707.1861 [nucl-ex] . [27] P. G. Blunden and I. Sick, Phys.Rev. C72 , 057601 (2005), arXiv:nucl-th/0508037 [nucl-th] .
- [28] P. Blunden, W. Melnitchouk, and J. Tjon, Phys.Rev. C72 , 034612 (2005), arXiv:nucl-th/0506039 [nucl-th] .
- [29] S. Kondratyuk, P. Blunden, W. Melnitchouk, and J. Tjon, Phys.Rev.Lett. 95 , 172503 (2005), arXiv:nucl-th/0506026 [nucl-th] .
Krzysztof M. Graczyk Cezary Juszczak
The Proton Radius from Bayesian Inference
The methods of Bayesian statistics are used to extract the value of the
proton radius from the elastic $ep$ scattering data in a model independent way.
To achieve that goal a large number of parametrizations (equivalent to neural
network schemes) are considered and ranked by their conditional probability
$P(\mathrm{parametrization}\,|\,\mathrm{data})$ instead of using the minimal
error criterion. As a result the most probable proton radii values
($r_E^p=0.899\pm 0.003$ fm, $r_M^p=0.879\pm 0.007$ fm) are obtained and
systematic error due to freedom in the choice of parametrization is estimated.
Correcting the data for the two photon exchange effect turns out to be
important to obtain the agreement between the $r_E^p$ and $r_M^p$ values. The
results disagree with recent muonic atom measurements.