powered by
dgpd(x, sigma, xi, u = 0, log.d = FALSE) pgpd(q, sigma, xi, u = 0, lower.tail = TRUE, log.p = FALSE) qgpd(p, sigma, xi, u = 0, lower.tail = TRUE, log.p = FALSE) rgpd(n, sigma, xi, u = 0)
x <- rgpd(1000, sigma=1, xi=.5) hist(x) x <- rgpd(1000, sigma=exp(rnorm(1000, 1, .25)), xi=rnorm(1000, .5, .2)) hist(x) plot(pgpd(x, sigma=1, xi=.5))
Run the code above in your browser using DataLab