Distribución de Pareto

Generación de aleatorios de pareto:

parametros k: mínimo g: gama (forma)

suppressPackageStartupMessages(library(distrEx))
k <- 5
g <- 3
u <- rpareto(10000, g, k)
head(u)

## [1] 0.57427 0.24758 1.46873 0.21410 0.01524 4.12027

dpareto(4, g, k)

## [1] 0.05716

ppareto(4, g, k)

## [1] 0.8285

qpareto(0.9, g, k)

## [1] 5.772

curve(dpareto(x, g, k), from = k, to = 7, lwd = 2, main = expression(f(x) == 
    frac(gamma * k^gamma, x^(gamma + 1))))

plot of chunk unnamed-chunk-1 

mean(u)

## [1] 2.535

sd(u)

## [1] 4.181

hist(u, freq = FALSE, xlim = c(5, 20), ylim = c(0, 0.1), breaks = 100, main = "Histograma de aleatorios de pareto u", 
    ylab = "Densidad", col = "blue")

plot of chunk unnamed-chunk-1