Iteració de Panjer

De testwiki
Salta a la navegació Salta a la cerca

La iteració de Panjer és un algorisme per calcular l'aproximació a la distribució de probabilitat d'una variable aleatòria composta S=i=1NXi. on N i Xi són variables aleatòries i de tipus especials. En casos més generals la distribució de S és una distribució composta. La recursió per als casos especials considerats va ser introduïda en un article[1] de Harry Panjer (professor emèrit de la Universitat de Waterloo).[2] És molt utilitzat en ciències actuarials (vegeu també risc sistèmic).

Preliminars

Estem interessats en la variable aleatòria composta S=i=1NXi on N i Xi compleixen les condicions prèvies següents:

Distribució de la mida de la seqüència

Suposem que Xi sigui i.i.d. i independent de N. A més, Xi han de ser distribuïts en una xarxa h0 amb ample de xarxa h>0.

fk=P[Xi=hk].

En pràctica actuarial, Xis'obté per discretització de la funció de densitat de la seqüència (superior, inferior ...).

Distribució de números de la seqüència

El número de la seqüència N és una variable aleatòria, que es diu que té una «distribució de números de la seqüència», i que pot prendre valors 0, 1, 2, ..., etc. Per a la recursió de Panjer, la distribució de probabilitat de N ha de ser membre de la «classe Panjer», també coneguda com la classe (a, b, 0) de les distribucions. Aquesta classe consta de totes les variables aleatòries que compleixen la següent relació:

P[N=k]=pk=(a+bk)pk1,k1.

per a alguns a i b que compleixen a+b0. El valor inicial p0 es determina de tal manera que k=0pk=1.

La recursió Panjer fa ús d'aquesta relació iterativa per especificar una manera recursiva de construir la distribució de probabilitat de S. En el següent WN(x)denota la funció generatriu de probabilitat de N (per a això, vegeu la taula de classe (a, b, 0) de les distribucions).

En el cas del número de la seqüència, tingueu en compte l'algoritme De Pril.[3] Aquest algoritme és adequat per calcular la distribució de la suma de n variables discretes aleatòries.[4]

Recursió

L'algorisme proporciona una recursió per calcular gk=P[S=hk].

El valor inicial és g0=WN(f0)amb els casos especials

g0=p0exp(f0b) if a=0,

i

g0=p0(1f0a)1+b/a for a0,

i es procedeix amb

gk=11f0aj=1k(a+bjk)fjgkj.

Exemple

El següent exemple mostra la densitat aproximada de S=i=1NXi on NNegBin(3.5,0.3) i XFrechet(1.7,1) amb ample de xarxa h = 0.04. (Vegeu distribució de Fréchet).

Com s'observa, pot sorgir un problema en la inicialització de la recursió. Guégan i Hassani (2009) han proposat una solució per abordar aquest tema.[5]

Referències

Plantilla:Referències

Enllaços externs

Plantilla:Autoritat