Draws z ~ N(ez, 1) truncated to (lb, ub) for each cell,
working entirely on the log-probability scale in the tail that
contains the truncation interval. This stays exact for arbitrarily
extreme ez / bounds (log-scale pnorm/qnorm are
accurate to |ez| of several hundred), so no probability clamp
is needed. Infinite bounds reduce to the
one-sided case automatically.
rtnorm_interval_logp(ez, lb, ub)vector of truncated-normal draws, same length as ez.
vector of conditional means.
lower bounds (recycled to length(ez)).
upper bounds (recycled to length(ez)).
Cassy Dorff, Shahryar Minhas, Tosin Salau