exp(dpareto2(1, min = 10, 3, 4, log = TRUE))
p <- (1:10)/10
ppareto2(qpareto2(p, min = 10, 2, 3), min = 10, 2, 3)
## variance
mpareto2(2, min = 10, 4, 1) - mpareto2(1, min = 10, 4, 1)^2
## case with shape - order > 0
levpareto2(10, min = 10, 3, scale = 1, order = 2)
## case with shape - order < 0
levpareto2(10, min = 10, 1.5, scale = 1, order = 2)
Run the code above in your browser using DataLab