Discrete Distribution Approximations
distcrete
takes a distribution and a set of parameters
and returns a list object with functions as elements. Each one is the
equivalent to the function calls we typically expect to be able to do in
R for a given distribution:
dnorm
pnorm
qnorm
rnorm
Each function created by distcrete
corresponds to the
first letter of the non-discrete equivalent.
set.seed(415)
<- distcrete::distcrete("gamma", 1, shape = 3, w = 0)
d0 $d(1:10)
d0#> [1] 0.243022187 0.253486335 0.185086776 0.113451286 0.062683215
#> [6] 0.032332641 0.015882196 0.007521773 0.003462799 0.001558522
$p(c(.1,.5))
d0#> [1] 0.09958372 0.19115317
$q(c(.1,.5))
d0#> [1] 0 1
$r(10)
d0#> [1] 2 3 1 2 12 2 2 4 5 3
You can use any distribution that conforms to the following expectations:
p[dist name]
availableq[dist name]
availableThese can be loaded from a package or created on the fly, but must
exist when the distcrete()
function is called.
You can install it from CRAN with:
install.packages("distcrete")
You can install distcrete
the most up to date version
from github with:
# install.packages("devtools")
::install_github("reconhub/distcrete") devtools
::test()
devtools#> Loading distcrete
#> Loading required package: testthat
#> Testing distcrete
#> distcrete: .....................................................................................................................................
#> utils: ...............
#>
#> DONE ======================================================================