Utility function wrapping up the d/p/q/r distribution functions
dpqrdist(dist, type = c("d", "p", "q", "r"), ...)a character description of a distribution, for example
"norm", "t", or "chisq"
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 ...
# 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