Find E[(X-mu)^k] - Normal Random Variable

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
6 replies · 3K views
motherh
Messages
27
Reaction score
0
Hi,

I'm having a bit of a problem with a probability question. The question is

Let X be a normal random variable with mean [itex]\mu[/itex] and variance [itex]\sigma^{2}[/itex]. Find E[(X -[itex]\mu[/itex])[itex]^{k}[/itex]] for all k = 1,2,...

I'm not really sure what to do and need some help to confirm how to approach the question. I've tried messing around with moment generating functions, using the binomial theorem to expand E[(X -[itex]\mu[/itex])[itex]^{k}[/itex]] and also have tried computing the expectation directly. Are any of these methods the correct way to progress?

Any help is much appreciated, thanks!
 
Physics news on Phys.org
Letting [itex]m_k= \mathbb E[(X-\mu)^k][/itex], we know [itex]m_1=0[/itex] by definition of [itex]\mu[/itex]. Next, maybe try to use integration by parts (and the exact functional form of the p.d.f. [itex]\phi[/itex]) to express [itex]m_k[/itex] ([itex]k=2,3,...[/itex]) in terms of [itex]m_1,...,m_{k-1}[/itex].
 
Hi, thanks for your help so far! I'm a little bit confused though. What do I integrate by parts exactly? And I'm a little bit thrown by "the exact functional form of the p.d.f. ϕ" - what do you mean? Apologies, I use ϕ as part of my notation for moment generating functions.
 
Anyone? I'm still really confused. Do I need a binomial expansion of (X-mu)^k in my answer somewhere?
 
Hey. I was using the notation [itex]\phi(z) = \frac1{\sqrt{2\pi}}e^{-\frac{z^2}2}[/itex], the probability density function of [itex]\mathcal N(0,1)[/itex], so that the the p.d.f. of a [itex]\mathcal N(\mu,\sigma^2)[/itex] random variable is [itex]f(x)=\phi\left(\frac{x-\mu}{\sigma}\right)[/itex].

Then, we can express the desired numbers as [tex]m_k = \int_{-\infty}^\infty (x-\mu)^k f(x)\text{ d}x.[/tex]
 
To clean up notation, let's notice that [itex]m_k = \sigma^k n_k[/itex], where [itex]n_k = \mathbb E \left[ \left( \frac{x-\mu}{\sigma} \right)^k \right][/itex]. So we can just work with the system [tex]n_k = \int_{-\infty}^\infty z^k\phi(z)\text{ d}z,[/tex] which I would suggest attacking with IBP.
 
Ah, thank you so much! I got the relation m_(k+2) = sigma^2 * (k+1) * m_k, does that sound right?

With m_1 = 0 it means m_k = 0 for odd k which does seem to make sense.