Given 𝑛 independent Normally distributed random variables 𝑋ᵢ ∼ 𝑁(𝜇ᵢ,𝜎²ᵢ) and 𝑛 real constants 𝑎ᵢ∈ℝ, I need to find an acceptable Normal approximation of the distribution of 𝑌 random variable (assuming Pr[𝑋ᵢ≤0]≈0, to avoid divisions by zero)

Y = ∑aᵢXᵢ / ∑Xᵢ

I thought to split 𝑌 into single components

Y = a₁X₁ / ∑Xᵢ + a₂X₂ / ∑Xᵢ + ... + aₙXₙ / ∑Xᵢ

Y = a₁Y₁ + a₂Y₂ + ... + aₙYₙ

where the distribution of each 𝑌ᵢ can be found noting that

Yᵢ = Xᵢ / (Xᵢ + ∑Xⱼ) , for j≠i

and that

1/Yᵢ = (Xᵢ + ∑Xⱼ) / Xᵢ = 1 + ∑Xⱼ / Xᵢ

so, calling ∑Xⱼ = Uᵢ we can say that Xᵢ and Uᵢ are independent and, according to Díaz-Francés et Al. 2012, a Normal approximation of 1/𝑌ᵢ can be the one in figure 1 and, considering 1 ~ N(1,0), the r.v. Yᵢ can be approximated to figure 2. Thus, approximation of each aᵢYᵢ is the one in figure 3.

But now... I'm stuck at the sum of aᵢYᵢ because, not being independent, I don't know how to approximate the variance of their sum.

Any advice? Any more straightforward or more efficient method?

More Max Pierini's questions See All
Similar questions and discussions