x <- rlomax(1e5, 5, 16)
xx <- seq(-100, 100, by = 0.001)
hist(x, 100, freq = FALSE)
lines(xx, dlomax(xx, 5, 16), col = "red")
hist(plomax(x, 5, 16))
plot(ecdf(x))
lines(xx, plomax(xx, 5, 16), col = "red", lwd = 2)
Run the code above in your browser using DataLab