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

Finding an expression for the standard deviation.

  1. Jan 4, 2010 #1
    The problem is to find a closed form for the expectation (mu) and standard deviation (sigma) of the probability density function underneath (Maple notation):

    f(t) = a*exp(a*t)*BesselK(0, 2*sqrt(exp(a*t)))/BesselK(1, 2)

    I was able to find an expression for mu::

    mu = Int(t*a*exp(a*t)*BesselK(0, 2*sqrt(exp(a*t)))/BesselK(1, 2), t = 0 .. infinity)

    Substitution with x = exp(a*t) leads to:

    mu = Int((1/(x*a))*ln(x)*x*BesselK(0, 2*sqrt(x))/BesselK(1, 2),x=1..infinity)

    Partial integration integration leads to:

    -ln(x)*sqrt(x)*BesselK(1, 2*sqrt(x))/(a*BesselK(1, 2))-Int((-sqrt(x)*BesselK(1,
    2*sqrt(x))/(a*BesselK(1, 2)))*(1/x),x)

    Integration leads to:

    -(ln(x)*sqrt(x)*BesselK(1, 2*sqrt(x))+BesselK(0, 2*sqrt(x)))/(a*BesselK(1, 2))

    Finally the following closed form can be obtained:

    mu = BesselK(0, 2)/(a*BesselK(1, 2))

    I was not able to find a closed form for sigma. I suspect, however, that the solution must be proportional to: BesselK(0, 2)/(a*BesselK(1, 2)).

    Who can help me?

    Ad van der Ven.

    PS. Mijnhomepage: http://www.socsci.ru.nl/~advdv/
  2. jcsd
  3. Jan 4, 2010 #2


    User Avatar
    Science Advisor

    Did you try to get the second moment?
    mu2 = Int(t2*a*exp(a*t)*BesselK(0, 2*sqrt(exp(a*t)))/BesselK(1, 2), t = 0 .. infinity)
  4. Jan 5, 2010 #3
    Yes, I did, but I was not able to find a solution.
  5. Jan 5, 2010 #4
    There is something wrong with your PDF f(x), since
    int(f(x), t = 0 .. infinity)
    does not equal one for any value of the variable a (thus f(x) is not a PDF).
  6. Jan 6, 2010 #5
    I am sorry to tell you, but according to Maple 12

    Int(a*exp(a*t)*BesselK(0, 2*sqrt(exp(a*t)))/BesselK(1, 2),t=0..infinity) = 1.

    Note, that f(t) is only defined for 0 < t.
  7. Jan 6, 2010 #6
    f(t) = a*exp(a*t)*BesselK(0, 2*sqrt(exp(a*t)))/BesselK(1, 2) is defined for 0 < t and for 0 < a.
  8. Jan 7, 2010 #7
    You are right, f(t) does indeed integrate to 1 for a>0 (on 0<t<inf). I was using an older version of Maple that does not think the integral of the PDF converges.

    Your calculation of the mean is correct (I have verified now using Mathematica). I have not been able to find an analytical solution for the variance, but you can easily calculate it numerically, e.g. in Maple by

    Int(x*x*exp(x)*BesselK(0, 2*exp(x/2)),x=0..infinity)/BesselK(1, 2)/a^2
    -(Int(x*exp(x)*BesselK(0, 2*exp(x/2)),x=0..infinity)/BesselK(1, 2)/a)^2

    ans = 0.3715365403*1/(a^2)

    Don't know if that is of any help... (in above, I have moved the variable "a" out of the integrals using variable substitution)

  9. Jan 7, 2010 #8
    Hi Emanuel,

    Thanks a lot for looking into the problem. I already kwew, that you can calculate it numerically using Maple. I also know for certain, that the solution for the standard deviation is proportional to the mean. However, I do not know the expression for the proportionality constant. At the same time I am also rather sure, that the proportionality constant is some function of Bessel functions of the second kind (BesselK(v, x)). But I am aware that this will not be much of a help.

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