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

CDF and PDF of order statistics

  1. May 28, 2012 #1
    Hi,

    I have K i.i.d. exponentially distributed random variables with mean unity. I need to find the CDF and PDF of the summation of the largest two random variables. How can I do that? The problem in this case is that the combinations are not independent.

    Thanks in advance
     
  2. jcsd
  3. May 28, 2012 #2

    AKG

    User Avatar
    Science Advisor
    Homework Helper

    EDIT: Oops, sorry, when I said it was clear that the [itex]X_i+X_j[/itex] are clearly iid I was mistaken. I'm not sure that they're not independent, but if they are independent, it certainly isn't clear.

    --------

    The CDF of this random variable, let's call it [itex]X[/itex], is given by:

    [tex]a \mapsto P(X < a)[/tex]

    What's the probability that the sum of the largest two of your exponentials is less than [itex]a[/itex]? Well, the sum of the largest two is less than [itex]a[/itex] iff the sum of every two is less than [itex]a[/itex]. In other words, if we let [itex]X_1, \dots, X_K[/itex] be your exponentials, then:

    [tex]P(X < a) = P(X_1 + X_2 < a\mbox{ and }X_1 + X_3 < a\mbox{ and } \dots \mbox{ and }X_{K-1} + X_K < a)[/tex]

    Since the [itex]X_i[/itex] are iid, [STRIKE]it's not hard to see that the [itex]X_i+X_j[/itex] are iid[/STRIKE]. As such, we can rewrite the above:

    [tex]P(X < a) = \prod_{1\leq i<j\leq K}P(X_i + X_j < a) = P(X_1 + X_2 < a)^{{K\choose 2}}[/tex]

    [itex]a \mapsto P(X_1 + X_2 < a)[/itex] is the CDF of a [itex]\Gamma (2,1)[/itex] distributed random variable.
     
    Last edited: May 29, 2012
  4. May 29, 2012 #3
    Interesting! But are the events X1+X2<a and X1+X3<a are independent?
     
  5. May 29, 2012 #4

    chiro

    User Avatar
    Science Advisor

    Hey S_David.

    The first thing is to get the distribution for the largest and second largest distribution. This can be done with order statistics.

    Once you have these distributions, then if they are of the same type you can use MGF's to find the result type of adding the two (usually this is a good idea because in many situations adding to distributions that are i.i.d results in the same distribution with different parameters).

    If the type is complex, use the convolution theorem to get the CDF and hence the PDF of the sum of the two variables. These variables have to be independent, but not necessarily identical distributed.

    For convolution algorithm and more depth:

    http://www.dartmouth.edu/~chance/teaching_aids/books_articles/probability_book/Chapter7.pdf

    Formula for order statistic given a known PDF can be found here:

    http://www.encyclopediaofmath.org/index.php/Order_statistic
     
  6. May 29, 2012 #5

    haruspex

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    CDF of highest of n = Fn(x)
    2nd highest < x if either highest < x or exactly one (any of the n) > x:
    CDF of 2nd highest = Fn(x) + (1-Fn(x))*n*Fn-1(x)
    = n*Fn-1(x) - (n-1)*Fn(x)
    CDF of sum of highest two = ∫y(n*Fn-1(x-y) - (n-1)*Fn(x-y)).dFn(y)
    = ∫y(n*Fn-1(x-y) - (n-1)*Fn(x-y))*n*Fn-1(y).dF(y)
     
  7. May 29, 2012 #6

    haruspex

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    On second thoughts, still an independence problem there.
    Try 2:
    P[2nd highest < x | highest = y > x] = (F(x)/F(y))n-1
    P[2nd highest < x | highest = y < x] = 1
    CDF of sum of highest two = Fn(x/2) + ∫y>x/2(F(x-y)/F(y))n-1.dFn(y)
    = Fn(x/2) + n∫y>x/2F(x-y)n-1.dF(y)
    Looks nicer at least.
     
    Last edited: May 29, 2012
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: CDF and PDF of order statistics
  1. PDF and CDF (Replies: 16)

  2. Order Statistics PDF (Replies: 2)

  3. Cdf and pdf (Replies: 4)

Loading...