Hi,
I am trying to simulate claim amounts and claim numbers over 10,000 simulations. I can figure out how to determine the number of claims but I want to then use that information to simulate the amount of each claim.
The code below simulates the number of claims for 10,000 observations. However, I want to go a step further and to simulate claim amounts for each similation. For example, if the Poisson generates 212 for the first observation I want to produce 212 claim amounts with a lognormal distribution across the columns? My lognormal parameters are mu = 9.1 and shape = 1.5.m
data poisson(keep = x); call streaminit (4321); lambda = 212; do i = 1 to 10000; x = rand("Poisson",lambda); output; end; run;
The resulting dataset would be a 10000 rows by approximately 300 colums.
I would also like to impose caps on the severity but I think I can do this.
Thanks
... View more