dist <- dist_burr(shape1 = c(1,1,1,2,3,0.5), shape2 = c(1,2,3,1,1,2))
dist
#> <distribution[6]>
#> [1] Burr12(1, 1, 1) Burr12(1, 2, 1) Burr12(1, 3, 1) Burr12(2, 1, 1)
#> [5] Burr12(3, 1, 1) Burr12(0.5, 2, 1)
mean(dist)
#> [1] Inf 1.570796 1.209200 1.000000 0.500000 Inf
variance(dist)
#> [1] NaN Inf 0.9562355 Inf 0.7500000 NaN
support(dist)
#> <support_region[6]>
#> [1] [0,Inf) (0,Inf) (0,Inf) [0,Inf) [0,Inf) (0,Inf)
generate(dist, 10)
#> [[1]]
#> [1] 0.66111892 0.52207968 2.03658634 0.02095592 0.39823740 0.14596066
#> [7] 0.01700044 3.57533998 0.50482281 1.56697206
#>
#> [[2]]
#> [1] 4.5507249 0.7880165 0.8190931 1.2074297 0.4062703 0.9652413 0.1454112
#> [8] 7.6006116 0.6963448 1.3013295
#>
#> [[3]]
#> [1] 0.4469956 0.7802311 0.7955033 0.6858941 0.9443517 1.4703917 0.9882107
#> [8] 1.1177725 0.7215172 0.8282467
#>
#> [[4]]
#> [1] 0.58828240 0.02089587 0.82983909 3.46316661 0.31740138 1.14222754
#> [7] 1.81878902 0.03243623 0.11714403 0.14855035
#>
#> [[5]]
#> [1] 0.23369863 0.22289518 1.18450878 0.37064711 0.79691686 0.13126608
#> [7] 0.13986661 0.01857069 0.72088796 0.01067822
#>
#> [[6]]
#> [1] 2.3819379 1.1680588 0.7120442 14.0559715 1.6133667 0.8459186
#> [7] 2.0676876 1.5087294 4.7981248 32.2154967
#>
density(dist, 2)
#> [1] 0.11111111 0.16000000 0.14814815 0.07407407 0.03703704 0.17888544
density(dist, 2, log = TRUE)
#> [1] -2.197225 -1.832581 -1.909543 -2.602690 -3.295837 -1.721010
cdf(dist, 4)
#> [1] 0.8000000 0.9411765 0.9846154 0.9600000 0.9920000 0.7574644
quantile(dist, 0.7)
#> [1] 2.3333333 1.5275252 1.3263524 0.8257419 0.4938016 3.1797973