The toolbox includes:
1. fast random number generators from the truncated univariate and multivariate student/normal distributions;
2. (Quasi-) Monte Carlo estimator of the cumulative distribution function of the multivariate student/normal;
3. accurate computation of the quantile function of the normal distribution in the extremes of its tails.
Reference:
Z. I. Botev (2017), The Normal Law Under Linear Restrictions: Simulation and Estimation via Minimax Tilting, Journal of the Royal Statistical Society, Series B, Volume 79, Part 1, pp. 1-24
Zdravko Botev (2021). Truncated Normal and Student's t-distribution toolbox (https://www.mathworks.com/matlabcentral/fileexchange/53796-truncated-normal-and-student-s-t-distribution-toolbox), MATLAB Central File Exchange. Retrieved .
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Create scripts with code, output, and formatted text in a single executable document.
Thanks for your sampling method and it really helped me a lot, but when I tried the scalar case, error always occurred saying that
Inner matrix dimensions must agree.
Error in cholperm (line 20)
tl=(l(I)-L(I,1:j-1)*z(1:j-1))./s;
Could you please take a look at it?
What can I do to sample from a truncated multivariate Gaussian with non-zero mean when the number of inequality constraints is greater than the dimensionality of the random vector (matrix A is not squared)?
Thanks for sharing. Best regards
Thanks for making this toolbox available!