Quelle est la distribution de la moyenne arrondie des variables aléatoires de Poisson?

20

Si j'ai des variables aléatoires qui sont des distributions de Poisson avec les paramètres , quelle est la distribution de (c'est-à-dire le plancher entier de la moyenne)?X1,X2,,Xnλ1,λ2,,λnY=i=1nXin

Une somme de Poissons est également Poisson, mais je ne suis pas assez confiant dans les statistiques pour déterminer s'il en est de même pour le cas ci-dessus.

Lubo Antonov
la source
@amoeba J'ai annulé votre modification du titre car il ne s'agit pas en fait d'un "arrondi". Le montage précédent de Cardinal, bien que pas tout à fait aussi précis, semble préférable car il est précis.
whuber
@whuber D'accord. J'hésitais lors de cette modification, mais j'ai décidé d'inclure le mot "arrondi" car actuellement le titre ne fait pas allusion à la principale difficulté ici (et est donc en quelque sorte trompeur). Le terme approprié devrait être "arrondi vers le bas", alors peut-être "Quelle est la distribution d'une moyenne de variables aléatoires de Poisson, arrondie vers le bas ?" - même si je dois admettre que cela semble un peu lourd.
amibe dit Réintégrer Monica
@amoeba D'autres modifications sont bien sûr les bienvenues!
whuber

Réponses:

27

Une généralisation de la question demande la distribution de Y=X/m lorsque la distribution de X est connue et appuyée sur les nombres naturels. (Dans la question, X a une distribution de Poisson du paramètre λ=λ1+λ2++λn et m=n .)

La distribution de Y est facilement déterminé par la distribution de mY , dont la fonction de génération probabilité (PGF) peut être déterminée en fonction de la PGF de X . Voici un aperçu de la dérivation.


Écrivez pour le pgf de , où (par définition) . est construit à partir de de telle sorte que son pgf, , soitp(x)=p0+p1x++pnxn+Xpn=Pr(X=n)mYXq

q(X)=(p0+p1++pm-1)+(pm+pm+1++p2m-1)Xm++(pnm+pnm+1++p(n+1)m-1)Xnm+.

Parce que cela converge absolument pour , nous pouvons réorganiser les termes en une somme de morceaux du formulaire|X|1

Dm,tp(x)=pt+pt+mxm++pt+nmxnm+

pour . La série de puissance des fonctions constituée de chaque terme de la série de commençant par le : on parle parfois de décimation de . Les recherches Google ne fournissent actuellement pas beaucoup d'informations utiles sur les décimations, donc pour être complet, voici une dérivation d'une formule.x t D m , t p m th p t tht=0,1,,m1xtDm,tpmthptthp

Soit n'importe quelle racine primitive de l'unité; par exemple, prenez . Il résulte alors de et quem th ω = exp ( 2 i π / m ) ω m = 1 m - 1 j = 0 ω j = 0ωmthω=exp(2iπ/m)ωm=1j=0m1ωj=0

xtDm,tp(x)=1mj=0m1ωtjp(x/ωj).

Pour voir cela, notez que l'opérateur est linéaire, il suffit donc de vérifier la formule sur la base . Appliquer le côté droit à donne { 1 , x , x 2 , , x n , } x nxtDm,t{1,x,x2,,xn,}xn

xtDm,t[xn]=1mj=0m1ωtjxnωnj=xnmj=0m1ω(tn)j.

Lorsque et diffèrent d'un multiple de , chaque terme de la somme est égal à et nous obtenons . Sinon, les termes parcourent les puissances de et ces sommes à zéro. D'où cet opérateur conserve toutes les puissances de congru à modulo et tue toutes les autres: c'est précisément la projection souhaitée.tnm1xnωtnxtm

Une formule pour suit facilement en changeant l'ordre de sommation et en reconnaissant l'une des sommes comme géométrique, l'écrivant ainsi sous forme fermée:q

q(x)=t=0m1(Dm,t[p])(x)=t=0m1xt1mj=0m1ωtjp(ωjx)=1mj=0m1p(ωjx)t=0m1(ωj/x)t=x(1xm)mj=0m1p(ωjX)X-ωj.

Par exemple, le pgf d'une distribution de Poisson du paramètre est . Avec , et le pgf de seraλp(X)=exp(λ(X-1))m=2ω=-12Oui

