Quelle est la distribution de la somme des variables aléatoires chi carré?

8

Quelle serait la distribution de l'équation suivante:

y=a2+2ad+d2

où et sont des variables aléatoires chi carré non centrales indépendantes avec degrés de liberté.ad2M

OBS.: Les RV générant à la fois et ont et , disons .adμ=0σ21σ2=c

Felipe Augusto de Figueiredo
la source
2
1. Comment sont et liées? 2. Les variables aléatoires du chi carré ont déjà une moyenne> 0 Pourquoi auriez-vous besoin de l'indiquer explicitement? (Ou essayez-vous de faire référence à un chi carré non central?)ad
Glen_b -Reinstate Monica
Je viens d'ajouter quelques informations supplémentaires à la question. Ce sont des RV chi carré non centraux car ils ont été générés par des variables aléatoires gaussiennes complexes symétriques circulaires non standard.
Felipe Augusto de Figueiredo
2M sont les degrés de liberté pour chacun des deux?
Alecos Papadopoulos
2
Felipe, dans votre question , vous dites et faire « ont » , mais maintenant dans votre dernier commentaire que vous ils déclarent ne pas avoir cette propriété. Lequel est-ce?? ad μ=0
whuber
2
Merci d'avoir essayé d'expliquer, mais je ne comprends toujours pas. Lorsque vous écrivez « et sont des variables aléatoires chi carré non centrales indépendantes», il semble que vous additionnez des carrés de variables aléatoires normales qui ont des moyennes non nulles, car c'est ainsi que les variables non centrales chi carré apparaissent généralement. Mais plus tard, vous écrivez «Les RV générant à la fois et ont », ce qui suggère que vous travaillez avec des variables centrales du chi carré. Je soupçonne que ce sont les incohérences qui ont incité le commentaire initial de @Glen_b. Pourriez-vous montrer explicitement ce etadadμ=0adsont?
whuber

Réponses:

7

Si sont indépendants, alors aura la . Puisque est non négatif, CDF de peut être trouvé en notantPar conséquent,a,dχ2M2X=a+dχ4M2XY=a2+2ad+d2=(a+d)2=X2

FY(y)=P(Yy)=P(X2y)=P(Xy)=FX(y).
fY(y)=12yfX(y)=122M+1Γ(2M)yM1ey/2.

Si et sont corrélés, les choses sont beaucoup plus complexes. Voir par exemple la fonction de distribution cumulative de NH Gordon et PF Ramig de la somme des variables aléatoires du chi carré corrélé (1983) pour une définition du chi carré multivarié et la distribution de sa somme.ad

Si vous avez affaire à un chi carré non central, donc ce qui précède ne sera plus valide. Ce message peut fournir un aperçu.μ2M

EDIT: Sur la base des nouvelles informations, il semble et soient formés en additionnant le rv normal avec une variance non unitaire. Rappelez-vous si alors . Depuis maintenant deux auront une distribution chi carré mise à l'échelle par , c'est-à-dire la distribution . Dans ce cas, sera distribué par . Par conséquent, pour nous avonsadZN(0,1)cZN(0,c)

a=ci=12MZi2=d,
a,dcΓ(M,2c)X=a+dΓ(2M,2c)Y=X2
fY(y)=12(2c)2MΓ(2M)yM1ey/2c.
Francis
la source
Comment mu entre-t-il? Est-il supposé être la moyenne d'une des variables khi-deux? Je soupçonne que cela n'a rien à voir avec le problème.
Michael R. Chernick
@MichaelChernick: signifie probablement que peut être un chi carré non central? a,d
Francis
Je suppose que vous pouvez faire cette hypothèse mais l'OP ne fait aucune connexion. Je pense que vous avez adopté la bonne approche, le non-central ne pouvait pas entrer dans ce problème. X est le carré d'un chi carré ici. Dans le cas de l'indépendance que vous avez utilisée ici, comment s'appelle cette distribution?
Michael R. Chernick
@MichaelChernick Je ne sais pas s'il y a un nom spécial associé à la distribution. "chi-tesseracté" peut-être?
Francis
a et sont des chi-carrés non centraux. d
Felipe Augusto de Figueiredo
4

Puisqu'un khi carré non central est une somme de RV indépendants, alors la somme de deux chi carrés non centraux indépendants est également un chi carré non central avec des paramètres la somme des paramètres correspondants de les deux composantes, (degrés de liberté), (paramètre de non centralité).X=a+bkx=ka+kbλx=λa+λb

Pour obtenir la fonction de distribution de son carré , on peut appliquer la "méthode CDF" (comme dans la réponse @francis),Y=X2

FY(y)=P(Yy)=P(X2y)=P(Xy)=FX(y)

et où

FX(x)=1Qkx/2(λx,x)

donc

FY(y)=1Qkx/2(λx,y1/4)

où est ici la fonction Q de Marcum .Q

Ce qui précède s'applique aux chi-carrés non centraux formés comme des sommes de normales quadratiques indépendantes chacune avec une variance unitaire mais une moyenne différente.

ADDENDA RÉPONDANT À LA MODIFICATION DE LA QUESTION

Si les RV de base sont , alors le carré de chacun est un voir https://stats.stackexchange.com/a/122864/28746 .N(0,c)Gamma(1/2,2c)

Donc, le rv et donc aussi (paramétrage de l'échelle de forme, et voir l'article de wikipedia pour le propriétés additives pour Gamma). aGamma(M,2c)bGamma(M,2c)X=a+bGamma(2M,2c)

Ensuite, on peut appliquer à nouveau la méthode CDF pour trouver le CDF du carréY=X2

Alecos Papadopoulos
la source
@FelipeAugustodeFigueiredo Désolé, je ne connais pas les VR complexes. Ma réponse a pris comme étant donné que nous partons de chi-carrés non centraux.
Alecos Papadopoulos
Et si les RV sont des variables aléatoires gaussiennes complexes symétriques circulaires avec et ? μ=0σ=cI
Felipe Augusto de Figueiredo
oublions les VR complexes. Et si les RV générant et sont des RV gaussiens avec et ? Tous les RV gaussiens ont la même variance, appelons-le . adμ=0σ1c
Felipe Augusto de Figueiredo
pourriez-vous s'il vous plaît m'aider avec la question suivante: stats.stackexchange.com/questions/253764/… . Tout indice serait très apprécié. Merci!
Felipe Augusto de Figueiredo
@FelipeAugustodeFigueiredo J'ai bien peur de n'avoir rien à offrir pour cette question.
Alecos Papadopoulos