somme des variables aléatoires du chi carré non central

21

J'ai besoin de trouver la distribution de la variable aléatoire

Y=i=1n(Xi)2
XiN(μi,σi2) et tous les Xi s sont indépendants. Je sais qu'il est possible de trouver d 'abord le produit de toutes les fonctions génératrices de moments pour Xi s, puis de retransformer pour obtenir la distribution de YCependant, je me demande s'il existe une forme générale pour Y comme le cas gaussien: nous savons que la somme de gaussien indépendant est encore gaussienne, et donc nous n'avons besoin que de connaître la moyenne et la variance sommées.

Que diriez-vous de tout ? Cette condition fera-t-elle une solution générale?σi2=σ2

piège
la source
1
En regardant le premier paragraphe ici , clairement la condition finale donne un chi carré non centrale échelle (diviser par par (le facteur d'échelle vous prendre à l'avant) et faire σ i = 1 dans Σ k i = 1 ( X i / σ i ) 2 ). La forme plus générale avec laquelle vous avez commencé ressemble à une combinaison linéaire ou à une moyenne pondérée échelonnée, avec des coefficients σ 2 i plutôt qu'une simple somme de carrés échelonnés ... et je crois que la distribution requise ne sera généralement pas requise. σ2σi=1i=1k(Xi/σi)2σi2
Glen_b -Reinstate Monica
Selon ce dont vous avez besoin, dans certains cas spécifiques, vous pouvez effectuer une convolution numérique ou une simulation.
Glen_b -Reinstate Monica
Ceci est généralisé par la distribution de la «somme pondérée des log chi-carrés au pouvoir». Mon package R sadistsfournit des fonctions «dpqr» approximatives pour ; cf github.com/shabbychef/sadistsY
shabbychef

Réponses:

17

Comme Glen_b l'a noté dans les commentaires, si les variances sont toutes les mêmes, vous vous retrouvez avec un chi carré non central à l'échelle.

Sinon, il existe un concept de distribution khi carré généralisée , c'est-à-dire pour x N ( μ , Σ ) et A fixe. Dans ce cas, vous avez le cas particulier de la diagonale Σ ( Σ i i = σ 2 i ), et A = I .xTAxxN(μ,Σ)AΣΣii=σi2A=I

Il y a eu du travail sur le calcul des choses avec cette distribution:

Vous pouvez également l'écrire sous la forme d'une combinaison linéaire de variables chi carré non centrales indépendantes , auquel cas:Y=i=1nσi2(Xi2σi2)

Bausch (2013) propose un algorithme plus efficace en termes de calcul pour la combinaison linéaire de khi deux centrés; son travail pourrait être extensible aux chi-carrés non centraux, et vous pourriez trouver des pointeurs intéressants dans la section de travail connexe.

Dougal
la source
2
Une comparaison des méthodes d'approximation se trouve dans Duchesne et al. 2010. Statistiques computationnelles et analyse des données, 54, 858–862. Les auteurs maintiennent le package R CompQuadForm avec des implémentations.
caracal
-10

Ce sera le chi carré de n degrés de liberté.

Ahmed
la source
6
Je crois que vous avez oublié que le peut être différent de zéro. Les commentaires à la question ainsi que la réponse existante sont informatifs. μi
whuber