Page 257 - riesgo2012

8.10. Ejercicios 247

return(aux)
}
}
#...............................
#Asignaci\’on en elarreglo"gc"y graficaci\’on.
#...............................
for (i in 1:R) {
gc[i] <- g(i)
}
#Nota:Se omiteen lagr\’aficaelvalorde ladensidaden cero "gc0".
barplot(gc,main="Funci?n de densidad de S",xlab="r", ylab="g(r)")
#################################################################
#Fin dec\’odigo
#################################################################

F´ormula de Panjer en R
(Caso Poisson)

####################################
# F\’ormula de Panjer en R v1.0 #
# [Caso Poisson] #
####################################
#
R<- 20# Valorm\’aximoparar eng(r)
#
#...............................
#c\’alculodep_k=P(N=k)(CasoPoisson)
#...............................
a<- 0
b<- 3.5 #lambda
p0 <- 2.7172^{-b}
p<- array(1:R, dim=c(R))
p[1] <- (a+b)*p0
for (k in 2:R) {
p[k] <- (a+b/k)*p[k-1]
}
#...............................
#c\’alculodef_r=P(Y=r),r>=1
#...............................
#