Given $X_i \sim \mathcal{N}(\mu_i,\sigma_i^2)$, for $i = 1,\dots,n$. How does one find the distribution of $D = \sum_{i=1}^n X_i^2$? In the case that all the standard deviations are the same (i.e. $\sigma_i = \sigma_1$ for all $i$), the random variable $D/\sigma_1^2$ has a noncentral chi-squared distribution. But when I don't have this simplifying assumption I don't know what to do. Is there a well known distribution for this?
Remember to vote up questions/answers you find interesting or helpful (requires 15 reputation points)
|
4
|
|
|||
|
1
|
All you could conceivably want to know about the subject (and many things you might not) are in Mathai + Provost, Quadratic Forms in Random Variables. |
||
|
|
You can accept an answer to one of your own questions by clicking the check mark next to it. This awards 15 reputation points to the person who answered and 2 reputation points to you.
|
1
|
I can write a formula, but I doubt you will like it. Set $$ \rho(x_1,\dotsc, x_n)=\frac{1}{(2\pi)^{\frac{n}{2}}\sigma_1\cdots\sigma_n} \exp\left(-\sum\frac{(x_i-\mu_i)^2}{2\sigma_i^2}\right) $$ denote the joint probability density of $(X_1,\dotsc, X_n)$ which I assumed to be independent. The probability density of $D$ is $$ d(t)=\frac{1}{2\sqrt{t}}\int_{|x|=\sqrt{t}} \rho(x) dS(x), $$ where $dS(x)$ denotes the area element on the sphere ${|x|=\sqrt{t}}$. |
|||
|
|
$D = sum X_i^2$not just the sum as written; the latter would follow a normal distribution. Could you clarify? – quid Feb 28 2012 at 17:29