rtrun {bayesm} | R Documentation |
rtrun
draws from a truncated univariate normal distribution.
rtrun(mu, sigma, a, b)
mu |
mean |
sigma |
standard deviation |
a |
lower bound |
b |
upper bound |
Note that due to the vectorization of the rnorm
and qnorm
commands in R, all arguments can be vectors of equal length. This makes the inverse CDF method the most efficient to use in R.
Draw (possibly a vector)
This routine is a utility routine that does not check the input arguments for proper dimensions and type.
**Note also that rtrun
is currently affected by the numerical accuracy of the inverse CDF function when trunctation points are far out in the tails of the distribution, where “far out” means |a - \mu| / \sigma > 6
and/or |b - \mu| / \sigma > 6
. A fix will be implemented in a future version of bayesm
.
Peter Rossi, Anderson School, UCLA, perossichi@gmail.com.
For further discussion, see Chapter 2, Bayesian Statistics and Marketing by Rossi, Allenby, and McCulloch.
set.seed(66)
rtrun(mu=c(rep(0,10)), sigma=c(rep(1,10)), a=c(rep(0,10)), b=c(rep(2,10)))