q(x)=x(1x2)2j=021p((1)jx)x(1)j=x1/x2(exp(λ(x1))x1+exp(λ(x1))x+1)=exp(λ)(sinh(λx)x+cosh(λx)).

Une utilisation de cette approche est de calculer les moments de et . La valeur de la dérivée du pgf évalué à est le moment factoriel . Le moment est une combinaison linéaire des premiers moments factoriels. En utilisant ces observations, nous trouvons, par exemple, que pour un Poisson distribué , sa moyenne (qui est le premier moment factoriel) est égale à , la moyenne de est égal à , et la moyenne de est égale àXmYkthx=1kthkthkXλ2(X/2)λ12+12e2λ3(X/3)λ1+e3λ/2(sin(3λ2)3+cos(3λ2)) :

Veux dire

Les moyennes pour sont affichées en bleu, rouge et jaune, respectivement, en fonction de : asymptotiquement, la moyenne diminue de par rapport à la moyenne de Poisson d'origine.m=1,2,3λ(m1)/2

Des formules similaires pour les variances peuvent être obtenues. (Ils deviennent désordonnés lorsque augmente et sont donc omis. Une chose qu'ils établissent définitivement est que lorsque aucun multiple de n'est Poisson: il n'a pas l'égalité caractéristique de moyenne et de variance) Voici un tracé des variances en fonction de pour :mm>1Yλm=1,2,3

Écarts

Il est intéressant de noter que pour des valeurs plus élevées de les variances augmentent . Intuitivement, cela est dû à deux phénomènes concurrents: la fonction de plancher consiste en fait à regrouper des groupes de valeurs qui étaient à l'origine distinctes; cela doit entraîner une diminution de la variance . En même temps, comme nous l'avons vu, les moyens changent aussi (car chaque bac est représenté par sa plus petite valeur); cela doit entraîner l'ajout d'un terme égal au carré de la différence de moyennes. L'augmentation de la variance pour les grands devient plus grande avec des valeurs plus élevées de .λλm

Le comportement de la variance de avec est étonnamment complexe. Terminons par une simulation rapide (in ) montrant ce qu'elle peut faire. Les graphiques montrent la différence entre la variance de et la variance de pour Poisson distribué avec différentes valeurs de allant de à . Dans tous les cas, les parcelles semblent avoir atteint leurs valeurs asymptotiques à droite.mYmRmX/mXXλ15000

set.seed(17)
par(mfrow=c(3,4))
temp <- sapply(c(1,2,5,10,20,50,100,200,500,1000,2000,5000), function(lambda) {
  x <- rpois(20000, lambda)
  v <- sapply(1:floor(lambda + 4*sqrt(lambda)), 
              function(m) var(floor(x/m)*m) - var(x))
  plot(v, type="l", xlab="", ylab="Increased variance", 
       main=toString(lambda), cex.main=.85, col="Blue", lwd=2)
})

Parcelles

whuber
la source
1
C'est une excellente réponse! Il me faudra probablement un peu de temps pour digérer :)
Lubo Antonov
1
et c'est pourquoi j'ai dit "Utiliser la fonction de plancher ... affecte aussi légèrement la variance mais de manière plus compliquée."
Henry
1
+1 Merci pour la réponse détaillée. Il existe certainement des façons compliquées dans lesquelles la fonction de plancher affecte la variance.
Dilip Sarwate
1
+1 pour la simulation en R avec du code --- c'est un très bel exemple d'utilisation sapply()pour la simulation. Merci.
Assad Ebrahim
1
@Roberto Merci. Cependant, la distinction entre " " et " ", étant purement une question de notation, est tout à fait triviale et n'a aucune importance mathématique ou statistique. xs
whuber
12

Comme le dit Michael Chernick, si les variables aléatoires individuelles sont indépendantes, la somme est Poisson avec le paramètre (moyenne et variance) que vous pourriez appeler .i=1nλiλ

La division par réduit la moyenne à et la variance sorte que la variance sera inférieure à la distribution de Poisson équivalente. Comme Michael le dit, toutes les valeurs ne seront pas des entiers.nλ/nλ/n2

