Is there any "easy" way to calculate fractional moments from Laplace transform. To be more specyfic let us consider the following example. Let $X$ be a positive random variable and $L(\theta) := E \exp (\theta X)$ be its Laplace transform. Of course it is easy to calculate $E X^n$ where $n$ is a natural number but what with e.g. $E X^{1/2}$.
If $X$ is positive the following works. Let $F(\theta)=E(e^{\theta X})$ be the Laplace transform. Given $s\in \mathbb{R}$, write $s=n\alpha$ with $n$ an integer and $\alpha >0$. Then $$E(X^{s}) = (1)^n\frac{1}{\Gamma(\alpha)} \int_0^\infty F^{(n)} (\theta) \theta^{\alpha1} d \theta$$ with $\Gamma$ the usual Gamma function, $$\Gamma(\alpha) = \int_0^\infty \theta^{\alpha 1} e^{\theta} d \theta.$$
Indeed by Tonelli's theorem and then a change of variable, $$\int_0^\infty F^{(n)}(\theta) \theta^{\alpha1} d \theta = (1)^n E \Bigl(\int_0^\infty X^n e^{\theta X} \theta^{\alpha1} d \theta \Bigr)= (1)^n E(X^{n\alpha}) \int_0^\infty \theta^{\alpha1} e^{\theta} d \theta,$$ and so long as $\alpha >0$ the integral on the right is convergent.

$\begingroup$ How do I choose the parameter 'Alpha' ? In the question 's'=1/2. Is 'n'=1, 'Alpha'=1/2 the only choice? Or can I use as well 'n'=2, 'Alpha'=3/2 with no difference ? $\endgroup$ Jan 8 '14 at 10:16

$\begingroup$ Yes. You can use $alpha =3/2$ and $n=2$. This amounts to integrating by parts once in the integral over $\theta$. $\endgroup$ Jan 8 '14 at 18:22

For your example of X^{1/2} you can evaluate the fractional half derivative of the Laplace transform (see for example the Wikipedia article on fractional calculus) at theta = 0.

$\begingroup$ I thought about this also. I have even tried to use it but it does not very handy to use. I was hoping for something more "natural". Anyway, I have another question. Let us now drop the assumption that X is positive but assume that the Laplace transform exists. How can one calculate $EX$? $\endgroup$ Nov 14 '09 at 18:07

$\begingroup$ If you write the Laplace transform of the probability measure of X as a sum of two pieces, the first corresponds to the integral from inf to 0 and the second from 0 to +inf, i.e., F(s) = F+(s) + F(s), then the Laplace transform of the probability measure of X will be F+(s)+F(s), from which you can derive the moments of X. $\endgroup$ Nov 16 '09 at 6:18

$\begingroup$ Yes, but the problem is that I have the Laplace transform of "whole" X that is F(s). How can I calculate F+ and F using F only? May be I missing something elementary, but I can not see how to do it in an easy way (of course in principle one can apply the inverse Laplace transform and so on but this is not easy usualy ). $\endgroup$ Nov 17 '09 at 14:31