I want to know how to treat propagation of errors in general.(adsbygoogle = window.adsbygoogle || []).push({});

When the transformation of variables is a transformation of $R^n$ to $R^n$ it

simply involves a jacobian:

[tex]g(\vec{y}) = f(x(\vec{y}))|J|[/tex]

With

[tex]J_{ij} = \frac{\partial x_i}{\partial y_j} [/tex]

(see http://pdg.lbl.gov/2005/reviews/probrpp.pdf for instance)

But there are also situation of $R^n$ to $R^m$ possible.

This is how far I got:

(Much of this can be found in http://arxiv.org/abs/hep-ex/0002056 appendix A)

The characteristic function (CF) is defined as:

[tex]\phi_X(t) = E[e^{itX}] = \int e^{itx} f(x) dx[/tex]

The inverse (FT):

[tex]f(x) = \frac{1}{2\pi} \int e^{-itX} \phi_X(t) dt[/tex]

If we have a function Y=g(X) the pdf according to Kendal and Stuart (1943) is:

[tex]\phi_Y(t) = \int e^{itg(x)} f(x) dx[/tex]

Taking the inverse, and some rewriting

[tex]f(y) = \frac{1}{2\pi} \int e^{-ity} \phi_Y(t) dt = \int f(x)\delta(y-g(x))dx[/tex]

In vector form

[tex]f(y) = \int f(\vec{x})\delta(y-g(\vec{x}))d\vec{x}[/tex]

Take for instance two independent variables taken from the same distribution:

[tex]Y = g(\vec{X}) = g(X_1, X_2) = X_1 + X_2[/tex]

Then the resulting pdf is with:

[tex]f(y) = \int f(x_1,x_2)\delta(y-x_1-x_2)dx_1 dx_2 = \int f(x_1) f(y - x_1 - u)\delta(u)dx_1 du = \int f(x_1) f(y - x_1) dx_1 [/tex]

Using:

[tex]u = y - x_1 - x_2 [/tex]

[tex]x_2 = y - x_1 - u [/tex]

[tex]dx_2 = du [/tex]

Which is a simply a convolution, which is a known result,

see for instance: http://en.wikipedia.org/wiki/Sum_of_normally_distributed_random_variables

Or for $R^1$ to $R^1$:

[tex]f(y) = \int f(x)\delta(y-g(x))dx = \int f(g^{-1}(y-u))\frac{1}{g'(x)}\delta(u)du = f(g^{-1}(y))\frac{1}{g'(x)} = f(x(y)) \frac{dx}{dy}[/tex]

Using:

[tex]u = y - g(x)[/tex]

[tex]du = \frac{dg(x)}{dx}dx[/tex]

[tex]dx = \frac{1}{g'(x)}du[/tex]

[tex]x = g^{-1}(y-u)[/tex]

Also a general result, just a change of variable.

But how to do this for the case of $R^n$ to $R^m$?, I don't know how to evaluate

the dirac delta function in general.

I did find on wikipedia (http://en.wikipedia.org/wiki/Dirac_delta_function)

[tex] \int_V f(\mathbf{r}) \, \delta(g(\mathbf{r})) \, d^nr = \int_{\partial V}\frac{f(\mathbf{r})}{|\mathbf{\nabla}g|}\,d^{n-1}r [/tex]

But I couldn't find a reference where this is explained.

But then again, if y is a vector, how to solve this?:

[tex]f(\vec{y}) = \int f(\vec{x})\vec{\delta}(\vec{y}-g(\vec{x}))d\vec{x}[/tex]

Anyone got some hints? Or am I going the wrong direction with this?

**Physics Forums - The Fusion of Science and Community**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Propagation of errors

Loading...

Similar Threads - Propagation errors | Date |
---|---|

I Error Propagation in Transcendental Equation | Nov 12, 2017 |

A Uncertainty Propagation of Complex Functions | Feb 26, 2017 |

A Uncertainty Propagation in Coupled Oscillator | Feb 22, 2017 |

I Propagating Error in both numerator and denominator of ratio | Jul 11, 2016 |

Error propagation for a sum of means | Feb 8, 2015 |

**Physics Forums - The Fusion of Science and Community**