Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Error propagation question

  1. Mar 23, 2008 #1

    I've been trying to reproduce the output of an analytical machine here at work by doing the calculation myself in Excel but I can't get the error to match. Perhaps I am propagating the error incorrectly...

    The calculation is z = 2^(x-y)

    The values are:
    x = 24.96
    y = 25.98
    and the std. devs are:

    I figured that the error in z should be:
    delta z = z * sqrt((2*deltax/x)^2 + (2*deltay/y)^2)
    which would give 0.010958

    However, the software/machine reports the error as 0.07
    I am not 100% sure how the software is calculating the error but it says it uses the standard deviations.

    Am I making a mistake in my formula for error propagation? (I haven't done this in a while, frankly).

    Any advice much appreciated.

  2. jcsd
  3. Mar 23, 2008 #2
    Remember that whatever variation is present in x and y will be exponentiated by the time it appears in Z. I am not sure about your error analysis method, but I prefer to go back to basics and do a simulation:

    Code (Text):

    x = 24.96;
    y = 25.98;

    sx = 0.265016;
    sy = 0.085049;

    (* this simulates a million measurements of X and Y
     by drawing from a normal distribution with the
     mean and standard deviation you have given *)

    distX = RandomReal[NormalDistribution[x, sx], 1000000];
    distY = RandomReal[NormalDistribution[y, sy], 1000000];

    distZ = 2^(distX - distY);


    I get a somewhat different answer from the instrument you are reading, so either the measurements are not normally distributed (i.e. the instrument is automatically correcting for a systematic bias) or else the standard deviations are not yet representative of the true population.

    I attached a bitmap graphic of the Histograms for all three variables in this simulation, although for visual clarity I reduced the size of the simulation from 1 million measurements to only 1 thousand. The main feature to notice is that the Z distribution is slightly skewed towards larger values, likely as a result of exponentiating x - y. Hopefully someone more knowledgeable about statistics can explain how to properly quantify the error in the Z distribution.

    Attached Files:

    Last edited: Mar 23, 2008
  4. Mar 23, 2008 #3


    User Avatar
    Science Advisor
    Homework Helper

    Well in my opinion the error should be given by: (this is the error formula I've always been taught)

    [tex]\sigma^2_z = (\frac{\partial z}{\partial x})^2\sigma^2_x + (\frac{\partial z}{\partial y})^2\sigma^2_y}[/tex]

    which gives an error of 0.095135
  5. Mar 24, 2008 #4
    Thanks to you both.
    I have to say that you both give very reasonable answers - and indeed they report similar errors. I don't know why the value reported by the machine's software is different (mine was just plain wrong). Indeed it may be that it takes so other invisible factor into account. I'll have to ask the manufacturer.
    Thanks again,
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook