Processus de mélange et de division

9

Sur la figure suivante à gauche, deux réalisations de processus ponctuels avec une densité (intensité) différente et sont mélangées en faisant correspondre le centre des zones d'appartenance pour construire un processus ponctuel au milieu avec une intensité . Ensuite, des points sélectionnés au hasard sous forme de deux ensembles extraits de celui-ci, comme illustré à droite. Questions: Est-ce que ? et ? Si deux à gauche étaient du Poisson PP, celui du milieu est-il un Poisson PP? Et les deux à droite?λ1λ2λ

λ=λ1+λ2λ=λ3+λ4

entrez la description de l'image ici

Développeur
la source
6
Les mots clés que vous recherchez sont la superposition et l' amincissement du processus de Poisson. La réponse, avec quelques réserves, est oui . Mais, une réponse affirmative dépend intimement (i) de l'indépendance des deux processus dans le premier cas et (ii) de la façon dont le fractionnement est effectué dans le second cas. :)
Cardinal
1
Merci pour les mots clés. Je vous serais reconnaissant de bien vouloir fournir une explication complète en guise de réponse. Pour (i) puisque les deux sont des poissons PP, ils sont indépendants (je pense). Pour (ii), disons qu'un sélecteur aléatoire de Poisson peut être poursuivi.
Développeur
3
Comme l'a dit le cardinal, l'indépendance des processus ponctuels est importante. Vous pouvez facilement définir deux processus de poisson dépendants dont la superposition ne serait pas un processus de poisson; par exemple: disons que les points du processus n ° 2 sont exactement les mêmes que ceux du processus n ° 1, juste décalés vers la droite d'une unité.
Karl
1
@Karl: J'aime l'essence de votre exemple, bien que le deuxième processus ne soit pas tout à fait un processus de Poisson puisque la probabilité d'une arrivée dans est nulle dans le deuxième cas. :)[0,1)
cardinal
3
@cardinal - Je pensais aux processus de Point sur le plan complet.
Karl

Réponses:

12

Pour répondre à cette question, nous avons besoin d'un peu de contexte et de notation. Dans la terminologie générale, désigne un processus ponctuel dans le plan, ce qui signifie que pour tout ensemble de Borel, , dans le plan, est une variable aléatoire de valeur entière (y compris ), qui compte le nombre des points . De plus, est une mesure pour chaque réalisation du procédé du point .NAN(A)+AAN(A)N

Le processus ponctuel est associé à la mesure d'espérance où l'attente est toujours bien définie, puisque , mais peut être . Il s'agit d'un exercice pour vérifier que est à nouveau une mesure. Pour éviter les problèmes techniques, supposons que , ce qui est également raisonnable si le processus ne vit vraiment que sur un ensemble borné tel que la boîte de la figure que l'OP a publiée. Elle implique que que pour tous .

Aμ(A):=E(N(A))
N(A)0+μμ(R2)<N(A)<A

Les définitions et observations suivantes suivent.

  • On dit que a une intensité si a une densité rapport à la mesure de Lebesgue, c'est-à-dire si N λμλ
    μ(A)=Aλ(x)dx.
  • Si et sont deux processus ponctuels, nous définissons la superposition comme la somme . Cela revient à superposer un motif de points au-dessus de l'autre.N1N2N1+N2
  • Si et sont deux processus ponctuels (indépendants ou non) d'intensités et alors la superposition a une intensité .N1N2λ1λ2λ1+λ2
  • Si et sont des processus de Poisson indépendants , la superposition est un processus de Poisson. Pour montrer cela, nous observons d'abord que est Poisson à partir des propriétés de convolution de la distribution de Poisson, puis que si sont disjoints alors sont indépendants car et sont indépendants et les processus de Poisson eux-mêmes. Ces deux propriétés caractérisent un processus de Poisson. N1N2N1(A)+N2(A)A1,,AnN1(A1)+N2(A1),,N1(An)+N2(An)N1N2

Résumé I: Nous avons montré que chaque fois qu'un processus ponctuel est une somme, ou superposition, de deux processus ponctuels avec des intensités, alors la superposition a pour intensité la somme des intensités. Si, par ailleurs, les processus sont indépendants de Poisson, la superposition est Poisson.

