Utility function wrapping up the d/p/q/r distribution functions

dpqrdist(dist, type = c("d", "p", "q", "r"), ...)

Arguments

dist

a character description of a distribution, for example "norm", "t", or "chisq"

type

one of "x", "p", "q", or "r"

...

additional arguments passed on to underlying distribution function. Note that one of d, p, q, or n must be a named argument in ...

Examples

# 3 random draws from N(1,2)
dpqrdist("norm", "r", n = 3, mean = 1, sd = 2)
#>            [,1]
#> [1,]  0.9701317
#> [2,] -3.8634415
#> [3,]  1.2331878
# These should all be the same
dpqrdist("norm", "d", x = 0) == dnorm(x = 0)
#> [1] TRUE
dpqrdist("norm", "p", q = 0, mean = 1, sd = 2) == pnorm(q = 0, mean = 1, sd = 2)
#> [1] TRUE
dpqrdist("norm", "q", p = 0.5, mean = 1, sd = 2) == qnorm(p = 0.5, mean = 1, sd = 2)
#> [1] TRUE