
RtoDPQ.d(r, e = getdistrOption("RtoDPQ.e"))
RtoDPQ.d
generates $10^e$ random numbers, by default $e =$
RtoDPQ.e
which are used to produce a density, cdf and
quantile function.
Of course, the results are usually not exact as they rely on random numbers.UnivariateDistribution-class
,
density
,
approxfun
,
ecdf
rn2 <- function(n){rnorm(n)^2}
x <- RtoDPQ(r = rn2, e = 4, n = 512)
# returns density, cumulative distribution and quantile function of
# squared standard normal distribution
x$dfun(4)
RtoDPQ(r = rn2, e = 5, n = 1024) # for a better result
rp2 <- function(n){rpois(n, lambda = 1)^2}
x <- RtoDPQ.d(r = rp2, e = 5)
# returns density, cumulative distribution and quantile function of
# squared Poisson distribution with parameter lambda=1
Run the code above in your browser using DataLab