Pour la partie restante de la question, nous supposons que comme pour tous les ensembles singleton . Ensuite, le processus ponctuel est appelé simple. Les processus de Poisson avec des intensités sont simples. Pour un processus ponctuel simple, il y a une représentation de comme c'est-à-dire comme une somme de mesures de Dirac aux points aléatoires. Si sont des variables aléatoires de Bernoulli, un amincissement aléatoire est le processus ponctuel simple Il est tout à fait clair qu'avec cela signifie que . Si nous IIDN({x})1{x}N

N=iδXi,
Zi{0,1}
N1=iZiδXi.
N2=i(1Zi)δXi
N=N1+N2amincissement aléatoire, ce qui signifie que les sont tous indépendants et identiquement distribués avec une probabilité de succès , disons, alors De là, Zip
N1(A)N(A)=nBin(n,p).
E(N1(A))=E(E(N1(A)N(A)))=E(N(A)p)=pμ(A).

Si est un processus de Poisson, il devrait être clair que pour puis sont à nouveau indépendants, et Cela montre que est un processus de Poisson. De même, est un processus de Poisson (avec mesure moyenneNA1,,AnN1(A1),,N1(An)

P(N1(A)=k)=n=kP(N1(A)=kN(A)=n)P(N(A)=n)=eμ(A)n=k(nk)pk(1p)nkμ(A)nn!=(pμ)kk!eμ(A)n=k((1p)μ(A))nk(nk)!=(pμ(A))kk!eμ(A)+(1p)μ(A)=epμ(A)(pμ(A))kk!.
N1N2(1p)μ). Il reste à montrer que et sont en fait indépendants. Nous prenons un coin ici et disons qu'il suffit en fait de montrer que et sont indépendants pour arbitraire , et cela découle de N1N2N1(A)N2(A)A
P(N1(A)=k,N2(A)=r)=P(N1(A)=k,N(A)=k+r)=P(N1(A)=kN(A)=k+r)P(N(A)=k+r)=eμ(A)(k+rk)pk(1p)rμ(A)k+r(k+r)!=epμ(A)(pμ(A))kk!e(1p)μ(A)((1p)μ(A))rr!=P(N1(A)=k)P(N2(A)=r).

Résumé II: Nous concluons que l'amincissement aléatoire iid avec probabilité de succès d'un processus ponctuel simple, , avec une intensité conduit à deux processus ponctuels simples, et , avec des intensités et , respectivement, et est la superposition de et . Si, de plus, est un processus de Poisson, alors et sont des processus de Poisson indépendants.pNλN1N2pλ(1p)λNN1N2NN1N2

Il est naturel de se demander si l'on pourrait éclaircir indépendamment sans supposer que les sont distribués de manière identique et obtenir des résultats similaires. C'est possible, mais un peu plus compliqué à formuler, car la distribution de doit alors être liée au manière ou d'une autre. Par exemple, pour une fonction donnée . Il est alors possible d'afficher le même résultat que ci-dessus mais avec l'intensité signifiant la fonction . Nous sautons la preuve. La meilleure référence mathématique générale couvrant les processus ponctuels spatiaux est Daley et Vere-JonesZiZiXiP(Zi=1N)=p(xi)ppλp(x)λ(x). Møller et Waagepetersen sont les deuxièmes plus proches des statistiques et des algorithmes de simulation .

NRH
la source
1
+1 La lecture de cette réponse est vraiment incroyable et utile. J'ai personnellement appris beaucoup de choses. C'est l'une des réponses les plus complètes que j'aie jamais reçues. J'apprécie vraiment cela.
Développeur
@Developer, merci. Heureux d'avoir pu vous aider.
NRH
C'est mieux qu'un manuel ...
Michael Mark
Merci pour votre réponse. Je pense que vous devez mentionner ici que pour les processus Point généraux, vous devez connaître l'intensité conditionnelle afin de pouvoir caractériser pleinement. Actuellement, ce que vous avez écrit pourrait être interprété comme quoi est constant. λ(t|Ht)λ
Sus20200
@ Sus20200, la densité conditionnelle, telle que vous l'écrivez, est utilisée pour les processus ponctuels temporels, tandis que la question concerne les processus ponctuels dans le plan sans ordre temporel. Sinon, je suis d'accord qu'il faut faire attention à distinguer l'intensité déterministe de l'intensité conditionnelle (ou stochastique). La première spécifie uniquement la mesure moyenne et non la distribution entière du processus ponctuel. Sauf pour un processus de Poisson, qui est complètement spécifié par sa mesure moyenne et donc son intensité. Notez que l'intensité n'est pas constante mais fonction des coordonnées spatiales.
NRH