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

Transformation of random variable

  1. Nov 11, 2011 #1

    I have a really tough time understanding this concept although this isn't anything more complex than a composition of functions.

    I have a done example from my book that I am trying to interpret.

    [itex] p_\xi (x)=\frac{1}{\pi} [/itex]

    for an interval between [0,∏]

    The transformation is given by law:

    [itex] \eta=sin(\xi) [/itex] interval from [0,∏]
    Find the PDF for η.

    Now I understand that sin(x) isn't monotonous on whole interval.

    So we have to divide this into 2 monotonous functions.

    From [0,pi/2] and from [pi/2, pi]

    And we use formula for that monotonous transformation 2 times.

    [PLAIN]http://pokit.org/get/d635e99a44c3c266a5c93e734d61af32.jpg [Broken]

    I understand everything up to the part where the final function is drawn.

    [PLAIN]http://pokit.org/get/b750c434c9d1e0ba259da0972158b514.jpg [Broken]

    This is pη graph which is given in my book.

    But I cannot wrap my head around it.

    First question:

    How can we just add those 2 monotonous transformations, aren't we losing pieces of information?

    And where is that second part, where transformation is monotonously decreasing?

    Any help with my intuition here?
    Last edited by a moderator: May 5, 2017
  2. jcsd
  3. Nov 11, 2011 #2

    Stephen Tashi

    User Avatar
    Science Advisor

    If that is supposed to be the graph of a probability density function, I don't understand it either. Is probability plotted on the horizontal axis? And why does the graph have y coordinates as large as 4.5?
    How did the book show the addition? ( The term is "monotonic" , not "monotonous", isn't it?)

    When you find the cumulative distribution y = F(x) it's clear you must add the probabilities of two disjoint intervals on the x-axis to get the values of the random variable that are equal or less than y.
  4. Nov 12, 2011 #3
    This is the final graph. Graph which the problem is asking for. The output function.

    This is the law of transformation:
    [PLAIN]http://pokit.org/get/76f8d6e186f55a03e47686496639aff5.jpg [Broken]

    ergo the sine wave from 0,pi

    and this is the input function:

    [PLAIN]http://pokit.org/get/72cd80aa4ae524097162ce20d8ebb105.jpg [Broken]

    ergo the uniformly distributed PDF
    Last edited by a moderator: May 5, 2017
  5. Nov 12, 2011 #4

    Stephen Tashi

    User Avatar
    Science Advisor

    The last two graphs make sense. The graph of [itex] p_{\eta}[/itex] doesn't make sense (to me). Are you taking these figures directly from your course materials?

    The figure showing the equation for [itex] p_\eta (y) [/itex] has "[itex] x = f_1^{-1}(x)[/itex]" instead of "[itex]x = f_1^{-1}(y)[/itex]", so it has errors also.

    I think the simplest way to understand why the "transformation" approach works is to first study the "distribution method", which finds the cumulative distribution. Have you studied that method?
  6. Nov 12, 2011 #5

    [PLAIN]http://pokit.org/get/efb0245c8d7ba024e4a8c3433dec0cb1.jpg [Broken]

    [PLAIN]http://pokit.org/get/90c6dc7a629e55a12b0e10c2e197f0f9.jpg [Broken]

    Here is how he got that.

    The fact that the y goes up to 1 is ok with me. But the, f(y) is not.
    Last edited by a moderator: May 5, 2017
  7. Nov 12, 2011 #6
    actually it does make sense. Because this is made out of 2 monotonous function, and the problem is ambiguity.

    The first part [0,pi/2] it goes from ~0.3 up until one and goes up to infinity. The second time it goes back from infinity to ~0.3. Add that up and will get ~0.6.

    I understand everything now. Thank you for effort anyway. You pointed out what I should be looking for. I did that in matlab, and everything adds up.
  8. Nov 12, 2011 #7
    Nope false alarm, Im clueless again. Matlab is showing the stuff we are seeing, when I am transforming these variables but I don't get it :D
  9. Nov 13, 2011 #8
    It's much easier if you look at the CDF: F(x) = P[X<=x] (which Stephen hinted at in post #2).

Share this great discussion with others via Reddit, Google+, Twitter, or Facebook