rtrunc {mc2d} | R Documentation |
Provides samples from classical R distributions and mc2d specific distributions truncated between linf and lsup.
rtrunc(distr=runif, n, linf=-Inf, lsup=Inf, ...)
distr |
A function providing random data or its name as character. The function 'rdistr' should have a 'qdistr' form (with argument 'p') and a 'pdistr' form (with argument 'q'). Example : 'rnorm' (has a 'qnorm' and a 'pnorm' form), 'rbeta', 'rbinom', 'rgamma', ... |
n |
The size of the sample. |
linf |
A vector of lower bounds. |
lsup |
A vector of upper bounds. |
... |
All arguments to be passed to pdistr and qdistr. |
The function 1) evaluates the p values corresponding to linf and lsup using pdistr; 2) samples n values using runif(n, min=pinf, max=psup), and 3) takes the n corresponding quantiles from the specified distribution using qdistr.
All distributions (but sample) implemented in the stats library could be used. The arguments in ... should be named. Do not use 'log' or 'log.p' or 'lower.tail'.
A vector of n values.
The inversion of the quantile function leads to time consuming functions for some distributions.
Regis Pouillot
rtrunc("rnorm", n=10, linf=0) range(rtrunc(rnorm, n=1000, linf=3, lsup=5, sd=10))