L'utilisation de la fonction de plancher réduit légèrement la moyenne, d'environ , et affecte légèrement la variance mais de manière plus compliquée. Bien que vous ayez des valeurs entières, la variance sera toujours sensiblement inférieure à la moyenne et vous aurez donc une distribution plus étroite que le Poisson.1212n

Henri
la source
merci, pas un résultat que je peux utiliser, mais au moins je sais maintenant :)
Lubo Antonov
Si les lambdas ne sont pas tous égaux, le résultat ne devrait-il pas ressembler davantage à un binôme négatif qu'à un Poisson (en ignorant pour l'instant la partie non entière)? Qu'est-ce que j'oublie ici?
gung - Rétablir Monica
2
@gung: Il vous manque le point que l'individu n'affecte que la distribution à travers leur somme et combien il y en a. Peu importe les valeurs particulières qu'ils prennent: donnera le même résultat que . λ 1 = 1 , λ 2 = 2 , λ 3 = 9 λ 1 = 4 , λ 2 = 4 , λ 3 = 4λiλ1=1,λ2=2,λ3=9λ1=4,λ2=4,λ3=4
Henry
10

La fonction de masse de probabilité de la moyenne de variables aléatoires de Poisson indépendantes peut être écrite explicitement, bien que la réponse ne vous aide pas beaucoup. Comme Michael Chernick l'a noté dans les commentaires sur sa propre réponse, la somme des variables aléatoires de Poisson indépendantes avec les paramètres respectifs est une variable aléatoire de Poisson avec le paramètre . Par conséquent, Ainsi, est une variable aléatoire prenant la valeur avec probabilitéi X i X i λ i λ = i λ i P { n i = 1 X i = k } = exp ( - λ ) λ kn iXiXiλiλ=iλi

P{i=1nXi=k}=exp(λ)λkk!,  k=0,1,2,,
Y^=n1i=1nXik/nexp(λ)λkk! . Notez que n'est pas une variable aléatoire à valeur entière (bien qu'elle prenne des valeurs rationnelles uniformément espacées). Il s'ensuit facilement que est une variable aléatoire à valeur entière prenant la valeur avec une probabilité Ce n'est pasY^Y=Y^m
P{Y=m}=P{1ni=1nXi=m}=exp(λ)i=0n1λmn+i(mn+i)!,  m=0,1,2,,
la fonction de masse de probabilité d'une variable aléatoire de Poisson. Les formules pour la moyenne et la variance peuvent être écrites en utilisant cette fonction de masse de probabilité, mais elles ne conduisent évidemment pas à de belles réponses simples en termes de et . Des valeurs approximatives peuvent être obtenues comme l'a souligné Henry.λn
Dilip Sarwate
la source
1 Il sont des formules fermées pour les moments de , bien que. Y
whuber
Merci pour la formulation rigoureuse! Y a-t-il une chance que vous aimeriez essayer les formules de moyenne et de variance?
Lubo Antonov
2
Peut-être que @whuber affichera un lien (ou une citation d'un livre ou d'un article de journal) où les formules fermées pour les moments peuvent être trouvées, ou écrira une réponse donnant les formules elles-mêmes, avec ou sans dérivation détaillée.
Dilip Sarwate
@Dilip Ma réclamation concernant les formules fermées ne reposait sur rien de publié, j'ai donc posté une réponse séparée indiquant ce que j'avais en tête et comment cela pourrait être utilisé pour comprendre cette situation.
whuber
3

Y ne sera pas Poisson. Notez que les variables aléatoires de Poisson prennent des valeurs entières non négatives. Une fois que vous avez divisé par une constante, vous créez une variable aléatoire qui peut avoir des valeurs non entières. Il aura toujours la forme du Poisson. C'est juste que les probabilités discrètes peuvent se produire à des points non entiers.

Michael R. Chernick
la source
Oui
@ lucas1024 Je ne pense pas, mais je ne suis pas sûr.
Michael R. Chernick
Xjen-1
@JDav La somme est Poisson avec le paramètre de taux égal à la somme des paramètres de taux individuels. Mais l'OP évolue de 1 / n et veut ensuite tronquer le entier juste en dessous de Y. Je ne sais pas exactement ce que cela fait à la distribution.
Michael R. Chernick
Mon commentaire précédent supposait l'indépendance.
Michael R. Chernick