[R] combining distributions

Hi all, I am trying to generate a series of random numbers. What I want to do is create n random Poisson numbers (easily done with the rpois function) and then use those generated numbers as the number of genenerated from another distribution (in this case, a function for the bounded Pareto

[R] creating a vector of sums

Hi, I am trying to create a function for a goodness-of-fit test for the Pareto Distribution for some loss data that I have. So far I have the following: function(X=OTOL) { n <- length(X)-1 #calculated the number of values (extra as 0 included) i <- 2:640 #values of i j <- 1:639 #values of i-1