MathOverflow is a question and answer site for professional mathematicians. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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?

share|cite|improve this question
4  
From context it seems you want $D = sum X_i^2$ not just the sum as written; the latter would follow a normal distribution. Could you clarify? – user9072 Feb 28 '12 at 17:29
    
Yes, I made a mistake originally, but Douglas Zare fixed it. Thanks Douglas! – Jonathan Feb 28 '12 at 22:17
1  
This is something that I don't remember well from several courses in which it was touched on. I'd start by googling "Satterthwaite". – Michael Hardy Feb 29 '12 at 2:31
up vote 4 down vote accepted

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.

share|cite|improve this answer

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}\}$.

share|cite|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.